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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

A Micromechanics Finite Element Model for Studying the Mechanical Behavior of Spray-On Foam Insulation (SOFI)

A micromechanics model has been constructed to study the mechanical behavior of spray-on foam insulation (SOFI) for the external tank. The model was constructed using finite elements representing the fundamental repeating unit of the SOFI microstructure. The details of the micromechanics model were based on cell observations and measured average cell dimensions discerned from photomicrographs. The unit cell model is an elongated Kelvin model (fourteen-sided polyhedron with 8 hexagonal and six quadrilateral faces), which will pack to a 100% density. The cell faces and cell edges are modeled using three-dimensional 20-node brick elements. Only one-eighth of the cell is modeled due to symmetry. By exercising the model and correlating the results with the macro-mechanical foam behavior obtained through material characterization testing, the intrinsic stiffness and Poisson s Ratio of the polymeric cell walls and edges are determined as a function of temperature. The model is then exercised to study the unique and complex temperature-dependent mechanical behavior as well as the fracture initiation and propagation at the microscopic unit cell level.

Ghosn, Louis J.↗
NASA NTRS1996

Finite Element Modeling of the Thermographic Inspection for Composite Materials

The performance of composite materials is dependent on the constituent materials selected, material structural geometry, and the fabrication process. Flaws can form in composite materials as a result of the fabrication process, handling in the manufacturing environment, and exposure in the service environment to anomalous activity. Often these flaws show no indication on the surface of the material while having the potential of substantially degrading the integrity of the composite structure. For this reason it is important to have available inspection techniques that can reliably detect sub-surface defects such as inter-ply disbonds, inter-ply cracks, porosity, and density changes caused by variations in fiber volume content. Many non-destructive evaluation techniques (NDE) are capable of detecting sub-surface flaws in composite materials. These include shearography, video image correlation, ultrasonic, acoustic emissions, and X-ray. The difficulty with most of these techniques is that they are time consuming and often difficult to apply to full scale structures. An NDE technique that appears to have the capability to quickly and easily detect flaws in composite structure is thermography. This technique uses heat to detect flaws. Heat is applied to the surface of a structure with the use of a heat lamp or heat gun. A thermographic camera is then pointed at the surface and records the surface temperature as the composite structure cools. Flaws in the material will cause the thermal-mechanical material response to change. Thus, the surface over an area where a flaw is present will cool differently than regions where flaws do not exist. This paper discusses the effort made to thermo-mechanically model the thermography process. First the material properties and physical parameters used in the model will be explained. This will be followed by a detailed discussion of the finite element model used. Finally, the result of the model will be summarized along with recommendations for future work.

Bucinell, Ronald B.↗
NASA NTRS1986

Thermodynamic evaluation of transonic compressor rotors using the finite volume approach

The finite volume explicit time marching method was refined and improved. Previously, extension had been made to the finite volume method to improve the accuracy of the calculation of total pressure in inviscid flow, extend the method to allow the calculation of laminar and turbulent boundary layers in internal flows, and improve the shock capturing properties of the method by introducing a Mach number dependent interpolation scheme for the pressure used in the calculating the density. The current work extends these developments by using the new pressure interpolation scheme in two dimensional viscous calculations, including a more complete description of the viscous stresses, introducing a criteria for the transverse upwind differencing which is a function of the ratio of transverse and streamwise mass fluxes, and allowing the calculation of internal flow where boundary layers are present on both walls of the duct. The manner in which the viscous stresses are evaluated in the nonorthogonal, nonuniform grid is detailed. The convergence is investigated and results for calculations of laminar flow in a converging duct are presented. Results for calculations of transonic flow in a converging-diverging nozzle are presented and the results are compared with Sajben's measurements and calculations by others.

Nicholson, S.↗
NASA NTRS2019

Warm Absorber Diagnostics of AGN Dynamics

