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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Solar cycle variation of interplanetary shocks, coronal mass ejections, and stream interactions observed at 0.7 AU

A survey of the Pioneer Venus Orbiter (PVO) magnetometer and plasma data from 1979-1980, shows that the occurrence frequency of interplanetary shocks, coronal mass ejections (CMEs) and stream interactions observed at 0.7 AU exhibits a solar cycle variation. As previously found at 1 AU, the observed number of both interplanetary shocks and CMEs peaks during solar maximum (approximately 16 and approximately 27 per year, respectively) and reaches a low during solar minimum (approximately 0 and approximately 7 per year, respectively), in phase with the variation in smoothed sunspot number. The number of stream interactions observed varies in the opposite manner, having a minimum during solar maximum (approximately 15 per year) and a maximum during solar minimum (approximately 34 per year). The percentage of CMEs and stream interactions producing interplanetary shocks also varies during the solar-cycle and exhibits interesting behavior during the declining phase. While the number of CMEs observed during this phase is decreasing, the percentage of CMEs producing interplanetary shocks reaches a maximum. Also, while the number of stream interactions observed is increasing, but has not reached maximum during the declining phase, the percentage of stream interactions producing interplanety shocks is at a maximum.

Lindsay, G. M.↗
NASA NTRS1993

Accurate upwind methods for the Euler equations

A new class of piecewise linear methods for the numerical solution of the one-dimensional Euler equations of gas dynamics is presented. These methods are uniformly second-order accurate, and can be considered as extensions of Godunov's scheme. With an appropriate definition of monotonicity preservation for the case of linear convection, it can be shown that they preserve monotonicity. Similar to Van Leer's MUSCL scheme, they consist of two key steps: a reconstruction step followed by an upwind step. For the reconstruction step, a monotonicity constraint that preserves uniform second-order accuracy is introduced. Computational efficiency is enhanced by devising a criterion that detects the 'smooth' part of the data where the constraint is redundant. The concept and coding of the constraint are simplified by the use of the median function. A slope steepening technique, which has no effect at smooth regions and can resolve a contact discontinuity in four cells, is described. As for the upwind step, existing and new methods are applied in a manner slightly different from those in the literature. These methods are derived by approximating the Euler equations via linearization and diagonalization. At a 'smooth' interface, Harten, Lax, and Van Leer's one intermediate state model is employed. A modification for this model that can resolve contact discontinuities is presented. Near a discontinuity, either this modified model or a more accurate one, namely, Roe's flux-difference splitting. is used. The current presentation of Roe's method, via the conceptually simple flux-vector splitting, not only establishes a connection between the two splittings, but also leads to an admissibility correction with no conditional statement, and an efficient approximation to Osher's approximate Riemann solver. These reconstruction and upwind steps result in schemes that are uniformly second-order accurate and economical at smooth regions, and yield high resolution at discontinuities.

Huynh, Hung T.↗
DOE OSTI2024

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability↗
NASA NTRS1986

Analysis of growth patterns during gravitropic curvature in roots of Zea mays by use of a computer-based video digitizer

A computer-based video digitizer system is described which allows automated tracking of markers placed on a plant surface. The system uses customized software to calculate relative growth rates at selected positions along the plant surface and to determine rates of gravitropic curvature based on the changing pattern of distribution of the surface markers. The system was used to study the time course of gravitropic curvature and changes in relative growth rate along the upper and lower surface of horizontally-oriented roots of maize (Zea mays L.). The growing region of the root was found to extend from about 1 mm behind the tip to approximately 6 mm behind the tip. In vertically-oriented roots the relative growth rate was maximal at about 2.5 mm behind the tip and declined smoothly on either side of the maximum. Curvature was initiated approximately 30 min after horizontal orientation with maximal (50 degrees) curvature being attained in 3 h. Analysis of surface extension patterns during the response indicated that curvature results from a reduction in growth rate along both the upper and lower surfaces with stronger reduction along the lower surface.

NASA Discipline Plant Biology↗
DOE OSTI2023

Customizable adaptive regularization techniques for B-spline modeling

B-spline models are a powerful way to represent scientific data sets with a functional approximation. However, these models can suffer from spurious oscillations when the data to be approximated are not uniformly distributed. Model regularization (i.e., smoothing) has traditionally been used to minimize these oscillations; unfortunately, it is sometimes impossible to sufficiently remove unwanted artifacts without smoothing away key features of the data set. In this article, we present a method of model regularization that preserves significant features of a data set while minimizing artificial oscillations. Our method varies the strength of a smoothing parameter throughout the domain automatically, removing artifacts in poorly-constrained regions while leaving other regions unchanged. Further, the proposed method selectively incorporates regularization terms based on first and second derivatives to maintain model accuracy while minimizing numerical artifacts. The behavior of our method is validated on a collection of two- and three-dimensional data sets produced by scientific simulations. In addition, a key tuning parameter is highlighted and the effects of this parameter are presented in detail. This paper is an extension of our previous conference paper at the 2022 International Conference on Computational Science (ICCS) (Lenz et al., 2022).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

