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 1,099 records · Page 61

Modeling of transient two-component flow using a four-point implicit method

The four-point, centered implicit scheme that is extensively used in open channel flow simulation is shown to be applicable to rapid and slow pressure transient problems in conduits with nearly single phase and two-phase flows. It is only necessary to choose the proper weighting factor value, theta, of the Courant number. For rapid pressure transients such as waterhammer, the implicit method can yield reasonable results with limited numerical dispersion and attenuation if theta is only slightly greater than the critical value of 0.5. For slower pressure gradients in single and two-phase flows, reasonable numerical solutions may be achieved for Courant number values as high as 20.

Wiggert, D. C.↗

Design and development of a facility for compressible dynamic stall studies of a rapidly pitching airfoil

A facility for the study of dynamic stall of an airfoil undergoing a transient ramp-type pitching motion is described. The facility can produce pitch rates of 3600 deg/sec to an angle of attack of 60 deg by using a specially designed hydraulic drive with feedback control. The ramp motion generator can also generate arbitrary motion of the airfoil and thus can simulate an arbitrary aircraft maneuver. A unique airfoil support system allows unobstructed flow visualization including the complete airfoil contour, thus permitting the use of nonintrusive optical diagnostic methods for flow measurement close to the surface as well as simultaneous far-field measurements. Schlieren pictures obtained during the study reveal the instantaneous density gradients associated with dynamic stall, even under conditions of very low Mach numbers.

Chandrasekhara, M. S.↗

Analytical estimates of radial segregation in Bridgman growth from low-level steady and periodic accelerations

Approximate analytical solutions are developed that are mathematically tractable but provide the elements of CFD models for estimating convective flows subjected to low-level accelerations. The estimates are applicable to dilute systems that have horizontal temperature gradients in the vertical Bridgman configuration. The Navier-Stokes momentum equation is solved in 1D and subsequently in 2D by the first-order perturbation method. The analysis used is the case of plane-front solidification of a dilute system in the vertical thermally stable configuration. The model is found to predict the degree of radial segregation within a factor of two for a range of material and processing parameters. The method provides detailed information on the effects of processing on solute distribution in the grown crystal which are of interest in the orbital experiments emphasizing the control of radial segregation.

Naumann, Robert J.↗

Three-dimensional aerodynamic shape optimization using discrete sensitivity analysis

An aerodynamic shape optimization procedure based on discrete sensitivity analysis is extended to treat three-dimensional geometries. The function of sensitivity analysis is to directly couple computational fluid dynamics (CFD) with numerical optimization techniques, which facilitates the construction of efficient direct-design methods. The development of a practical three-dimensional design procedures entails many challenges, such as: (1) the demand for significant efficiency improvements over current design methods; (2) a general and flexible three-dimensional surface representation; and (3) the efficient solution of very large systems of linear algebraic equations. It is demonstrated that each of these challenges is overcome by: (1) employing fully implicit (Newton) methods for the CFD analyses; (2) adopting a Bezier-Bernstein polynomial parameterization of two- and three-dimensional surfaces; and (3) using preconditioned conjugate gradient-like linear system solvers. Whereas each of these extensions independently yields an improvement in computational efficiency, the combined effect of implementing all the extensions simultaneously results in a significant factor of 50 decrease in computational time and a factor of eight reduction in memory over the most efficient design strategies in current use. The new aerodynamic shape optimization procedure is demonstrated in the design of both two- and three-dimensional inviscid aerodynamic problems including a two-dimensional supersonic internal/external nozzle, two-dimensional transonic airfoils (resulting in supercritical shapes), three-dimensional transport wings, and three-dimensional supersonic delta wings. Each design application results in realistic and useful optimized shapes.

Burgreen, Gregory W.↗

The Thin Oil Film Equation

A thin film of oil on a surface responds primarily to the wall shear stress generated on that surface by a three-dimensional flow. The oil film is also subject to wall pressure gradients, surface tension effects and gravity. The partial differential equation governing the oil film flow is shown to be related to Burgers' equation. Analytical and numerical methods for solving the thin oil film equation are presented. A direct numerical solver is developed where the wall shear stress variation on the surface is known and which solves for the oil film thickness spatial and time variation on the surface. An inverse numerical solver is also developed where the oil film thickness spatial variation over the surface at two discrete times is known and which solves for the wall shear stress variation over the test surface. A One-Time-Level inverse solver is also demonstrated. The inverse numerical solver provides a mathematically rigorous basis for an improved form of a wall shear stress instrument suitable for application to complex three-dimensional flows. To demonstrate the complexity of flows for which these oil film methods are now suitable, extensive examination is accomplished for these analytical and numerical methods as applied to a thin oil film in the vicinity of a three-dimensional saddle of separation.

