Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quasi-static approximation”

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 73 records · Page 4

On the theory of thermal runaway in the hydrogen envelope of a white dwarf

An evolutionary sequence was calculated for a white dwarf on which a hydrogen rich envelope is assumed to increase with time. The stellar models were computed by integrating numerically the system of stellar structure equations by means of the Henyey method. The accretion of matter was assumed to be quasi-static. Near the point of the chemical discontinuity, a temperature maximum is shown to arise, and hydrogen ignites; a thin shell energy source is formed which is thermally unstable. The resulting thermal runaway was investigated numerically and by means of a simple approximation based on the energy balance in the hydrogen burning shell source.

Redkoborodyy, Y. N.↗

CFD Simulations of Supersonic Highly Swirling Flow Exiting a Turbine Vane Row Compared with Experimental Observations

Marshall Space Flight Center (MSFC) has been heavily involved in developing the J2-X engine. The Center has been testing a Work Horse Gas Generator (WHGG) to supply gas products to J2-X turbine components at realistic flight-like operating conditions. Three-dimensional time accurate CFD simulations and analytical fluid analysis have been performed to support WHGG tests at MSFC. The general purpose CFD program LOCI/Chem was utilized to simulate flow of products from the WHGG through a turbine manifold, a stationary row of turbine vanes, into a Can and orifice assembly used to control the back pressure at the turbine vane row and finally through an aspirator plate and flame bucket. Simulations showed that supersonic swirling flow downstream of the turbine imparted a much higher pressure on the Can wall than expected for a non-swirling flow. This result was verified by developing an analytical model that predicts wall pressure due to swirling flow. The CFD simulations predicted that the higher downstream pressure would cause the pressure drop across the nozzle row to be approximately half the value of the test objective. With CFD support, a redesign of the Can orifice and aspirator plate was performed. WHGG experimental results and observations compared well with pre-test and post-test CFD simulations. CFD simulations for both quasi-static and transient test conditions correctly predicted the pressure environment downstream of the turbine row and the behavior of the gas generator product plume as it exited the WHGG test article, impacted the flame bucket and interacted with the external environment.

West, Jeff S.↗

Dynamic Characterization of Galfenol (Fe81.6Ga18.4)

Galfenol has the potential to transform the smart materials industry by allowing for the development of multifunctional, load-bearing devices. One of the primary technical challenges faced by this development is the very limited experimental data on Galfenol's frequency-dependent response to dynamic stress, which is critically important for the design of such devices. This report details a novel and precise characterization of the constitutive behavior of polycrystalline Galfenol (Fe81.6Ga18.4) under quasi-static (1 Hz) and dynamic (4 to 1000 Hz) stress loadings. Mechanical loads are applied using a high-frequency load frame. Quasi-static minor and major hysteresis loop measurements of magnetic flux density and strain are presented for constant electromagnet currents (0 to 1.1 A) and constant magnetic fields 0 to 14 kA/m (0 to 180 Oe). The dynamic stress amplitude for minor and major loops is 2.88 and 31.4 MPa (418 and 4550 psi), respectively. Quasi-static material properties closely match published values for similar Galfenol materials. Quasi-static actuation responses are also measured and compared to quasi-static sensing responses; the high degree of reversibility seen in the comparison is consistent with published measurements and modeling results. Dynamic major and minor loops are measured for dynamic stresses up to 31 MPa (4496 psi) and 1 kHz, and the bias condition resulting in maximum, quasi-static sensitivity. Eddy current effects are quantified by considering solid and laminated Galfenol rods. Three key sources of error in the dynamic measurements are accounted for: (1) electromagnetic noise in strain signals due to Galfenol's magnetic response, (2) error in load signals due to the inertial force of fixturing, and (3) phase misalignment between signals due to conditioning electronics. For dynamic characterization, strain error is kept below 1.2 percent of full scale by wiring two collocated gauges in series (noise cancellation) and through leadwire weaving. Inertial force error is kept below 0.41 percent by measuring the dynamic force in the specimen using a nearly collocated piezoelectric load washer. The phase response of all conditioning electronics is explicitly measured and corrected for. In general, as frequency is increased, the sensing response becomes more linear because of an increase in eddy currents. As frequency increases above approximately 100 Hz, the elbow in the strain-versus-stress response disappears as the active (soft) regime stiffens toward the passive (hard) regime. Under constant-field conditions, the loss factors of the solid rod peak between 200 and 600 Hz, rather than exhibiting a monotonic increase. Compared to the solid rod, the laminated rod exhibits much slower increases in hysteresis with frequency, and its quasi-static behavior extends to higher frequencies. The elastic modulus of the laminated rod decreases between 100 and 300 Hz; this trend is currently unexplained.

