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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Internal gravity waves in the solar atmosphere. II - Effects of radiative damping

In the solar photosphere, temperature fluctuations associated with acoustic-gravity waves may be rapidly smoothed by the transfer mechanism of radiation between hotter and cooler regions. The present investigation of the radiative effects on internal gravity waves takes into account the parameterization of the radiative energy, employing the Newtonian cooling approximation. A linear analysis of the propagation of internal gravity waves is carried out in a model of the solar atmosphere which is taken to be homogeneous in the horizontal coordinates. Linear wave properties both with and without radiative cooling are summarized, and the variation with height of energy fluxes and of nonlinearities in the waves is discussed. Attention is given to the significance of the obtained results in terms of energy balance in the chromosphere and in relation to spectral line observations.

Mihalas, B. W.↗
NASA NTRS1994

Nonlinear instability of accelerating shock waves with application to supernovae

We consider the stability of a planar accelerating shock front in an exponential atmosphere, a situation that can be described by a self-similar solution. A previous linear instability analysis showed that there are three regimes, depending on the wavelength of the perturbation along the direction of the shock front: at long wavelengths the shock is unstable, at intermediate wavelengths it is overstable, and at short wavelenghts it is stable, where the characteristic length is set by the initial density scale height. We have carried out numerical simulations of such an accelerating shock front and have confirmed the results of the semianalytic linear analysis in the linear regime. In the nonlinear regime, the evolution again depends on the three wavelength domains. At long wavelengths, the instability continuously grows but the flow remains smooth, and at short wavelengths, the flow is stable. The intermediate wavelength, overstable regime shows more complex evolution. The growing linear oscillations saturate soon after entering the nonlinear regime and the oscillations continue with the same approximate period. The shock front develops moving points of intersection, which generate weak shock fronts and density and pressure structure in the immediate postshock flow. The density fluctuations, with a contrast of a factor of 2-3, become frozen into the downstream flow. The growth of the instability is slow, so substantial initial perturbations are needed; thesemay be present in Type II supernovae with red supergiant progenitor stars which have outer convective regions. The clumping in the outer supernova atmosphere may affect spectral line formation and may play a role in the formation of fast knots.

Luo, Ding↗
NASA NTRS2003

Light Scattering by Gaussian Particles: A Solution with Finite-Difference Time Domain Technique

The understanding of single-scattering properties of complex ice crystals has significance in atmospheric radiative transfer and remote-sensing applications. In this work, light scattering by irregularly shaped Gaussian ice crystals is studied with the finite-difference time-domain (FDTD) technique. For given sample particle shapes and size parameters in the resonance region, the scattering phase matrices and asymmetry factors are calculated. It is found that the deformation of the particle surface can significantly smooth the scattering phase functions and slightly reduce the asymmetry factors. The polarization properties of irregular ice crystals are also significantly different from those of spherical cloud particles. These FDTD results could provide a reference for approximate light-scattering models developed for irregular particle shapes and can have potential applications in developing a much simpler practical light scattering model for ice clouds angular-distribution models and for remote sensing of ice clouds and aerosols using polarized light. (copyright) 2003 Elsevier Science Ltd. All rights reserved.

Sun, W.↗
DOE OSTI2021

A minimally invasive, efficient method for propagation of full-field uncertainty in solid dynamics

In this work, we present a minimally invasive method for forward propagation of material property uncertainty to full-field quantities of interest in solid dynamics. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known a priori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the evolution equations based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the driving term of the evolution equations. The implementation is model independent, storage and communication efficient, and straightforward. We demonstrate its efficiency, accuracy, scaling with dimension of the parameter space, and convergence in distribution with two problems: a quasi-one-dimensional bar impact, and a two material notched plate impact. For the bar impact problem, we provide an analytical solution to PDF of the solution fields for method validation. With the notched plate problem, we also demonstrate good parallel efficiency and scaling of the method.

42 ENGINEERING↗
DOE OSTI2020

Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains

In this work, we propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modeling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved domains with piecewise smooth curved boundaries that may present corners. The solution method we present is based on the hybridizable discontinuous Galerkin method and sidesteps the need for geometry-conforming triangulations thanks to a transfer technique that allows to approximate the solution using only a polygonal subset as computational domain. Moreover, the solver features automatic mesh refinement driven by a residual-based a posteriori error estimator. As the mesh is locally refined, the computational domain is automatically updated in order to always maintain the distance between the actual boundary and the computational boundary of the order of the local mesh diameter. Numerical evidence is presented of the suitability of the estimator as an approximate error measure for physically relevant equilibria with pressure pedestals, internal transport barriers, and current holes on realistic geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Superhorizon isocurvature fluctuations relax tensions

Here, we present a new class of models that have potential to alleviate tensions present in the cosmological data today. We postulate about variation in the sound horizon scale on superhorizon scales, i.e., on scales that are larger than that of the present observable low-redshift Universe (≳ 1Gpc) while at the same time smaller than the largest scales probed by the cosmic microwave background (CMB) (≲ 10Gpc). In this scenario, CMB peaks are naturally smoothed as preferred by the Planck data, while at the same time the low-redshift baryon acoustic oscillation calibration is partially decoupled from the CMB. Taking superhorizon variations in baryon fraction as an example and using approximate modeling, we find improvement in the best fit Planck power spectrum model Δχ 2 ~ 6 for 1 extra degree of freedom with the relevant extension parameter 10 3 σ b = 2.10 ± 0.60, implying about 10% variations in baryon fraction across the Universe. At the same time, S 8 drops by about 1 sigma, easing tension with weak lensing surveys. We find that the combination of Planck 2018 data, eBOSS BAO data, and Riess et al. distance ladder Hubble parameter determination produce a satisfactory fit in our model if we allow for a phantom dark energy equation of state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Turbulent Processes that Influence Boundary-Layer Cloud Structure (Final Technical Report)

Many models of the atmosphere simulate vertical profiles of observable fields, such as temperature, moisture, and vertical velocity. To assess the realism of such model output, it is important to evaluate the modeled profiles against observations. However, many extant evaluations consist of comparing time averaged profiles from the model and the observations. Such time averages can smooth out unrealistic features that appear in instantaneous profiles. In this way, the time averages may obscure certain model errors from detection. For instance, many models misrepresent turbulent processes, leading to profiles that are too smoothed out or too noisy. E.g., if a model has coarse vertical resolution and overly strong diffusion, then the model may be unable to produce observed sharp peaks. On the other hand, a model may have artificial approximations that produce spurious jumps in the vertical. Much can be learned by browsing through instantaneous snapshots of model output. But the volume of output produced by a global model is far too great to rely purely on visual inspection. In addition, we need a way to statistically characterize the profiles, so that anomalies can be detected more conveniently.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1977

The relationship between mean anomaly block sizes and spherical harmonic representations

The frequently used rule specifying the relationship between a mean gravity anomaly in a block whose side length is theta degrees and a spherical harmonic representation of these data to degree l-bar is examined in light of the smoothing parameter used by Pellinen (1966). It is found that if the smoothing parameter is not considered, mean anomalies computed from potential coefficients can be in error by about 30% of the rms anomaly value. It is suggested that the above mentioned rule should be considered only a crude approximation.

Rapp, R. H.↗
NASA NTRS1979

Capturing and tracking performance of the horizontal guidance and control systems of the terminal configured vehicle

A twin-jet commercial transport equipped with digital navigation, guidance, and control systems and advanced electronic display system was used for airborne operational research. The results of flight tests which evaluated a second-order horizontal-path guidance control law and the autopilot and airplane response are presented. This evaluation was accomplished through analysis of recorded flight data and through pilot opinion of the airplane maneuvers during automatic path-capture scenarios and path tracking. Four different path captures were flown at a ground speed of approximately 160 knots and repeated at approximately 300 knots. Path tracking tracking error was measured in terms of cross track error along two paths: one pat at cruise speeds and one at airport-terminal-area speeds. The path tracking accuracy and the smoothness of the airplane maneuvers were judged satisfactory for high speeds; however, at lower speeds the control law design should be improved so that tracking will be more accurate for operations in the airport terminal area.

