Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite density”

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 361 records · Page 20

Numerical solution of transonic full stream function equations in conservation form

The stream function equation in conservation form is solved iteratively based on the artificial compressibility method. The density is not a unique function of the mass flux. In order to avoid the ambiguity near the sonic line, the density is updated in terms of the velocity, which is obtained through a simple integration of a first order equation step by step in the flow field. Iteration algorithms and finite difference approximations are discussed and numerical results of both conservative and nonconservative calculations are presented.

Hafez, M. M.↗

Transonic flutter calculations using the Euler equations

In transonic flutter problems where shock motion plays an important part, it is believed that accurate predictions of the flutter boundaries will require the use of codes based on the Euler equations. Only Euler codes can obtain the correct shock location and shock strength, and the crucially important shock excursion amplitude and phase lag. The present study is based on the finite volume scheme developed by Jameson and Venkatakrishnan for the 2-D unsteady Euler equations. The equations are solved in integral form on a moving grid. The variable are pressure, density, Cartesian velocity components, and total energy.

Bendiksen, Oddvar O.↗

Computing Deformations Of Rubbery Materials

Better use made of experimental data in finite-element computations. New formulation of constitutive equations of rubbery, nonlinearly elastic material enables finite-element analysis of boundary-value stress-and-strain problems involving arbitrary shapes and loads. In development of formulation, principal stretches used as arguments of strain-energy-density function.

Peng, Steven T. J.↗

Energetic particle signatures of satellites and rings in Neptune's magnetosphere

The cosmic ray system on Voyager 2 found a trapped radiation environment in Neptune's inner magnetosphere which is controlled primarily by absorption at the rings and satellite surfaces. The intensity of electrons with kinetic energies approximately greater than 1 MeV shows particularly strong and narrow signatures associated with absorption by the satellite 1989N1 at an orbital radius of 4.75 Neptune radii. Closer to the planet are several signatures of the inner satellites and rings. Absorption limits the intensity of the inner radiation belt sufficiently for the maximum intensity to occur outside the orbit of 1989N1 at a magnetic L shell of about 7. Radial profiles of the electron phase space density show that electrons diffuse inward from a source in the outer magnetosphere. Many of the inward-diffusing electrons are absorbed upon reaching a satellite orbital radius, but the finite absorption efficiency allows some of the electrons to pass by unaffected. The locations of the satellite and ring signatures also provide constraints on the nondipolar components of the planetary magnetic field.

Selesnick, R. S.↗

Effects of Subscale Size and Shape on Global Energy Dissipation in a Multiscale Model of a Fiber-Reinforced Composite Exhibiting Post-Peak Strain Softening Using Abaqus and FEAMAC

A mesh objective crack band model is implemented in the generalized method of cells (GMC) micromechanics model to predict failure of a composite repeating unit cell (RUC). The micromechanics calculations are achieved using the MAC/GMC core engine within the ImMAC suite of micromechanics codes, developed at the NASA Glenn Research Center. The microscale RUC is linked to a macroscale Abaqus/Standard finite element model using the FEAMAC multiscale framework (included in the ImMAC suite). The effects of the relationship between the characteristic length of the finite element and the size of the microscale RUC on the total energy dissipation of the multiscale model are investigated. A simple 2-D composite square subjected to uniaxial tension is used to demonstrate the effects of scaling the dimensions of the RUC such that the length of the sides of the RUC are equal to the characteristic length of the finite element. These results are compared to simulations where the size of the RUC is fixed, independent of the element size. Simulations are carried out for a variety of mesh densities and element shapes, including square and triangular. Results indicate that a consistent size and shape must be used to yield preserve energy dissipation across the scales.

Pineda, Evan, J.↗

Prediction of electron density in the boundary layer on entry vehicles with ablation

A rather general computer code has been developed for the numerical solution of the laminar boundary layer and thin shock layer (stagnation point only) equations for a multi-component gas mixture with finite reaction rates. The purpose of this paper is to indicate the capabilities of this computer program by presenting solutions for pure air flows and flows with carbon ablation. At the body surface the effects of oxidation and sublimation of carbon are taken into account with the following species included in the gas model: 0 2 , N 2 , O, N, NO, NO + , C 1 , C 2 , C 3 , CN, CO, and CO 2 . Results for the boundary layer along a 10° half-angle hyperboloid with ablation of carbon are compared to pure air results. Shock layer results for several wall temperatures and altitudes are presented and compared to pure air results. Other types of problems that can be solved with this computer code are indicated.

F. G. Blottner↗

East–West Proton Flux Anisotropy Observed with the PAMELA Mission

