Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable geometry”

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 793 records · Page 44

Two-and-a-half-dimensional magnetohydrodynamic turbulence

The homogeneous turbulence for which fluctuating magnetic fields and velocity fields are independent of one spatial coordinate but still possess all three components are studied to form a generalized two-and-half dimensional geometry. The integral of the z-vector potential and the magnetic helicity are shown to be ideal invariants, and the basic dynamical variables and the equations of uniform-density incompressible magnetohydrodynamics are defined. The possibility of simultaneous inverse cascades is considered, and Kolmogoroff dimensional analysis is employed to infer omnidirectional inverse cascade spectra. Implications of a selective decay hypothesis, where a tendency exists in the initial value problem for the ideal invariants directly cascaded to higher wavenumbers to be selectively dissipated, and the possibility of a dynamo action in the considered geometry, are examined.

Montgomery, D.↗

A review of tidal energy—Resource, feedbacks, and environmental interactions

The ocean contains a variety of renewable energy resources, little of which has been exploited. Here, we review both tidal range and tidal stream energy, with a focus on the resource, feedbacks, and environmental interactions. The review covers a wide range of timescales of relevance to tidal energy, from fortnightly (spring-neap) and semi-diurnal variability, down to array, and device-scale turbulence. When simulating the regional tidal energy resource, and to assess environmental impacts, it is necessary to account for feedbacks between the tidal array and the resource itself. We critically review various methods for simulating energy extraction, from insights gained through theoretical studies of “tidal fences” in idealized channels, to realistic three-dimensional model studies with complex geometry and arrays of turbines represented by momentum sinks and additional turbulence due to the presence of rotors and support structures. We discuss how variability can be reduced by developing multiple (aggregated) sites with a consideration of the enhanced phase diversity offered by exploiting less energetic tidal currents. This leads to future research questions that have not yet been explored in depth at first-generation tidal sites in relatively sheltered channels (e.g., the interaction of waves with currents). Such enhanced understanding of real sea conditions, including the effects of wind and waves, leads to our other identified primary future research direction—reduced uncertainties in turbulence predictions, including the development of realistic models that simulate the interaction between ambient turbulence and the turbulence resulting from multiple wakes, and changes to system-wide hydrodynamics, water quality, and sedimentation.

Environmental impacts↗

Seasonal and Dayurnal Planetary Albedo Variability from Six Years of DSCOVR EPIC Data

Deep Space Climate Observatory (DSCOVR) measurements of Earth’s reflected solar radiation from the Lissajous orbital position near the Lagrangian L1 point provide continuous monitoring of the Earth’s sunlit hemisphere. CERES-based angle models were used to convert the near-hourly reflected radiances of the EPIC images into a climate-style planetary albedo data-point over the sunlit hemisphere. Integration over the sunlit hemispheres averages out the meteorological weather noise, but retains the seasonal and planetary-scale variability. As the Earth rotates, this generates variations in the Earth’s planetary albedo that are precisely aligned in longitude, which constitutes the dayurnal cycle. This dayurnal variability in the planetary albedo arises from planetary-scale changes in cloud radiative properties that can be directly compared to similarly sampled climate GCM output data. Six years of EPIC data have been analyzed, showing characteristic patterns in the seasonal and dayurnal variability of the Earth’s planetary albedo. Much of the seasonal change in planetary albedo is associated with the changing DSCOVR viewing geometry and the change in solar declination. But throughout the year, the highest planetary albedos are observed over the Central Asia (Iraq) longitude, while the lowest planetary albedos occur over the Central Pacific longitude. For these longitudes, the relative seasonal changes in the planetary albedo are slowly varying and anti-correlated. Dayurnal amplitude maxima tend to occur during the July-September timeframe, with April-May and December exhibit distinct minima in the dayurnal amplitude. West Africa and the West Pacific longitudes likewise exhibit anti-correlated seasonal variability, while they also undergo anti-correlated short period oscillations. On the other hand, in nearby longitudes, there are short period spikes in planetary albedo of a few-days duration, as well as longer period oscillations that may range from a week to several months, that generally tend to be correlated.

