Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hyperbolic”

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,117 records · Page 62

A Fluorescent Aerogel for Capture and Identification of Interplanetary and Interstellar Dust

Aerogels are extremely low-density solids whose superiority as capturing media for hypervelocity (v > 0.5 km/s) grains has been well established. A prominent example is the use of silica aerogel as the collecting medium for cometary and interstellar grains on NASA s Stardust mission. Aerogel collectors have been deployed in low-earth orbit, but severe background from anthropogenic orbital debris has so far prevented the identification of more than a handful of interplanetary particles. No interstellar particles have been identified so far. Since they are on hyperbolic orbits, extraterrestrial particles are faster than orbital debris, so could in principle be identified on that basis, but existing aerogels give little information on impact velocity. With this in mind, we have developed a novel calorimetric aerogel which passively records the kinetic energy of captured hypervelocity particles.

Dominguez, Gerardo↗

Mission Steering Profiles of Outer Planetary Orbiters Using Radioisotope Electric Propulsion

Radioisotope Electric Propulsion (REP) has the potential to enable small spacecraft to orbit outer planetary targets with trip times comparable to flyby missions. The ability to transition from a flyby to an orbiter mission lies in the availability of continuous low power electric propulsion along the entire trajectory. The electric propulsion system s role is to add and remove energy from the spacecraft s trajectory to bring it in and out of a heliocentric hyperbolic escape trajectory for the outermost target bodies. Energy is added and the trajectory is reshaped to rendezvous with the closer-in target bodies. Sample REP trajectories will be presented for missions ranging for distances from Jupiter orbit to the Pluto-Kuiper Belt.

Fiehler, Douglas↗

Macroscopic Models of Radiative Transfer as Applied to Computation of the Radiation Field in the Solar Atmosphere

