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 433 records · Page 24

Algorithm for Unsteady Potential Flow About Airfoils

Implicit finite-difference scheme efficiently computes unsteady potential flow about airfoils. Formulation uses density and velocity potential as dependent variables. Conservation form is retained to assure that shock wave location and speed are computed correctly. Scheme fills need for method to calculate efficiently unsteady potential flow about airfoils and to predict flutter and other unsteady aeroelastic phenomena in transonic flow regimes.

Chipman, R.↗

Systematic generation of nonlinear discretized dynamic equilibrium equations of spinning cantilevers

General nonlinear discretized governing equations of motion of spinning elastic solids and structures are adjusted for the case of a spinning cantilever with initial geometric imperfections. Consideration is given to second degree nonlinearities in the strain-displacement and velocity-displacement relationships. Parameters of the discretization are developed to include the type and number of the coordinate functions used in the admissible trial solution in order to unify the discretization approaches associated with stationarity principles. The coordinate functions comprise both sets of continuous and piecewise continuous functions employed in the Rayleigh-Ritz and the finite element methods, respectively. Coefficient matrices are provided which contain the energy density expressions and which are adaptable to computer programming.

El-Essawi, M.↗

Interpretation of the geoid

The superposition of the first satellite geoid determined by Iszak upon Ootilla's geoid was based on surface gravity determinations. Good agreement was observed except over the Pacific area of the globe. The poor agreement over the Pacific was interpreted as the result of inadequate observations there. Many geoids were determined from satellite observations, including Doppler measurements. It is found that the geoid is the result of density differences in the mantle maintained since the primeval Earth by its finite strength. Various models based on this assumption are developed.

Runcorn, S. K.↗

Involute composite design evaluation using global design sensitivity derivatives

An optimization capability for involute structures has been developed. Its key feature is the use of global material geometry variables which are so chosen that all combinations of design variables within a set of lower and upper bounds correspond to manufacturable designs. A further advantage of global variables is that their number does not increase with increasing mesh density. The accuracy of the sensitivity derivatives has been verified both through finite difference tests and through the successful use of the derivatives by an optimizer. The state of the art in composite design today is still marked by point design algorithms linked together using ad hoc methods not directly related to a manufacturing procedure. The global design sensitivity approach presented here for involutes can be applied to filament wound shells and other composite constructions using material form features peculiar to each construction. The present involute optimization technology is being applied to the Space Shuttle SRM nozzle boot ring redesigns by PDA Engineering.

Hart, J. K.↗

A continuum method for modeling surface tension

In the novel method presented for modeling the effects of surface tension on fluid motion, the interfaces between fluids with different, color-represented properties are finite-thickness transition regions across which the color varies continuously. A force density proportional to the surface curvature of constant color is defined at each point in the transition region; this force-density is normalized in such a way that the conventional description of surface tension on an interface is recovered when the ratio of local transition-reion thickness to local curvature radius approaches zero. The properties of the method are illustrated by computational results for 2D flows.

Brackbill, J. U.↗

Simulation of Low-density Nozzle Plumes in Non-zero Ambient Pressures

The direct simulation Monte-Carlo (DSMC) method was applied to the analysis of low-density nitrogen plumes exhausting from a small converging-diverging nozzle into finite ambient pressures. Two cases were considered that simulated actual test conditions in a vacuum facility. The numerical simulations readily captured the complicated flow structure of the overexpanded plumes adjusting to the finite ambient pressures, including Mach disks and barrel shaped shocks. The numerical simulations compared well to experimental data of Rothe.

Chung, Chan-Hong↗

NASA Lewis Launch Collision Probability Model Developed and Analyzed

