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 973 records · Page 54

Film stability in a vertical rotating tube with a core-gas flow.

The linear hydrodynamic stability of a thin-liquid layer flowing along the inside wall of a vertical tube rotating about its axis in the presence of a core-gas flow is examined. The stability problem is formulated under the conditions that the liquid film is thin, the density and viscosity ratios of gas to liquid are small and the relative (axial) pressure gradient in the gas is of the same order as gravity. The resulting eigenvalue problem is first solved by a perturbation method appropriate to axisymmetric long-wave disturbances. The damped nature (to within the thin-film and other approximations made) of the nonaxisymmetric and short-wave disturbances is noted. In view of the limitations on a truncated perturbation solution when the disturbance wavenumber is not small, an initial value method using digital computer is presented. Stability characteristics of neutral, growing, and damped modes are presented showing the influences of rotation, surface tension, and the core-gas flow. Energy balance in a neutral mode is also illustrated.

Sarma, G. S. R.↗

On laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for a number of separated flows and the accuracy of the method is verified.

Klineberg, J. M.↗

Geophysical, geochemical, and geological investigations of the Dunes geothermal system, Imperial Valley, California

The Dunes anomaly is a water-dominated geothermal system in the alluvium of the Salton Trough, lacking any surface expression. It was discovered by shallow-temperature gradient measurements. A 612-meter-deep test well encountered several temperature-gradient reversals, with a maximum of 105 C at 114 meters. The program involves surface geophysics, including electrical, gravity, and seismic methods, down-hole geophysics and petrophysics of core samples, isotopic and chemical studies of water samples, and petrological and geochemical studies of the cores and cuttings. The aim is (1) to determine the source and temperature history of the brines, (2) to understand the interaction between the brines and rocks, and (3) to determine the areal extent, nature, origin, and history of the geothermal system. These studies are designed to provide better definition of exploration targets for hidden geothermal anomalies and to contribute to improved techniques of exploration and resource assessment.

Elders, W. A.↗

Repeated-cascade theory of strong turbulence in a magnetized plasma

A two-dimensional Navier-Stokes equation of vorticity in fluid turbulence is used to model drift turbulence in a plasma with a strong constant magnetic field and a constant mean density gradient. The nonlinear eddy diffusivity is described by a time-integrated Lagrangian correlation of velocities, and the repeated-cascade method is employed to choose the rank accounting for nearest-neighbor interactions, to calculate the Lagrangian correlation, and to close the correlation hierarchy. As a result, the diffusivity becomes dependent on the plasma's induced diffusion and is represented by a memory chain that is cut off by similarity and inertial randomization. Spectral laws relating the kinetic-energy spectrum to the -5, -5/2, -3, and -11 powers of wavenumber are derived for the velocity subranges of production, approach to inertia, inertia, and dissipation, respectively. It is found that the diffusivity is proportional to some inverse power of the magnetic field, that power being 1, 2/3, 5/6, and 2, respectively, for the four velocity subranges.

Tchen, C. M.↗

Practical aspects of silicon float zone crystal growth

The float zone (FZ) growth method for silicon crystallization has advantages and disadvantages when compared to the other principal growth method, the Czochralski (CZ) technique. CZ growth is attractive for the following reasons: large diameters, high single crystal yields, low thermal gradients, low operator skill requirements, efficient power utilization, and low feed material costs. Some areas where FZ growth is advantageous include high purity, axial resistivity uniformity (on a macroscale), visibility of the growth region, low consumable material costs, and high growth rates.

Ciszek, T. T.↗

An algorithm for obtaining bang-bang control laws

A new computational method for solving control problems, when it is desired that the control action be bang-bang is presented. A gradient-based algorithm is developed where the switching times are updated to minimize the final state missed distance. The algorithm explains how a change in a switching time propagates through the other switching times and affects the final state. The algorithm is applicable to linear and nonlinear problems, including multiinput systems.

Wen, J.↗