Dynamic sensing↗

Linear evolution of current sheets in sheared force-free magnetic fields with discontinuous connectivity

Thin current sheets arising in tenuous, magnetized solar coronal plasmas may constitute an important mechanism for energy buildups and subsequent energy releases; they could arise from the continuous-and-random motion of magnetic footprints associated with photospheric velocity fields. A model is presented for study of the quasi-static evolution of current sheets due to shearing of the footpoints, in a highly idealized geometry that incorporates an abrupt jump in field-line connectivity. The model highlights that formation of thin current layers and allows large shearing motions prior to violation of the linear approximation. Excess energy comparable to that released by solar flares can be stored in a sheared field.

Wolfson, Richard↗

The Compressive Behavior of Isocyanate-crosslinked Silica Aerogel at High Strain Rates

Aerogels are low-density, highly nano-porous materials. Their engineering applications are limited due to their brittleness and hydrophilicity. Recently, a strong lightweight crosslinked silica aerogel has been developed by encapsulating the skeletal framework of amine-modified silica aerogels with polyureas derived by isocyanate. The mesoporous structure of the underlying silica framework is preserved through conformal polymer coating, and the thermal conductivity remains low. Characterization has been conducted on the thermal, physical properties and the mechanical properties under quasi-static loading conditions. In this paper, we present results on the dynamic compressive behavior of the crosslinked silica aerogel (CSA) using a split Hopkinson pressure bar (SHPB). A new tubing pulse shaper was employed to help reach the dynamic stress equilibrium and constant strain rate. The stress-strain relationship was determined at high strain rates within 114-4386/s. The effects of strain rate, density, specimen thickness and water absorption on the dynamic behavior of the CSA were investigated through a series of dynamic experiments. The Young's moduli (or 0.2% offset compressive yield strengths) at a strain rate approx.350/s were determined as 10.96/2.08, 159.5/6.75, 192.2/7.68, 304.6/11.46, 407.0/20.91 and 640.5/30.47 MPa for CSA with densities 0.205, 0.454, 0.492, 0.551,0.628 and 0.731 g/cu cm, respectively. The deformation and failure behaviors of a native silica aerogel with density (0.472 g/cu cm ), approximately the same as a typical CSA sample were observed with a high speed digital camera. Digital image correlation technique was used to determine the surface strains through a series of images acquired using high speed photography. The relative uniform axial deformation indicated that localized compaction did not occur at a compressive strain level of approx.17%, suggesting most likely failure mechanism at high strain rate to be different from that under quasi-static loading condition. The Poisson s ratio was determined to be 0.162 in nonlinear regime under high strain rates. CSA samples failed generally by splitting, but were much more ductile than native silica aerogels.

Luo, H.↗

Numerical solution of three-dimensional unsteady transonic flow over wings including inviscid/viscous interactions

