SEARCH · Engineering Papers

Results for “dissipative dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1978

Modeling structural damping for solids having distinct shear and dilatational loss factors

For steady state time harmonic problems (rigid format 8), the NASTRAN program as currently configured treats internal structural damping through the introduction of a single structural element damping coefficient that typically is viewed as the ratio of the complex to real modulus of elasticity. For problems dealing with two or three dimensional dynamic linear viscoelasticity (e.g. a Kelvin-Voigt viscoelastic model), the present NASTRAN capability cannot directly handle this situation wherein two independent damping coefficients are required to properly model the dissipation phenomenon. A technique is presented whereby the user can adapt the standard versions of NASTRAN (without resorting to either DMAP and/or FORTRAN coding changes) for the purpose of treating this class of problem.

Kalinowski, A. J.↗

NASA NTRS1979

Gravitational potential energy of the earth - A spherical harmonic approach

A spherical harmonic equation for the gravitational potential energy of the earth is derived for an arbitrary density distribution by conceptually bringing in mass-elements from infinity and building up the earth shell upon spherical shell. The zeroth degree term in the spherical harmonic expansion agrees with the usual expression for the energy of a radial density distribution. The second degree terms give a maximum nonhydrostatic energy in the crust and mantle of -2.77 x 10 to the 29th ergs, an order of magnitude below McKenzie's (1966) estimate. McKenzie's result stems from mathematical error. Our figure is almost identical with Kaula's (1963) estimate of the minimum shear strain energy in the mantle, a not unexpected result on the basis of the virial theorem. If the earth is assumed to be a homogeneous viscous oblate spheroid relaxing to an equilibrium shape, then a lower limit to the mantle viscosity of 1.3 x 10 to the 20th P is found by assuming that the total geothermal flux is due to viscous dissipation of energy. This number is almost six orders of magnitude below MacDonald's (1966) estimate of the viscosity and removes his objection to convection. If the nonequilibrium figure is dynamically maintained by the earth acting as a heat engine at 1% efficiency, then the viscosity is 10 to the 22nd P, a number preferred by Cathles (1975) and Peltier and Andrew (1976) as the viscosity of the mantle.

Rubincam, D. P.↗

NASA NTRS1990

Three-dimensional flux-split Euler schemes involving unstructured dynamic meshes

Improved algorithms for the solution of the 3-D time dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time integration scheme using a multistage Runge-Kutta procedure or an implicit time integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the ONERA M6 wing to demonstrate applications of the new Euler solvers. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

Batina, John T.↗

NASA NTRS1993

Turbulence modeling for high speed compressible flows

The following grant objectives were delineated in the proposal to NASA: to offer course work in computational fluid dynamics (CFD) and related areas to enable mechanical engineering students at North Carolina A&T State University (N.C. A&TSU) to pursue M.S. studies in CFD, and to enable students and faculty to engage in research in high speed compressible flows. Since no CFD-related activity existed at N.C. A&TSU before the start of the NASA grant period, training of students in the CFD area and initiation of research in high speed compressible flows were proposed as the key aspects of the project. To that end, graduate level courses in CFD, boundary layer theory, and fluid dynamics were offered. This effort included initiating a CFD course for graduate students. Also, research work was performed on studying compressibility effects in high speed flows. Specifically, a modified compressible dissipation model, which included a fourth order turbulent Mach number term, was incorporated into the SPARK code and verified for the air-air mixing layer case. The results obtained for this case were compared with a wide variety of experimental data to discern the trends in the mixing layer growth rates with varying convective Mach numbers. Comparison of the predictions of the study with the results of several analytical models was also carried out. The details of the research study are described in the publication entitled 'Compressibility Effects in Modeling Turbulent High Speed Mixing Layers,' which is attached to this report.

Chandra, Suresh↗

NASA NTRS2002

Design, Fabrication and Testing of a Crushable Energy Absorber for a Passive Earth Entry Vehicle

A conceptual study was performed to investigate the impact response of a crushable energy absorber for a passive Earth entry vehicle. The spherical energy-absorbing concept consisted of a foam-filled composite cellular structure capable of omni-directional impact-load attenuation as well as penetration resistance. Five composite cellular samples of hemispherical geometry were fabricated and tested dynamically with impact speeds varying from 30 to 42 meters per second. Theoretical crush load predictions were obtained with the aid of a generalized theory which accounts for the energy dissipated during the folding deformation of the cell-walls. Excellent correlation was obtained between theoretical predictions and experimental tests on characteristic cell-web intersections. Good correlation of theory with experiment was also found to exist for the more complex spherical cellular structures. All preliminary design requirements were met by the cellular structure concept, which exhibited a near-ideal sustained crush-load and approximately 90% crush stroke.

