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 289 records · Page 16

Bow-tie antennas on a dielectric half-space - Theory and experiment

A new formulation is discussed for the rigous calculation of the radiation pattern of a bow-tie antenna of finite length and infinitesimal thickness, placed on a lossless dielectric substrate. The analysis is based on a representation of the current density on the metal surface of the antenna as a sum of an imposed (quasistatic) term and a set of current modes with unknown amplitudes. Free-space fields that are expressed in terms of continuous spectra of symmetrized plane waves are matched to the current modes using the method of moments. The resulting set of equations are solved for the unknown current amplitudes. The calculations show that for increasing bow length the antenna impedance spirals rapidly to a value predicted by transmission line theory. The theory also shows that the E-plane pattern of a two-wavelength, 60-deg bow-tie antenna is dominated by low-loss current modes propagating at the dielectric wavenumber. As the bow tie narrows, the loss of the modes increases, and the dominant wavenumber tends to the quasi-static value. Pattern measurements made at 94 GHz are shown to agree well with theoretical predictions. Measurements for a long-wire antenna, a linear array of bow-tie elements, and a log-periodic antenna are also presented.

Compton, Richard C.↗

Resonant excitation of internal gravity waves in cluster cooling flows

Orbital motions of galaxies in X-ray clusters will resonantly excite internal gravity waves (g-waves) that propagate in the plasma. The resonance corresponds to a match between the local Brunt-Vaisala oscillation frequency in the plasma and the appropriate Fourier component of the galactic gravitational potential. Radiated internal waves have an inward group velocity that carries them to the central region of the flow where they become tightly wrapped, geometrically amplified, and eventually dissipated. Waves transport energy and angular momentum, but probably not in amounts large enough to grossly alter the structure of a cooling flow, should one be present. All physical flow variables are finite at the resonance, and all energy and angular momentum deposited by a galaxy is carried inward by the waves. Nonlinear density and velocity fluctuations are likely to result from wave amplitude growth in the central regions of the flow, and a connection with the formation of emission-line filaments is a possibility.

Balbus, Steven A.↗

Variational finite element-tensor formulation for the large deflection random vibration of composite plates

A numerical integration routine is derived from a set of unified single step integration algorithms using a weighted satisfaction of the equilibrium equation governing the large deflection random response of laminated composite plates. The equilibrium equation is derived using a constant matrix large deflection finite-element formulation. In-plane inertia terms are considered in the formulation, however, rotary inertia terms are assumed negligible. Probability density, spectral density and autocorrelation functions of the maximum displacement and strain responses are presented for three acoustic excitation levels. Classical thin plate boundary conditions and pseudo white noise excitation are used in this investigation.

Robinson, Jay H.↗

A Kernel-Free Particle-Finite Element Method for Hypervelocity Impact Simulation

An improved hybrid particle-finite element method has been developed for the simulation of hypervelocity impact problems. Unlike alternative methods, the revised formulation computes the density without reference to any kernel or interpolation functions, for either the density or the rate of dilatation. This simplifies the state space model and leads to a significant reduction in computational cost. The improved method introduces internal energy variables as generalized coordinates in a new formulation of the thermomechanical Lagrange equations. Example problems show good agreement with exact solutions in one dimension and good agreement with experimental data in a three dimensional simulation.

Park, Young-Keun↗

Low and High Frequency Models of Response Statistics of a Cylindrical Orthogrid Vehicle Panel to Acoustic Excitation