Planetary Albedo↗

Impact of lead on an axisymmetric, single bead blown powder DED overhung geometry

Metal additive manufacturing can be utilized for the near net-shape manufacture of components in a layer wise technique. Traditionally, the buildup process is conducted in the direction of gravity or bottom to top in the vertical orientation; however, with the availability of commercial systems with additive heads or fixturing that can index from vertical, and the need to manufacture larger parts without modifying available systems, the limitation of manufacturing from bottom to top is removed enabling the buildup of larger components printed at an angle from vertical. The effect of gravity on the melt pool and as-deposited component quality when printing off vertical is unknown in literature, especially for axisymmetric components. This understanding is critical for the advancement of manufacturing components of increased size in 4 and 5-axis additive systems. This investigation utilizes blown powder directed energy deposition to evaluate the change in as-printed geometry when the start point is altered in relation to gravity while the part rotates to manufacture an axisymmetric component. The objectives of this work are to determine the impact of the deposition location on the geometric variability on axisymmetric components. This investigation tests the hypothesis that differences in layer height exist due to a change in catchment efficiency when the deposition location is moved to a tangent of the round geometry due to a change in the melt pool dynamics. It was found the ideal deposition location when printing at 26.6-degrees from vertical was at −90-degrees from the top center point of the component, along the tangent, where gravity was pushing the melt pool down, but the rotation of the part was pulling the deposited material towards the top center of the component. This work provides an understanding of layer height stability and catchment efficiency to guide print orientation strategy for high-aspect ratio components. It was found the −90-degree lead deposition had the best layer height stability at 0.6 % as compared to the top-center at 6.3 % and +90-degree deposition location at 9.2 % for the nominal programmed layer height. The change in layer height also effected the final diameter the greatest for the top center deposition location where the diameter diverged by 2.5 % during the overbuild condition and converged by 0.6 %. In conclusion, this finding will increase the manufacturing efficiency of axisymmetric components by increasing the passive stability of the printing process for the successful manufacture of parts without the need for perfectly calibrated manufacturing parameters.

36 MATERIALS SCIENCE↗

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗

Melting rate correlation with batch properties and melter operating conditions during conversion of nuclear waste melter feeds to glasses

The rate of conversion of nuclear waste melter feed to glass is affected by the selection of melter feed materials and by melter design and operation. The melting rate correlation (MRC) is an equation that relates the glass production rate with two types of variables: (1) feed and melt properties: conversion heat, cold-cap bottom temperature, and glass melt viscosity; and (2) melter design and operation parameters: melter geometry, melter operating temperature, and gas bubbling rate. The MRC shows good agreement for an extended melting-rate data set of high-level waste (HLW) melter feeds and a data set generated for low-activity waste (LAW) melter feeds. Laboratory observation of heated melter feed samples is often used to assess the cold-cap bottom temperature of HLW melter feeds (moderately foaming feeds), but this technique appears inadequate for LAW melter feeds (vigorously foaming feeds). For LAW feeds, an adequate assessment of the cold-cap bottom temperature was achieved using evolved gas analysis, which allows identification of the collapse of primary foam for oxidized feeds. This assessment shows that the cold-cap bottom temperature for vigorously foaming LAW feeds is higher than that for moderately foaming HLW feeds. When the results of MRC are compared, LAW feeds are generally less sensitive to the bubbling rate and melt viscosity, and more sensitive to the cold-cap bottom temperature than HLW feeds. The MRC qualifies as a promising tool to support the selection of melter feed materials and melter operating conditions, which is determined from expensive independent scaled melter experiments, and sophisticated mathematical models.

Lee, Seung Min↗

Comparison of Two Algebraic Methods for Curve/curve Intersection