Kellas, Sotiris↗

NASA NTRS2013

Evaluation of Nanomaterial Approaches to Damping in Epoxy Resin and Carbon Fiber/Epoxy Composite Structures by Dynamic Mechanical Analysis

Vibration mitigation in composite structures has been demonstrated through widely varying methods which include both active and passive damping. Recently, nanomaterials have been investigated as a viable approach to composite vibration damping due to the large surface available to generate energy dissipation through friction. This work evaluates the influence of dispersed nanoparticles on the damping ratio of an epoxy matrix. Limited benefit was observed through dispersion methods, however nanoparticle application as a coating resulting in up to a three-fold increase in damping.

Miller, G.↗

NASA NTRS1986

The effects of space radiation on a chemically modified graphite-epoxy composite material

The effects of the space environment on the engineering properties and chemistry of a chemically modified T300/934 graphite-epoxy composite system are characterized. The material was subjected to 1.0 x 10 to the 10th power rads of 1.0 MeV electron irradiation under vacuum to simulate 30 years in geosynchronous earth orbit. Monotonic tension tests were performed at room temperature (75 F/24 C) and elevated temperature (250 F/121 C) on 4-ply unidirectional laminates. From these tests, inplane engineering and strength properties (E sub 1, E sub 2, Nu sub 12, G sub 12, X sub T, Y sub T) were determined. Cyclic tests were also performed to characterize energy dissipation changes due to irradiation and elevated temperature. Large diameter graphite fibers were tested to determine the effects of radiation on their stiffness and strength. No significant changes were observed. Dynamic-mechanical analysis demonstrated that the glass transition temperature was reduced by 50 F(28 C) after irradiation. Thermomechanical analysis showed the occurrence of volatile products generated upon heating of the irradiated material. The chemical modification of the epoxy did not aid in producing a material which was more radiation resistant than the standard T300/934 graphite-epoxy system. Irradiation was found to cause crosslinking and chain scission in the polymer. The latter produced low molecular weight products which plasticize the material at elevated temperatures and cause apparent material stiffening at low stresses at room temperature.

Reed, S. M.↗

NASA NTRS1990

Implicit flux-split Euler schemes for unsteady aerodynamic analysis involving unstructured dynamic meshes

Improved algorithm for the solution of the time-dependent Euler equations are presented for unsteady aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves an implicit time-integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. For example, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady and unsteady flow results are presented for the NACA 0012 airfoil to demonstrate applications of the new Euler solvers. The unsteady results were obtained for the airfoil pitching harmonically about the quarter chord. The resulting instantaneous pressure distributions and lift and moment coefficients during a cycle of motion compare well with experimental data. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

Batina, John T.↗

NASA NTRS1991

Linear prediction filters for linear and nonlinear modeled geomagnetic activity

It is shown that the Faraday loop analog model of geomagnetic activity exhibits both the directly driven and loading-unloading magnetospheric responses to solar wind input. It is further shown that the directly driven component is a linear response to loading while the loading-unloading response is nonlinear. Linear prediction filters which relate model input to output are discussed. By either allowing or suppressing the loading-unloading model response filters that relate to nonlinear or linear dynamics, respectively, have been computed. Filters that described the directly driven response are finite ranged; they asymptote to zero with increasing lag on a time scale that is fixed by the dissipation rate of the model. Filters that describe the nonlinear total response are infinite ranged; they asymptote with increasing lag to large amplitude periodic oscillations. Some implications of these infinite ranged filters are discussed.

Klimas, A. J.↗

NASA NTRS1991

Nonlinear symmetric instability and intraseasonal oscillations in the tropical atmosphere

Symmetric inertial instability (SII) is studied as a mechanism for stratospheric and tropospheric phenomena in the equatorial regions. The linear and nonlinear dynamics of SII in a two-layer, zonally symmetric model on an equatorial beta plane, in the presence of a basic flow with horizontal and vertical shear, with and without dissipative effects are investigated. Linear symmetric instabilities are, in accordance with previously published results, purely exponential, that is, nonoscillatory. Nonlinear SII, studied here for the first time on a planetary scale, can produce finite-amplitude oscillatory behavior, periodic or chaotic. The period of oscillations in the inviscid case depends on the initial data. In the presence of dissipative effects, all solutions tend to a limit cycle or to a strange attractor. The dominant period in this case, over a wide range of parameters and whether vertical shear is present or not, is in the intraseasonal, 20-30-day range. It appears therefore that nonlinear SII might be a contributing mechanism to low-frequency oscillations in the tropical atmosphere.

