SEARCH · Engineering Papers

Results for “Smoothing methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1975

Application of Fourier transforms for microwave radiometric inversions

Existing microwave radiometer technology now provides a suitable method for remote determination of the ocean surface's absolute brightness temperature. To extract the brightness temperature of the water from the antenna temperature, an unstable Fredholm integral equation of the first kind is solved. Fourier transform techniques are used to invert the integral after it is placed into a cross correlation form. Application and verification of the methods to a two-dimensional modeling of a laboratory wave tank system are included. The instability of the ill-posed Fredholm equation is examined and a restoration procedure is included which smooths the resulting oscillations. With the recent availability and advances of fast Fourier transform (FFT) techniques, the method presented becomes very attractive in the evaluation of large quantities of data.

Holmes, J. J.↗

NASA NTRS2019

Efficient computation of optimal low thrust gravity perturbed orbit transfers

We have developed a new method for solving low-thrust fuel-optimal orbit transfer problems in the vicinity of a large body (planet or asteroid), considering a high-fidelity spherical harmonic gravity model. The algorithm is formulated via the indirect optimization method, leading to a two-point boundary value problem (TPBVP). We make use of a hyperbolic tangent smoothing law for performing continuation on the thrust magnitude to reduce the sharpness of the control switches in early iterations and thus promote convergence. The TPBVP is solved using the method of particular solutions (MPS) shooting method and Picard-Chebyshev numerical integration. Application of Picard-Chebyshev integration affords an avenue for increased efficiency that is not available with step-by-step integrators. We demonstrate that computing the particular solutions with only a low-fidelity force model greatly increases the efficiency of the algorithm while ultimately achieving near machine precision accuracy. A salient feature of the MPS is that it is parallelizable, and thus further speedups are available. It is also shown that, for near Earth orbits and over a small number of en-route revolutions around the Earth, only the zonal perturbation terms are required in the costate equations to obtain a solution that is accurate to machine precision and optimal to engineering precision. The proposed framework can be used for trajectory design around small asteroids and also for orbit debris rendezvous and removal tasks.

Junkins, John L.↗

NASA NTRS1988

Convergence speeding up in the calculation of the viscous flow about an airfoil

A finite volume method to solve the three dimensional Navier-Stokes equations was developed. It is based on a cell-vertex scheme with central differences and explicit Runge-Kutta time steps. A good convergence for a stationary solution was obtained by the use of local time steps, implicit smoothing of the residues, a multigrid algorithm, and a carefully controlled artificial dissipative term. The method is illustrated by results for transonic profiles and airfoils. The method allows a routine solution of the Navier-Stokes equations.

Radespiel, R.↗

NASA NTRS2018

NASA ACC High Energy Dynamic Impact Methodology and Outcomes

High fidelity analysis methods known as Progressive Damage Analysis (PDA) methods, capable of reliably predicting the onset and progression of damage in composite materials, are being developed for HEDI event simulation. Four state-of-the-art PDA methods are being investigated: LS-DYNA MAT162, LS-DYNA MAT261, Smooth-Particle Galerkin (SPG), and EMU Peridynamics. Validation of these PDA models aims to follow a building block approach, starting with coupon- and component-level development. As part of the overall effort, material characterization testing was performed to bridge gaps between existing experimental data and the material property inputs required to predict ballistic impact behavior at the component-level. The material models were updated with the results of the coupon-level testing and were then used to predict panel behavior and damage in ballistic impact events. To assess the accuracy of the PDA methods, these pre-test predictions were compared to the results of ballistic impact testing.

Hunziker, Kenneth J.↗

NASA NTRS1990

Heat transfer measurements from a NACA 0012 airfoil in flight and in the NASA Lewis icing research tunnel