Most geometric modeling systems use either polynomial or rational functions to represent geometry. In such systems most computational problems can be formulated as systems of polynomials in one or more variables. Classical elimination theory can be used to solve such systems. Here Cayley's method of elimination is summarized and it is shown how it can best be used to solve the curve/curve intersection problem. Cayley's method was found to be a more straightforward approach. Furthermore, it is computationally simpler, since the elements of the Cayley matrix are one variable instead of two variable polynomials. Researchers implemented and tested both methods and found Cayley's to be more efficient. Six pairs of curves, representing mixtures of lines, circles, and cubic arcs were used. Several examples had multiple intersection points. For all six cases Cayley's required less CPU time than the other method. The average time ratio of method 1 to method 2 was 3.13:1, the least difference was 2.33:1, and the most dramatic was 6.25:1. Both of the above methods can be extended to solve the surface/surface intersection problem.

Demontaudouin, Y.↗

A multistage mesh generator for solving the average-passage equation system

One means of numerically simulating the 3-D flow field within a multistage turbomachine is through the solution of the average-passage equation system. One requirement of a current algorithm used to solve this system of equations has been the ability to generate multiple blade row meshes which satisfy specific geometrical constraints. In addition to meeting this criterion, one desires a mesh generation code which requires minimal user input, utilizes variable mesh control parameters, generates diagnostics helpful to the user, and possesses the capability to handle widely varying geometries. A mesh generation code with these features was written and has been used in solving the inviscid form of the average-passing equation system for both ducted and unducted multiple blade row geometries. This paper serves as a user reference guide, with a description of the mesh generation algorithm, a sample input file, and examples of typical meshes generated.

Mulac, Richard A.↗

Two-scale model of a geomagnetic field variation

The effect of the vertical scale is investigated by considering a simple kinematic two-scale model of fluid flow inducing a variable magnetic field. Depending on the time constant, the induced magnetic field displays a variety of behaviors and geometries. In the high-frequency case, for example, a strong magnetic field tangential to the core mantle boundary, and hidden in the Delta layer, can be generated. A detailed computation and description of this magnetic field are presented. Some possible features of the secular variation of the actual geomagnetic field are discussed in the light of the model proposed here.

Braginsky, S. I.↗

Multi-grid domain decomposition approach for solution of Navier-Stokes equations in primitive variable form

The new multigrid (or adaptive) pseudospectral element method was carried out for the solution of incompressible flow in terms of primitive variable formulation. The desired features of the proposed method include the following: (1) the ability to treat complex geometry; (2) high resolution adapted in the interesting areas; (3) requires minimal working space; and (4) effective in a multiprocessing environment. The approach for flow problems, complex geometry or not, is to first divide the computational domain into a number of fine-grid and coarse-grid subdomains with the inter-overlapping area. Next, it is necessary to implement the Schwarz alternating procedure (SAP) to exchange the data among subdomains, where the coarse-grid correction is used to remove the high frequency error that occurs when the data interpolation from the fine-grid subdomain to the coarse-grid subdomain is conducted. The strategy behind the coarse-grid correction is to adopt the operator of the divergence of the velocity field, which intrinsically links the pressure equation, into this process. The solution of each subdomain can be efficiently solved by the direct (or iterative) eigenfunction expansion technique with the least storage requirement, i.e. O(N(exp 3)) in 3-D and O(N(exp 2)) in 2-D. Numerical results of both driven cavity and jet flow will be presented in the paper to account for the versatility of the proposed method.

Ku, Hwar-Ching↗

Independent pixel and Monte Carlo estimates of stratocumulus albedo