Zhao, J.-X.↗

NASA NTRS2002

Efficient Low Dissipative High Order Schemes for Multiscale MHD Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes for the compressible Euler and Navier-Stokes equations has been developed and verified by the authors and collaborators. These schemes are suitable for the problems in question. Basically, the scheme consists of sixth-order or higher non-dissipative spatial difference operators as the base scheme. To control the amount of numerical dissipation, multiresolution wavelets are used as sensors to adaptively limit the amount and to aid the selection and/or blending of the appropriate types of numerical dissipation to be used. Magnetohydrodynamics (MHD) waves play a key role in drag reduction in highly maneuverable high speed combat aircraft, in space weather forecasting, and in the understanding of the dynamics of the evolution of our solar system and the main sequence stars. Although there exist a few well-studied second and third-order high-resolution shock-capturing schemes for the MHD in the literature, these schemes are too diffusive and not practical for turbulence/combustion MHD flows. On the other hand, extension of higher than third-order high-resolution schemes to the MHD system of equations is not straightforward. Unlike the hydrodynamic equations, the inviscid MHD system is non-strictly hyperbolic with non-convex fluxes. The wave structures and shock types are different from their hydrodynamic counterparts. Many of the non-traditional hydrodynamic shocks are not fully understood. Consequently, reliable and highly accurate numerical schemes for multiscale MHD equations pose a great challenge to algorithm development. In addition, controlling the numerical error of the divergence free condition of the magnetic fields for high order methods has been a stumbling block. Lower order methods are not practical for the astrophysical problems in question. We propose to extend our hydrodynamics schemes to the MHD equations with several desired properties over commonly used MHD schemes.

Sjoegreen, Bjoern↗

NASA NTRS2010

Global Response to Local Ionospheric Mass Ejection

We revisit a reported "Ionospheric Mass Ejection" using prior event observations to guide a global simulation of local ionospheric outflows, global magnetospheric circulation, and plasma sheet pressurization, and comparing our results with the observed global response. Our simulation framework is based on test particle motions in the Lyon-Fedder-Mobarry (LFM) global circulation model electromagnetic fields. The inner magnetosphere is simulated with the Comprehensive Ring Current Model (CRCM) of Fok and Wolf, driven by the transpolar potential developed by the LFM magnetosphere, and includes an embedded plasmaspheric simulation. Global circulation is stimulated using the observed solar wind conditions for the period 24-25 Sept 1998. This period begins with the arrival of a Coronal Mass Ejection, initially with northward, but later with southward interplanetary magnetic field. Test particles are launched from the ionosphere with fluxes specified by local empirical relationships of outflow to electrodynamic and particle precipitation imposed by the MIlD simulation. Particles are tracked until they are lost from the system downstream or into the atmosphere, using the full equations of motion. Results are compared with the observed ring current and a simulation of polar and auroral wind outflows driven globally by solar wind dynamic pressure. We find good quantitative agreement with the observed ring current, and reasonable qualitative agreement with earlier simulation results, suggesting that the solar wind driven global simulation generates realistic energy dissipation in the ionosphere and that the Strangeway relations provide a realistic local outflow description.

Moore, T. E.↗

NASA NTRS2012

Vibration Propagation of Gear Dynamics in a Gear-Bearing-Housing System Using Mathematical Modeling and Finite Element Analysis

Vibration and noise caused by gear dynamics at the meshing teeth propagate through power transmission components to the surrounding environment. This study is devoted to developing computational tools to investigate the vibro-acoustic propagation of gear dynamics through a gearbox using different bearings. Detailed finite element/contact mechanics and boundary element models of the gear/bearing/housing system are established to compute the system vibration and noise propagation. Both vibration and acoustic models are validated by experiments including the vibration modal testing and sound field measurements. The effectiveness of each bearing type to disrupt vibration propagation is speed-dependent. Housing plays an important role in noise radiation .It, however, has limited effects on gear dynamics. Bearings are critical components in drivetrains. Accurate modeling of rolling element bearings is essential to assess vibration and noise of drivetrain systems. This study also seeks to fully describe the vibro-acoustic propagation of gear dynamics through a power-transmission system using rolling element and fluid film wave bearings. Fluid film wave bearings, which have higher damping than rolling element bearings, could offer an energy dissipation mechanism that reduces the gearbox noise. The effectiveness of each bearing type to disrupt vibration propagation in explored using multi-body computational models. These models include gears, shafts, rolling element and fluid film wave bearings, and the housing. Radiated noise is mapped from the gearbox surface to surrounding environment. The effectiveness of rolling element and fluid film wave bearings in breaking the vibro-acoustic propagation path from the gear to the housing is investigated.