Knox, C. E.↗
NASA NTRS1986

Radiation pattern analysis of the tapered slot antenna, appendix 1

A theoretical model for the tapered slot antenna is presented. The model is valid for any smooth taper of the slot. The problem is solved by modeling the slot electric field and using the half plane Green's function to compute the far fields. The aperture field is obtained by affecting a stepped approximation to the continuous taper and utilizing the uniform wide slot line data. The uniform wide slot line is solved by the spectral Galerkin's technique and closed form experssions are developed for the slot wave length and the slot characteristic impedance. Numerous comparisons with measurement are made to demonstrate the versatility of the model in treating an arbitrary slot taper.

Janaswamy, Ramakrishna↗
NASA NTRS1992

Development of a 3-D upwind PNS code for chemically reacting hypersonic flowfields

Two new parabolized Navier-Stokes (PNS) codes were developed to compute the three-dimensional, viscous, chemically reacting flow of air around hypersonic vehicles such as the National Aero-Space Plane (NASP). The first code (TONIC) solves the gas dynamic and species conservation equations in a fully coupled manner using an implicit, approximately-factored, central-difference algorithm. This code was upgraded to include shock fitting and the capability of computing the flow around complex body shapes. The revised TONIC code was validated by computing the chemically-reacting (M(sub infinity) = 25.3) flow around a 10 deg half-angle cone at various angles of attack and the Ames All-Body model at 0 deg angle of attack. The results of these calculations were in good agreement with the results from the UPS code. One of the major drawbacks of the TONIC code is that the central-differencing of fluxes across interior flowfield discontinuities tends to introduce errors into the solution in the form of local flow property oscillations. The second code (UPS), originally developed for a perfect gas, has been extended to permit either perfect gas, equilibrium air, or nonequilibrium air computations. The code solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that was modified to account for real gas effects. The dissipation term associated with this algorithm is sufficiently adaptive to flow conditions that, even when attempting to capture very strong shock waves, no additional smoothing is required. For nonequilibrium calculations, the code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. This code was used to calculate the hypersonic, laminar flow of chemically reacting air over cones at various angles of attack. In addition, the flow around the McDonnel Douglas generic option blended-wing-body was computed and comparisons were made between the perfect gas, equilibrium air, and the nonequilibrium air results.

Tannehill, J. C.↗
NASA NTRS1996

A dynamic subgrid-scale model for LES of the G-equation