A numerical procedure is presented for computing the unsteady transonic flow field about three dimensional swept wings undergoing general time dependent motion. The outer inviscid portion of the flow is assumed to be governed by the modified unsteady transonic small disturbance potential equation which is integrated in the time domain by means of an efficient alternating direction implicit approximate factorization algorithm. Gross dominant effects of the shock boundary layer interaction are accounted for by a simple empirically defined model. Viscous flow regions adjacent to the wing surface and in the trailing wake are described by a set of integral equations appropriate for compressible turbulent shear layers. The two dimensional boundary layer equations are applied quasi-statically stripwise across the span. Coupling with the outer inviscid flow is implemented through use of the displacement thickness concept within the limitations of small disturbance theory. Validity of the assumptions underlying the method is established by comparison with experimental data for the flow about a high aspect ratio transport wing having an advanced airfoil section.

Rizzetta, D. P.↗

Spectral and Temporal Characteristics of X-Ray-Bright Stars in the Pleiades

We follow up our deep ROSAT imaging survey of the Pleiades (Stauffer et al. 1994) with an analysis of the spectral and temporal characteristics of the X-ray-bright stars in the Pleiades. Raymond & Smith (1977) one and two-temperature models have been used to fit the position-sensitive proportional counter (PSPC) pulse-height spectra of the dozen or so brightest sources associated with late-type Pleiades members. The best-fit temperatures suggest hot coronal temperatures for K, M, and rapidly rotating G stars, and cooler temperatures for F and slowly rotating G stars. In order to probe the many less X-ray-luminous stars, we have generated composite spectra by combining net counts from all Pleiades members according to spectral type and rotational velocity. Model fits to the composite spectra confirm the trend seen in the individual spectral fits. Particularly interesting is the apparent dependence of coronal temperature on L(sub x)/L(sub bol). A hardness-ratio analysis also confirms some of these trends. The PSPC data have also revealed a dozen or so strong X-ray flares with peak X-ray luminosities in excess of approx. 10(exp 30) ergs/sec. We have modeled the brightest of these flares with a simple quasi-static cooling loop model. The peak temperature and emission measure and the inferred electron density and plasma volume suggest a very large scale flaring event. The PSPC data were collected over a period of approx. 18 months, allowing us to search for source variability on timescales ranging from less than a day (in the case of flares) to more than a year between individual exposures. On approximately year-long timescales, roughly 25% of the late-type stars are variable. Since the Pleiades was also intensively monitored by the imaging instruments on the Einstein Observatory, we have examined X-ray luminosity variations on the 10 yr timescale between Einstein and ROSAT and find that up to 40% of the late-type stars are X-ray variable. Since there is only marginal evidence for increased variability on decade-long timescales, the variability observed on long and short timescales may have a common physical origin.

Gagne, Marc↗

Evolving force-free magnetic fields. III - States of nonequilibrium and the preflare stage

The paper considers whether a neighboring magnetostatic equilibrium exists to allow a magnetic field initially in a force-free configuration to accommodate any imposed weak pressure. The following problem is treated. The foot points of the field are fixed and the plasma is frozen into the field lines under the approximation of infinite electrical conductivity. A weak pressure is introduced. It is determined infinitesimal plasma displacements exist to adjust the field lines to a new equilibrium without changing the field line connectivity. The analysis is carried out for the bipolar force-free fields forming one of two evolutionary sequences modeling the development of the preflare stage. It was found that for the force-free field corresponding to the quasi-static stage of evolution, the neighboring magnetostatic equilibrium always exists and the imposed gas pressure can be accommodated with a slight departure of the field from being exactly force free.

Low, B. C.↗

Deflection prediction for serial manipulators

A method of deflection analysis for a general serial manipulator which considers end-effector, gravity, and inertial loads is described. The conceptualized manipulator can have many degrees of freedom and be of the most general geometry. The analysis returns end-effector error based on a quasi-static approach. Thus, small oscillations (due to vibrations) are ignored and the inertial loads are considered applied loads. A kinematic analysis based on geometric influence coefficients is used to obtaining resultant joint loads. A flexibility analysis sufficient to describe general link deformation is performed to yield the local link deformations. These local values are transformed to end-effector deflections through the use of the kinematic influence coefficients as a first order approximation. End-effector deflection due to flexibilities in the joint actuators is addressed. A description of the generalized end-effector spring is presented.