Parker, Robert G.↗

NASA NTRS2014

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗

NASA NTRS1981

On the role of the lower hybrid drift instability in substorm dynamics

Recent studies of the lower hybrid drift instability have shed new light on the role of this mode in field-reversed plasmas. For substorm magnetotail conditions it is found that the lower hybrid drift instability can penetrate to the neutral line and can dissipate magnetic energy at a rate of approximately 4 x 10 to the 17th erg/s. Thus this instability is capable of playing a major role in the onset of substorms and providing resistivity for reconnection processes in the context of the neutral line substorm model.

Huba, J. D.↗

NASA NTRS1998

EOS Interdisciplinary Investigation: Observational and Modeling Studies of Radiative, Chemical, and Dynamical Interactions in the Earth's Atmosphere

A 3-dimensional model of atmospheric dynamics and photochemistry has been developed from the primitive equations in isentropic coordinates. The model extends from the upper troposphere through the middle atmosphere and is driven by observed dynamical and chemical structure in the troposphere. The model's formulation is entirely spectral: Horizontal structure is represented in terms of vector Hough functions, which explicitly partition the motion into its rotational and divergent components. This formulation enables the model's computational performance to be increased dramatically by filtering high-frequency gravity waves, which do not affect PV conservation. Vertical structure is represented in terms of eigenfunctions that follow from the primitive equations in isentropic coordinates. The model's fully spectral formulation enables scale-selective dissipation, necessary for numerical stability, to be applied at 6th order-in all 3 coordinate directions. This feature leaves all but the shortest vertical scales undamped and, consequently, allows potential vorticity and chemical tracers to be conserved quite accurately, distinctly better than in the model's layered counterpart. These dynamical features are complemented by a basic but fairly complete treatment of gas phase photochemistry, which accounts for some 48 chemical species, diurnally-varying SW absorption by O2 and O3, and LW cooling calculated via a detailed band calculation of radiative transfer. Driven by observed tropospheric structure, the 3D model reproduces observed structure in the middle atmosphere, inclusive of transport by the Brewer-Dobson circulation and accompanying chemical variations. Calculated residual motion is consistent with diabatic cooling rates, poleward transport, and descent in the polar night derived from satellite measurements of chemical tracers made from LIMS, SAGE, and UARS. The model also reproduces the observed structure of chemical species like O3, HNO3, and ClO. At the same time, it provides a respectable 3D simulation of the wintertime increase of midlatitude total ozone, as well as its geographical structure observed by Nimbus-7 TOMS. The 3D model has been used to explore the origin of diabatic upwelling in the tropical lower stratosphere, which has been related to planetary wave absorption in the winter- time stratosphere, as well as to influences from tropospheric convection at the tropical tropopause.

Salby, Murry L.↗

NASA NTRS1989

Small-scale structure and turbulence observed in MAP/WINE)

During MAP/WINE small scale structure and turbulence in the mesosphere and lower thermosphere was studied in situ by rocket-borne instruments as well as from the ground by remote sensing techniques. The eight salvoes launched during the campaign resulted in a wealth of information on the dynamical structure of these regions. The experimental results are reviewed and their interpretation is discussed in terms of gravity waves and turbulence. It is shown that eddy diffusion coefficients and turbulent energy dissipation rates may be derived from the in situ measurements in a consistent manner. The observations are also shown to be consistent with the hypothesis that turbulence can be created by a process of gravity wave saturation.

Blix, T. A.↗

NASA NTRS1988

Unsteady transonic airfoil computation using implicit Euler scheme on body-fixed grid

The unsteady Euler equations have been derived for the flow relative motion with respect to a frame of reference that is rigidly attached to the moving airfoil. The grid is generated once by an elliptic solver without a need for dynamic grid computation. An implicit factored finite-volume scheme has been developed and implemented through a fully vectorized computer program. Implicit second-order and explicit second and fourth-order dissipations are added to the scheme. The boundary conditions are explicitly satisfied. The scheme is applied to steady and unsteady transonic airfoil flows and the results are in good agreement with the experimental data. For forced harmonic airfoil motions, periodic solutions are achieved within the third cycle of oscillation.

Kandil, Osama A.↗