Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing technique”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Unsteady transonic Euler solutions using finite elements

A finite element solution of the unsteady Euler equations is presented and demonstrated for 2D airfoil configurations oscillating in transonic flows. Computations are performed by spatially discretizing the conservation equations using the Galerkin weighted residual method and then employing a multistage Runge-Kutta scheme to march forward in time. A mesh deformation scheme has been developed to efficiently move interior points in a smooth fashion as the airfoil undergoes rigid body pitch and plunge motion. Both steady and unsteady results are presented, and a comparison is made with solutions obtained using finite-volume techniques. The effects of using either a lumped or consistent mass matrix are presented; the finite element method provides an accurate solution for unsteady transonic flows about isolated airfoils.

Davis, Gary A.↗
NASA NTRS1994

Experimental studies of transpiration cooling with shock interaction in hypersonic flow, part B

This report describes the result of experimental studies conducted to examine the effects of the impingement of an oblique shock on the flowfield and surface characteristics of a transpiration-cooled wall in turbulent hypersonic flow. The principal objective of this work was to determine whether the interaction between the oblique shock and the low-momentum region of the transpiration-cooled boundary layer created a highly distorted flowfield and resulted in a significant reduction in the cooling effectiveness of the transpiration-cooled surface. As a part of this program, we also sought to determine the effectiveness of transpiration cooling with nitrogen and helium injectants for a wide range of blowing rates under constant-pressure conditions in the absence of shock interaction. This experimental program was conducted in the Calspan 48-Inch Shock Tunnel at nominal Mach numbers of 6 and 8, for a Reynolds number of 7.5 x 10(exp 6). For these test conditions, we obtained fully turbulent boundary layers upstream of the interaction regions over the transpiration-cooled segment of the flat plate. The experimental program was conducted in two phases. In the first phase, we examined the effects of mass-addition level and coolant properties on the cooling effectiveness of transpiration-cooled surfaces in the absence of shock interaction. In the second phase of the program, we examined the effects of oblique shock impingement on the flowfield and surface characteristics of a transpiration-cooled surface. The studies were conducted for a range of shock strengths with nitrogen and helium coolants to examine how the distribution of heat transfer and pressure and the characteristics of the flowfield in the interaction region varied with shock strength and the level of mass addition from the transpiration-cooled section of the model. The effects of the distribution of the blowing rate along the interaction regions were also examined for a range of blowing rates through the transpiration-cooled panels. The regions of shockwave/boundary layer interaction examined in these studies were induced by oblique shocks generated with a sharp, flat plate, inclined to the freestream at angles of 5 degrees, 7.5 degrees, and 10 degrees. It was found that, in the absence of an incident shock, transpiration cooling was a very effective method for reducing both the heat transfer and the skin friction loads on the surface. The helium coolant was found to be significantly more effective than nitrogen, because of its low molecular weight and high specific heat. The studies of shock-wave/transpiration-cooled surface interaction demonstrated that the interaction region between the incident shock and the low-momentum transpiration-cooled boundary layer did not result in a significant increase in the size of attached or separated interaction regions, and did not result in significant flowfield distortions above the interaction region. The increase in heating downstream of the shock-impingement point could easily be reduced to the values without shock impingement by a relatively small increase in the transpiration cooling in this region. Surprisingly, this increase in cooling rate did not result in a significant increase in size of the region ahead of the incident shock or create a significantly enlarged interaction region with a resultant increase in the distortion level in the inviscid flow. Thus, transpiration cooling appears to be a very effective technique to cool the internal surfaces of scramjet engines, where shocks in the engine would induce large local increases in wall heating and create viscous/inviscid interactions that could significantly disturb the smooth flow through the combustor. However, if hydrogen is used as the coolant, burning upstream of shock impingement might result in localized hot spots. Clearly, further research is needed in this area.

Holden, Michael S.↗
NASA NTRS2003

Development of Numerical Codes for Modeling Electromagnetic Behavior at High Frequencies Near Large Objects