Warm absorbers and related phenomena are some of the observable manifestations of outflows or winds from active galactic nuclei (AGNs). Warm absorbers are common in low-luminosity AGNs. They have been extensively studied observationally and are well described by simple phenomenological models. However, major open questions remain. What is the driving mechanism? What is the density and geometrical distribution? How much associated fully ionized gas is there? What is the relation to the quasi-relativistic "ultrafast outflows"? In this paper we present synthetic spectra for the observable properties of warm absorber flows and associated quantities. We use ab initio dynamical models, i.e., solutions of the equations of motion for gas in finite difference form. The models employ various plausible assumptions for the origin of the warm absorber gas and the physical mechanisms affecting its motion. The synthetic spectra are presented as an observational test of these models. In this way we explore various scenarios for warm absorber dynamics. We show that observed spectra place certain requirements on the geometrical distribution of the warm absorber gas, and that not all dynamical scenarios are equally successful at producing spectra similar to what is observed.

Kallman, Timothy R.↗
NASA NTRS1990

Laboratory experiments on stability and entrainment of oceanic stratocumulus. Part 2: Entrainment experiment

A stratified interface is stable to the buoyancy reversal instability for surprisingly large values of D (buoyancy reversal parameter). A new instability mechanism is proposed, which considers the mixing process at the interface. For the type of density curves studied here, under strong perturbations, the mixed parcel must have a buoyancy reversal comparable to the initial stratification before the interface is unstable. This is in accord with a simple model of the interface mixing process, as well as aircraft observations of long-live marine stratocumulus clouds. These clouds' remarkable longevity in the face of finite D indicates that they can be stable (Hanson, 1984; Albrecht et al., 1985; Siems et al., 1989). It is suggested that buoyancy reversal as well as the disturbance must be large for Cloudtop Entrainment Stability. The effect of buoyancy reversal (evaporative cooling) does not always enhance the entrainment rate over that in the inert case, but it may be negligible if Ri (Richardson number) is large (Ri is larger than 50) and D is small (D is smaller than 0.5). This work may shed some light on the fundamental mechanism of the breakup process of the subtropical stratocumulus clouds into tradewind cumulus. These results may also be related to the instability in the Weddell Sea off of Antarctica.

Shy, Shenqyang S.↗
NASA NTRS1993

Coupled structural, thermal, phase-change and electromagnetic analysis for superconductors, volume 2

Two families of parametrized mixed variational principles for linear electromagnetodynamics are constructed. The first family is applicable when the current density distribution is known a priori. Its six independent fields are magnetic intensity and flux density, magnetic potential, electric intensity and flux density and electric potential. Through appropriate specialization of parameters the first principle reduces to more conventional principles proposed in the literature. The second family is appropriate when the current density distribution and a conjugate Lagrange multiplier field are adjoined, giving a total of eight independently varied fields. In this case it is shown that a conventional variational principle exists only in the time-independent (static) case. Several static functionals with reduced number of varied fields are presented. The application of one of these principles to construct finite elements with current prediction capabilities is illustrated with a numerical example.

Felippa, Carlos A.↗
NASA NTRS2002

Influence of Cooling Channel Geometry on the Thermal Response in Silicon Nitride Plates Studied

Engine manufacturers are continually attempting to improve the performance and efficiency of internal combustion engines. Usually they raise the operating temperature or reduce the cooling air requirement for the hot section turbine components. However, the success of these attempts depends on finding materials that are lightweight, are strong, and can withstand high temperatures. Ceramics are among the top candidate materials considered for such harsh applications. They hold low-density, high-temperature strength, and thermal conductivity, and they are undergoing investigation as potential materials for replacing nickel-base alloys and superalloys that are currently used for engine hot-section components. Ceramic structures can withstand higher operating temperatures and a harsh combustion environment. In addition, their low densities relative to metals help reduce component mass. The long-term objectives of the High Temperature Propulsion Components (HOTPC) Project are to develop manufacturing technology, thermal and environmental barrier coatings (TBC/EBC), and the analytical modeling capability to predict thermomechanical stresses in minimally cooled silicon nitride turbine nozzle vanes under simulated engine conditions. Two- and three-dimensional finite element analyses with TBC were conducted at the NASA Glenn Research Center. Nondestructive evaluation was used to determine processing defects. The study included conducting preliminary parametric analytical runs of heat transfer and stress analyses under steady-state conditions to demonstrate the feasibility of using cooled Si3N4 parts for turbine applications. The influence of cooling-channel shapes (such as circular, square, and ascending-order cooling channels) on cooling efficiency and thermal stresses was investigated. Temperature distributions were generated for all cases considered under both cooling and no-cooling conditions, with air being the cooling medium. The table shows the magnitude of the maximum and minimum temperature obtained for the plates under air cooling. Each channel's cross-sectional shape delivered a different temperature; however, the two-dimensional analyses for circular and square or equal-side rectangular holes produced close results. Moreover, the model of the panel with ascending order cooling channels experienced the lowest temperature. A difference of near 260 C was found among the three cooling-hole configurations investigated. The ascending-order cooling channels arrangement showed superior performance by attaining the lowest temperature (1077 C) in comparison to the circular (1379 C) and square (1343 C) channels for the same cooling-hole size. This indicates that the panel with ascending-order cooling channels is the most suitable configuration regardless of the complexity involved in its manufacture. More details pertaining to this study are reported.