Computational technique for compressible vortex flows using the integral equation solution

The steady full-potential equation is written in the form of Poisson's equation, and the solution for the velocity field is expressed in terms of an integral equation. The integral solution consists of two surface integrals and one volume integral. The solution is obtained through successive iteration cycles. Each cycle of iteration consists of two sub-cycles, an inner cycle for wake relaxation and an out cycle for the strength of the source distribution integrals representing the flow compressibility. The density gradients in the source distribution is computed by using a type-differencing scheme of the Murman-Cole type. The method is applied to delta wings and the numerical examples show that a curved shock is captured on the wing suction side beneath the leading edge vortex sheet. Recently, a modified version of the scheme was applied to rectangular wings. In this modified scheme, the surface integral terms were computed by using a bilinear distribution of vorticity on triangular vortex panels which represent the wing and its wake. The results were compared with the available experimental data and they are in good agreement.

Kandil, Osama A.↗

Recent experiences using finite-element-based structural optimization

Structural optimization has been available to the structural analysis community as a tool for many years. The popular use of displacement method finite-element techniques to analyze linearly elastic structures has resulted in an ability to calculate the weight and constraint gradients inexpensively for numerical optimization of structures. Here, recent experiences in the investigation and use of structural optimization are discussed. In particular, experience with the commercially available ADS/NASOPT code is addressed. An overview of the ADS/NASOPT procedure and how it was implemented is given. Two example problems are also discussed.

Paul, B. K.↗

Natural convection between concentric spheres

A finite-difference solution for steady natural convective flow in a concentric spherical annulus with isothermal walls has been obtained. The stream function-vorticity formulation of the equations of motion for the unsteady axisymmetric flow is used; interest lying in the final steady solution. Forward differences are used for the time derivatives and second-order central differences for the space derivatives. The alternating direction implicit method is used for solution of the discretization equations. Local one-dimensional grid adaptation is used to resolve the steep gradients in some regions of the flow at large Rayleigh numbers. The break-up into multi-cellular flow is found at high Rayleigh numbers for air and water, and at significantly low Rayleigh numbers for liquid metals. Excellent agreement with previous experimental and numerical data is obtained.

Garg, Vijay K.↗

Some Features of Artificially Thickened Fully Developed Turbulent Boundary Layers with Zero Pressure Gradient

Report gives an account of an investigation conducted to determine the feasibility of artificially thickening a turbulent boundary layer on a flat plate. A description is given of several methods used to thicken artificially the boundary layer. It is shown that it is possible to do substantial thickening and obtain a fully developed turbulent boundary layer, which is free from any distortions introduced by the thickening process, and, as such, is a suitable medium for fundamental research.

Klebanoff, P S↗

Comparing the Performance of Two Dynamic Load Distribution Methods

Parallel processing of symbolic computations on a message-passing multi-processor presents one challenge: To effectively utilize the available processors, the load must be distributed uniformly to all the processors. However, the structure of these computations cannot be predicted in advance. go, static scheduling methods are not applicable. In this paper, we compare the performance of two dynamic, distributed load balancing methods with extensive simulation studies. The two schemes are: the Contracting Within a Neighborhood (CWN) scheme proposed by us, and the Gradient Model proposed by Lin and Keller. We conclude that although simpler, the CWN is significantly more effective at distributing the work than the Gradient model.

Kale, L. V.↗

Frontal Polymerization in Microgravity

Frontal polymerization systems, with their inherent large thermal and compositional gradients, are greatly affected by buoyancy-driven convection. Sounding rocket experiments allowed the preparation of benchmark materials and demonstrated that methods to suppress the Rayleigh-Taylor instability in ground-based research did not significantly affect the molecular weight of the polymer. Experiments under weightlessness show clearly that bubbles produced during the reaction interact very differently than under 1 g.

Pojman, John A.↗

Flame Synthesis Of Single-Walled Carbon Nanotubes And Nanofibers