A study into the problem of determining electromagnetic solutions at high frequencies for problems involving complex geometries, large sizes and multiple sources (e.g. antennas) has been initiated. Typical applications include the behavior of antennas (and radiators) installed on complex conducting structures (e.g. ships, aircrafts, etc..) with strong interactions between antennas, the radiation patterns, and electromagnetic signals is of great interest for electromagnetic compatibility control. This includes the overall performance evaluation and control of all on-board radiating systems, electromagnetic interference, and personnel radiation hazards. Electromagnetic computational capability exists at NASA LaRC, and many of the codes developed are based on the Moment Method (MM). However, the MM is computationally intensive, and this places a limit on the size of objects and structures that can be modeled. Here, two approaches are proposed: (i) a current-based hybrid scheme that combines the MM with Physical optics, and (ii) an Alternating Direction Implicit-Finite Difference Time Domain (ADI-FDTD) method. The essence of a hybrid technique is to split the overall scattering surface(s) into two regions: (a) a MM zone (MMZ) which can be used over any part of the given geometry, but is most essential over irregular and "non-smooth" geometries, and (b) a PO sub-region (POSR). Currents induced on the scattering and reflecting surfaces can then be computed in two ways depending on whether the region belonged to the MMZ or was part of the POSR. For the MMZ, the current calculations proceed in terms of basis functions with undetermined coefficients (as in the usual MM method), and the answer obtained by solving a system of linear equations. Over the POSR, conduction is obtained as a superposition of two contributions: (i) currents due to the incident magnetic field, and (ii) currents produced by the mutual induction from conduction within the MMZ. This effectively leads to a reduction in the size of linear equations from N to N - Npo with N being the total number of segments for the entire surface and Npo the number of segments over the POSR. The scheme would be appropriate for relatively large, flat surfaces, and at high frequencies. The ADI-FDTD scheme provides for both transient and steady state analyses. The restrictive Courant-Friedrich-Levy (CFL) condition on the time-step is removed, and so large time steps can be chosen even though the spatial grids are small. This report includes the problem definition, a detailed discussion of both the numerical techniques, and numerical implementations for simple surface geometries. Numerical solutions have been derived for a few simple situations.

Joshi, R. P.↗
NASA NTRS1974

Renal parameter estimates in unrestrained dogs

A mathematical formulation has been developed to describe the hemodynamic parameters of a conceptualized kidney model. The model was developed by considering regional pressure drops and regional storage capacities within the renal vasculature. Estimation of renal artery compliance, pre- and postglomerular resistance, and glomerular filtration pressure is feasible by considering mean levels and time derivatives of abdominal aortic pressure and renal artery flow. Changes in the smooth muscle tone of the renal vessels induced by exogenous angiotensin amide, acetylcholine, and by the anaesthetic agent halothane were estimated by use of the model. By employing totally implanted telemetry, the technique was applied on unrestrained dogs to measure renal resistive and compliant parameters while the dogs were being subjected to obedience training, to avoidance reaction, and to unrestrained caging.

Rader, R. D.↗
NASA NTRS1982

Long term frequency stability analysis of the GPS NAVSTAR 6 Cesium clock

Time domain measurements, taken between the NAVSTAR 6 Spacecraft Vehicle (SV) and the Vandenberg Global Positioning System (GPS) Monitor Site, by a pseudo random noise receiver, were collected over an extended period of time and analyzed to estimate the long term frequency stability of the NAVSTAR 6 onboard frequency standard, referenced to the Vandenberg MS frequency standard. The technique employed separates the clock offset from the composite signal by first applying corrections for equipment delays, ionospheric delay, tropospheric delay, Earth rotation and the relativistic effect. The data are edited and smoothed using the predicted SV ephemeris to calculate the geometric delay. Then all available passes from each of the four GPS monitor stations, are collected at 1-week intervals and used to calculate the NAVSTAR orbital elements. The procedure is then completed by subtracting the corrections and the geometric delay, using the final orbital elements, from the composite signal, thus leaving the clock offset and random error.

Mccaskill, T. B.↗
NASA NTRS1989

Thrust chamber thermal barrier coating techniques

Methods for applying thermal barrier coatings to the hot-gas side wall of rocket thrust chambers in order to significantly reduce the heat transfer in high heat flux regions has been the focus of technology efforts for many years. A successful technique developed by NASA-Lewis that starts with the coating on a mandrel and then builds the thrust chamber around it by electroforming appropriate materials is described. This results in a smooth coating with exceptional adherence, as was demonstrated in hot fire rig tests. The low cycle fatigue life of chambers with coatings applied in this manner was increased dramatically compared to uncoated chambers.

Quentmeyer, Richard J.↗
NASA NTRS1996

Multilevel Algorithm for Atmospheric Data Assimilation

A multiscale algorithm for the problem of optimal statistical interpolation of observed data has been developed. This problem includes the calculation of the vector of the 'analyzed' (best estimated) atmosphere flow field w(sup a) by the formula: w(sup a) = w(sup f) + P(sup f) H(sup T) y, where the quantity y is defined by the equation (H P(sup f) H(sup T) + R)y = w(sup o) - H w(sup f), using the given model forecast first guess w(sup f) and the vector of observations w(sup o); H is an interpolation operator from the regular grid to the observation network, P(sup f) is the forecast error covariance matrix, and R is the observation error covariance matrix. At this initial stage the case of univariate analysis of single level radiosonde height data is considered. The matrix R is assumed to be diagonal, and the matrix P(sup f) is assumed to be given by the formula P(sub ij)(sup f) = sigma(sub i)(sup f) mu(sub ij) sigma(sub j)(sub f), where mu(sub ij) is a smooth, decreasing function of the distance between the i-th and the j-th points. In this paper we describe a multiscale iterative process based on a multiresolution, simultaneous displacement technique and a localized variational calculation of iteration parameters.