There are nearly 10,000 tracked objects orbiting the earth. These objects encompass manned objects, active and decommissioned satellites, spent rocket bodies, and debris. They range from a few centimeters across to the size of the MIR space station. Anytime a new satellite is launched, the launch vehicle with its payload attached passes through an area of space in which these objects orbit. Although the population density of these objects is low, there always is a small but finite probability of collision between the launch vehicle and one or more of these space objects. Even though the probability of collision is very low, for some payloads even this small risk is unacceptable. To mitigate the small risk of collision associated with launching at an arbitrary time within the daily launch window, NASA performs a prelaunch mission assurance Collision Avoidance Analysis (or COLA). For the COLA of the Cassini spacecraft, the NASA Lewis Research Center conducted an in-house development and analysis of a model for launch collision probability. The model allows a minimum clearance criteria to be used with the COLA analysis to ensure an acceptably low probability of collision. If, for any given liftoff time, the nominal launch vehicle trajectory would pass a space object with less than the minimum required clearance, launch would not be attempted at that time. The model assumes that the nominal positions of the orbiting objects and of the launch vehicle can be predicted as a function of time, and therefore, that any tracked object that comes within close proximity of the launch vehicle can be identified. For any such pair, these nominal positions can be used to calculate a nominal miss distance. The actual miss distances may differ substantially from the nominal miss distance, due, in part, to the statistical uncertainty of the knowledge of the objects positions. The model further assumes that these position uncertainties can be described with position covariance matrices. With these, and some additional simplifying assumptions, a closed-form solution for the probability of collision is obtained. This solution provides clear insights into how each of the independent parameters affects the probability of collision. It shows that for a given maximum probability of collision and prior knowledge of the objects position uncertainties and sizes, only knowledge of the nominal closest approach distance is required to make the launch/no launch decision. The model was completed and used for the mission assurance COLA analysis for the Cassini spacecraft, which was launched on a Titan IV/Centaur rocket on October 15, 1997. Although the model was specifically developed for the Cassini mission, it is clearly applicable for other launches as well. The effect of COLA closures on the launch window is shown. The bar represents the entire 140-min launch window on October 15, 1997; the blackened areas represent the loss of launch opportunities due to the potential for a collision with an orbiting object.

Bollenbacher, Gary↗

Development and Characterization of Embedded Sensory Particles Using Multi-Scale 3D Digital Image Correlation

A method for detecting fatigue cracks has been explored at NASA Langley Research Center. Microscopic NiTi shape memory alloy (sensory) particles were embedded in a 7050 aluminum alloy matrix to detect the presence of fatigue cracks. Cracks exhibit an elevated stress field near their tip inducing a martensitic phase transformation in nearby sensory particles. Detectable levels of acoustic energy are emitted upon particle phase transformation such that the existence and location of fatigue cracks can be detected. To test this concept, a fatigue crack was grown in a mode-I single-edge notch fatigue crack growth specimen containing sensory particles. As the crack approached the sensory particles, measurements of particle strain, matrix-particle debonding, and phase transformation behavior of the sensory particles were performed. Full-field deformation measurements were performed using a novel multi-scale optical 3D digital image correlation (DIC) system. This information will be used in a finite element-based study to determine optimal sensory material behavior and density.

Cornell, Stephen R.↗

Microscopic-macroscopic level densities for low excitation energies

Level density ρ(E,Q) is derived within the micro-macroscopic approximation (MMA) for a system of strongly interacting Fermi particles with the energy E and additional integrals of motion Q , in line with several topics of the universal and fruitful activity of A. S. Davydov. Within the extended Thomas Fermi and semiclassical periodic orbit theory beyond the Fermi-gas saddle-point method, we obtain ρ ∝ I ν (S)/S ν , where I ν (S) is the modified Bessel function of the entropy S . For small shell-structure contribution, one finds ν = κ/2 + 1, where κ is the number of additional integrals of motion. This integer number is a dimension of Q, Q = { N, Z , …} for the case of two-component atomic nuclei, where N and Z are the numbers of neutrons and protons, respectively. For much larger shell structure contributions, one obtains ν = κ /2 + 2. The MMA level density ρ reaches the well-known Fermi gas asymptote for large excitation energies and the finite micro-canonical combinatoric limit for low excitation energies. Further, the additional integrals of motion can also be the projection of the angular momentum of a nuclear system for nuclear rotations of deformed nuclei, number of excitons for collective dynamics, and so on. Fitting the MMA total level density ρ( E , Q) for a set of the integrals of motion Q = { N, Z }, to experimental data on a long nuclear isotope chain for low excitation energies, one obtains the results for the inverse level-density parameter K , which differs significantly from those of neutron resonances due to shell, isotopic asymmetry, and pairing effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Finite Element Modeling of Thermal Cycling Induced Microcracking in Carbon/Epoxy Triaxial Braided Composites