Brown, James L.↗

Investigation of VEGGIE Root Mat

VEGGIE is a plant growth facility that utilizes the phenomenon of capillary action as its primary watering system. A cloth made of Meta Aramid fiber, known as Nomex is used to wick water up from a reservoir to the bottom of the plants roots. This root mat system is intended to be low maintenance with no moving parts and requires minimal crew interface time. Unfortunately, the water wicking rates are inconsistent throughout the plant life cycle, thus causing plants to die. Over-wicking of water occurs toward the beginning of the cycle, while under-wicking occurs toward the middle. This inconsistency of wicking has become a major issue, drastically inhibiting plant growth. The primary objective is to determine the root cause of the inconsistent wicking through experimental testing. Suspect causes for the capillary water column to break include: a vacuum effect due to a negative pressure gradient in the water reservoir, contamination of material due to minerals in water and back wash from plant fertilizer, induced air bubbles while using syringe refill method, and material limitations of Nomex's ability to absorb and retain water. Experimental testing will be conducted to systematically determine the cause of under and over-wicking. Pressure gages will be used to determine pressure drop during the course of the plant life cycle and during the water refill process. A debubbler device will be connected to a root mat in order to equalize pressure inside the reservoir. Moisture and evaporation tests will simultaneously be implemented to observe moisture content and wicking rates over the course of a plant cycle. Water retention tests will be performed using strips of Nomex to determine materials wicking rates, porosity, and absorptivity. Through these experimental tests, we will have a better understanding of material properties of Nomex, as well as determine the root cause of water column breakage. With consistent test results, a forward plan can be achieved to resolve the issue and give valuable insight for the next generation of VEGGIE.

Subbiah, Arun M.↗

Wind Tunnel Testing of High Efficiency Low Power (HELP) Actuation for Active Flow Control

A High Efficiency Low Power (HELP) methodology has been successfully developed and tested on an active flow control (AFC) version of the 10%-scale high-lift Common Research Model (CRM-HL) at the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel (14x22). The AFC variant of the CRM-HL, designated as CRM-SHL-AFC, was integrated with modular HELP actuator cartridges on the shoulder of its highly deflected (≥50°) simplehinged flaps. A reference conventional CRM-HL configuration equipped with Fowler flaps and a nacelle chine was chosen to provide the targeted lift enhancement goals and for comparison with the results of the CRM-SHL-AFC. The current data are presented with the Transonic Wall Interference Correction System (TWICS) method applied. The HELP actuators, which use a combination of unsteady sweeping jets and steady discrete jets in tandem, were designed to overcome strong adverse pressure gradients, while minimizing the pneumatic power usage. The power coefficient (C(sub π)), which takes account of both supply air pressure and mass flow usage for the AFC actuators, is a useful parameter for judging the actuators’ performance efficiency and perhaps for scaling-up of the AFC system. Full HELP actuation coverage with a constant spanwise nozzle pressure ratio (NPR) was able to achieve the lift enhancement goals for the entire lift curve and was found to be the most effective AFC case. For the most effective case, mass flow rates greater than ~0.91 lbm/s and NPR values greater than ~1.8 (corresponding to C(sub π) = 0.18) are needed to achieve the lift enhancement goal at α = 9°, whereas mass flow rates greater than ~1.23 lbm/s and NPR values greater than ~2.3 (corresponding to C(sub π) = 0.3) are needed to achieve the lift enhancement goal at α = 17°. Surface pressure distributions indicate that the HELP actuation on the flap shoulder increased the suction pressures (and flow circulation) globally in both the streamwise and the spanwise directions, and thereby enhanced the lift over the entire high-lift system. The CRM-SHLAFC configuration equipped with HELP actuation was able to match or exceed the lift performance of the reference conventional CRM-HL, thus meeting the objective of the research.

John C Lin↗

Results of correlations for transition location on a clean-up glove installed on an F-14 aircraft and design studies for a laminar glove for the X-29 aircraft accounting for spanwise pressure gradient

