Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient 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.

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

At least 595 records · Page 33

Analysis of acoustic radiation in a jet flow environment

The analysis of sound fields from arbitrary source distributions in terms of Legendre and spherical Hankel functions is well known. The purpose of this paper is to extend this classical method of analysis to environments such as jet flows where flow and flow gradients are inherently present. The wave-equation governing the radiation of sound in such an environment is derived. The steady state flow and flow gradients in the axial and transverse directions appear as coefficients in the terms of the wave-equation. A semi-numerical method is used to solve the wave-equation in terms of modified spherical harmonics yielding the phase velocities and the directivities of an infinite set of modes. The directivity of each mode is obtained in terms of modified Legendre functions by numerical integration. Some results of these directivity and phase-velocity calculations are presented for a limited number of frequency and flow parameters. Both convective and shear refraction are shown to be important.

Mungur, P.↗

Aeroacoustics of Turbulent High-Speed Jets

Aeroacoustic noise generation in a supersonic round jet is studied to understand in particular the effect of turbulence structure on the noise without numerically compromising the turbulence itself. This means that direct numerical simulations (DNS's) are needed. In order to use DNS at high enough Reynolds numbers to get sufficient turbulence structure we have decided to solve the temporal jet problem, using periodicity in the direction of the jet axis. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. Therefore in order to answer some questions about the turbulence we will partially compromise the overall structure of the jet. The first section of chapter 1 describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. In the second section we present preliminary work done using a TVD numerical scheme on a CM5. This work is only two-dimensional (plane) but shows very interesting results, including weak shock waves. However this is a nonviscous computation and the method resolves the shocks by adding extra numerical dissipation where the gradients are large. One wonders whether the extra dissipation would influence small turbulent structures like small intense vortices. The second chapter is an extensive discussion of preliminary numerical work using the spectral method to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which are solved in O(N) steps. A very accurate highly resolved DNS of a turbulent jet flow is expected.

Rao, Ram Mohan↗

Aeroacoustic Flow Phenomena Accurately Captured by New Computational Fluid Dynamics Method

One of the challenges in the computational fluid dynamics area is the accurate calculation of aeroacoustic phenomena, especially in the presence of shock waves. One such phenomenon is "transonic resonance," where an unsteady shock wave at the throat of a convergent-divergent nozzle results in the emission of acoustic tones. The space-time Conservation-Element and Solution-Element (CE/SE) method developed at the NASA Glenn Research Center can faithfully capture the shock waves, their unsteady motion, and the generated acoustic tones. The CE/SE method is a revolutionary new approach to the numerical modeling of physical phenomena where features with steep gradients (e.g., shock waves, phase transition, etc.) must coexist with those having weaker variations. The CE/SE method does not require the complex interpolation procedures (that allow for the possibility of a shock between grid cells) used by many other methods to transfer information between grid cells. These interpolation procedures can add too much numerical dissipation to the solution process. Thus, while shocks are resolved, weaker waves, such as acoustic waves, are washed out.

Blech, Richard A.↗

Interchange Method in Compressible Magnetized Couette Flow: Magnetorotational and Magnetoconvective Instabilities

We obtain the general forms of the axisymmetric stability criteria in a magnetized compressible Couette flow using an energy variational principle, the so-called interchange or Chandrasekhar s met hod, which we applied successfully in the incompressible case. This formulation accounts for the simultaneous presence of gravity, rotation, a toroidal magnetic field, a weak axial magnetic field, entropy gradients, and density gradients in the initial equilibrium state. The power of the method lies in its simplicity which allows us to derive extremely compact and physically clear expressions for the relevant stability criteria despite the inclusion of so many physical effects. In the implementation of the method, all the applicable conservation laws are explicitly taken into account during the variations of a quantity with dimensions of energy which we call the free energy function. As in the incompressible case, the presence of an axial field invalidates the conservation laws of angular momentum and azimuthal magnetic flux and introduces instead isorotation and axial current conservation along field lines. Our results are therefore markedly different depending on whether an axial magnetic field is present, and generalize in two simple expressions all previously known, partial stability criteria for the appearance of magnetorotational instability. Furthermore, the coupling between magnetic tension and buoyancy and its influence to the dynamics of nonhomoentropic magnetized flows becomes quite clear from our results. In the limits of plane-parallel atmospheres and homoentropic flows, our formulation easily recovers the stability criteria for suppression of convective and Parker instabilities, as well as some related special cases studied over 40 years ago by Newcomb and Tserkovnikov via laborious variational techniques.

Christodoulou, Dimitris M.↗

A Preconditioning Method for Shape Optimization Governed by the Euler Equations

We consider a classical aerodynamic shape optimization problem subject to the compressible Euler flow equations. The gradient of the cost functional with respect to the shape variables is derived with the adjoint method at the continuous level. The Hessian (second order derivative of the cost functional with respect to the shape variables) is approximated also at the continuous level, as first introduced by Arian and Ta'asan (1996). The approximation of the Hessian is used to approximate the Newton step which is essential to accelerate the numerical solution of the optimization problem. The design space is discretized in the maximum dimension, i.e., the location of each point on the intersection of the computational mesh with the airfoil is taken to be an independent design variable. We give numerical examples for 86 design variables in two different flow speeds and achieve an order of magnitude reduction in the cost functional at a computational effort of a full solution of the analysis partial differential equation (PDE).

Arian, Eyal↗

Effect of temperature gradient on the optical quality of mercurous chloride crystals

Single crystals of mercurous chloride were grown at temperature gradients of 8, 11 and 17 K/cm by the physical vapor transport method. The optical quality of these crystals was evaluated by measuring bulk scattering and inhomogeneity of refractive index by birefringence interferometry. It was observed that a high temperature gradient at the solid-vapor interface induced thermal stresses and crystals showed higher scattering and irregular fringes.

Singh, N. B.↗

Two-Flux Method for Transient Radiative Transfer in a Semitransparent Layer

The two-flux method was used to obtain transient solutions for a plane layer including internal reflections and scattering. The layer was initially at uniform temperature, and was heated or cooled by external radiation and convection. The two-flux equations were examined as a means for evaluating the radiative flux gradient in the transient energy equation. Comparisons of transient temperature distributions using the two-flux method were made with results where the radiative flux gradient was evaluated from the exact radiative transfer equations. Good agreement was obtained for optical thicknesses from 0.5 to 5 and for refractive indices of 1 and 2. Illustrative results obtained with the two-flux method demonstrate the effect of isotropic scattering coupled with changing the refractive index. For small absorption with large scattering the maximum layer temperature is increased when the refractive index is increased. For larger absorption the effect is opposite, and the maximum temperature decreases with increased refractive index .

Siegel, Robert↗

Horizon Detection for Mars Surface Operations

A new image processing method for automated horizon detection is presented in this paper. This method was developed for the Mars Exploration Rovers mission’s surface operations to aid in determining orbiter visibility periods for downlink data, and is now being considered for all Mars surface missions and for a variety of applications. Starting with images received from Mars rovers or landers, the horizon pixel location is found by (1) obtaining the gradient images using the Sobel operator, (2) calculating the preliminary horizon locations by optimizing the maximum difference between sky and ground regions in the gradient domain, and (3) applying a multi-variable thresholding method. The results are analyzed using a variety of images from the Mars rovers, as an example one dataset’s results were categorized as 94.2% Good, 4.11% Okay, and 1.65% Poor. The pixel location corresponding to the found horizon in an image is converted to azimuth and elevation values using the camera model. The elevation of the terrain surrounding the rover or lander directly affects the duration orbiters are visible for downlink data volume, and when the sun is visible for image quality or solar energy consumption. Having an automated process to detect the elevation of the terrain allows the operation team to better predict the outcome of the plan, reducing heath and safety risks and allowing new operational limits to be defined.

Oij, Stephanie Lynn↗

Front Delineation and Tracking with Multiple Underwater Vehicles

This work describes a method for detecting and tracking ocean fronts using multiple autonomous underwater vehicles. Multiple vehicles — equally-spaced along the expected frontal boundary — complete near parallel transects orthogonal to the front. Lateral gradients are used to determine the location of the front crossing from each individual vehicle transect by detecting a change in the observed water property. Adaptive control of the vehicles ensure they remain perpendicular to the estimated front boundary as it evolves over time. This method was demonstrated in and around Monterey Bay, California in May of 2017. We compare the front detection method to previously used methods. We introduce a metric in order to evaluate the adaptive control techniques presented. We show the capability of this method for repeated sampling across a dynamic two-dimensional ocean front using short-range Iver AUVs. This method extends to tracking gradients of different properties using a variety of vehicles.

Chavez, Francisco P.↗

Planning and Execution for Front Delineation and Tracking with Multiple Underwater Vehicles

This work describes a planning architecture for a heterogeneous fleet of marine assets as well as a method for detecting and tracking ocean fronts using multiple autonomous underwater vehicles. Multiple vehicles — equally-spaced along the expected frontal boundary — complete near parallel transects orthogonal to the front. Lateral gradients are used to determine the location of the front crossing from each individual vehicle transect by detecting a change in the observed water property. Adaptive control of the vehicles ensure they remain perpendicular to the estimated frontal boundary as it evolves over time. This method was demonstrated in several experiment periods totaling weeks, in and around Monterey Bay, California in May and June of 2017. We discuss the challenges associated with the implementation of the planning system. We show the capability of this method for repeated sampling across a dynamic two-dimensional ocean front using a fleet of three types of platforms: short-range Iver AUVs, Tethys-Class Long-Range AUVs, and Seagliders. This method extends to tracking gradients of different properties using a variety of vehicles.

Chavez, Francisco P.↗

Optical experiments on thermophoretically augmented submicron particle deposition from 'dusty' high temperature gas flows

A real-time laser reflectivity method and Pt ribbon targets are used to obtain experimental data on the initial deposition rate of MgO(s) particles of approximately 700 nm diameter from otherwise clean combustion products as a function of target temperature (about 950-1450 K) and mainstream gas temperature (about 1500-1600 K). These preliminary data are used to demonstrate the dominant role of thermophoresis (particle drift down a temperature gradient) and to assess the utility of recently developed theoretical methods for predicting and correlating the temperature dependence of thermophoretically augmented convective-diffusion 'dust' deposition rates from flowing hot gases.

Rosner, Daniel E.↗

A Theoretical and Experimental Investigation of the Influence of Temperature Gradients on the Deformation and Burst Speeds of Rotating Disks

The purposes of this investigation were to evaluate the influence of temperature gradients and to test the validity of a recently developed method of calculating plastic flow in disks by comparing the calculated results with experimental observations. Short-time spin tests on parallel-sided, 10-inch-diameter disks were conducted under conditions that subject the disks to a range of temperatures from 70 to 1440 F. The agreement between the theoretical and experimental results was good over the range of temperature conditions investigated. Temperature gradients produced little reduction in the burst speed of the disks which had a high ductility..

LOADS AND STRESSES, STRUCTURAL↗

System identification computational considerations.

System identification methods may be accurately characterized as inverse computational problems. Three approaches to the computational problem are considered. Equation error methods are based on regression analysis techniques. Attention is given to the maximum conditional likelihood estimate and the maximum unconditional likelihood (Bayesian) estimate. Gradient dependent algorithms for solving the minimization problems are discussed, taking into account Gaussian methods, quadratic model methods, and the homogeneous function model method of Jacobson and Oksman. Some advanced methods are also examined.

Bowles, R. L.↗

Feedback control for unsteady flow and its application to the stochastic Burgers equation

The study applies mathematical methods of control theory to the problem of control of fluid flow with the long-range objective of developing effective methods for the control of turbulent flows. Model problems are employed through the formalism and language of control theory to present the procedure of how to cast the problem of controlling turbulence into a problem in optimal control theory. Methods of calculus of variations through the adjoint state and gradient algorithms are used to present a suboptimal control and feedback procedure for stationary and time-dependent problems. Two types of controls are investigated: distributed and boundary controls. Several cases of both controls are numerically simulated to investigate the performances of the control algorithm. Most cases considered show significant reductions of the costs to be minimized. The dependence of the control algorithm on the time-descretization method is discussed.

Choi, Haecheon↗

Pressure-gradient effects on hypersonic turbulent skin friction and boundary-layer profiles.

Local skin friction, total-temperature profiles, and pitot-pressure profiles were measured on the wall of a Mach-7.4 wind tunnel. The wall to adiabatic wall temperature ratio was varied from 0.3 to 0.5. Boundary-layer characteristics were compared with those predicted by a finite-difference method. Local skin friction was predicted to within 15%. Pressure-gradient effects on the temperature and Mach number distributions and the shape factor (displacement thickness/momentum thickness) were underpredicted, but the velocity distributions were closely predicted.

Hopkins, E. J.↗

Multigrid, Fractional-Step Computation Of Flow

Speed of computer code solving three-dimensional Navier-Stokes equations of flow of incompressible fluid by fractional-step method increased significantly by use of multigrid procedures. In method, equations solved on general nonorthogonal curvilinear coordinate grid, using volume fluxes. At each time step, computations performed in two fractional steps. In first step, equations of conservation of momentum solved by use of gradient of pressure from previous time step via explicit approximate-factorization method, yielding approximate flow field that does not satisfy equation of conservation of mass. In second step, discrete Poisson-like equation with Neumann-type boundary conditions, formed by combining equations of conservation of momentum and mass, solved iteratively.

Kwak, Dochan↗

Errors incurred in profile reconstruction and methods for increasing inversion accuracies for occultation type measurements

A method for augmenting the occultation experiment is described for slightly refractive media. This method which permits separation of the components of the gradient of refractivity, appears applicable to most of the planets for a major portion of their atmospheres and ionospheres. The analytic theory is given, and the results of numerical tests with a radially and angularly varying model of an ionosphere are discussed.

Gross, S. H.↗

Calculations of the pressure distribution on axisymmetric boattails including effects of viscous interactions and exhaust jets in subsonic flow

A method of calculating the pressure distributions on boattails is proposed. This method accounts for viscous effects including the presence of a separated region for base flows by combining an inviscid analysis with a boundary layer analysis in an iterative calculation. Details of the reversed flow region are not considered. Some preliminary results have been obtained for boattails at subsonic free stream Mach number with turbulent boundary layers separating at the boattail base. In some cases convergence could not be obtained using the present computer program. It is possible, in principle, to extend this method to the calculation of boattail flows with pressure gradient induced separation on the boattail.

Rom, J.↗