We present a study of the east–west anisotropy of trapped-proton fluxes in low-Earth orbit based on the measurements of the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics(PAMELA)experiment. The differential intensities of eastward- and westward-traveling protons detected in the South Atlantic Anomaly region were estimated as a function of equatorial pitch angle and drift shell, for six energy bins between80 MeV and 2 GeV. We found that, as a consequence of the strong atmospheric gradient coupled with the large gyro radius in this energy range, the intensities of eastward fluxes exceed those of westward fluxes by a factor of∼10–20. However, the reported directional asymmetry also depends on the sign of the local flux gradient, resulting in more intense westward fluxes beyond the radial distances where the inner belt peaks. PAMELA observations can be used to improve the description of the near-Earth radiation environment at lowest altitudes and highest trapping energies, where current theoretical and empirical models are affected by the largest uncertainties. Unified Astronomy Thesaurus concepts: Cosmic rays(329);Van Allen radiation belts(1758)1. Introduction Low-altitude inner-belt protons are strongly influenced by the density distribution of Earth’s atmosphere, mostly through interactions with its neutral constituents, which induce significant flux anisotropies. In prim is, the atmospheric loss cone results in a steep pitch-angle distribution, which becomes narrower for lower drift shells. A further, azimuthal anisotropy originates from finite gyro radius effects at proton energies in excess of a few tens of MeV(Haerendel1962; Lenchek & Singer1962). In fact, for a given spacecraft position, protons with the same pitch angle but different gyro phase—the azimuth angle associated with the gyration motion—have their guiding centers on different drift shells. In particular, protons from the west and from the east gyrate around magnetic field lines located at higher and lower altitudes, respectively. The guiding-center separationΔhincreaseswith increasing energy, so that protons moving eastward will encounter progressively lower drift-averaged densities, thus experiencing less atmospheric absorption; the opposite situation will occur for protons traveling westward, resulting in an east–west asymmetry of flux intensities whenΔhbecomes comparable to or larger than the flux scale height(Garmire1963; Heckman &Nakano1963).The trapped-flux anisotropy is a relevant aspect of the modeling of the low Earth orbit(LEO)radiation environment, given the significant engineering implications, especially for The Astrophysical Journal,919:114(6pp), 2021 October 1https://doi.org/10.3847/1538-4357/ac1677© 2021. The American Astronomical Society.

A. Bruno↗

Nuclear gamma rays from stellar flares

Solar flare observations are consistent with the phenomenological description that a loop of magnetic flux is convected to the surface of the star and twisted. The resulting inductive current parallel to the field is dissipated at an enhanced rate throughout the field volume by current limiting instabilities. The steady state balance between joule heating and thermal conduction along the field lines of force to the denser, cooler surface establishes a temperature distribution. The expansion of heated and ionized surface layers leads to a pressure balance and hence predictable density and X-ray emission measure. The current limitation instabilities result observationally in the parallel current being transferred to run-away ions that reach a kinetic energy of some finite fraction of the inductive potential drop. The nuclear excitation gamma rays produced by such a run-away ion current are calculated for a white dwarf flare.

Colgate, S. A.↗

On equilibrium states of heated self-gravitating gas clouds cooling by conduction in an external gravitational field

Exact analytic solutions are presented for equilibrium states of a self-gravitating one-dimensional cloud of gas, embedded in an external gravitational field due to a plane of 'stars' being heated at a rate proportional to the local gas density, and cooling by thermal conduction. It is found that the general topology of the solutions is such that the gas density has a minimum on the plane of 'stars', rising to an infinite but integrable peak away from the plane so that the total mass of gas in the cloud is finite. The results may be of interest in investigations of interstellar molecular clouds and of filamentary structures in supernova remnants as well as in the modeling of gas distributions around 'cocoon' protostars.

Lerche, I.↗

Numerical Investigation of the Micromechanics of Composite Fracture

A generalized program for the two-dimensional analysis of static crack growth problems was developed and is being tested. In addition to a detailed examination and development of fracture criteria needed to track cracks in composites at the micromechanical level, the program employs fully automatic finite element mesh generation. This is an entirely new approach in program structure that is being taken to allow for the automatic tracking of discrete crack growth. The maximum strain energy density criterion (T-criterion) was found to be a good alternative to the minimum strain energy density criterion (S-criterion).

Source record↗

Computation of hypersonic flow fields

The present discussion of the theory, solution methods, and development status of chemically reactive flow CFD codes pertaining to the flow field around a hypersonic vehicle encompasses the formulation of multicomponent viscous equations, grid-generation techniques, and finite-difference algorithms. The Navier-Stokes equations presented focus on such particularities of high speed vehicles as their aerodynamic configuration, lee and wake flow, bow shock, and chemistry and low density effects. Issues for which further development is desirable encompass boundary treatments, grid quality, shock capturing, and the decoupling of chemistry from flow variables.

Li, Chien-Peng↗

Modified Mode-I Cracked Sandwich Beam (CSB) Fracture Test