Local heat transfer coefficients from a smooth and roughened NACA 0012 airfoil were measured using a steady state heat flux method. Heat transfer measurements on the specially constructed 0.533 meter chord airfoil were made both in flight on the NASA Lewis Twin Otter Research Aircraft and in the NASA Lewis Icing Research Tunnel (IRT). Roughness was obtained by the attachment of small, 2 mm diameter, hemispheres of uniform size to the airfoil surface in four distinct patterns. The flight data was taken for the smooth and roughened airfoil at various Reynolds numbers based on chord in the range of 1.24x10(exp 6) to 2.50x10(exp 6) and at various angles of attack up to 4 degrees. During these flight tests the free stream velocity turbulence intensity was found to be very low (less than 0.1 percent). The wind tunnel data was taken in the Reynolds number range of 1.20x10(exp 6) to 4.52x10(exp 6) and at angles of attack from -4 degrees to +8 degrees. The turbulence intensity in the IRT was 0.5 to 0.7 percent with the cloud making spray off. Results for both the flight and tunnel tests are presented as Frossling number based on chord versus position on the airfoil surface for various roughnesses and angle of attack. A table of power law curve fits of Nusselt number as a function of Reynolds number is also provided. The higher level of turbulence in the IRT versus flight had little effect on heat transfer for the lower Reynolds numbers but caused a moderate increase in heat transfer at the higher Reynolds numbers. Turning on the cloud making spray air in the IRT did not alter the heat transfer. Roughness generally increased the heat transfer by locally disturbing the boundary layer flow. Finally, the present data was not only compared with previous airfoil data where applicable, but also with leading edge cylinder and flat plate heat transfer values which are often used to estimate airfoil heat transfer in computer codes.

Poinsatte, Philip E.↗

NASA NTRS1991

A Runge-Kutta discontinuous finite element method for high speed flows

A Runge-Kutta discontinuous finite element method is developed for hyperbolic systems of conservation laws in two space variables. The discontinuous Galerkin spatial approximation to the conservation laws results in a system of ordinary differential equations which are marched in time using Runge-Kutta methods. Numerical results for the two-dimensional Burger's equation show that the method is (p+1)-order accurate in time and space, where p is the degree of the polynomial approximation of the solution within an element and is capable of capturing shocks over a single element without oscillations. Results for this problem also show that the accuracy of the solution in smooth regions is unaffected by the local projection and that the accuracy in smooth regions increases as p increases. Numerical results for the Euler equations show that the method captures shocks without oscillations and with higher resolution than a first-order scheme.

Bey, Kim S.↗

NASA NTRS2007

Solar Cycle 24 and the Solar Dynamo

We will discuss the polar field precursor method for solar activity prediction, which predicts cycle 24 will be significantly lower than recent activity cycles, and some new ideas rejuvenating Babcock's shallow surface dynamo. The polar field precursor method is based on Babcock and Leighton's dynamo models wherein the polar field at solar minimum plays a major role in generating the next cycle's toroidal field and sunspots. Thus, by examining the polar fields of the Sun near solar minimum, a forecast for the next cycle's activity is obtained. With the current low value for the Sun's polar fields, this method predicts solar cycle 24 will be one of the lowest in recent times, with smoothed F10.7 radio flux values peaking near 135 plus or minus 35 (2 sigma), in the 2012-2013 timeframe (equivalent to smoothed Rz near 80 plus or minus 35 [2 sigma]). One may have to consider solar activity as far back as the early 20th century to find a cycle of comparable magnitude. We discuss unusual behavior in the Sun's polar fields that support this prediction. Normally, the solar precursor method is consistent with the geomagnetic precursor method, wherein geomagnetic variations are thought to be a good measure of the Sun's polar field strength. Because of the unusual polar field, the Earth does not appear to be currently bathed in the Sun's extended polar field (the interplanetary field), hence negating the primal cause behind the geomagnetic precursor technique. We also discuss how percolation may support Babcock's original shallow solar dynamo. In this process ephemeral regions from the solar magnetic carpet, guided by shallow surface fields, may collect to form pores and sunspots.

Pesnell, W. D.↗

NASA NTRS1989

Friction factor data for flat plate tests of smooth and honeycomb surfaces

Friction factors for honeycomb surfaces were measured with a flat plate tester. The flat plate test apparatus was described and a method was discussed for determining the friction factor experimentally. The friction factor model was developed for the flat plate test based on the Fanno Line Flow. The comparisons of the friction factor were plotted for smooth surfaces and six-honeycomb surfaces with three-clearances, 6.9 bar to 17.9 bar range of inlet pressures, and 5,000 to 100,000 range of the Reynolds number. The optimum geometries for the maximum friction factor were found as a function of cell width to cell depth and cell width to clearance ratios.

Ha, Tae Woong↗

NASA NTRS1992

Friction-factor data for flat-plate tests of smooth and honeycomb surfaces

Ha, T. W.↗

NASA NTRS1988

