Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “differential 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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

A Preliminary Assessment of LANDSAT-4 Thematic Mapper Data

The geometric errors for image to map rectification of one Thematic Mapper (TM) scene of an area near Windsor, Ontario were studied. The scene had been produced on computer compatible tape by NASA and contained radiometric and system corrections for geometric distortions. The geometric properties of TM photographic imagery permitted very good fitting to map detail using simple scaling techniques to localized areas and, using simple scaling, the overall geometry remained with 500 meters or 0.4 mm at the image scale of 1:1,141,600. An affine transformation, permitting differential scaling, slightly improves the fit to about 400 meters or 0.35 mm at image scale. The imagery shows promise of having the needed additional resolution and spectral discrimination to provide map revision information in urban-rural areas where the MSS sensor is considered inadequate. The late-season prairie image, however, did not hold such promise, and judgement must be reserved until images are acquired at other seasons in this particular geographic area.

Goodenough, D. G.↗
NASA NTRS1986

Analysis of the harmonic gyrotron traveling wave amplifier

A single nonlinear differential equation is derived for describing the spatial evolution of the wave field in a gyrotron amplifier on the cylindrical guide geometry. The equation is then used to determine the efficiencies and the optimized interaction lengths for gyrotron operated with TE sub N1 modes at the N-th cyclotron harmonic. The equation offers a significant advantage in simplifying numerical work in comparison with the brute force particle simulation approach. The efficiencies of the harmonic operations by using a hollow beam decrease quickly with increasing harmonic number, but are enhanced by using an axis-encircling beam instead.

Kuo, S. P.↗
NASA NTRS2002

Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation

The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.

Chi, X.↗
NASA NTRS1975

Stress analysis and stress intensity factors for finite geometry solids containing rectangular surface cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement field in a finite geometry bar containing a variable depth rectangular surface crack under extensionally applied uniform loading. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. Using the obtained displacement field, normal stresses and the stress intensity factor variation along the crack periphery are calculated for different crack depth to bar thickness ratios. Crack opening displacements and stress intensity factors are also obtained for a through-thickness, center cracked bar with variable thickness. The reported results show a considerable potential for using this method in calculating stress intensity factors for commonly encountered surface crack geometries in finite solids.

Gyekenyesi, J. P.↗
NASA NTRS1977

Stress analysis and stress-intensity factors for finite geometry solids containing rectangular surface cracks

The line method of analysis is applied to the Navier-Cauchy equations of elastic equilibrium to calculate the displacement field in a finite geometry bar containing a variable depth rectangular surface crack under extensionally applied uniform loading. The application of this method to these equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. Using the obtained displacement field, normal stresses, and the stress-intensity factor variation along the crack periphery are calculated for different crack depth to bar thickness ratios. Crack opening displacements and stress-intensity factors are also obtained for a through-thickness, center-cracked bar with variable thickness. The reported results show a considerable potential for using this method in calculating stress-intensity factors for commonly encountered surface crack geometries in finite solids

Gyekenyesi, J. P.↗
NASA NTRS1996

Fluid Creep Effects on Near-Wall Solute Transport for Non-Isothermal Ampoules

There is a growing practical and theoretical interest in developing accurate macroscopic modelling for flows arising in chemical or physical vapor transport (VT) crystal growth experiments, including those conducted in reduced gravity environments. Rosner was the first person to point out that previously neglected rarefield gas dynamics phenomena (Stefan and bouyancy-driven flows) become rather important sources of convection. In particular, the combination of rarefaction and strong gradients of temperature (and/or concentration) tangential to the side-walls of the ampoule induces convective flows known as thermal (and concentration) 'creep' respectively. His order-of-magnitude estimates revealed that thermal creep effects can be non-negligible even at normal gravitational levels. On the macroscopic level, the bulk fluid mechanics can be adequately described by the familiar macroscopic equations as long as the boundary conditions are modified to account for the integrated effect of kinetic boundary layers adjacent to solid boundaries. Motivated by the growing importance of these phenomena, we have embarked on a series of computational studies to elucidate these fundamental creep-induced effects for a rarefied gas in simple, two-dimensional confined geometries. However, unlike previous related studies, we resort to a microscopic description of the gas, mathematically expressed by the Boltzmann integro-differential equation. We employ the direct simulation Monte Carlo (DSMC) method of Bird, the theoretical foundations and several practical applications. In the case of thermally induced flows, the no-time counter method of Bird is used, as implemented for a hard-sphere gas. The scheme has been also extended to account for realistic molecular interaction models, an extension necessary if the diffusion physics underlying concentration creep are to be captured.