Five composite sandwich panels were fabricated using vacuum assisted resin transfer molding (VARTM). Four of these panels had E-glass/vinylester facesheets and one had carbon/epoxy facesheets. The sandwich panels had different density PVC foam cores. The four E-glass panels had core densities of 80, 100, 130, 200 kg/cu m. The sandwich with carbon/epoxy 3 facesheets had a core with density of 100 kg/cu m. Fracture tests were conducted using a modified Cracked Sandwich Beam (CSB) test configuration. Load displacement curves were obtained for loading and unloading of the specimens during crack growth. Various increments of crack growth were monitored. Critical Strain Energy Release Rates (SERR) were determined from the tests using the area method. The critical values of SERR can be considered the fracture toughness of the sandwich material. The fracture toughness ranged 367 J/sq m to 1350 J/sq m over the range of core densities. These results are compared to the Mode-I fracture toughness of the PVC foam core materials and values obtained for foam-cored sandwiches using the TSD specimen. Finite-element analyses (FEA) were performed for the test configuration and Strain Energy Release Rates were calculated using the Virtual Crack Closure Technique (VCCT). The SERR values determined from the FEA were scaled to the fracture loads, or critical loads, obtained from the modified CSB tests. These critical loads were in close agreement with the test values.

Smith, S. A.↗

Numerical Simulation of Illumination and Thermal Conditions at the Lunar Poles Using LOLA DTMs

We are interested in illumination conditions and the temperature distribution within the upper two meters of regolith near the lunar poles. Here, areas exist receiving almost constant illumination near areas in permanent shadow, which were identified as potential exploration sites for future missions. For our study a numerical simulation of the illumination and thermal environment for lunar near-polar regions is needed. Our study is based on high-resolution, twenty meters per pixel and 400 x 400 km large polar Digital Terrain Models (DTMs), which were derived from Lunar Orbiter Laser Altimeter (LOLA) data. Illumination conditions were simulated by synthetically illuminating the LOLA DTMs using the horizon method considering the Sun as an extended source. We model polar illumination for the central 50 x 50 km subset and use it as an input at each time-step (2 h) to evaluate the heating of the lunar surface and subsequent conduction in the sub-surface. At surface level we balance the incoming insolation with the subsurface conduction and radiation into space, whereas in the sub-surface we consider conduction with an additional constant radiogenic heat source at the bottom of our two-meter layer. Density is modeled as depth-dependent, the specific heat parameter as temperature-dependent and the thermal conductivity as depth- and temperature-dependent. We implemented a fully implicit finite-volume method in space and backward Euler scheme in time to solve the one-dimensional heat equation at each pixel in our 50 x 50 km DTM. Due to the non-linear dependencies of the parameters mentioned above, Newton's method is employed as the non-linear solver together with the Gauss-Seidel method as the iterative linear solver in each Newton iteration. The software is written in OpenCL and runs in parallel on the GPU cores, which allows for fast computation of large areas and long time scales.

Glaser, P.↗

Are BL Lac-type objects nearby black holes

It is pointed out that isolated black holes accreting interstellar gas can account for the characteristic properties of the Lacertids. Emission spectra for various interstellar gas densities and black hole masses are compared with the data plotted by Strittmatter et al. (1972) for the BL Lac-type objects. Rough estimates indicate that there may indeed be a finite number of stellar-mass black holes close to the earth as required by the theory. If it is determined that the BL Lac-type objects lie outside of the galactic disk a black hole accretion model may still apply if certain conditions are satisfied.

Shapiro, S. L.↗

Conservative implicit schemes for the full potential equation applied to transonic flows

Implicit approximate factorization techniques (AF) were investigated for the solution of matrix equations resulting from finite difference approximations to the full potential equation in conservation form. For transonic flows, an artificial viscosity, required to maintain stability in supersonic regions, was introduced by an upwind bias of the density. Two implicit AF procedures are presented and their convergence performance is compared with that of the standard transonic solution procedure, successive line overrelaxation (SLOR). Subcritical and supercritical test cases are considered. The results indicate that the AF schemes are substantially faster than SLOR.

Holst, T. L.↗

Fast, conservative schemes for the full potential equation applied to transonic flows

Implicit approximate factorization techniques (AF) are investigated for the solution of matrix equations resulting from finite-difference approximations to the full potential equation in conservation form. For transonic flows, an artificial viscosity, required to maintain stability in supersonic regions, is introduced by an upwind bias of the density. Two implicit AF procedures are presented, and their convergence performance is compared with that of the standard transonic solution procedure: successive line overrelaxation (SLOR). Subcritical and supercritical test cases are considered. Results indicate a substantial improvement in convergence rate for AF schemes relative to SLOR.

Holst, T. L.↗

Flow aerodynamics modeling of an MHD swirl combustor - Calculations and experimental verification

The paper describes a computer code for calculating the flow dynamics of a constant-density flow in the second-stage trumpet shaped nozzle section of a two-stage MHD swirl combustor for application to a disk generator. The primitive pressure-velocity variable, finite-difference computer code has been developed for the computation of inert nonreacting turbulent swirling flows in an axisymmetric MHD model swirl combustor. The method and program involve a staggered grid system for axial and radial velocities, and a line relaxation technique for the efficient solution of the equations. The code produces as output the flow field map of the nondimensional stream function, axial and swirl velocity. It was found that the best location for seed injection to obtain a uniform distribution at the combustor exit is in the central location for seed injected at the entrance to the second stage combustor.

Gupta, A. K.↗