Carbon nanotubes are widely sought for a variety of applications including gas storage, intercalation media, catalyst support and composite reinforcing material [1]. Each of these applications will require large scale quantities of CNTs. A second consideration is that some of these applications may require redispersal of the collected CNTs and attachment to a support structure. If the CNTs could be synthesized directly upon the support to be used in the end application, a tremendous savings in post-synthesis processing could be realized. Therein we have pursued both aerosol and supported catalyst synthesis of CNTs. Given space limitations, only the aerosol portion of the work is outlined here though results from both thrusts will be presented during the talk. Aerosol methods of SWNT, MWNT or nanofiber synthesis hold promise of large-scale production to supply the tonnage quantities these applications will require. Aerosol methods may potentially permit control of the catalyst particle size, offer continuous processing, provide highest product purity and most importantly, are scaleable. Only via economy of scale will the cost of CNTs be sufficient to realize the large-scale structural and power applications on both earth and in space. Present aerosol methods for SWNT synthesis include laser ablation of composite metalgraphite targets or thermal decomposition/pyrolysis of a sublimed or vaporized organometallic [2]. Both approaches, conducted within a high temperature furnace, have produced single-walled nanotubes (SWNTs). The former method requires sophisticated hardware and is inherently limited by the energy deposition that can be realized using pulsed laser light. The latter method, using expensive organometallics is difficult to control for SWNT synthesis given a range of gasparticle mixing conditions along variable temperature gradients; multi-walled nanotubes (MWNTs) are a far more likely end products. Both approaches require large energy expenditures and produce CNTs at prohibitive costs, around $500 per gram. Moreover these approaches do not possess demonstrated scalability. In contrast to these approaches, flame synthesis can be a very energy efficient, low-cost process [3]; a portion of the fuel serves as the heating source while the remainder serves as reactant. Moreover, flame systems are geometrically versatile as illustrated by innumerable boiler and furnace designs. Addressing scalability, flame systems are commercially used for producing megatonnage quantities of carbon black [4]. Although it presents a complex chemically reacting flow, a flame also offers many variables for control, e.g. temperature, chemical environment and residence times [5]. Despite these advantages, there are challenges to scaling flame synthesis as well.

Wal, Randy L. Vander↗

Optimal Control within the Context of Multidisciplinary Design, Analysis, and Optimization

Multidisciplinary design, analysis and optimization involves modeling the interactions of complex systems across a variety of disciplines. The optimization of such systems can be a computationally expensive exercise with multiple levels of nested nonlinear solvers running under an optimizer.The application of optimal control in project development often involves performing trajectory optimization for fixed vehicle designs or parametric sweeps across some key vehicle properties.This information is then relayed to the subsystem design teams who update their designs and relay some bulk characteristics back to the trajectory optimization procedure.This iteration is then repeated until the design closes.However, with increasing interest in more tightly coupled systems, such as electric and hybrid-electric aircraft propulsion and boundary layer ingestion, this process is prone to ignore subtle coupling between vehicle subsystem designs and vehicle operation on a given mission.Integrating trajectory optimization into a tightly coupled multidisciplinary design procedure can be computationally prohibitive, depending on the complexity of the subsystem analyses and the optimal control technique applied.To address these issues a new optimal control software tool, Dymos, has been developed.Dymos is built upon NASA's OpenMDAO software and can leverage its capabilities to efficiently compute gradients for the optimization and optimize complex models in parallel on distributed memory systems.This report provides some explanation into the numerical methods employed in Dymos and provides several use cases that demonstrate its performance on traditional optimal control problems and improvements ino techniques have been used extensively in recent decades to solve a variety of optimal control problems, typically in the form of aerospace vehicle trajectory optimization.

pseudospectral↗

3D Propagation of Relativistic Solar Protons through Interplanetary Space