Fresonke, D. A.↗

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗

Self-similar evolution of magnetized plasmas. I - Quasi-static solution

The concept of linear expansion suggested by Wei-Hong (1989 and 1990), describes the self-similar evolution of a magnetic structure. Linear expansion can be represented by a single function which connects the evolving physical parameters of the plasma with their initial values in explicit forms. A general self-similar dynamic equation, therefore, is derived. As the first step toward more general consideration, the quasi-static solution is investigated in this paper. It is shown that a gamma = 4/3 polytrope may evolve through consecutive equilibria if its magnetic field expands self-similarly. The change of the energy everywhere inside the plasma equals the work done by the internal plasma pressure and magnetic field for the expansion. For the special case of an expanding force-free magnetic field, the self-similar expansion is a clean expansion. No free magnetic energy is left anywhere inside the magnetic structure. The approximation in quasi-state modeling of a pressure confined magnetized plasmoid is analyzed.

Yang, Wei-Hong↗

Effective permittivity of dielectric mixtures

General mixing formulas are derived for discrete scatterers immersed in a host medium. The inclusion particles are assumed to be ellipsoidal. The electric field inside the scatterers is determined by quasi-static analysis, assuming the diameter of the inclusion particles to be much smaller than one wavelength. The results are applicable to general multiphase mixtures, and the scattering ellipsoids of the different phases can have different sizes and arbitrary ellipticity distribution and axis orientation, i.e., the mixture may be isotropic or anisotropic. The resulting mixing formula is nonlinear and is suitable for iterative solutions. The formula contains a quantity called the apparent permittivity, and with different choices of this quantity, the result leads to the generalized Lorentz-Lorenz formula, the generalized Polder-van Santen formula, and the generalized coherent potential-quasicrystalline approximation formula. The results are applied to calculating the complex effective permittivity of dry and wet snow, and sea ice.

Sihvola, Ari H.↗

Hydrodynamics of coalescing binary neutron stars: Ellipsoidal treatment

We employ an approximate treatment of dissipative hydrodynamics in three dimensions to study the coalescence of binary neutron stars driven by the emission of gravitational waves. The stars are modeled as compressible ellipsoids obeying a polytropic equation of state; all internal fluid velocities are assumed to be linear functions of the coordinates. The hydrodynamics equations then reduce to a set of coupled ordinary differential equations for the evolution of the principal axes of the ellipsoids, the internal velocity parameters, and the binary orbital parameters. Gravitational radiation reaction and viscous dissipation are both incorporated. We set up exact initial binary equilibrium configurations and follow the transition from the quasi-static, secular decay of the orbit at large separation to the rapid dynamical evolution of the configurations just prior to contact. A hydrodynamical instability resulting from tidal interactions significantly accelerates the coalescence at small separation, leading to appreciable radial infall velocity and tidal lag angles near contact. This behavior is reflected in the gravitational waveforms and may be observable by gravitational wave detectors under construction. In cases where the neutron stars have spins which are not aligned with the orbital angular momentum, the spin-induced quadrupole moment can lead to precession of the orbital plane and therefore modulation of the gravitational wave amplitude even at large orbital radius. However, the amplitude of the modulation is small for typical neutron star binaries with spins much smaller than the orbital angular momentum.

Lai, Dong↗

Parameters Impacting Columnated Granular Soil Pneumatic Seal Performance