In this paper a new three dimensional half-moment model for radiative transfer is presented for a gray medium. It describes the evolution of the zeroth and first directional half moments of the radiative intensity. The closure is provided, similarly to Dubroca and Klar, by the maximum entropy concept. This work generalizes that model to three dimensions. The model presented here (the derivation being done in Ripoll and Wray, called the M(sup 1/2)(sub 1) model, is a hyperbolic system consisting of a total of eight equations in three dimensions, four equations for each direction. Each half model has the classical form of a macroscopic moment model in which the pressure tensor is constructed from the well-known Eddington tensor with a particular Eddington factor. Moreover, different source and border terms occur. The latter introduce couplings between the macroscopic and microscopic quantities and between the + and - streams, through the intensity in the plane perpendicular to the flux. The main theoretical application of the half moment model, treated in this paper, is its reduction to a full moment model, called M(sup +)(sub 1), for the particular but important case of a hot, opaque source radiating in a cold transparent (or semi-transparent) medium for very specific applications, such as stellar interiors or atmospheres, or combustion problems. The structure of the paper is as follows. In section 2, the model M(sup 1/2)(sub 1) is presented. In section 3, for the particular case of a hot, opaque source radiating into a cold medium, the half moment model is reduced to the M(sup +)(sub 1) model. In section 4, we first solve a simple and academic problem to validate the models, followed by a simplified solar atmosphere.

Ripoll, J.-F.↗

Combustion Of Moving Droplets And Of Droplets Suspended Within A Convective Environment: Transient Numerical Results

The problem considered is that of a single-component liquid fuel (n-heptane) droplet undergoing evaporation and combustion in a hot, convective, low pressure, zero-gravity environment of infinite expanse. For a moving droplet, the relative velocity (U(sub infinity)) between the droplet and freestream is subject to change due to the influence of the drag force on the droplet. For a suspended droplet, the relative velocity is kept constant. The governing equations for the gas-phase and the liquid-phase consist of the unsteady, axisymmetric equations of mass, momentum, species (gas-phase only) and energy conservation. Interfacial conservation equations are employed to couple the two phases. Variable properties are used in the gas- and liquid-phase. Multicomponent diffusion in the gas-phase is accounted for by solving the Stefan-Maxwell equations for the species diffusion velocities. A one-step overall reaction is used to model the combustion. The governing equations are discretized using the finite volume and SIMPLEC methods. A colocated grid is adopted. Hyperbolic tangent stretching functions are used to concentrate grid points near the fore and aft lines of symmetry and at the droplet surface in both the gas- and liquid-phase. The discretization equations are solved using the ADI method with the TDMA used on each line of the two alternating directions. Iterations are performed within each time-step until convergence is achieved. The grid spacing, size of the computational domain and time-step were tested to ensure that all solutions are independent of these parameters. A detailed discussion of the numerical model is given.

Gogos, George↗

An Investigation of Wave Propagations in Discontinuous Galerkin Method

Analysis of the discontinuous Galerkin method has been carried out for one- and two-dimensional system of hyperbolic equations. Analytical, as well as numerical, properties of wave propagation in a DGM scheme are derived and verified with direct numerical simulations. In addition to a systematic examination of the dissipation and dispersion errors, behaviours of a DG scheme at an interface of two different grid topologies are also studied. Under the same framework, a quantitative discrete analysis of various artificial boundary conditions is also conducted. Progress has been made in numerical boundary condition treatment that is closely related to the application of DGM in aeroacoustics problems. Finally, Fourier analysis of DGM for the Convective diffusion equation has also be studied in connection with the application of DG schemes for the Navier-Stokes equations. This research has resulted in five(5) publications, plus one additional manuscript in preparation, four(4) conference presentations, and three(3) departmental seminars, as summarized in part II. Abstracts of papers are given in part 111 of this report.

Hu, Fang Q.↗

Variable Step Integration Coupled with the Method of Characteristics Solution for Water-Hammer Analysis, A Case Study

One-dimensional water-hammer modeling involves the solution of two coupled non-linear hyperbolic partial differential equations (PDEs). These equations result from applying the principles of conservation of mass and momentum to flow through a pipe, and usually the assumption that the speed at which pressure waves propagate through the pipe is constant. In order to solve these equations for the interested quantities (i.e. pressures and flow rates), they must first be converted to a system of ordinary differential equations (ODEs) by either approximating the spatial derivative terms with numerical techniques or using the Method of Characteristics (MOC). The MOC approach is ideal in that no numerical approximation errors are introduced in converting the original system of PDEs into an equivalent system of ODEs. Unfortunately this resulting system of ODEs is bound by a time step constraint so that when integrating the equations the solution can only be obtained at fixed time intervals. If the fluid system to be modeled also contains dynamic components (i.e. components that are best modeled by a system of ODEs), it may be necessary to take extremely small time steps during certain points of the model simulation in order to achieve stability and/or accuracy in the solution. Coupled together, the fixed time step constraint invoked by the MOC, and the occasional need for extremely small time steps in order to obtain stability and/or accuracy, can greatly increase simulation run times. As one solution to this problem, a method for combining variable step integration (VSI) algorithms with the MOC was developed for modeling water-hammer in systems with highly dynamic components. A case study is presented in which reverse flow through a dual-flapper check valve introduces a water-hammer event. The predicted pressure responses upstream of the check-valve are compared with test data.

Turpin, Jason B.↗

An assay for intermolecular exchange of alpha crystallin

An affinity column of alpha crystallin linked to cyanogen bromide-activated Sepharose was developed to study the exchange of alpha subunits. Alpha crystallin bound to the Sepharose-alpha complex was dissociated with 8 mol/l urea, followed by quantitation using high-performance reverse-phase liquid chromatography. The time course of binding at 37 degrees C showed a hyperbolic binding pattern reaching equilibrium between 6-18 hr. Under these conditions, binding of beta and gamma crystallins to the same matrix was less than 10% of the alpha values, as was binding of alpha to glycine-coupled Sepharose. This assay was used to demonstrate changes in the subunit exchange of alpha crystallins present in high molecular weight versus lower molecular weight aggregates of the human lens. These results show that this binding procedure was a specific reproducible assay that might be used to study intermolecular interactions of the alpha crystallins.

NASA Discipline Cell Biology↗

Pressure Distribution Measurements at High Speed and Oblique Incidence of Flow

The present report contains the results of a series of observations obtained for a wing of symmetrical profile for different angles of yaw. The shock tunnel with 0.4m x 0.4m cross section, of the Institute for Gas Dynamics at L.F.A. - Braunchweig was available for this work. The profile used was a 9 percent, thick hyperbolic profile with maximum thickness 40 percent aft, which was calculated by the method of F. Ringleb (2). Since the conformal transformation of this profile is known, the theoretical pressure distribution could be determined exactly for the case of incompressible flow. Then by the use of the Prandtl rule one may extend the comparison of theory and experiment to the case of higher velocity of incident flow.

Lippisch, A.↗

Methods of Stress Calculation in Rotating Disks

The paper describes nethods of computing the stresses in disks of a given profile as well as methods of choosing the disk profiles for a given stress distribution for turhines, turbo blowers, and so forth. A new method of in tegrating the differential equations of Stodola leads to a simplification of the computation for disks of hyperbolic profile.

Tumarkin, S.↗

Linear Stability Analysis of Gravitational Effects on a Low-Density Gas Jet Injected into a High-Density Medium

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The Briggs-Bers criterion was combined with the spatio-temporal stability analysis to determine the nature of the absolute instability of the jet whether absolutely or convectively unstable. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the absolute instability of the jet were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be.

Lawson, Anthony L.↗

Effects of Mean Flow Profiles on Instability of a Low-Density Gas Jet Injected into a High-Density Gas

The objective of this study was to investigate the effects of the mean flow profiles on the instability characteristics in the near-injector region of low-density gas jets injected into high-density ambient gas mediums. To achieve this, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round gas jet injected vertically upwards into a high-density ambient gas were performed by assuming three different sets of mean velocity and density profiles. The flow was assumed to be isothermal and locally parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The first set of mean velocity and density profiles assumed were those used by Monkewitz and Sohn for investigating absolute instability in hot jets. The second set of velocity and density profiles assumed for this study were the ones used by Lawson. And the third set of mean profiles included a parabolic velocity profile and a hyperbolic tangent density profile. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results for each set of mean profiles were delineated. Additional information is included in the original extended abstract.

Vedantam, Nanda Kishore↗

Instability Analysis of a Low-Density Gas Jet Injected into a High-Density Gas

The objective of this study was to determine the effects of buoyancy on the absolute instability of low-density gas jets injected into high-density gas mediums. Most of the existing analyses of low-density gas jets injected into a high-density ambient have been carried out neglecting effects of gravity. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis and a spatio-temporal stability analysis of a low-density round jet injected into a high-density ambient gas were performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. It was found that the presence of variable density within the shear layer resulted in an increase in the temporal amplification rate of the disturbances and an increase in the range of unstable frequencies, accompanied by a reduction in the phase velocities of the disturbances. Also, the temporal growth rates of the disturbances were increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity. The spatio-temporal stability analysis was performed to determine the nature of the absolute instability of the jet. The roles of the density ratio, Froude number, Schmidt number, and the lateral shift between the density and velocity profiles on the jet s absolute instability were determined. Comparisons of the results with previous experimental studies show good agreement when the effects of these variables are combined together. Thus, the combination of these variables determines how absolutely unstable the jet will be. Experiments were carried out to observe the qualitative differences between a round low-density gas jet injected into a high-density gas (helium jet injected into air) and a round constant density jet (air jet injected into air). Flow visualizations and velocity measurements in the near-injector region of the helium jet show more mixing and spreading of the helium jet than the air jet. The vortex structures develop and contribute to the jet spreading causing the helium jet to oscillate.

Lawson, Anthony Layiwola↗

Linear Temporal Stability Analysis of a Low-Density Round Gas Jet Injected into a High-Density Gas

It has been observed in previous experimental studies that round helium jets injected into air display a repetitive structure for a long distance, somewhat similar to the buoyancy-induced flickering observed in diffusion flames. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis of a round helium jet injected into air was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. The temporal growth rates of the disturbances increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity.

Lawson, Anthony L.↗

A Nonlinear Viscoelastic Model for Ceramics at High Temperatures

High-temperature creep behavior of ceramics is characterized by nonlinear time-dependent responses, asymmetric behavior in tension and compression, and nucleation and coalescence of voids leading to creep rupture. Moreover, creep rupture experiments show considerable scatter or randomness in fatigue lives of nominally equal specimens. To capture the nonlinear, asymmetric time-dependent behavior, the standard linear viscoelastic solid model is modified. Nonlinearity and asymmetry are introduced in the volumetric components by using a nonlinear function similar to a hyperbolic sine function but modified to model asymmetry. The nonlinear viscoelastic model is implemented in an ABAQUS user material subroutine. To model the random formation and coalescence of voids, each element is assigned a failure strain sampled from a lognormal distribution. An element is deleted when its volumetric strain exceeds its failure strain. Element deletion has been implemented within ABAQUS. Temporal increases in strains produce a sequential loss of elements (a model for void nucleation and growth), which in turn leads to failure. Nonlinear viscoelastic model parameters are determined from uniaxial tensile and compressive creep experiments on silicon nitride. The model is then used to predict the deformation of four-point bending and ball-on-ring specimens. Simulation is used to predict statistical moments of creep rupture lives. Numerical simulation results compare well with results of experiments of four-point bending specimens. The analytical model is intended to be used to predict the creep rupture lives of ceramic parts in arbitrary stress conditions.

Powers, Lynn M.↗

A Study of Multigrid Preconditioners Using Eigensystem Analysis

The convergence properties of numerical schemes for partial differential equations are studied by examining the eigensystem of the discrete operator. This method of analysis is very general, and allows the effects of boundary conditions and grid nonuniformities to be examined directly. Algorithms for the Laplace equation and a two equation model hyperbolic system are examined.

Roberts, Thomas W.↗

Mechanical Properties and Microstructure Investigation of Lead Free Solder

While the electronics industry appears to be focusing on Sn-Ag-Cu as the alloy of choice for lead free electronics assembly, ,the exact composition varies by geographic region, supplier and user. Add to that dissolved copper and silver from the printed circuit board traces and surface finish, and there can be significant variation in the final solder joint composition. A systematic study of the mechanical and microstructural properties of Sn-Ag-Cu alloys with Ag varying from 2wt% to 4wt% and Cu varying from 0.5wt% to lSwt%, was undertaken in this research study. Different sample preparation techniques (water quenched, oil quenched and water quenched followed by reflow) were explored and the resulting microstructure compared to that of a typical reflowed lead free chip scale package (CSP) solder joint. Tensile properties (modulus, 0.2% yield strength and the ultimate tensile strength) and creep behavior of selected alloy compositions (Sn-4Ag-1 X u , Sn-4Ag-OSCu, Sn- 2Ag-1 X u , Sn-2Ag-OSCu, Sn-3.5Ag-O.SCu) were determined for three conditions: as- cast; aged for 100 hours at 125OC; and aged for 250 hours at 125OC. There was no significant difference in Young's Modulus as a function of alloy composition. After an initial decrease in modulus after 100 hours at 125"C, there was an insignificant change with further aging. The distribution of 0.2% strain yield stress and ultimate tensile strength as a function of alloy composition was more significant and decreased with aging time and temperature. The microstructures of these alloys were examined using light and scanning electron microscopy (LM and SEM) respectively and SEM based energy dispersive x-ray spectroscopy (EDS). Fracture surface and cross-section analysis were performed on the specimens after creep testing. The creep testing results and the effect of high temperature aging on mechanical properties is presented for the oil quenched samples. In general the microstructure of oil quenched specimen exhibited a eutectic region of Sn with moderately dispersed Ag3Sn intermetallic, surrounded by a dendritic Sn-rich phase. The SEM images of the fracture surface indicated the presence of a tough shear surface at the initial cavity break area and a break line in the middle of specimen along the failure direction. A hyperbolic-sine creep model was adopted and used to fit the creep experiment data. The effect on the mechanical properties by adding the quaternary element bismuth to the Sn-3.5Ag-0.8Cu alloy was measured and compared with the mechanical properties of the ternary alloys. The results of this research study provide necessary data for the modeling of solder joint reliability for a range of Sn-Ag-Cu compositions and a baseline for evaluating the effects of subsequent quaternary additions.

Wang, Qing↗

A Robust Absorbing Boundary Condition for Compressible Flows

An absorbing non-reflecting boundary condition (NRBC) for practical computations in fluid dynamics and aeroacoustics is presented with theoretical proof. This paper is a continuation and improvement of a previous paper by the author. The absorbing NRBC technique is based on a first principle of non reflecting, which contains the essential physics that a plane wave solution of the Euler equations remains intact across the boundary. The technique is theoretically shown to work for a large class of finite volume approaches. When combined with the hyperbolic conservation laws, the NRBC is simple, robust and truly multi-dimensional; no additional implementation is needed except the prescribed physical boundary conditions. Several numerical examples in multi-dimensional spaces using two different finite volume schemes are illustrated to demonstrate its robustness in practical computations. Limitations and remedies of the technique are also discussed.

Loh, Ching Y.↗

Baseline Computational Fluid Dynamics Methodology for Longitudinal-Mode Liquid-Propellant Rocket Combustion Instability

A computational method for the analysis of longitudinal-mode liquid rocket combustion instability has been developed based on the unsteady, quasi-one-dimensional Euler equations where the combustion process source terms were introduced through the incorporation of a two-zone, linearized representation: (1) A two-parameter collapsed combustion zone at the injector face, and (2) a two-parameter distributed combustion zone based on a Lagrangian treatment of the propellant spray. The unsteady Euler equations in inhomogeneous form retain full hyperbolicity and are integrated implicitly in time using second-order, high-resolution, characteristic-based, flux-differencing spatial discretization with Roe-averaging of the Jacobian matrix. This method was initially validated against an analytical solution for nonreacting, isentropic duct acoustics with specified admittances at the inflow and outflow boundaries. For small amplitude perturbations, numerical predictions for the amplification coefficient and oscillation period were found to compare favorably with predictions from linearized small-disturbance theory as long as the grid exceeded a critical density (100 nodes/wavelength). The numerical methodology was then exercised on a generic combustor configuration using both collapsed and distributed combustion zone models with a short nozzle admittance approximation for the outflow boundary. In these cases, the response parameters were varied to determine stability limits defining resonant coupling onset.

Litchford, R. J.↗