Abdul-Aziz, Ali↗
NASA NTRS2011

On Calculation Methods and Results for Straight Cylindrical Roller Bearing Deflection, Stiffness, and Stress

The purpose of this study was to assess some calculation methods for quantifying the relationships of bearing geometry, material properties, load, deflection, stiffness, and stress. The scope of the work was limited to two-dimensional modeling of straight cylindrical roller bearings. Preparations for studies of dynamic response of bearings with damaged surfaces motivated this work. Studies were selected to exercise and build confidence in the numerical tools. Three calculation methods were used in this work. Two of the methods were numerical solutions of the Hertz contact approach. The third method used was a combined finite element surface integral method. Example calculations were done for a single roller loaded between an inner and outer raceway for code verification. Next, a bearing with 13 rollers and all-steel construction was used as an example to do additional code verification, including an assessment of the leading order of accuracy of the finite element and surface integral method. Results from that study show that the method is at least first-order accurate. Those results also show that the contact grid refinement has a more significant influence on precision as compared to the finite element grid refinement. To explore the influence of material properties, the 13-roller bearing was modeled as made from Nitinol 60, a material with very different properties from steel and showing some potential for bearing applications. The codes were exercised to compare contact areas and stress levels for steel and Nitinol 60 bearings operating at equivalent power density. As a step toward modeling the dynamic response of bearings having surface damage, static analyses were completed to simulate a bearing with a spall or similar damage.

Krantz, Timothy L.↗
NASA NTRS2004

Physics of Galaxy Clusters and How it Affects Cosmological Tests

The main activities in 2004 were focused on completion of the new 400 square degrees ROSAT PSPC survey for distant galaxy clusters. We observed and reduced optical spectra for all X-ray candidates and now we have complete identification for a statistically complete sample of distant 283 clusters. The papers describing the cluster catalog and first science results are in preparation and will be submitted in early 2005. We also completed a project to measure temperature and density profiles at large radii using Chandra observations of a 11 well exposed low-redshift clusters. We were able to demonstrate that the density, temperature, and total mass profiles are self-similar at large radii. This analysis has led to significant improvements in determination of the cluster baryon fraction as well as cosmologically important scaling relations, such as Mtot-T. The paper describing these results is submitted to ApJ in November, 2004. We continued to study evolution of the cluster scaling relations at high redshifts using Chandra and XMM data. We developed code for image and spectral deconvolution of the XMM observations. This code was used to reconstruct the distribution of baryons and total mass from observations of distant clusters which suffer from the finite size of the XMM PSF. This study allowed us to derive a high-redshift relation between cluster temperature and mass and compare it with the local relation obtained. The paper describing the first results is submitted to the ApJ. However, the project is still on-going as more distant cluster observations enter XMh4 and Chandra public data archives. We continued our work on improving techniques for accurate measurements of the cluster mass function and obtaining cosmological constraints from such observations. We published (ApJ, 601, 610) a study in which we derived the baryon mass function for a complete sample of low-redshift clusters. These papers argued that it was an excellent proxy for the total mass function. The baryon mass function can be used to constrain the amplitude and slope of the density fluctuation power spectrum on cluster scales. This method does not use observational determinations of the total mass and thus bypasses major uncertainties in the traditional analyses based on the X-ray temperature function. We derived the measurements for the amplitude of density perturbations, sigma8 = 0.72 +/- 0.04 and the shape parameter Omega*h = 0.l3 +/- 0.07, in good agreement with a number of independent methods.