The ability of a column of loose granular soil to form a pneumatic seal was investigated by varying the diameter of the soil column, the effective column height, and the level of compaction in the soil. The soil column diameter was tested at three levels using pipes with inner diameters measuring 5.08, 10.16, and 15.24 centimeters (2, 4, and 6 inches). Soil was filled in each pipe to form 15, 30, and 45 centimeter (6, 12, and 18 inch) tall soil columns. Each diameter/height configuration was also tested at three levels of soil compaction, compared by calculating the bulk density of the soil with mass and volume measurements. GRC-1a simulant was used, with approximate low/medium/high bulk densities of 1.6, 1.75, and 1.9 g/cc achieved with a combination of vibration and tamping. The order of tests for a given column diameter was randomized and repeated three times. With the top of the soil open to atmosphere at room temperature, compressed air was injected through a small diffuser at the column base with several small downward-facing holes. The number and size of these holes was scaled such that a constant total inlet orifice area to column cross sectional area ratio was maintained for each column diameter. Inlet air pressure was slowly increased via a precision regulator to preserve quasi-static equilibrium in the soil column to minimize the impact of dynamics. Air pressure was increased all the way through the static and bubbling regimes until slugging or turbulent behavior was observed in the soil to ensure that the entire static regime had been captured during data collection. A three-factor, three-level analysis of variance (ANOVA) statistical analysis was performed on the resulting data to determine the extent to which each physical parameter impacted the soil column seal performance. It was concluded that there is statistically significant evidence that column height, and the interaction between column height and diameter impact soil seal performance. In all other cases there was insufficient data to identify a statistically significant causal relationship. Additionally, plots were generated comparing experimental data to the predictive formula developed by Ogino et al. for fluidized beds in 1993. Because this model was developed for industrial spouted fluidized beds, the accuracy of its output prior to fluidization in the static seal ‘edge case’ is unknown, especially considering in this application the working gas was diffused across the column base rather than being injected through a spout. Further, the fidelity of the model had not yet been tested with lunar soil simulants. Plotting the Ogino et al. model alongside test data allows for a more intuitive sense of the impact of test parameters on soil seal performance, as well as providing a quick means to further tune this predictive model for more accurate use with static seals across granular lunar soil simulants. The model provided by Ogino et al. was further tuned using test data to determine the degree to which a spouted bed model could be applied to a slightly modified set of testing conditions: lunar soil simulants and a more diffuse, homogeneous application of pneumatic pressure. A Matlab script was created to test different values for the leading coefficient and exponents in Ogino’s formula. The script swept preset ranges, then iterated with higher resolutions across narrower ranges to converge on the optimal value for each parameter. The resulting adjusted model was compared to the original, as well as test data with noticeable improvements across the entire test domain.

Jack Stewart↗

Modeling of High-Strain-Rate Deformation, Fracture, and Impact Behavior of Advanced Gas Turbine Engine Materials at Low and Elevated Temperatures

Gamma titanium aluminides have received considerable attention over the last decade. These alloys are known to have low density, good high temperature strength retention, and good oxidation and corrosion resistance. However, poor ductility and low fracture toughness have been the key limiting factors in the full utilization of these alloys. More recently, Gamma-met PX has been developed by GKSS, Germany. These alloys have been observed to have superior strengths at elevated temperatures and quasi-static deformation rates and good oxidation resistance at elevated temperatures when compared with other gamma titanium aluminides. The present paper discusses results of a study to understand dynamic response of gamma-met PX in uniaxial compression. The experiments were conducted by using a modified split Hopkinson pressure bar between room temperature and 900 C and strain rates of up to 3500 per second. The Gamma met PX alloy showed superior strength when compared to nickel based superalloys and other gamma titanium aluminides at all test temperatures. It also showed strain and strain-rate hardening at all levels of strain rates and temperatures and without yield anomaly up to 900 C. After approximately 600 C, thermal softening is observed at all strain rates with the rate of thermal softening increasing dramatically between 800 and 900 C. However, these flow stress levels are comparatively higher in Gamma met PX than those observed for other TiAl alloys.

Shazly, Mostafa↗