Monte Carlo radiative transfer methods are employed here to estimate the plane-parallel albedo bias for marine stratocumulus clouds. This is the bias in estimates of the mesoscale-average albedo, which arises from the assumption that cloud liquid water is uniformly distributed. The authors compare such estimates with those based on a more realistic distribution generated from a fractal model of marine stratocumulus clouds belonging to the class of 'bounded cascade' models. In this model the cloud top and base are fixed, so that all variations in cloud shape are ignored. The model generates random variations in liquid water along a single horizontal direction, forming fractal cloud streets while conserving the total liquid water in the cloud field. The model reproduces the mean, variance, and skewness of the vertically integrated cloud liquid water, as well as its observed wavenumber spectrum, which is approximately a power law. The Monte Carlo method keeps track of the three-dimensional paths solar photons take through the cloud field, using a vectorized implementation of a direct technique. The simplifications in the cloud field studied here allow the computations to be accelerated. The Monte Carlo results are compared to those of the independent pixel approximation, which neglects net horizontal photon transport. Differences between the Monte Carlo and independent pixel estimates of the mesoscale-average albedo are on the order of 1% for conservative scattering, while the plane-parallel bias itself is an order of magnitude larger. As cloud absorption increases, the independent pixel approximation agrees even more closely with the Monte Carlo estimates. This result holds for a wide range of sun angles and aspect ratios. Thus, horizontal photon transport can be safely neglected in estimates of the area-average flux for such cloud models. This result relies on the rapid falloff of the wavenumber spectrum of stratocumulus, which ensures that the smaller-scale variability, where the radiative transfer is more three-dimensional, contributes less to the plane-parallel albedo bias than the larger scales, which are more variable. The lack of significant three-dimensional effects also relies on the assumption of a relatively simple geometry. Even with these assumptions, the independent pixel approximation is accurate only for fluxes averaged over large horizontal areas, many photon mean free paths in diameter, and not for local radiance values, which depend strongly on the interaction between neighboring cloud elements.

Cahalan, Robert F.↗

Non-Coalescence Effects in Microgravity

Non-coalescence of two bodies of the same liquid and the suppression of contact between liquid drops and solid surfaces is being studied through a pair of parallel investigations being conducted at the Georgia Institute of Technology and the Microgravity Research and Support (MARS) Center in Naples, Italy. Both non-coalescence and contact suppression are achieved by exploiting the mechanism of thermocapillary convection to drive a lubricating film of surrounding gas (air) into the space between the two liquid free surfaces (non-coalescence) or between the drop free surface and the solid (contact suppression). Earlier experiments performed included flow-visualization experiments in both axisymmetric and (nearly) two-dimensional geometries and quantitative measurements of film thickness in the contact-suppression case in both geometries. Work done in the second year has focused on obtaining quantitative results relating to the effects of variable air pressure, development of analytical and numerical models of non-coalescing droplets and to pursuing potential applications of these self-lubricated systems.

Neitzel, G. Paul↗

The Oxidation Kinetics of Continuous Carbon Fibers in a Cracked Ceramic Matrix Composite

Experimental observations and results suggest two primary regimes as a function of temperature, i.e., diffusion and reaction controlled kinetics. Thermogravimetric analysis of carbon fiber in flowing oxygen gave an activation energy of 64.1 kJ/mol in the temperature range of 500 to 600 C and an apparent activation energy of 7.6 kJ/mol for temperatures from 600 to 1400 C. When C/SiC composite material was unstressed, matrix effects at temperatures from 900 to 1400 C protected the internal fibers. When under stress, self-protection was not observed. Increasing the stress from 10 to 25 ksi caused a 67 to 82 percent reduction in times to failure at temperatures from 750 to 1500 C. Based on experimental results, observation, and theory, a finite difference model was developed, which simulates the diffusion of oxygen into a matrix crack that is bridged by carbon fibers. The model allows the influence of important variables on oxidation kinetics to be studied systematically, i.e., temperature, reaction rate constant, diffusion coefficient, environment, and sample geometry.