Papadopoulos, Dimitrios↗
NASA NTRS1980

Gigahertz analog modulation and differential delay of GaAlAs lasers - Temperature and current behavior

Gigahertz analog modulation characteristics of broad-area commercially available GaAlAs lasers have been investigated as a function of temperature and current in the vicinity of the upper frequency limit, where the resonance phenomena occur. The optimum temperature for small-signal amplitude modulation was found to be around -15 C for our particular broad-stripe geometry double-heterostructure laser. The Q was found to increase by a factor of 2 and the bandwidth by about 2%; the external quantum efficiency was maximized in this range. The optimum dc current bias was about 2% above the threshold current. Differential delays have also been measured down to a few picosecond accuracy by a unique phase-angle measurement method using a vector voltmeter. Some of the temperature effects observed may be related to mode changes and multimode and superradiance behavior.

Eng, S. T.↗
NASA NTRS1995

Zero side force volute development

Collector scrolls on high performance centrifugal pumps are currently designed with methods which are based on very approximate flowfield models. Such design practices result in some volute configurations causing excessive side loads even at design flowrates. The purpose of this study was to develop and verify computational design tools which may be used to optimize volute configurations with respect to avoiding excessive loads on the bearings. The new design methodology consisted of a volute grid generation module and a computational fluid dynamics (CFD) module to describe the volute geometry and predict the radial forces for a given flow condition, respectively. Initially, the CFD module was used to predict the impeller and the volute flowfields simultaneously; however, the required computation time was found to be excessive for parametric design studies. A second computational procedure was developed which utilized an analytical impeller flowfield model and an ordinary differential equation to describe the impeller/volute coupling obtained from the literature, Adkins & Brennen (1988). The second procedure resulted in 20 to 30 fold increase in computational speed for an analysis. The volute design analysis was validated by postulating a volute geometry, constructing a volute to this configuration, and measuring the steady radial forces over a range of flow coefficients. Excellent agreement between model predictions and observed pump operation prove the computational impeller/volute pump model to be a valuable design tool. Further applications are recommended to fully establish the benefits of this new methodology.

Anderson, P. G.↗
NASA NTRS1993

On the large eddy simulation of turbulent flows in complex geometry

Application of the method of Large Eddy Simulation (LES) to a turbulent flow consists of three separate steps. First, a filtering operation is performed on the Navier-Stokes equations to remove the small spatial scales. The resulting equations that describe the space time evolution of the 'large eddies' contain the subgrid-scale (sgs) stress tensor that describes the effect of the unresolved small scales on the resolved scales. The second step is the replacement of the sgs stress tensor by some expression involving the large scales - this is the problem of 'subgrid-scale modeling'. The final step is the numerical simulation of the resulting 'closed' equations for the large scale fields on a grid small enough to resolve the smallest of the large eddies, but still much larger than the fine scale structures at the Kolmogorov length. In dividing a turbulent flow field into 'large' and 'small' eddies, one presumes that a cut-off length delta can be sensibly chosen such that all fluctuations on a scale larger than delta are 'large eddies' and the remainder constitute the 'small scale' fluctuations. Typically, delta would be a length scale characterizing the smallest structures of interest in the flow. In an inhomogeneous flow, the 'sensible choice' for delta may vary significantly over the flow domain. For example, in a wall bounded turbulent flow, most statistical averages of interest vary much more rapidly with position near the wall than far away from it. Further, there are dynamically important organized structures near the wall on a scale much smaller than the boundary layer thickness. Therefore, the minimum size of eddies that need to be resolved is smaller near the wall. In general, for the LES of inhomogeneous flows, the width of the filtering kernel delta must be considered to be a function of position. If a filtering operation with a nonuniform filter width is performed on the Navier-Stokes equations, one does not in general get the standard large eddy equations. The complication is caused by the fact that a filtering operation with a nonuniform filter width in general does not commute with the operation of differentiation. This is one of the issues that we have looked at in detail as it is basic to any attempt at applying LES to complex geometry flows. Our principal findings are summarized.

Ghosal, Sandip↗
NASA NTRS1977

Gravitational scattering of zero-rest-mass plane waves

The Feyman-diagram technique is used to calculate the differential cross sections for the scattering of zero-rest-mass plane waves of spin 0, 1, and 2 by linearized Schwarzschild and Kerr geometries in the long-wavelength weak-field limit. It is found that the polarization of right (or left) circularly polarized electromagnetic waves is unaffected by the scattering process (i.e., helicity is conserved) and that the two helicity (polarization) states of the photon are scattered differently by the Kerr geometry. This coupling between the photon helicity and the angular momentum of the scatterer also leads to a partial polarization of unpolarized incident light. For gravitational waves, on the other hand, there is neither helicity conservation nor helicity-dependent scattering; the angular momentum of the scatterer has no polarizing effect on incident unpolarized gravitational waves.