Results of correlative and design studies for transition location, laminar and turbulent boundary-layer parameters, and wake drag for forward swept and aft swept wings are presented. These studies were performed with the use of an improved integral-type boundary-layer and transition-prediction methods. Theoretical predictions were compared with flight measurements at subsonic and transonic flow conditions for the variable aft swept wing F-14 aircraft for which experimental pressure distributions, transition locations, and turbulent boundary-layer velocity profiles were measured. Flight data were available at three spanwise stations for several values of sweep, freestream unit Reynolds number, Mach numbers, and lift coefficients. Theory/experiment correlations indicate excellent agreement for both transition location and turbulent boundary-layer parameters. The results of parametric studies performed during the design of a laminar glove for the forward swept wing X-29 aircraft are also presented. These studies include the effects of a spanwise pressure gradient on transition location and wake drag for several values of freestream Reynolds numbers at a freestream Mach number of 0.9.

Goradia, S. H.↗

Lunar soil properties and soil mechanics. Flow in porous media under rarefied gas conditions. Research phase: Fluid conductivity of lunar surface materials

Theoretical and experimental research on fluid conductivity of lunar surface materials is summarized. Theoretical methods were developed for the analysis of transitional and free-molecular flows, and for analysis of lunar permeability probe data in general. Experimental studies of rarefied flows under conditions of a large pressure gradient show flows in the continuum regime to be responsible for the largest portion of the pressure drop between source and sink for one dimensional flow, provided the entrance Knudsen number is sufficiently small. The concept of local similarity leading to a universal nondimensional function of Knudsen number was shown to have approximate validity; flows in all regimes may be described in terms of an area fraction and a single length parameter. Synthetic porous media prepared from glass beads exhibited flow behavior similar in many regards to that of a natural sandstone; studies using artificial stones with known pore configurations may lead to new insight concerning the structure of natural materials. The experimental method involving the use of segmented specimens of large permeability is shown to be fruitful.

Hurlbut, F. C.↗

Cause and Cure - Deterioration in Accuracy of CFD Simulations With Use of High-Aspect-Ratio Triangular Tetrahedral Grids

Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD re-searchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions and also cause numerical instability. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where triangular/tetrahedral elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identities the reason behind the difficulties in use of such high-aspect ratio triangular/tetrahedral elements is presented here. As will be shown, it turns out that the degree of accuracy deterioration of gradient computation involving a triangular element is hinged on the value of its shape factor Gamma def = sq sin Alpha1 + sq sin Alpha2 + sq sin Alpha3, where Alpha1; Alpha2 and Alpha3 are the internal angles of the element. In fact, it is shown that the degree of accuracy deterioration increases monotonically as the value of Gamma decreases monotonically from its maximal value 9/4 (attained by an equilateral triangle only) to a value much less than 1 (associated with a highly obtuse triangle). By taking advantage of the fact that a high-aspect ratio triangle is not necessarily highly obtuse, and in fact it can have a shape factor whose value is close to the maximal value 9/4, a potential solution to avoid accuracy deterioration of gradient computation associated with a high-aspect ratio triangular grid is given. Also a brief discussion on the extension of the current mathematical framework to the tetrahedral-grid case along with some of the practical results of this extension is also provided. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, the effectiveness of the gradient evaluation procedures within the CESE framework (that have their basis on the analysis presented here) to produce accurate and stable results on such high-aspect ratio meshes is also showcased.

CESE Method↗

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.

Seybert, C. D.↗

Magnetocaloric Pumping of Liquid Oxygen

The field-induced force density on a magnetic fluid is proportional to the magnetic susceptibility times the gradient of the magnetic field squared. The direction of the force is towards increasing magnetic field (positive gradient). Applying a magnetic field to a magnetic fluid will result in a force from all directions towards the location of peak field. Since the magnetic field is conservative and there are no magnetic monopoles, the net field-induced force on any fluid of constant susceptibly will be zero. The only manner to obtain a nonzero net field-induced force is to vary the susceptibility of the fluid. At the gas/liquid interface of liquid oxygen, the susceptibility varies drastically, and the exploitation of the resultant large net forces. An alternative method of varying the magnetic susceptibility is to vary the temperature of the fluid. The magnetic susceptibility of paramagnetic liquid oxygen obeys the Curie-Weiss law: it is inversely proportional to temperature. By applying a temperature gradient in the presence of a symmetric magnetic field, a nonzero net force results. Much of the theory of the so-called Magnetocaloric Effect has previously been developed for and applied to ferromagnetic fluids, or ferrofluids, but is readily applied to paramagnetic liquid oxygen.