Context. Solar Energetic Particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce Ground Level Enhancements (GLEs). The traditional approach to interpreting and modelling GLE observations assumes particle propagation only parallel to the magnetic field lines of interplanetary space, i.e. it is spatially 1D. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the ~100 MeV-1 GeV range at high spectral resolution. Aims. We model the transport of GLE-energy solar protons using a 3D approach, to assess the effect of the Heliospheric Current Sheet (HCS) and drifts associated to the gradient and curvature of the Parker spiral. We derive 1 AU observables and compare the simulation results with data from PAMELA. Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used to compare with PAMELA data. Results. Spatial patterns of 1 AU crossings and the average number of crossings are strongly influenced by 3D effects, with significant differences between periods of A+ and A- polarities. The decay time constant of 1 AU intensity profiles varies depending on the position of the observer and is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, resulting in a rollover in the spectrum.

Solar energetic particles↗

3D propagation of relativistic solar protons through interplanetary space

Context.Solar energetic particles (SEPs) with energy in the GeV range can propagate to Earth from their acceleration region near the Sun and produce ground level enhancements (GLEs). The traditional approach to interpreting and modelling GLE observations assumes particle propagation which is only parallel to the magnetic field lines of interplanetary space, that is, spatially 1D propagation. Recent measurements by PAMELA have characterised SEP properties at 1 AU for the∼100 MeV–1 GeV range at high spectral resolution. Aims. We model the transport of GLE-energy solar protons using a 3D approach to assess the effect of the heliospheric current sheet(HCS) and drifts associated to the gradient and curvature of the Parker spiral. We derive 1 AU observables and compare the simulation results with data from PAMELA. Methods. We use a 3D test particle model including a HCS. Monoenergetic populations are studied first to obtain a qualitative picture of propagation patterns and numbers of crossings of the 1 AU sphere. Simulations for power law injection are used to derive intensity profiles and fluence spectra at 1 AU. A simulation for a specific event, GLE 71, is used for comparison purposes with PAMELA data. Results. Spatial patterns of 1 AU crossings and the average number of crossings per particle are strongly influenced by 3D effects, with significant differences between periods of A+ and A−polarities. The decay time constant of 1 AU intensity profiles varies depending on the position of the observer and it is not a simple function of the mean free path as in 1D models. Energy dependent leakage from the injection flux tube is particularly important for GLE energy particles, resulting in a rollover in the spectrum.

S. Dalla↗

Backshell Thermal Protection Materials for Additive Manufacturing

Leveraging AM technology focused on the Fused Deposition Modeling (FDM) process could automate the heat shield manufacturing process and allow it to reduce cost and build time and improve the design to 3D print and process TPS materials directly on the structure in monolith segments. Furthermore, AM allows the production of complicated three-dimensional geometry and density gradient layers without custom molds or special tools. NASA invests in additive manufacturing (AM) to enable new mission architectures and design methods; AM is suited for backshell thermal protection materials that require many complex closeouts that lead to integration challenges. Current work will review backshell thermal protection material formulation, development, and characterization.3D printer was used to 3D print arc-jet test coupons. The coupons were tested at the NASA Ames arc jet facility, and the results show promising properties compared with the heritage materials.

Tane Boghozian↗

Acoustic propagation in rigid ducts with blockage

Acoustic levitation has been suggested for moving nonmagnetic material in furnaces for heat processing in space experiments. Basically, acoustic standing waves under resonant conditions are excited in the cavity of the furnace while the material blockage is located at a pressure node and thus at a maximum gradient. The position of the blockage is controlled by displacing the node as a result of frequency change. The present investigation is concerned with the effect of blockage on the longitudinal and transverse resonances of a cylindrical cavity, taking into account the results of a one-dimensional and three-dimensional (3-D) analysis. Based on a Green's function surface element method, 3-D analysis is tested experimentally and proved to be accurate over a wide range of geometric parameters and boundary shapes. The shift in resonance depends on the change in pressure gradient and duct shortening caused by the blockage.

El-Raheb, M.↗