Oliversen, Ronald J.↗
NASA NTRS1977

Wind enhanced planetary escape - Collisional modifications

Effects of collisions and finite winds characteristic of a highly perturbed atmosphere on the thermal escape of terrestrial hydrogen and helium are investigated using a Monte Carlo approach. The limiting cases of vertical and horizontal winds are considered, and the relaxation layer between the collisionless exosphere and the collision-dominated thermosphere is modeled as a plane-parallel slab of given column density, depth, and atmospheric density. For both gases, the upwardly injected flux at the base of the relaxation layer is compared with the returning downward flux distribution at the same location; the technique is also applied to the atmosphere of Titan. The results show that inclusion of collisions in the escape model for terrestrial hydrogen with winds effectively throttles the escape process, that collisional throttling is negligible for helium when the exobase temperature is at least 5000 K, and that the escape of a planetary-atmosphere constituent depends on the ratio of its gravitational and kinetic energies as well as on the ratio of its mass to that of the background gas.

Curtis, S. A.↗
NASA NTRS1990

CRIT II electric, magnetic, and density measurements within an ionizing neutral stream

Measurements from rocket-borne sensors inside a high-velocity neutral barium beam show a-factor-of-six increase in plasma density in a moving ionizing front. This region was colocated with intense fluctuating electric fields at frequencies well under the lower hybrid frequency for a barium plasma. Large quasi-dc electric and magnetic field fluctuations were also detected with a large component of the current and the electric field parallel to B(0). An Alfven wave with a finite electric field component parallel to the geomagnetic field was observed to propagate along B(0), where it was detected by an instrumented subpayload.

Swenson, C. M.↗
NASA NTRS1989

Navier-Stokes analysis of cold scramjet-afterbody flows

The progress of two efforts in coding solutions of Navier-Stokes equations is summarized. The first effort concerns a 3-D space marching parabolized Navier-Stokes (PNS) code being modified to compute the supersonic mixing flow through an internal/external expansion nozzle with multicomponent gases. The 3-D PNS equations, coupled with a set of species continuity equations, are solved using an implicit finite difference scheme. The completed work is summarized and includes code modifications for four chemical species, computing the flow upstream of the upper cowl for a theoretical air mixture, developing an initial plane solution for the inner nozzle region, and computing the flow inside the nozzle for both a N2/O2 mixture and a Freon-12/Ar mixture, and plotting density-pressure contours for the inner nozzle region. The second effort concerns a full Navier-Stokes code. The species continuity equations account for the diffusion of multiple gases. This 3-D explicit afterbody code has the ability to use high order numerical integration schemes such as the 4th order MacCormack, and the Gottlieb-MacCormack schemes. Changes to the work are listed and include, but are not limited to: (1) internal/external flow capability; (2) new treatments of the cowl wall boundary conditions and relaxed computations around the cowl region and cowl tip; (3) the entering of the thermodynamic and transport properties of Freon-12, Ar, O, and N; (4) modification to the Baldwin-Lomax turbulence model to account for turbulent eddies generated by cowl walls inside and external to the nozzle; and (5) adopting a relaxation formula to account for the turbulence in the mixing shear layer.

Baysal, Oktay↗
NASA NTRS1989

Cyclotron line features from neutron star atmospheres

Cyclotron line features from the polar cap plasmas of accreting magnetized neutron stars are modeled. The polar cap environment is modeled by a homogeneous isothermal plane-parallel slab threaded by a uniform superstrong magnetic field oriented orthogonal to the slab. Two physical cyclotron photon source functions are used. For conventional stopping depths, the emergent line spectra are fairly independent of the source functions. For hot atmospheres the details of the line transfer using relativistic kinematics differ substantially from those using nonrelativistic kinematics. Aside from finite recoil, both a finite natural linewidth and conventional stopping depths suppress the escape of photons blueward of line center and potential line features tend to be smeared out into the low-energy continuum. The formation of a line features is most favorable in low-field high-temperature atmospheres where cyclotron photon production is confined to a layer of small column density.