Brandt, Achi↗
NASA NTRS1980

Theoretical considerations of soil retention

The performance of solar energy devices is adversely affected by surface soiling, and generally, the loss of performance increases with increases in the quantity of soil retained on their surfaces. To minimize performance losses caused by soiling, solar devices should not only be deployed in low soiling geographical areas, but employ surfaces or surfacing materials having low affinity for soil retention, maximum susceptibility to be naturally cleaned by wind, rain and snow, and to be readily cleanable by simple and inexpensive maintenance cleaning techniques. This article describes known and postulated mechanisms of soil retention on surfaces, and infers from these mechanisms that low soiling and easily cleanable surfaces should have low surface energy, and be hard, smooth, hydrophobic and chemically clean of sticky materials and water soluble salts.

Cuddihy, E. F.↗
NASA NTRS1982

High spatial resolution Mg/Al maps of the western Crisium and Sulpicius Gallus regions

High spatial resolution Mg/Al ratio maps of the western Crisium and Sulpicius Gallus regions of the moon are presented. The data is from the X-ray fluorescence experiment and the image enhancement technique in the Laplacian subtraction method using a special least-squares version of the Laplacian to reduce noise amplification. In the highlands region west of Mare Crisium several relatively small patches of smooth material have high local Mg/Al ratio similar to values found in mare sites, suggesting volcanism in the highlands. In the same highland region there were other smooth areas with no high Mg/Al local values and they are probably Cayley Formation material produced by impact mass wasting. The Sulpicius Gallus region has variable Mg/Al ratios. In this region there are several high Mg/Al ratio spots, two of which occur at the highland-mare interface. Another high Mg/Al ratio area corresponds to the Sulpicius Gallus Rima I region. The high Mg/Al ratio material in the Sulpicius Gallus region is probably pyroclastic.

Schonfeld, E.↗
NASA NTRS1988

An intelligent interface for satellite operations: Your Orbit Determination Assistant (YODA)

An intelligent interface is often characterized by the ability to adapt evaluation criteria as the environment and user goals change. Some factors that impact these adaptations are redefinition of task goals and, hence, user requirements; time criticality; and system status. To implement adaptations affected by these factors, a new set of capabilities must be incorporated into the human-computer interface design. These capabilities include: (1) dynamic update and removal of control states based on user inputs, (2) generation and removal of logical dependencies as change occurs, (3) uniform and smooth interfacing to numerous processes, databases, and expert systems, and (4) unobtrusive on-line assistance to users of concepts were applied and incorporated into a human-computer interface using artificial intelligence techniques to create a prototype expert system, Your Orbit Determination Assistant (YODA). YODA is a smart interface that supports, in real teime, orbit analysts who must determine the location of a satellite during the station acquisition phase of a mission. Also described is the integration of four knowledge sources required to support the orbit determination assistant: orbital mechanics, spacecraft specifications, characteristics of the mission support software, and orbit analyst experience. This initial effort is continuing with expansion of YODA's capabilities, including evaluation of results of the orbit determination task.

Schur, Anne↗
NASA NTRS1992

An efficient and robust grid optimization algorithm

The development of an efficient and robust grid optimization is presented. This algorithm is developed by combining the best characteristics of algebraic, elliptic, and hyperbolic grid generation techniques. This development is based on the following observations and evaluations: (1) algebraic systems are fast and economical; (2) precise spacing control is always achieved; (3) grid generation by elliptic systems is always smooth; and (4) the hyperbolic system preserves the orthogonality at the solid boundary and the point distribution in the field. Computational examples representing practical internal flow configurations are presented to demonstrate the algorithm.

Soni, Bharat K.↗
NASA NTRS1992

Gravitational influences on the behavior of confined diffusion flames

Closed form, mathematical representations for the shapes of laminar diffusion flames stabilized on two-dimensional slot burners are derived for arbitrary gravitational levels. The solution technique employs the Schvab-Zeldovich conserved species variable. Effects of gravity on the velocity field are incorporated through the momentum equation; these effects, in turn, influence the species field. The results reduce smoothly to the classical Burke-Schumann solution as gravity approaches zero. It is shown that both flame height and width increase with decreasing gravity levels. Several interesting cases are considered explicitly; these include effects of negative and extremely high values of gravity. In the latter case, it is shown that the flame collapses onto the burner surface. The analysis finds application in the modeling of steady and unsteady gas and solid propellant diffusion flames.