The microcrack distribution and mass change in PR520/T700s and 3502/T700s carbon/epoxy braided composites exposed to thermal cycling was evaluated experimentally. Acoustic emission was utilized to record the crack initiation and propagation under cyclic thermal loading between -55 C and 120 C. Transverse microcrack morphology was investigated using X-ray Computed Tomography. Different performance of two kinds of composites was discovered and analyzed. Based on the observations of microcrack formation, a meso-mechanical finite element model was developed to obtain the resultant mechanical properties. The simulation results exhibited a decrease in strength and stiffness with increasing crack density. Strength and stiffness reduction versus crack densities in different orientations were compared. The changes of global mechanical behavior in both axial and transverse loading conditions were studied. Keywords: Thermal cycles; Microcrack; Finite Element Model; Braided Composite

Zhang, Chao↗

Finite-momentum instability of a dynamical axion insulator

Due to the chiral anomaly, Weyl semimetals can exhibit a signature topological magnetoelectric response known as an axion term which is determined by the microscopic band structure. In the presence of strong interactions Weyl fermions may form a chiral condensate, with the intrinsic dynamics and fluctuations of the associated condensate phase producing a dynamical contribution to the axion response. Here we show that an imbalance in the density of right- and left-handed electrons drives an instability of the chiral condensate towards finite momentum and leads to strong fluctuations in the axion response. We derive a long-wavelength theory of Lifschitz type governing the dynamics of the Goldstone mode and use this to characterize its associated spatial fluctuations, which manifest as an inhomogeneous anomalous Hall effect. In conclusion, we show that these fluctuations produce signatures in inelastic light-scattering experiments across a broad spectrum of frequencies, and can be used to determine the structure factor for the axionic collective mode.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Coronal closed structures. IV - Hydrodynamical stability and response to heating perturbations

The response of magnetically confined atmospheres to perturbations in the temperature and density distribution, and the local heating rate by means of a one-dimensional time-dependent hydrodynamical code, which incorporates the full energy, momentum and mass conservation equations is studied. These studies extend the linear instability analysis of Habbal and Rosner (1979) into the finite-amplitude regime, and generalize the confined atmosphere models of Serio et al., to the time-dependent domain. The results show that closed coronal atmospheres are stable against finite-amplitude perturbations if the chromospheric response is taken into account; and observed correlated increases in coronal density and temperature can only be achieved under quiescent conditions by increasing the heat deposition rate relatively more in the chromosphere than in the corona.

Peres, G.↗

Shell-structure and asymmetry effects in level densities

Level density [Formula: see text] is derived for a nuclear system with a given energy [Formula: see text], neutron [Formula: see text], and proton [Formula: see text] particle numbers, within the semiclassical extended Thomas–Fermi and periodic-orbit theory beyond the Fermi-gas saddle-point method. We obtain [Formula: see text], where [Formula: see text] is the modified Bessel function of the entropy [Formula: see text], and [Formula: see text] is related to the number of integrals of motion, except for the energy [Formula: see text]. For small shell structure contribution one obtains within the micro–macroscopic approximation (MMA) the value of [Formula: see text] for [Formula: see text]. In the opposite case of much larger shell structure contributions one finds a larger value of [Formula: see text]. The MMA level density [Formula: see text] reaches the well-known Fermi gas asymptote for large excitation energies, and the finite micro-canonical limit for low excitation energies. Fitting the MMA [Formula: see text] to experimental data on a long isotope chain for low excitation energies, due mainly to the shell effects, one obtains results for the inverse level density parameter [Formula: see text], which differs significantly from that of neutron resonances.

Physics↗

Explaining an unusual electromigration behavior—A comprehensive experimental and theoretical analysis using finite element method