De Logi, W. K.↗
NASA NTRS2015

Improvements to the Tendon-Actuated Lightweight In-Space MANipulator (TALISMAN)

Devices for manipulating and precisely placing payloads are critical for efficient space operations including berthing of spacecraft, in-space assembly, construction and repair. Key to the success of many NASA space activities has been the availability of long-reach crane-like devices such as the Shuttle Remote Manipulation System (SRMS) and the Space Station Remote Manipulation System (SSRMS). These devices have been used for many operations including berthing visiting spacecraft to the International Space Station, deployment of spacecraft, space station assembly, astronaut positioning, payload transfer, and spacecraft inspection prior to atmospheric re-entry. Retiring the Space Transportation System has led to the removal of the SRMS from consideration for in-space missions, thus creating a capability gap. Recognizing this gap, work was initiated at NASA on a new architecture for long-reach space manipulators. Most current devices are constructed by joining revolute joints with carbon composite tubes, with the joints accounting for the majority of the device mass. For example in the case of the SRMS, the entire device mass is 410 kg (904 lbm); the joint structure, motors, gear train, cabling, etc., accounts for the majority of the system mass because the carbon composite tubes mass is 46 kg (101 lbm). An alternate space manipulator concept, the Tendon-Actuated Lightweight In-Space MANipulator (TALISMAN) was created to address deficiencies in the current state-of-the-art in long-reach manipulators. The antagonistic tendon actuated joint architecture allows the motors actuating the joint to be removed from the joint axis, which simplifies the joint design while simultaneously providing mechanical advantage for the motors. The improved mechanical advantage, in turn, reduces the size and power requirements for the motor and gear train. This paper will describe recent architectural improvements to the TALISMAN design that: 1) improve the operational robustness of the system by enabling maneuvers not originally possible by varying the TALISMAN geometry; 2) enable efficient active antagonistic control of a joint while sharing cable between antagonistic tension networks; and 3) uses a unique arrangement of differential capstans to reduce motor torque requirements by an order of magnitude. The paper will also summarize recent efforts to enable autonomous deployment of a TALISMAN including the deployment concept of operations and associated hardware system design. The deployment forces are provided by the same motor systems that are used for articulation, thus reducing the mass associated with the deployment system. The deployment approach is being tested on a TALISMAN prototype which is designed to provide the same operational performance as a shuttle-class manipulator. The prototype has been fabricated and is operational in a new facility at NASA Langley Research Center that has a large area (15.2 m by 21.3 m [50 ft by 70 ft]) air-bearing floor.

Doggett, William R.↗
NASA NTRS1996

Multigrid Methods for EHL Problems

In many bearings and contacts, forces are transmitted through thin continuous fluid films which separate two contacting elements. Objects in contact are normally subjected to friction and wear which can be reduced effectively by using lubricants. If the lubricant film is sufficiently thin to prevent the opposing solids from coming into contact and carries the entire load, then we have hydrodynamic lubrication, where the lubricant film is determined by the motion and geometry of the solids. However, for loaded contacts of low geometrical conformity, such as gears, rolling contact bearings and cams, this is not the case due to high pressures and this is referred to as Elasto-Hydrodynamic Lubrication (EHL) In EHL, elastic deformation of the contacting elements and the increase in fluid viscosity with pressure are very significant and cannot be ignored. Since the deformation results in changing the geometry of the lubricating film, which in turn determines the pressure distribution, an EHL mathematical model must simultaneously satisfy the complex elasticity (integral) and the Reynolds lubrication (differential) equations. The nonlinear and coupled nature of the two equations makes numerical calculations computationally intensive. This is especially true for highly loaded problems found in practice. One novel feature of these problems is that the solution may exhibit sharp pressure spikes in the outlet region. To this date both finite element and finite difference methods have been used to solve EHL problems with perhaps greater emphasis on the use of the finite difference approach. In both cases, a major computational difficulty is ensuring convergence of the nonlinear equations solver to a steady state solution. Two successful methods for achieving this are direct iteration and multigrid methods. Direct iteration methods (e.g Gauss Seidel) have long been used in conjunction with finite difference discretizations on regular meshes. Perhaps one of the best examples of the application of such methods is the recent Effective Influence Method of Dowson and Wang. Multigrid methods have also been used with great success by Venner and Venner and Lubrecht with a good summary being given by Venner. As both these finite difference discretization based approaches appear to provide an efficient way of solving EHL problems, it is important to understand their relative merits. This paper is a first attempt at providing such an understanding in the context of EHL point contact problem, (contact of two spheres), in which the contact zone is a point and an ellipse or circle for unloaded and loaded dry contacts respectively. Since the film thickness and the contact width are generally small compared to the local radius of curvature of the two surfaces, the reduced geometry of the surfaces in the contact area can be accurately approximated to the contact between a paraboloid and a flat surface. The layout of the remainder of this paper is as follows. In section 2 we introduce the form of the equations to be solved. The Effective Influence Newton Method is described in Section 3 while Section 4 describes the Multigrid method to be used. Sections 5 and 6 describe the test problems to be used in the comparison between the two methods and compare the performance of the two methods. Section 7 concludes the paper with an argument of the two methods and suggests some future research directions.