Hegde, U.↗
NASA NTRS1992

Investigation of advanced counterrotation blade configuration concepts for high speed turboprop systems. Task 4: Advanced fan section aerodynamic analysis

The purpose of this study is the development of a three-dimensional Euler/Navier-Stokes flow analysis for fan section/engine geometries containing multiple blade rows and multiple spanwise flow splitters. An existing procedure developed by Dr. J. J. Adamczyk and associates and the NASA Lewis Research Center was modified to accept multiple spanwise splitter geometries and simulate engine core conditions. The procedure was also modified to allow coarse parallelization of the solution algorithm. This document is a final report outlining the development and techniques used in the procedure. The numerical solution is based upon a finite volume technique with a four stage Runge-Kutta time marching procedure. Numerical dissipation is used to gain solution stability but is reduced in viscous dominated flow regions. Local time stepping and implicit residual smoothing are used to increase the rate of convergence. Multiple blade row solutions are based upon the average-passage system of equations. The numerical solutions are performed on an H-type grid system, with meshes being generated by the system (TIGG3D) developed earlier under this contract. The grid generation scheme meets the average-passage requirement of maintaining a common axisymmetric mesh for each blade row grid. The analysis was run on several geometry configurations ranging from one to five blade rows and from one to four radial flow splitters. Pure internal flow solutions were obtained as well as solutions with flow about the cowl/nacelle and various engine core flow conditions. The efficiency of the solution procedure was shown to be the same as the original analysis.

Crook, Andrew J.↗
NASA NTRS1993

Field by field hybrid upwind splitting methods

A new and general approach to upwind splitting is presented. The design principle combines the robustness of flux vector splitting schemes in the capture of nonlinear waves and the accuracy of some flux difference splitting schemes in the resolution of linear waves. The new schemes are derived following a general hybridization technique performed directly at the basic level of the field by field decomposition involved in FDS methods. The scheme does not use a spatial switch to be tuned up according to the local smoothness of the approximate solution.

Coquel, Frederic↗
NASA NTRS1993

Simulations of particle acceleration in parallel shocks: Direct comparison between Monte Carlo and one-dimensional hybrid codes

We have made a direct comparison between two different computer simulations of a plane, parallel, collisionless shock including particle acceleration to energies typical of those of diffuse ions observed at the earth bow shock. Despite the fact that the one-dimensional hybrid and Monte Carlo techniques employ entirely different algorithms, they give surprisingly close agreement in the overall shapes of the complete distribution functions for protons as well as heavier ions. Both methods show that energetic ions emerge smoothly from the background thermal plasma with approximately the same relative injection rate and that the fraction of the incoming plasma's energy flux that is converted into downstream enthalpy flux of the accelerated population (i.e., the acceleration efficiency) is similar in the two cases. The fraction of the downstream proton distribution made up of superthermal particles is quite large, with at least 10% of the energy flux going into protons with energies above 10 keV. In addition, an upstream precursor, produced by backstreaming energetic particles, is present in both shocks, although the Monte Carlo precursor is considerably longer than that produced in the hybrid shock. These results offer convincing evidence that, at least in these ways, the two simulations are consistent in their description of parallel shock structure and particle acceleration, and they lay the groundwork for development of shock models employing a combination of both methods.

Ellison, Donald C.↗
NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗
NASA NTRS2019

Adjoint-Based Mesh Adaptation and Shape Optimization for Simulations with Propulsion

We demonstrate a well-posed formulation of permeable boundary conditions and mass- flow-rate functionals for adjoint-based mesh refinement and shape optimization governed by the steady Euler equations. The boundary conditions are used to model propulsion- system effects of inlets and nozzles. A two-shock diffuser with an analytic solution is used to verify the implementation. Numerical examples show that the adjoint solution is smooth at the boundary, indicating that the discretization is adjoint consistent when exit pressure is specified at subsonic outflow, and stagnation temperature and pressure at subsonic inflow. The results focus on improving simulation techniques for low-boom aircraft analysis and design. By including mass-flow-rate outputs, we obtain reliable estimates of engine flow rates concurrently with nearfield pressure signatures without increasing simulation cost. We also demonstrate the importance of mass-flow-rate constraints in shape optimization by examining trade-offs between maximizing performance of a shrouded supersonic nozzle and minimizing shocks in its nearfield.

Nemec, Marian↗