In metallic interconnects, it is generally assumed that electromigration (EM) failure location is independent of the applied electrical current and always occurs at the highest-current-density area. Our experiments show otherwise. We designed an Al interconnect that alters its failure location by only varying the applied current density. The failure occurs near the high for a current above 2 × 10 7 A/cm 2 , but at a location with 59% of the maximum for lower current densities. Thermoreflectance thermal imaging is employed to gather time-dependent high-resolution spatial temperature distributions of the Al interconnect during EM. More importantly, we propose a computationally inexpensive 2D finite element method that tracks EM evolution in time and matches well with the observations from different experimental conditions. A detailed analysis covering the major driving forces of EM is carried out to understand the complex physics behind EM. The atomic depletion rate contributed by each force is quantitatively studied. By examining the results from every tested experimental condition, the model reveals that the temperature gradient is the key reason causing atomic depletion near the failure location. Graphical illustrations and qualitative analysis are provided to intuitively show the key findings of our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling of Fillets in Thin-Walled Structures for Dynamic Analysis

Relatively new developments in manufacturing methods have made it possible to produce machined parts with a wall thickness of less than 0.010". While the parts are being machined. Fillets of relatively large radius to thickness ratios are created by the milling tool. This paper discusses an accurate new technique for finite-element modeling of fillets for dynamic and stiffness analysis using efficient plate elements rather than more computationally intensive, high density solid meshes, A simple filleted cantilever beam of 0.040" thickness, 1.6" length, 1" depth, and 0.250" radius in the corner was modeled using both solid elements and several different plate element geometries. The finite element results were then compared with static and modal testing of a machined sample. The highest-scoring plate element technique uses a bridge of elements that is constructed through the tangent point of the Fillet radius such that the thickness of the element matches tile volume of the fillet. This model produces errors of less than 4.6 percent for static loading and less than 7.4 percent for modal analysis. This simple plate element technique will prove critical for efficient, timely, and accurate dynamic analysis of complex thin-walled structures.

Seugling, Richard M.↗

Finite element simulation of thunderstorm electrodynamics in the proximity of the storm

Observations of electric fields, Maxwell current density, and air conductivity over thunderstorms were presented. The measurements were obtained using electric field mils and conductivity probes installed on a U2 aircraft as the aircraft passed approximately directly over an active thunderstorm at an altitude of 18 to 20 km. Accurate electrical observations of this type are rare and provide important information to those involved in numerically modeling a thunderstorm. A preliminary set of computer simulations based on this data were conducted and are described. The simulations show good agreement with measurements and are used to infer the thundercloud's charging current and amount of charge exchanged per flash.

Baginski, Michael Edward↗

Modeling Space Radiation Induced Bone Changes in Rat Femurs through Finite Element Analysis

As the duration of manned missions outside of the Earth’s protective shielding increase, astronauts are at risk for exposure to space radiation. Various organ systems may be damaged due to exposure. This study investigates the bonestrength changes using finite element modeling of Long Evans rats (n=85) subjected to graded, head-only proton (0, 10, 25, and 100 cGy, 150 MeV/n) and 28silicon (0, 10, 25, and 50 cGy, 300 MeV/n) radiation. The strength of the femoral neck will be examined due its clinical relevance to hip fractures. It has been shown in previous studies that bone mineral density was not reduced at the site of fracture. These findings question whether measurements of bone mineral density may be used to assess risk of hip fracture. The mechanisms leading to the irregular relationship between bone density and strength are still uncertain within literature and investigated to greater extent in clinical applications. Finite element analysis within this study simulated physiological loading of the femoral neck. No significant changes in femoral neck strength were found across doses of proton or 28silicon head-only radiation. Future work includes performing mechanical testing of the bone samples. Moving from mouse to larger animal models may also provide the increased lifespan for assessing the long-term outcomes of radiation exposure.

Dale Johnson↗

COMGEN: A computer program for generating finite element models of composite materials at the micro level

COMGEN (Composite Model Generator) is an interactive FORTRAN program which can be used to create a wide variety of finite element models of continuous fiber composite materials at the micro level. It quickly generates batch or session files to be submitted to the finite element pre- and postprocessor PATRAN based on a few simple user inputs such as fiber diameter and percent fiber volume fraction of the composite to be analyzed. In addition, various mesh densities, boundary conditions, and loads can be assigned easily to the models within COMGEN. PATRAN uses a session file to generate finite element models and their associated loads which can then be translated to virtually any finite element analysis code such as NASTRAN or MARC.

Melis, Matthew E.↗