Turbulent combustion is a difficult subject as it must deal with all of the issues found in both turbulence and combustion. (We consider only premixed flames in this paper, but some of the ideas can be applied to the non-premixed case.) As in many other fields, there are two limiting cases that are easier to deal with than the general case. These are the situations in which the chemical time scale is either much shorter or much longer than the time scale associated with the turbulence. We deal with the former case. In this limit, the flame is thin compared to the turbulence length scales and can be idealized as an infinitely thin sheet. This is commonly called the flamelet regime; it has been the subject of many papers and the basis for many models (see, e.g., Linan & Williams 1993). In the flamelet model, the local flame structure is assumed to be identical to the laminar flame structure; thus the flame propagates normal to itself at the laminar flame speed, S(sub L). This allows the use of simple approximations. For example, one expects the rate of consumption of fuel to be proportional to the area of the flame surface. This idea allowed Damkohler (1940) to propose that the wrinkled flame could be replaced by a smooth one which travels at the turbulent flame speed, S(sub T), defined by S(sub T)/S(sub L) = A(sub L) /A(sub P) where A(sub L) is the total flame surface area and AP is the area projected onto the mean direction of propagation. This relation can be expected to be valid when the flame structure is modified only slightly by the turbulence. More recent approaches have attempted to relate the turbulent flame speed to turbulence intensity, u(sub '), which presumably, characterizes the wrinkling of the flame.

Bourlioux, A.↗
NASA NTRS1997

Near-Surface Geologic Units Exposed Along Ares Vallis and in Adjacent Areas: A Potential Source of Sediment at the Mars Pathfinder Landing Site

A sequence of layers, bright and dark, is exposed on the walls of canyons, impact craters and mesas throughout the Ares Vallis region, Chryse Planitia, and Xanthe Terra, Mars. Four layers can be seen: two pairs of alternating dark and bright albedo. The upper dark layer forms the top surface of many walls and mesas. The upper dark-bright pair was stripped as a unit from many streamlined mesas and from the walls of Ares Valles, leaving a bench at the top of the lower dark layer, approximately 250 m below the highland surface on streamlined islands and on the walls of Ares Vallis itself. Along Ares Vallis, the scarp between the highlands surface and this bench is commonly angular in plan view (not smoothly curving), suggesting that erosion of the upper dark-bright pair of layers controlled by planes of weakness, like fractures or joints. These near-surface layers in the Ares Vallis area have similar thicknesses, colors, and resistances to erosion to layers exposed near the tops of walls in Valles Marineris (Treiman et al.) and may represent the same pedogenic hardpan units. From this correlation, and from analogies with hardpans on Earth, the light-color layers may be cemented by calcite or gypsum. The dark layers are likely cemented by an iron-bearing mineral. Mars Pathfinder instruments should permit recognition and useful analyses of hardpan fragments, provided that clean uncoated surfaces are accessible. Even in hardpan-cemented materials, it should be possible to determine the broad types of lithologies in the Martian highlands. However, detailed geochemical modeling of highland rocks and soils may be compromised by the presence of hardpan cement minerals.

Treiman, Allan H.↗
NASA NTRS2003

The Importance of Chaos and Lenticulae on Europa for the JIMO Mission

The Galileo Solid State Imaging (SSI) experiment provided high-resolution images of Europa's surface allowing identification of surface features barely distinguishable at Voyager's resolution. SSI revealed the visible pitting on Europa's surface to be due to large disrupted features, chaos, and smaller sub-circular patches, lenticulae. Chaos features contain a hummocky matrix material and commonly contain dislocated blocks of ridged plains. Lenticulae are morphologically interrelated and can be divided into three classes: domes, spots, and micro-chaos. Domes are broad, upwarped features that generally do not disrupt the texture of the ridged plains. Spots are areas of low albedo that are generally smooth in texture compared to other units. Micro-chaos are disrupted features with a hummocky matrix material, resembling that observed within chaos regions. Chaos and lenticulae are ubiquitous in the SSI regional map observations, which average approximately 200 meters per pixel (m/pxl) in resolution, and appear in several of the ultra-high resolution, i.e., better than 50 m/pxl, images of Europa as well. SSI also provided a number of multi-spectral observations of chaos and lenticulae. Using this dataset we have undertaken a thorough study of the morphology, size, spacing, stratigraphy, and color of chaos and lenticulae to determine their properties and evaluate models of their formation. Geological mapping indicates that chaos and micro-chaos have a similar internal morphology of in-situ degradation suggesting that a similar process was operating during their formation. The size distribution denotes a dominant size of 4-8 km in diameter for features containing hummocky material (i.e., chaos and micro-chaos). Results indicate a dominant spacing of 15 - 36 km apart. Chaos and lenticulae are generally among the youngest features stratigraphically observed on the surface, suggesting a recent change in resurfacing style. Also, the reddish non-icy materials on Europa's surface have high concentrations in many chaos and lenticulae features. Nonetheless, a complete global map of the distribution of chaos and lenticulae is not possible with the SSI dataset. Only <20% of the surface has been imaged at 200 m/pxl or better resolution, mostly of the near-equatorial regions. Color and ultra-high-res images have much less surface coverage. Thus we suggest that full global imaging of Europa at 200 m/pxl or better resolution, preferably in multi-spectral wavelengths, should be a high priority for the JIMO mission.

Spaun, Nicole A.↗
NASA NTRS2006

Effects of Laser and Shot Peening on Fatigue Crack Growth in Friction Stir Welds

The effects of laser, and shot peening on the fatigue life of Friction Stir Welds (FSW) have been investigated. The surface roughness resulting from various peening techniques was assessed, and the fracture surfaces microstructure was characterized. Laser peening resulted in an increase in fatigue life approximately 60%, while shot peening resulted in 10% increase when compared to the unpeened material. The surface roughness of shot peening was significantly higher compared to the base material, while specimens processed with laser peening were relatively smooth.

Hatamleh, Omar↗
DOE OSTI2026

Score-based deterministic density sampling

We propose a deterministic sampling framework using Score-Based Transport Modeling for sampling an unnormalized target density π given only its score ∇ log π. Our method approximates the Wasserstein gradient flow on KL($f_t$∥π) by learning the time-varying score ∇ log $f_t$ on the fly using score matching. While having the same marginal distribution as Langevin dynamics, our method produces smooth deterministic trajectories, resulting in monotone noise-free convergence. We prove that our method dissipates relative entropy at the same rate as the exact gradient flow, provided sufficient training. Numerical experiments validate our theoretical findings: our method converges at the optimal rate, has smooth trajectories, and is often more sample efficient than its stochastic counterpart. Experiments on high-dimensional image data show that our method produces high-quality generations in as few as 15 steps and exhibits natural exploratory behavior. The memory and runtime scale linearly in the sample size.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

Quantum mechanical streamlines. I - Square potential barrier

Exact numerical calculations are made for scattering of quantum mechanical particles hitting a square two-dimensional potential barrier (an exact analog of the Goos-Haenchen optical experiments). Quantum mechanical streamlines are plotted and found to be smooth and continuous, to have continuous first derivatives even through the classical forbidden region, and to form quantized vortices around each of the nodal points. A comparison is made between the present numerical calculations and the stationary wave approximation, and good agreement is found between both the Goos-Haenchen shifts and the reflection coefficients. The time-independent Schroedinger equation for real wavefunctions is reduced to solving a nonlinear first-order partial differential equation, leading to a generalization of the Prager-Hirschfelder perturbation scheme. Implications of the hydrodynamical formulation of quantum mechanics are discussed, and cases are cited where quantum and classical mechanical motions are identical.

Hirschfelder, J. O.↗
NASA NTRS1994

Correlated environmental corrections in TOPEX/POSEIDON, with a note on ionospheric accuracy

Estimates of the effectiveness of an altimetric correction, and interpretation of sea level variability as a response to atmospheric forcing, both depend upon assuming that residual errors in altimetric corrections are uncorrelated among themselves and with residual sea level, or knowing the correlations. Not surprisingly, many corrections are highly correlated since they involve atmospheric properties and the ocean surface's response to them. The full corrections (including their geographically varying time mean values), show correlations between electromagnetic bias (mostly the height of wind waves) and either atmospheric pressure or water vapor of -40%, and between atmospheric pressure and water vapor of 28%. In the more commonly used collinear differences (after removal of the geographically varying time mean), atmospheric pressure and wave height show a -30% correlation, atmospheric pressure and water vapor a -10% correlation, both pressure and water vapor a 7% correlation with residual sea level, and a bit surprisingly, ionospheric electron content and wave height a 15% correlation. Only the ocean tide is totally uncorrelated with other corrections or residual sea level. The effectiveness of three ionospheric corrections (TOPEX dual-frequency, a smoothed version of the TOPEX dual-frequency, and Doppler orbitography and radiopositioning integrated by satellite (DORIS) is also evaluated in terms of their reduction in variance of residual sea level. Smooth (90-200 km along-track) versions of the dual-frequency altimeter ionosphere perform best both globally and within 20 deg in latitude from the equator. The noise variance in the 1/s TOPEX inospheric samples is approximately (11 mm) squared, about the same as noise in the DORIS-based correction; however, the latter has its error over scales of order 10(exp 3) km. Within 20 deg of the equator, the DORIS-based correction adds (14 mm) squared to the residual sea level variance.

Zlotnicki, V.↗