Wang, John C. L.↗
NASA NTRS1990

Vertical structure of accretion disks - A simplified analytical model

A simplified model of the vertical structure of accretion disks is derived. Analytical expressions for the temperature and density structure, which represent a generalization of the gray model long known in the theory of classical stellar atmospheres, are presented. The formalism naturally explains similarities and differences between the structure of a disk and a stellar atmosphere. In particular, the influence of viscous dissipation and external irradiation of the disk by the central star, as well as of the finite optical thickness of the disk, may be easily accounted for and explained by the present model.

Hubeny, I.↗
NASA NTRS1981

The field and plasma configuration of a filament overlying a solar bipolar magnetic region

This paper presents an analytic model for a finite-size straight filament suspended horizontally in a steady state over a bipolar magnetic region. The equations of magnetostatic equilibrium are integrated exactly. The solution obtained illustrates the roles played by the electric current, magnetic field, pressure, and plasma weight in the balance of force everywhere in space. A specific example of a filament of diameter 50,000 km, with a density two orders of magnitude over the corona and supported by a magnetic field of about 4 gauss is included. The filament temperature can take values ranging from a small fraction to a few times the coronal temperature, depending on the internal electric current of the filament. To produce a cool filament, such as the quiescent prominence, the solution is required to have an internal field with a strong component along the filament, giving rise to helical structures. A hot filament such as the X-ray coronal loop can be produced as a twisted magnetic flux tube embedded in a strong background field aligned parallel to the filament and having lower density and temperature. The basic steps of construction can be used to develop models more realistic than the ones presented for their analytic simplicity.

Low, B. C.↗
NASA NTRS1987

Variational approach to probabilistic finite elements

Probabilistic finite element method (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties, and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

Belytschko, T.↗
NASA NTRS1989

Axisymmetric shapes and stability of charged drops in an external electric field

A highly conducting charged drop that is surrounded by a fluid insulator of another density can be levitated by suitably applying a uniform electric field. Axisymmetric equilibrium shapes and stability of the levitated drop are found by solving simultaneously the augmented Young-Laplace equation for surface shape and the Laplace equation for the elecric field, together with constraints of fixed drop volume, charge, and center of mass. The means are a method of subdomains, finite element basis functions, and Galerkin's method of weighted residuals, all facilitated by a large-scale computer. Shape families of fixed charge are treated systematically by first-order continuation. Previous analyses by Abbas et al. in 1967 and Abbas and Latham in 1969, in which the shapes of levitated drops are approximated as spheroids, are corrected. The new analysis shows that drops charged to less than the Rayleigh limit lose shape stability at turning points, with respect to external field strength, and that the instability seen in experiments of Doyle et al. in 1964 and others is not a bifurcation to a family of two-lobed shapes, but rather is a related imperfect bifurcation.

Basaran, O. A.↗
NASA NTRS1990

Mutual coupling between circular apertures on an infinite conducting ground plane and radiating into a finite width slab

The problem of electromagnetic coupling between two horns is of interest for the Microwave Reflectometer Ionization Sensor (MRIS) that will be used in the Aeroassist Flight Experiment (AFE). Laboratory measurements of mutual coupling between conical horns (using a flat metallic reflector to simulate a critically dense plasma outside) have shown a strong dependence on the finite dimensions of the shuttle tile over the apertures. Since both, the dielectric tile and the plasma outside the tile reflect microwaves, a study should be done to isolate the two mechanisms so that the MRIS reentry flight data can be interpreted correctly. Once the coupling due to the tile itself is determined then the location of the critial electron number density layers can be determined. As a first attempt to tackle this problem the Geometrical Theory of Diffraction was used to modify the existing solution to mutual coupling between apertures with infinite dielectric sheets. By using the equivalent current method, aperture theory to determine the radiated fields inside the dielectric tiles, and ray tracing the contributions to mutual coupling were determined. Results from two cases with different tile thicknesses have indicated that the main contribution to mutual coupling is due to diffraction from the bottom and top (back and front) wedges.

Christodoulou, Christos↗