System Identification for Nonlinear Control Using Neural Networks

An approach to incorporating artificial neural networks in nonlinear, adaptive control systems is described. The controller contains three principal elements: a nonlinear inverse dynamic control law whose coefficients depend on a comprehensive model of the plant, a neural network that models system dynamics, and a state estimator whose outputs drive the control law and train the neural network. Attention is focused on the system identification task, which combines an extended Kalman filter with generalized spline function approximation. Continual learning is possible during normal operation, without taking the system off line for specialized training. Nonlinear inverse dynamic control requires smooth derivatives as well as function estimates, imposing stringent goals on the approximating technique.

Stengel, Robert F.↗
DOE OSTI2021

Toward efficient elastocaloric systems: Predicting material thermal properties with high fidelity

A critical need to accurately model thermal behaviors of materials that exhibit strong elastocaloric effects, including predicting the effects themselves at varying stresses and temperatures, has been addressed using a simple and versatile physics-based approach. The key factor leading to the high precision is approximating the underlying elastic phase transition as a smooth modification of lattice entropy between coexisting phases. Once the phase transformation entropy is modeled to match experimentally measured strain as a function of temperature and applied stress, estimating the heat capacity, entropy, and isothermal entropy and adiabatic temperature changes in temperature-stress coordinates becomes straightforward. This approach provides insight into how thermal properties of elastocaloric materials vary through the transition based on strain measurements that are simple to perform and interpret. In addition to aiding in the rapid evaluation of new and existing elastocaloric materials, this advance is expected to prove invaluable for accurate heat transfer modeling aimed at designing efficient regenerative elastocaloric cooling devices.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Drag and heat transfer on surfaces with small longitudinal fins

An investigation is conducted regarding the drag characteristics of longitudinally ribbed surfaces using various configurations which attempt to confine the turbulent wall bursts to regions of small transverse extent in the initial birth region. The ratio of heat transfer increase to drag increase for the high area ratio rib configurations is determined and the heat transfer efficiency of the present surfaces is compared with data found in the literature. The drag data obtained indicate local 'wetted area skin friction' reductions but due to the increased surface area the planform drag is increased in most instances. There is no indication of drag reduction with the rectangular ribbed models that had sharp leading and fin tip edges. Heat transfer data indicate that several of the triangular grooved models have efficiencies approximately 10 percent higher than that of a smooth flat plate.

Walsh, M. J.↗
NASA NTRS1982

A grid spacing control technique for algebraic grid generation methods

A technique which controls the spacing of grid points in algebraically defined coordinate transformations is described. The technique is based on the generation of control functions which map a uniformly distributed computational grid onto parametric variables defining the physical grid. The control functions are smoothed cubic splines. Sets of control points are input for each coordinate directions to outline the control functions. Smoothed cubic spline functions are then generated to approximate the input data. The technique works best in an interactive graphics environment where control inputs and grid displays are nearly instantaneous. The technique is illustrated with the two-boundary grid generation algorithm.

Smith, R. E.↗
NASA NTRS1984

Far-infrared photometry of compact extragalactic sources: OJ 187 and BL Lac

The 50 and 100 micron emissions of OJ 287 were detected and upper limits for BL Lac were obtained. These first measurements of two BL Lac objects in the far-infrared show them to be similar to the few quasars previously observed in the far-infrared. In particular, there is no evidence for significant dust emission, and the lambda approximately 100 micron flux density fits on a smooth line joining the near-infrared and millimeter continuum fluxes. The implications of the results for models of the sources are discussed briefly.

Harvey, P. M.↗
NASA NTRS1984

Far-infrared photometry of compact extragalactic sources - OJ 287 and BL Lacertae

The 50 and 100 micron emissions of OJ 287 were detected and upper limits for BL Lac were obtained. These first measurements of two BL Lac objects in the far-infrared show them to be similar to the few quasars previously observed in the far-infrared. In particular, there is no evidence for significant dust emission, and the lambda approximately 100 micron flux density fits on a smooth line joining the near-infrared and millimeter continuum fluxes. The implications of the results for models of the sources are discussed briefly.

Harvey, P. M.↗
NASA NTRS1994

Transonic flight test of a laminar flow leading edge with surface excrescences

A flight experiment, conducted at NASA Dryden Flight Research Center, investigated the effects of surface excrescences, specifically gaps and steps, on boundary-layer transition in the vicinity of a leading edge at transonic flight conditions. A natural laminar flow leading-edge model was designed for this experiment with a spanwise slot manufactured into the leading-edge model to simulate gaps and steps like those present at skin joints of small transonic aircraft wings. The leading-edge model was flown with the flight test fixture, a low-aspect ratio fin mounted beneath an F-104G aircraft. Test points were obtained over a unit Reynolds number range of 1.5 to 2.5 million/ft and a Mach number range of 0.5 to 0.8. Results for a smooth surface showed that laminar flow extended to approximately 12 in. behind the leading edge at Mach number 0.7 over a unit Reynolds number range of 1.5 to 2.0 million/ft. The maximum size of the gap-and-step configuration over which laminar flow was maintained consisted of two 0.06-in. gaps with a 0.02-in. step at a unit Reynolds number of 1.5 million/ft.