Nurgat, Elyas↗
NASA NTRS2015

Aircraft Measurements of BrO, IO, Glyoxal, NO2, H2O, O2-O2 and Aerosol Extinction Profiles in the Tropics: Comparison with Aircraft-/Ship-Based in Situ and Lidar Measurements

Tropospheric chemistry of halogens and organic carbon over tropical oceans modifies ozone and atmospheric aerosols, yet atmospheric models remain largely untested for lack of vertically resolved measurements of bromine monoxide (BrO), iodine monoxide (IO) and small oxygenated hydrocarbons like glyoxal (CHOCHO) in the tropical troposphere. BrO, IO, glyoxal, nitrogen dioxide (NO2), water vapor (H2O) and O2-O2 collision complexes (O4/ were measured by the University of Colorado Airborne Multi-AXis Differential Optical Absorption Spectroscopy (CU AMAXDOAS) instrument, aerosol extinction by high spectral resolution lidar (HSRL), in situ aerosol size distributions by an ultra high sensitivity aerosol spectrometer (UHSAS) and in situ H2O by vertical-cavity surface-emitting laser (VCSEL) hygrometer. Data are presented from two research flights (RF12, RF17) aboard the National Science Foundation/ National Center for Atmospheric Research Gulfstream V aircraft over the tropical Eastern Pacific Ocean (tEPO) as part of the "Tropical Ocean tRoposphere Exchange of Reactive halogens and Oxygenated hydrocarbons" (TORERO) project (January/February 2012). We assess the accuracy of O4 slant column density (SCD) measurements in the presence and absence of aerosols. Our O4-inferred aerosol extinction profiles at 477 nm agree within 6% with HSRL in the boundary layer and closely resemble the renormalized profile shape of Mie calculations constrained by UHSAS at low (sub-Rayleigh) aerosol extinction in the free troposphere. CU AMAX-DOAS provides a flexible choice of geometry, which we exploit to minimize the SCD in the reference spectrum (SCDREF, maximize signal-to-noise ratio) and to test the robustness of BrO, IO and glyoxal differential SCDs. The RF12 case study was conducted in pristine marine and free tropospheric air. The RF17 case study was conducted above the NOAA RV Ka'imimoana (TORERO cruise, KA-12-01) and provides independent validation data from ship-based in situ cavity-enhanced DOAS and MAX-DOAS. Inside the marine boundary layer (MBL) no BrO was detected (smaller than 0.5 pptv), and 0.2-0.55 pptv IO and 32-36 pptv glyoxal were observed. The near-surface concentrations agree within 30% (IO) and 10% (glyoxal) between ship and air-craft. The BrO concentration strongly increased with altitude to 3.0 pptv at 14.5 km (RF12, 9.1 to 8.6 deg N; 101.2 to 97.4 deg W). At 14.5 km, 5-10 pptv NO2 agree with model predictions and demonstrate good control over separating tropospheric from stratospheric absorbers (NO2 and BrO). Our profile retrievals have 12-20 degrees of freedom (DoF) and up to 500m vertical resolution. The tropospheric BrO vertical column density (VCD) was 1.5 x 10(exp 13) molec cm(exp -2) (RF12) and at least 0.5 x 10(exp 13) molec cm(exp -2) (RF17, 0- 10 km, lower limit). Tropospheric IO VCDs correspond to 2.1 x 10(exp 12) molec cm(exp -2) (RF12) and 2.5 x 10(exp 12) molec cm(exp -2) (RF17) and glyoxal VCDs of 2.6 x 10(exp 14) molec cm(exp -2) (RF12) and 2.7 x 10(exp 14) molec cm(exp -2) (RF17). Surprisingly, essentially all BrO as well as the dominant IO and glyoxal VCD fraction was located above 2 km (IO: 58 plus or minus 5 %, 0.1-0.2 pptv; glyoxal: 52 plus or minus 5 %, 3-20 pptv). To our knowledge there are no previous vertically resolved measurements of BrO and glyoxal from aircraft in the tropical free troposphere. The atmospheric implications are briefly discussed. Future studies are necessary to better understand the sources and impacts of free tropospheric halogens and oxygenated hydrocarbons on tropospheric ozone, aerosols, mercury oxidation and the oxidation capacity of the atmosphere.

Volkamer, R.↗
NASA NTRS1980

Numerical study of baroclinic instability

The effect of a power law gravity field on baroclinic instability is examined with emphasis on the case of inverse fifth power gravity, since this is the power law produced when terrestrial gravity is simulated in spherical geometry by electrostatic means. Growth rates of unstable normal modes were obtained as a function of parameters of the problem by solving a second order differential equation numerically. Results are compared with those from an earlier study where gravity was a constant. The conclusion is that, over the range of parameter space explored here, there is no significant change in the character of theoretical regime diagrams if the vertically averaged gravity is used as a parameter.

Geisler, J. E.↗
NASA NTRS1986

Numerical prediction of cold turbulent flow in combustor configurations with different centerbody flame holders

A numerical study of cold turbulent flow in a combustor containing a centerbody flame holder is presented. The axisymmetric Navier-Stokes equations incorporating a k-epsilon turbulence model were solved in a nonorthogonal curvilinear coordinate system. The finite volume method, applied to a staggered grid system, was used to discretize the differential equations. These finite differenced equations were then solved iteratively utilizing the SIMPLE algorithm described by Patankar (1980). Solutions were obtained for two combustor duct geometries with various centerbody flame holders which included a disk, a cone, and a sphere. The extent of mixing due to these bodies was evaluated. A comparison with previously obtained experimental data yielded moderately good agreement.

Yung, C. N.↗
NASA NTRS1990

Experimental feasibility of investigating acoustic waves in Couette flow with entropy and pressure gradients

The feasibility is discussed for an experimental program for studying the behavior of acoustic wave propagation in the presence of strong gradients of pressure, temperature, and flow. Theory suggests that gradients effects can be experimentally observed as resonant frequency shifts and mode shape changes in a waveguide. A convenient experimental geometry for such experiments is the annular region between two co-rotating cylinders. Radial temperature gradients in a spinning annulus can be generated by differentially heating the two cylinders via electromagnetic induction. Radial pressure gradients can be controlled by varying the cylinder spin rates. Present technology appears adequate to construct an apparatus to allow independent control of temperature and pressure gradients. A complicating feature of a more advanced experiment, involving flow gradients, is the requirement for independently controlled cylinder spin rates. Also, the boundary condition at annulus terminations must be such that flow gradients are minimally disturbed. The design and construction of an advanced apparatus to include flow gradients will require additional technology development.

Parrott, Tony L.↗
NASA NTRS1998

Surface Micromachined Silicon Carbide Accelerometers for Gas Turbine Applications

A finite-element analysis of possible silicon carbide (SIC) folded-beam, lateral-resonating accelerometers is presented. Results include stiffness coefficients, acceleration sensitivities, resonant frequency versus temperature, and proof-mass displacements due to centripetal acceleration of a blade-mounted sensor. The surface micromachined devices, which are similar to the Analog Devices Inc., (Norwood, MA) air-bag crash detector, are etched from 2-pm thick, 3C-SiC films grown at 1600 K using atmospheric pressure chemical vapor deposition (APCVD). The substrate is a 500 gm-thick, (100) silicon wafer. Polysilicon or silicon dioxide is used as a sacrificial layer. The finite element analysis includes temperature-dependent properties, shape change due to volume expansion, and thermal stress caused by differential thermal expansion of the materials. The finite-element results are compared to experimental results for a SiC device of similar, but not identical, geometry. Along with changes in mechanical design, blade-mounted sensors would require on-chip circuitry to cancel displacements due to centripetal acceleration and improve sensitivity and bandwidth. These findings may result in better accelerometer designs for this application.

DeAnna, Russell G.↗
NASA NTRS1976

Pseudopotentials of Estabrook and Wahlquist, the geometry of solitons, and the theory of connections

The prolongation structure of Wahlquist and Estabrook is interpreted as a connection. In this way, some geometric insight might be provided for the description of those nonlinear partial differential equations which admit soliton solutions. A new geometric property - linked to the existence of an SL(2,R) connection - is proved for the solutions of the Korteweg-de Vries equation.

Hermann, R.↗