Vortex methods for separated flows

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented in an elementary fashion and includes the relationship with traditional point-vortex studies, the convergence to smooth solutions of the Euler equations, and the essential differences between two- and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. The overlap with the excellent review articles available is kept to a minimum and more emphasis is placed on the area of expertise, namely two-dimensional flows around bluff bodies. When solid walls are present, complete mathematical models are not available and a more heuristic attitude must be adopted. The imposition of inviscid and viscous boundary conditions without conformal mappings or image vortices and the creation of vorticity along solid walls are examined in detail. Methods for boundary-layer treatment and the question of the Kutta condition are discussed. Practical aspects and tips helpful in creating a method that really works are explained. The topics include the robustness of the method and the assessment of accuracy, vortex-core profiles, timemarching schemes, numerical dissipation, and efficient programming. Calculations of flows past streamlined or bluff bodies are used as examples when appropriate.

Spalart, Philippe R.↗

NASA NTRS1997

Direct Numerical Simulation of Incompressible Pipe Flow Using a B-Spline Spectral Method

A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Loulou, Patrick↗

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.↗

NASA NTRS2010

Alternative Modal Basis Selection Procedures for Nonlinear Random Response Simulation

Three procedures to guide selection of an efficient modal basis in a nonlinear random response analysis are examined. One method is based only on proper orthogonal decomposition, while the other two additionally involve smooth orthogonal decomposition. Acoustic random response problems are employed to assess the performance of the three modal basis selection approaches. A thermally post-buckled beam exhibiting snap-through behavior, a shallowly curved arch in the auto-parametric response regime and a plate structure are used as numerical test articles. The results of the three reduced-order analyses are compared with the results of the computationally taxing simulation in the physical degrees of freedom. For the cases considered, all three methods are shown to produce modal bases resulting in accurate and computationally efficient reduced-order nonlinear simulations.

Przekop, Adam↗

NASA NTRS1972

Locally nonsimilar solutions for radiating shock layer about smooth axisymmetric bodies

Locally nonsimilar solutions for nongray radiating shock layers about smooth axisymmetric bodies have been obtained based on a newly developed approximate method. Good agreement is found with numerical solutions for inviscid cases (both radiating and nonradiating) and with series solutions for the radiating viscous case. For the inviscid case, the effect of radiative cooling is to destroy the entropy layer; at a distance far from the stagnation point, the shock layer is nearly isothermal. For the viscous case, the radiative wall flux approaches that of the inviscid case at a distance far downstream of the stagnation point. The method can also treat surface mass injection.

Chou, Y. S.↗

NASA NTRS1973

Remote sensing of the turbulence characteristics of a planetary atmosphere by radio occultation of a space probe.

The purpose of this paper is to analyze the effects of small-scale turbulence on radio waves propagating through a planetary atmosphere. The analysis provides a technique for inferring the turbulence characteristics of a planetary atmosphere from the radio signals received from a spacecraft as it is occulted by the planet. The planetary turbulence is assumed to be localized and smoothly varying, with the structure constant varying exponentially with altitude. Rytov's method is used to derive the variance of log-amplitude and phase fluctuations of a wave propagating through the atmosphere.

Woo, R.↗

NASA NTRS1981

A procedure for designing forebodies with constraints on cross-section shape and axial area distribution

A method is described for designing a forebody with cross sections which vary smoothly from an initial prescribed nose shape to a different prescribed base shape in such a way that the cross-section areas conform to a preassigned axial area distribution. It is shown that these conditions can be satisfied with a remaining degree of freedon, which can be used to accomplish a modest amount of geometric or pressure tailoring of the forebody. An example is provided which involves modifying the pressure distribution along a given meridian line of the forebody.

Barger, R. L.↗

NASA NTRS1983

Grid generation in three dimensions by Poisson equations with control of cell size and skewness at boundary surfaces

An algorithm for generating computational grids about arbitrary three-dimensional bodies is developed. The elliptic partial differential equation (PDE) approach developed by Steger and Sorenson and used in the NASA computer program GRAPE is extended from two to three dimensions. Forcing functions which are found automatically by the algorithm give the user the ability to control mesh cell size and skewness at boundary surfaces. This algorithm, as is typical of PDE grid generators, gives smooth grid lines and spacing in the interior of the grid. The method is applied to a rectilinear wind-tunnel case and to two body shapes in spherical coordinates.

Sorenson, R. L.↗

NASA NTRS1986

Flexible-Rotor Balancing Demonstration

Report describes method for balancing high-speed rotors at relatively low speeds and discusses demonstration of method on laboratory test rig. Method ensures rotor brought up to speeds well over 20,000 r/min smoothly, without excessive vibration amplitude at critical speeds or at operating speed.

Giordano, J.↗