This presentation further develops the orthogrid vehicle panel work. Employed Hybrid Module capabilities to assess both low/mid frequency and high frequency models in the VA One simulation environment. The response estimates from three modeling approaches are compared to ground test measurements. Detailed Finite Element Model of the Test Article -Expect to capture both the global panel modes and the local pocket mode response, but at a considerable analysis expense (time & resources). A Composite Layered Construction equivalent global stiffness approximation using SEA -Expect to capture response of the global panel modes only. An SEA approximation using the Periodic Subsystem Formulation. A finite element model of a single periodic cell is used to derive the vibroacoustic properties of the entire periodic structure (modal density, radiation efficiency, etc. Expect to capture response at various locations on the panel (on the skin and on the ribs) with less analysis expense

Smith, Andrew↗

Gravitational Lenses and the Structure and Evolution of Galaxies

The grant has supported the completion of 16 papers and 4 conference proceedings to date. During the first year of the project we completed five papers, each of which represents a new direction in the theory and interpretation of gravitational lenses. In the first paper, "The Importance of Einstein Rings", we developed the first theory for the formation and structure of the Einstein rings formed by lensing extended sources like the host galaxies of quasar and radio sources. We applied the theory to three lenses with lensed host galaxies. For the time delay lens PG 1115+080 we found that the structure of the Einstein ring ruled out models of the gravitational potential which permitted a large Hubble constant (70 km/s Mpc). In the second paper, :Cusped Mass Models Of Gravitational Lenses", we introduced a new class of lens models where the central density is characterized by a cusp ( rho proportional to tau(sup -gamma), 1 less than gamma less than 2) as in most modern models and theories of galaxies rather than a finite core radius. In the third paper, "Global Probes of the Impact of Baryons on Dark Matter Halos", we made the first globally consistent models for the separation distribution of gravitational lenses including both galaxy and cluster lenses. We show that the key physics for the origin of the sharp separation cutoff in the separation distribution near 3 arc sec is the effect of the cooling baryons in galaxies on the density structure of the system.

Kochanek, Christopher↗

Predicting SLS Launch Environment using a Novel Multiphase Formulation

Powerful acoustic waves generated during ignition of launch vehicles may be dangerous to the vehicle, its payload, or the surrounding structures. The water-based Ignition Overpressure and Sound Suppression (IOP/SS) system at Kennedy Space Center’s (KSC) Launch Complex 39B (LC-39B) will be used to protect the Space Launch System (SLS) from the acoustic vibrations generated during launch. The IOP/SS system uses enormous amounts of water to dampen and attenuate these sound waves. To better understand the launch environment risks and to study the effectiveness of the IOP/SS system it is desirable to have time-accurate unsteady simulations of the vehicle ignition with water-based sound suppression. This paper presents results obtained with a novel, high-order accurate, and robust numerical method designed for simulating compressible multiphase flows. A positivity-preserving finite difference scheme is utilized which is formally high-order accurate and also provably robust. Robustness is critical due to the extreme nature of the flow which exhibits highly nonlinear shock and rarefaction waves interacting with liquid-gas interfaces with density ratios of the order of 1000:1. Furthermore, the high-order accuracy (and the high resolution property) is desirable for predicting wave phenomena like IOP waves since the signal can be resolved accurately and propagated long distances with fewer grid points. This finite-difference method was developed using NASA’s Launch, Ascent, and Vehicle Aerodynamics (LAVA) Cartesian immersed boundary framework. We present a validation case by applying our solver to the SLS Scale Model Acoustic Test (SMAT). The SLS SMAT is a well-instrumented 5% scale model test meant to represent the SLS at NASA KSC’s LC-39B pad. Scale IOP tests were performed with and without the sound suppression water and included many sensors which recorded the pressure waves produced during ignition. For this validation case we conduct two simulations, likewise with and without sound suppression water, and compare the SLS SMAT pressure sensor signals with our numerical signals at identical locations. Following this validation case we present a study of the SLS launch environment to examine engineering safety concerns about the mobile launch pad. Engineers at KSC redesigned the main flame deflector at LC-39B anticipating the increased loads from the SLS and to repair damage from prior Shuttle missions. This deflector redesign made use of surface pressure and temperature data from LAVA full-scale SLS simulations without the sound suppression system. The engineers were questioning the possibility of increased pressure loads on the underside of the mobile launcher due to the water in the flame trench. Based on the results established in our simulations of the SLS SMAT, we performed updated calculations for SLS at LC-39B with and without water systems active to assess the readiness of the launch pad for Artemis I launch. Our results show that the IOP/SS system is effective at reducing the overpressure signal and overall sound pressure levels felt by the vehicle and additionally that the pressure loads experienced by the mobile launcher (ML) during engine startup is not increased by the presence of water.

EGS↗

An analysis of the role of drift waves in equatorial spread F

An account is given of results of rocket measurements of the wave number spectrum of equatorial spread F irregularities, with emphasis on wavelengths less than 100 m. The measurements were made from two sounding rockets launched from Peru as part of Project Condor. The Condor density fluctuation spectra display a break at a wavelength near 100 m, identical to that found in the PLUMEX experiment (Kelley et al., 1982). The Condor data also confirm a subrange in which the density and the wave potential obey the Boltzmann relation - a strong indication of the presence of low-frequency electrostatic waves with finite wavelength parallel to the magnetic field, perhaps low-frequency drift waves as proposed by Kelley et al. The Condor data are also consistent with the previous conjecture that drift waves only exist above 300 km altitude. To investigate the difference in spectra observed over two altitude ranges, the data must be fitted to a form for the power spectrum taken from Keskinen and Ossakow (1981). The fitted spectrum, along with empirically determined growth and dissipation rates, is used to calculate the energy pumped into the spectrum at long wavelengths as well as the energy dissipated at shorter wavelengths. It is found that the energy is balanced by classical collisional effects in the low-altitude case, but energy balance in the high-altitude case requires an enhanced dissipation of about 500 times that due to classical diffusion. The model is consistent with, but does not uniquely imply, an inverse cascade of drift wave turbulence in equatorial spread F.

Labelle, J.↗

Linear mixing rule in screened binary ionic mixtures

The validity of the linear mixing rule is examined for the following two cases (1) when the response of the electron gas is taken into account in the effective ionic interaction and (2) when finite-temperature effects are included in the dielectric response of the electrons, i.e., when the ions interact with both temperature- and density-dependent screened Coulomb potentials. It is found that the linear mixing rule remains valid when the electron response is taken into account in the interionic potential at any density, even though the departure from linearity can reach a few percent for the asymmetric mixtures in the region of weak degeneracy for the electron gas. A physical explanation of this behavior is proposed which is based on a simple additional length scale.

Chabrier, G.↗

Self-organized criticality in sandpiles - Nature of the critical phenomenon

The scaling behavior of sandpile models is investigated analytically. First, it is shown that sandpile models contain a set of domain walls, referred to as troughs, which bound regions that can experience avalanches. It is further shown that the dynamics of the troughs is governed by a simple set of rules involving birth, death, and coalescence events. A simple trough model is then introduced, and it is proved that the model has a phase transition with the density of the troughs as an order parameter and that, in the thermodynamic limit, the trough density goes to zero at the transition point. Finally, it is shown that the observed scaling behavior is a consequence of finite-size effects.

Carlson, J. M.↗

A six-degree-of-freedom guidance and control analysis of Mars aerocapture

A six-degree-of-freedom (6DOF) simulation is developed to investigate the control and guidance issues of a Mars aerobraking vehicle. The guidance algorithm used is a predictor-corrector guidance formulation designed to control the exit orbital apoapsis and wedge angle using bank-angle modulation. Major features of this predictor-corrector guidance algorithm include: (1) integration of the 3DOF equations of motion within an inner-loop simulation; (2) load-relief logic; (3) finite roll rates; and (4) an aerodynamic feedback multiplier. The algorithm is capable of successfully guiding the vehicle through combinations of atmospheric density dispersions, aerodynamic mispredictions, and off-nominal atmospheric interface conditions. This study demonstrated that the addition of vehicle dynamics to the Mars aerobraking simulation does not significantly impact mission feasibility. That is, a robust control system design coupled with an adaptive guidance algorithm can assure mission success in the presence of numerous off-nominal conditions.

Powell, Richard W.↗

Performance and state-space analyses of systems using Petri nets

The goal of any modeling methodology is to develop a mathematical description of a system that is accurate in its representation and also permits analysis of structural and/or performance properties. Inherently, trade-offs exist between the level detail in the model and the ease with which analysis can be performed. Petri nets (PN's), a highly graphical modeling methodology for Discrete Event Dynamic Systems, permit representation of shared resources, finite capacities, conflict, synchronization, concurrency, and timing between state changes. By restricting the state transition time delays to the family of exponential density functions, Markov chain analysis of performance problems is possible. One major drawback of PN's is the tendency for the state-space to grow rapidly (exponential complexity) compared to increases in the PN constructs. It is the state space, or the Markov chain obtained from it, that is needed in the solution of many problems. The theory of state-space size estimation for PN's is introduced. The problem of state-space size estimation is defined, its complexities are examined, and estimation algorithms are developed. Both top-down and bottom-up approaches are pursued, and the advantages and disadvantages of each are described. Additionally, the author's research in non-exponential transition modeling for PN's is discussed. An algorithm for approximating non-exponential transitions is developed. Since only basic PN constructs are used in the approximation, theory already developed for PN's remains applicable. Comparison to results from entropy theory show the transition performance is close to the theoretic optimum. Inclusion of non-exponential transition approximations improves performance results at the expense of increased state-space size. The state-space size estimation theory provides insight and algorithms for evaluating this trade-off.

Watson, James Francis, III↗

The Dynamics of Miscible Interfaces: A Space Flight Experiment

Experiments as well as accompanying simulations are described that serve in preparation of a space flight experiment to study the dynamics of miscible interfaces. The investigation specifically addresses the importance of both nonsolenoidal effects as well as nonconventional Korteweg stresses in flows that give rise to steep but finite concentration gradients. The investigation focuses on the flow in which a less viscous fluid displaces one of higher viscosity and different density within a narrow capillary tube. The fluids are miscible in all proportions. An intruding finger forms that occupies a fraction of the total tube diameter. Depending on the flow conditions, as expressed by the Peclet number, a dimensionless viscosity ratio, and a gravity parameter, this fraction can vary between approximately 0.9 and 0.2. For large Pe values, a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the specific flow conditions, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip. For small Pe values, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). Comparisons between experiments and numerical simulations based on the 'conventional' assumption of solenoidal velocity fields and without Korteweg stresses yield poor agreement as far as the Pe value is concerned that distinguishes these two regimes. As one possibility, we attribute this lack of agreement to the disregard of these terms. An attempt is made to use scaling arguments in order to evaluate the importance of the Korteweg stresses and of the assumption of solenoidality. While these effects should be strongest in absolute terms when steep concentration fronts exist, i.e., at large Pe, they may be relatively most important at lower values of Pe. We subsequently compare these conventional simulations to more complete simulations that account for nonvanishing divergence as well as Korteweg stresses. While the exact value of the relevant stress coefficients are not known, ballpark numbers do exist, and their use in the simulations indicates that these stresses may indeed be important. We plan to evaluate these issues in detail by means of comparing a space experiment with corresponding simulations, in order to extract more accurate Korteweg stress coefficients, and to confirm or deny the importance of such stresses.

Maxworthy, Tony↗

Compression After Impact on Honeycomb Core Sandwich Panels with Thin Facesheets, Part 2: Analysis

A two part research study has been completed on the topic of compression after impact (CAI) of thin facesheet honeycomb core sandwich panels. The research has focused on both experiments and analysis in an effort to establish and validate a new understanding of the damage tolerance of these materials. Part 2, the subject of the current paper, is focused on the analysis, which corresponds to the CAI testings described in Part 1. Of interest, are sandwich panels, with aerospace applications, which consist of very thin, woven S2-fiberglass (with MTM45-1 epoxy) facesheets adhered to a Nomex honeycomb core. Two sets of materials, which were identical with the exception of the density of the honeycomb core, were tested in Part 1. The results highlighted the need for analysis methods which taken into account multiple failure modes. A finite element model (FEM) is developed here, in Part 2. A commercial implementation of the Multicontinuum Failure Theory (MCT) for progressive failure analysis (PFA) in composite laminates, Helius:MCT, is included in this model. The inclusion of PFA in the present model provided a new, unique ability to account for multiple failure modes. In addition, significant impact damage detail is included in the model. A sensitivity study, used to assess the effect of each damage parameter on overall analysis results, is included in an appendix. Analysis results are compared to the experimental results for each of the 32 CAI sandwich panel specimens tested to failure. The failure of each specimen is predicted using the high-fidelity, physicsbased analysis model developed here, and the results highlight key improvements in the understanding of honeycomb core sandwich panel CAI failure. Finally, a parametric study highlights the strength benefits compared to mass penalty for various core densities.

Mcquigg, Thomas D.↗

Microwave signature of thick-target electron beams in solar flares

The steady-state behavior of a flux of nonthermal electrons injected into a fully ionized thick target is examined. Owing to the (inverse square) energy dependence of the Coulomb collisional cross section, it is found that injected electron distributions that are monotonically decreasing functions of electron energy develop at finite depths into distributions that have 'humps' in velocity space; the electron energy corresponding to the hump correlates with the overlying particle column density to the target. This results in a two-stream unstable situation. The distribution is constantly being relaxed by quasi-linear relaxation and re-created by collisions; in this way a steady nonthermal level of Langmuir plasma waves is created, and these waves in turn produce microwave plasma radiation with a typical flux of 3 x 10 to the -16th erg/sq cm-sec. This flux can be enhanced by a factor of up to 100 by a high level of low-frequency (such as ion-acoustic) turbulence, which prevents quasi-linear relaxation for a sufficient fraction of the path length.

Emslie, A. G.↗

Influence of Compression and Shear on the Strength of Composite Laminates With Z-Pinned Reinforcement

The influence of compression and shear loads on the strength of composite laminates with z-pins is evaluated parametrically using a 2D Finite Element Code (FLASH) based on Cosserat couple stress theory. Meshes were generated for three unique combinations of z-pin diameter and density. A laminated plate theory analysis was performed on several layups to determine the bi-axial stresses in the zero degree plies. These stresses, in turn, were used to determine the magnitude of the relative load steps prescribed in the FLASH analyses. Results indicated that increasing pin density was more detrimental to in-plane compression strength than increasing pin diameter. Compression strengths of lamina without z-pins agreed well with a closed form expression derived by Budiansky and Fleck. FLASH results for lamina with z-pins were consistent with the closed form results, and FLASH results without z-pins, if the initial fiber waviness due to z-pin insertion was added to the fiber waviness in the material to yield a total misalignment. Addition of 10% shear to the compression loading significantly reduced the lamina strength compared to pure compression loading. Addition of 50% shear to the compression indicated shear yielding rather than kink band formation as the likely failure mode. Two different stiffener reinforced skin configurations with z-pins, one quasi-isotropic and one orthotropic, were also analyzed. Six unique loading cases ranging from pure compression to compression plus 50% shear were analyzed assuming material fiber waviness misalignment angles of 0, 1, and 2 degrees. Compression strength decreased with increased shear loading for both configurations, with the quasi-isotropic configuration yielding lower strengths than the orthotropic configuration.

OBrien, T. Kevin↗

Interstellar Optics

We discuss the effects of finite source size on the diffraction pattern produced by scattering in a thin screen, particularly as applied to radio-wave scattering, by density fluctuations in the interstellar plasma.

Interstellar Optics Scattering techniques telescop↗

Influence of Compression and Shear on the Strength of Composite Laminates with Z-Pinned Reinforcement

The influence of compression and shear loads on the strength of composite laminates with z-pins is evaluated parametrically using a 2D Finite Element Code (FLASH) based on Cosserat couple stress theory. Meshes were generated for three unique combinations of z-pin diameter and density. A laminated plate theory analysis was performed on several layups to determine the bi-axial stresses in the zero degree plies. These stresses, in turn, were used to determine the magnitude of the relative load steps prescribed in the FLASH analyses. Results indicated that increasing pin density was more detrimental to in-plane compression strength than increasing pin diameter. Compression strengths of lamina without z-pins agreed well with a closed form expression derived by Budiansky and Fleck. FLASH results for lamina with z-pins were consistent with the closed form results, and FLASH results without z-pins, if the initial fiber waviness due to z-pin insertion was added to the fiber waviness in the material to yield a total misalignment. Addition of 10% shear to the compression loading significantly reduced the lamina strength compared to pure compression loading. Addition of 50% shear to the compression indicated shear yielding rather than kink band formation as the likely failure mode. Two different stiffener reinforced skin configurations with z-pins, one quaiisotropic and one orthotropic, were also analyzed. Six unique loading cases ranging from pure compression to compression plus 50% shear were analyzed assuming material fiber waviness misalignment angles of 0, 1, and 2 degrees. Compression strength decreased with increased shear loading for both configurations, with the quasi-isotropic configuration yielding lower strengths than the orthotropic configuration.

O'Brien, T. Kevin↗