Zuniga, Fanny A.↗
NASA NTRS2004

Representation of Ice Geometry by Parametric Functions: Construction of Approximating NURBS Curves and Quantification of Ice Roughness--Year 1: Approximating NURBS Curves

Software was developed to construct approximating NURBS curves for iced airfoil geometries. Users specify a tolerance that determines the extent to which the approximating curve follows the rough ice. The user can therefore smooth the ice geometry in a controlled manner, thereby enabling the generation of grids suitable for numerical aerodynamic simulations. Ultimately, this ability to smooth the ice geometry will permit studies of the effects of smoothing upon the aerodynamics of iced airfoils. The software was applied to several different types of iced airfoil data collected in the Icing Research Tunnel at NASA Glenn Research Center, and in all cases was found to efficiently generate suitable approximating NURBS curves. This method is an improvement over the current "control point formulation" of Smaggice (v.1.2). In this report, we present the relevant theory of approximating NURBS curves and discuss typical results of the software.

Dill, Loren H.↗
NASA NTRS2012

Low-Cost, Rugged High-Vacuum System

A need exists for miniaturized, rugged, low-cost high-vacuum systems. Recent advances in sensor technology have led to the development of very small mass spectrometer detectors as well as other analytical instruments such as scanning electron microscopes. However, the vacuum systems to support these sensors remain large, heavy, and power-hungry. To meet this need, a miniaturized vacuum system was developed based on a very small, rugged, and inexpensive-to-manufacture molecular drag pump (MDP). The MDP is enabled by a miniature, very-high-speed (200,000 rpm), rugged, low-power, brushless DC motor optimized for wide temperature operation and long life. The key advantages of the pump are reduced cost and improved ruggedness compared to other mechanical hig-hvacuum pumps. The machining of the rotor and stators is very simple compared to that necessary to fabricate rotor and stator blades for other pump designs. Also, the symmetry of the rotor is such that dynamic balancing of the rotor will likely not be necessary. Finally, the number of parts in the unit is cut by nearly a factor of three over competing designs. The new pump forms the heart of a complete vacuum system optimized to support analytical instruments in terrestrial applications and on spacecraft and planetary landers. The MDP achieves high vacuum coupled to a ruggedized diaphragm rough pump. Instead of the relatively complicated rotor and stator blades used in turbomolecular pumps, the rotor in the MDP consists of a simple, smooth cylinder of aluminum. This will turn at approximately 200,000 rpm inside an outer stator housing. The pump stator comprises a cylindrical aluminum housing with one or more specially designed grooves that serve as flow channels. To minimize the length of the pump, the gas is forced down the flow channels of the outer stator to the base of the pump. The gas is then turned and pulled toward the top through a second set of channels cut into an inner stator housing that surrounds the motor. The compressed gas then flows down channels in the motor housing to the exhaust port of the pump. The exhaust port of the pump is connected to a commercially available diaphragm or scroll pump.

Sorensen, Paul↗
NASA NTRS2012

X-Ray Spectra from MHD Simulations of Accreting Black Holes

We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, nonrotating black hole. For the first time, we are able to explain from first principles in a self-consistent way the X-ray spectra observed from stellar-mass black holes, including a thermal peak, Compton reflection hump, power-law tail, and broad iron line. Varying only the mass accretion rate, we are able to reproduce the low/hard, steep power-law, and thermal-dominant states seen in most galactic black hole sources. The temperature in the corona is T(sub e) ~ 10 keV in a boundary layer near the disk and rises smoothly to T(sub e) greater than or approximately 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to approximately equal to 6M as the accretion rate decreases, we find that the shape of the Fe Ka line is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations.

Schnittman, Jeremy D.↗
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRSDate not supplied

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗
NASA NTRS1997

Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws

In these lecture notes we describe the construction, analysis, and application of ENO (Essentially Non-Oscillatory) and WENO (Weighted Essentially Non-Oscillatory) schemes for hyperbolic conservation laws and related Hamilton- Jacobi equations. ENO and WENO schemes are high order accurate finite difference schemes designed for problems with piecewise smooth solutions containing discontinuities. The key idea lies at the approximation level, where a nonlinear adaptive procedure is used to automatically choose the locally smoothest stencil, hence avoiding crossing discontinuities in the interpolation procedure as much as possible. ENO and WENO schemes have been quite successful in applications, especially for problems containing both shocks and complicated smooth solution structures, such as compressible turbulence simulations and aeroacoustics. These lecture notes are basically self-contained. It is our hope that with these notes and with the help of the quoted references, the reader can understand the algorithms and code them up for applications.

Shu, Chi-Wang↗