Immer, Christopher↗

Electron density profiles from ionograms - Comparisons with rocket profiles.

From a series of ten rocket flights at Wallops Island conducted by the University of Illinois, detailed electron density profiles are available in the altitude range 60-120 km. Concurrent ionograms from a nearby conventional ionosonde of good sensitivity and precision, have been reduced to N(h) profiles by the NOAA methods. Special attention is given to the needs for accurate ionosonde calibration, proper ionogram interpretation, and corrections for underlying ionization, 'valley' ionization, and effects of lateral-gradients. Proper corrections in the present cases yield profiles by ionogram inversion which agree to within a few percent in density, and a few hundred meters or less in height, with the rocket results.

Wright, J. W.↗

Structural optimization by methods of feasible directions.

A general design algorithm based on methods of feasible directions is presented. Zoutendijk's method of feasible directions is first presented as applied to structural design. This method is modified to improve numerical stability of the design process and is then further modified to deal efficiently with infeasible designs. The algorithm requires the analytic gradient of the objective function and the constraint functions which are active at a given stage in the design process. Gradient information is not required for nonactive constraints. Complex constraint functions may be ignored in the initial design stages because violation of these constraints is efficiently overcome later in the design process. The algorithm is demonstrated with elastic design of redundant trusses.

Vanderplaats, G. N.↗

Entry heating and thermal protection

The book provides a documentation of progress in entry heating and thermal protection as related to two very challenging and very dissimilar applications, viz., the Jovian atmospheric probe as a key element of the Galileo mission, and the Space Shuttle. The Jovian atmospheric probe presents what might be the ultimate challenge in terms of the severity of entry. Success of the Galileo probe depends upon a design based on computational analysis largely unverified by experiment. Topics of interest include methods for predicting off-stagnation-point flowfields for planetary entry probes, modeling of radiative heating in base region of Jovian entry probe, gas heating with pressure gradients for the Space Shuttle, external tank laminar and rarefied aeroheating, and a study of the aerothermal entry environment for the Galileo probe.

Olstad, W. B.↗

A vectorized, finite-volume, adaptive grid algorithm applied to planetary entry problems

An adaptive grid, finite-volume method has been applied to problems in planetary entry for computing complete flowfields. The adaption algorithm is implicit in nature and is keyed to resolve user specified gradients. The finite-volume algorithm is explicit, utilizing a maximum time step advancement at each grid point to accelerate convergence to the steady state. The present version of the code is for the laminar flow of a perfect gas. The role of the adaption algorithm in resolving various features of blunt body/wake flow for planetary entry conditions is emphasized.

Gnoffo, P. A.↗

Finite element flow-thermal-structural analysis of aerodynamically heated leading edges

Hypersonic vehicles operate in a hostile aerothermal environment which has a significant impact on their aerothermostructural performance. Significant coupling occurs between the aerodynamic flow field, structural heat transfer, and structural response creating a multidisciplinary interaction. A long term goal of the Aerothermal Loads Branch at the NASA Langley Research Center is to develop a computational capability for integrated fluid, thermal and structural analysis of aerodynamically heated structures. The integrated analysis capability includes the coupling between the fluid and the structure which occurs primarily through the thermal response of the structure, because: (1) the surface temperature affects the external flow by changing the amount of energy absorbed by the structure, and (2) the temperature gradients in the structure result in structural deformations which alter the flow field and attendant surface pressures and heating rates. In the integrated analysis, a finite element method is used to solve: (1) the Navier-Stokes equations for the flow solution, (2) the energy equation of the structure for the temperature response, and (3) the equilibrium equations of the structure for the structural deformation and stresses.

Dechaumphai, Pramote↗

Effect of spatial resolution on apparent sensitivity to initial conditions of a decaying flow as it becomes turbulent

Grids with 32(exp 3), 64(exp 3), and 128(exp 3) points are used in numerical solutions for a decaying flow. The sensitivity of initially neighboring solutions to small changes in initial conditions increases as the spatial resolution improves. A fourth-order finite-difference method is used for the solutions with 32(exp 3) and 64(exp 3) grid points, and a pseudospectral method is used for 128(exp 3) grid points. The latter solutions appear to be rather well-resolved, in spite of the formation of steep velocity gradients in the flow.

Deissler, Robert G.↗