A spall and diffraction study of nanosecond pressure release across the iron ε-α phase boundary

The extreme response of polycrystalline iron at high pressures and high strain rates is revealed by means of high-power laser pulses. The compression portion of the pulse coupled with x-ray diffraction identifies the expected body-centered cubic (α) to hexagonal close packed (ε) displacive transformation. Upon release, observation shows that the complete reverse transformation takes approximately 8 ns and that the structure returns to its initial microstructural configuration, in a reversible transformation path. This is in good agreement with molecular dynamics (MD) simulations which predict an inverse dependence between transformation time and strain rate. The grain size is reduced from μm to nm range during compression and begins increasing back to the original grain size on decompression. The kinetics of the transition is dictated by heterogenous nucleation as it follows the Johnson-Mehl-Avrami-Kolmogorov equation with the appropriate time exponent of ~1. This is confirmed by MD simulations which also identify profuse twinning and dislocation generation. Here, the tensile pulse generated upon reflection at the free surface is captured by time-resolved free surface velocity measurements from which a peak tensile stress of 7 GPa is obtained, in stark contrast with its quasi-static value of ~200 MPa. At these strain rates, the strength of grain interiors, which is determined by twinning and slip exceeds the strength of the boundaries, and failure initiates preferentially in the latter.

36 MATERIALS SCIENCE↗

ARTEMIS Observations of Plasma Waves in Laminar and Perturbed Interplanetary Shocks

The ‘Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun’ (ARTEMIS) mission provides a unique opportunity to study the structure of interplanetary shocks and the associated generation of plasma waves with frequencies between ~50-8000 Hz due to its long duration electric and magnetic field burst waveform captures. We compare wave properties and occurrence rates at 11 quasi-perpendicular interplanetary shocks with burst data within 10 minutes (~3200 proton gyroradii upstream, ~1900 downstream) of the shock ramp. A perturbed shock is defined as possessing a large amplitude whistler precursor in the quasi-static magnetic field with an amplitude greater than ⅓ the difference between the upstream and downstream average magnetic field magnitudes; laminar shocks lack these large precursors and have a smooth, step function-like transition. In addition to wave modes previously observed, including ion acoustic, whistler, and electrostatic solitary waves, waves in the ion acoustic frequency range that show rapid temporal frequency change are common. Three shocks had burst captures in the ramp; of these, the two laminar shocks with burst data in the ramp contained a wide range of large amplitude wave modes in the ramp whereas the one perturbed shock contained no such waves. Thus, energy dissipation through wave-particle interactions is more prominent in these two laminar shocks than in the perturbed shock. Based on observations from all 11 shocks, the The wave occurrence rates for laminar shocks are higher in the transition region, especially the ramp, than downstream. In contrast, perturbed shocks have approximately 2-3 times the wave occurrence rate downstream than laminar shocks.

L. A. Davis↗

Aeromechanical stability of a full-scale hingeless rotor in hover

The aeroelastic stability in hover of a 9.8-m-diam, hingeless helicopter rotor system is evaluated. Rotor inplane damping data were obtained for rotor operation between 350 and 425 rpm for thrust coefficients between 0.0 and 0.12. At constant rotor rotational speed, a minimum inplane damping level was obtained at a thrust coefficient of approximately 0.02. At constant rotor thrust, a minimum in rotor inplane damping was measured at 400 rpm. Correlation with a comprehensive aeroelastic stability program is found to be good. The influence of different aerodynamic inflow models on predicted damping levels is shown. Best correlation with experimental data is obtained using a dynamic inflow model rather than static or quasi-static inflow models. Comparison with other full-scale, hingeless rotor stability data in hover is presented. Also, direct comparisons are made with data from a full-scale bearingless main rotor test performed on the same general-purpose test apparatus. Although the bearingless rotor is more highly damped at design tip speed and 1-g thrust operation, greater sensitivity to operating conditions is shown.

Peterson, R. L.↗