Halbig, Michael C.↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called 'structured' (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called "structured" (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗

A Multidisciplinary Tool for Systems Analysis of Planetary Entry, Descent, and Landing (SAPE)

SAPE is a Python-based multidisciplinary analysis tool for systems analysis of planetary entry, descent, and landing (EDL) for Venus, Earth, Mars, Jupiter, Saturn, Uranus, Neptune, and Titan. The purpose of SAPE is to provide a variable-fidelity capability for conceptual and preliminary analysis within the same framework. SAPE includes the following analysis modules: geometry, trajectory, aerodynamics, aerothermal, thermal protection system, and structural sizing. SAPE uses the Python language-a platform-independent open-source software for integration and for the user interface. The development has relied heavily on the object-oriented programming capabilities that are available in Python. Modules are provided to interface with commercial and government off-the-shelf software components (e.g., thermal protection systems and finite-element analysis). SAPE runs on Microsoft Windows and Apple Mac OS X and has been partially tested on Linux.

Samareh, Jamshid A.↗

Software for Refining or Coarsening Computational Grids

A computer program performs calculations for refinement or coarsening of computational grids of the type called structured (signifying that they are geometrically regular and/or are specified by relatively simple algebraic expressions). This program is designed to facilitate analysis of the numerical effects of changing structured grids utilized in computational fluid dynamics (CFD) software. Unlike prior grid-refinement and -coarsening programs, this program is not limited to doubling or halving: the user can specify any refinement or coarsening ratio, which can have a noninteger value. In addition to this ratio, the program accepts, as input, a grid file and the associated restart file, which is basically a file containing the most recent iteration of flow-field variables computed on the grid. The program then refines or coarsens the grid as specified, while maintaining the geometry and the stretching characteristics of the original grid. The program can interpolate from the input restart file to create a restart file for the refined or coarsened grid. The program provides a graphical user interface that facilitates the entry of input data for the grid-generation and restart-interpolation routines.

Daines, Russell↗

In-Situ Nondestructive Evaluation of Kevlar(Registered Trademark)and Carbon Fiber Reinforced Composite Micromechanics for Improved Composite Overwrapped Pressure Vessel Health Monitoring

NASA has been faced with recertification and life extension issues for epoxy-impregnated Kevlar 49 (K/Ep) and carbon (C/Ep) composite overwrapped pressure vessels (COPVs) used in various systems on the Space Shuttle and International Space Station, respectively. Each COPV has varying criticality, damage and repair histories, time at pressure, and pressure cycles. COPVs are of particular concern due to the insidious and catastrophic burst-before-leak failure mode caused by stress rupture (SR) of the composite overwrap. SR life has been defined [1] as the minimum time during which the composite maintains structural integrity considering the combined effects of stress level(s), time at stress level(s), and associated environment. SR has none of the features of predictability associated with metal pressure vessels, such as crack geometry, growth rate and size, or other features that lend themselves to nondestructive evaluation (NDE). In essence, the variability or surprise factor associated with SR cannot be eliminated. C/Ep COPVs are also susceptible to impact damage that can lead to reduced burst pressure even when the amount of damage to the COPV is below the visual detection threshold [2], thus necessitating implementation of a mechanical damage control plan [1]. Last, COPVs can also fail prematurely due to material or design noncompliance. In each case (SR, impact or noncompliance), out-of-family behavior is expected leading to a higher probability of failure at a given stress, hence, greater uncertainty in performance. For these reasons, NASA has been actively engaged in research to develop NDE methods that can be used during post-manufacture qualification, in-service inspection, and in-situ structural health monitoring. Acoustic emission (AE) is one of the more promising NDE techniques for detecting and monitoring, in real-time, the strain energy release and corresponding stress-wave propagation produced by actively growing flaws and defects in composite materials [3,4,5,6,7,8]. To gain further insight into the mechanisms responsible for composite rupture, broadband modal acoustic emission analysis was used. Also, since AE data reduction proved to be very time consuming, specialized data reduction software was written to automate the process.

Waller, Jess↗