Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic timing analysis”

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 1,315 records · Page 73

An analysis code for the Rapid Engineering Estimation of Momentum and Energy Losses (REMEL)

Nonideal behavior has traditionally been modeled by defining efficiency (a comparison between actual and isentropic processes), and subsequent specification by empirical or heuristic methods. With the increasing complexity of aeropropulsion system designs, the reliability of these more traditional methods is uncertain. Computational fluid dynamics (CFD) and experimental methods can provide this information but are expensive in terms of human resources, cost, and time. This report discusses an alternative to empirical and CFD methods by applying classical analytical techniques and a simplified flow model to provide rapid engineering estimates of these losses based on steady, quasi-one-dimensional governing equations including viscous and heat transfer terms (estimated by Reynold's analogy). A preliminary verification of REMEL has been compared with full Navier-Stokes (FNS) and CFD boundary layer computations for several high-speed inlet and forebody designs. Current methods compare quite well with more complex method results and solutions compare very well with simple degenerate and asymptotic results such as Fanno flow, isentropic variable area flow, and a newly developed, combined variable area duct with friction flow solution. These solution comparisons may offer an alternative to transitional and CFD-intense methods for the rapid estimation of viscous and heat transfer losses in aeropropulsion systems.

Dechant, Lawrence J.↗

Application of OpenFOAM to Plume Impingement in Space Environments

After 30 years of continuous human presence in low-earth orbit, NASA is returning to the moon and eventually will go to Mars. Travelling beyond low earth orbit requires NASA to learn how humans can live in Deep Space environments – beyond the protection of Earth’s magnetosphere and at distances from Earth that prevent a quick return in case of trouble. To this end, NASA is constructing the Lunar Gateway, an ISS-like space station to be put in orbit around the moon to act as a home base for Lunar exploration for NASA astronauts. The Gateway Lunar outpost will be built incrementally, via modules which will arrive at separate times and dock to the existing structure. The incremental addition of Gateway modules, and the docking of visiting vehicles, is achieved via a sequence of firings from the approaching body’s onboard reaction control system (RCS) thrusters to achieve the required approach trajectory. The typical hypergolic chemical RCS thrusters work by firing hot gases to produce adverse thrust and the needed change in velocity to safely finish the docking process. The exhaust gas from the RCS thrusters form plumes that expand into the vacuum of space and can impinge onto the outer surfaces of the Lunar Gateway, causing unwanted forces and moments, heat loads, sediment deposition, and in extreme cases, even surface erosion - all mechanisms that can damage the Lunar Gateway and must be minimized. Both permanent and visiting modules will have this RCS thruster exhaust impingement problem. This research aims to establish existing OpenFOAM solvers as a methodology for improving simulation techniques of rocket exhaust plume impingement in space environments. The flow structure of a plume in a space environment is complex; a plume that originates from a hypergolic chemical RCS thruster and expands into a vacuum will experience several regimes of rarefication. This range includes the continuum flow in the rocket nozzle through the fully rarefied free molecular flow further from the nozzle. The flow physics is different at these two extremes, and as such, the simulation approach for plumes is generally divided into a traditional computational fluid dynamics (CFD) simulation in and near the nozzle which is coupled to a subsequent direct simulation Monte Carlo (DSMC) simulation. At this time, the scope of this research is developing, verifying, and validating a method using existing solvers in the OpenFOAM framework for performing coupled CFD/DSMC calculations to determine the extent of plume impingement loading on generic space structures. This presentation will detail code-to-code comparisons between the hyStrath dsmcFoam+ solver, developed using OpenFOAM and available as open-source, and NASA’s in-house DSMC Analysis Code (DAC). Comparisons to several open-source publication findings using DAC [3,4] are presented, and advantages of using an OpenFOAM based solver are also discussed. The presentation concludes with a discussion of future work, and a plan for coupling the dsmcFoam+ solver with CFD simulations of chemical rocket engines for unified coupled plume simulation.

DSMC↗

Human Exploration Missions Study Launch Window from Earth Orbit

The determination of orbital launch window characteristics is of major importance in the analysis of human interplanetary missions and systems. The orbital launch window characteristics are directly involved in the selection of mission trajectories, the development of orbit operational concepts, and the design of orbital launch systems. The orbital launch window problem arises because of the dynamic nature of the relative geometry between outgoing (departure) asymptote of the hyperbolic escape trajectory and the earth parking orbit. The orientation of the escape hyperbola asymptotic relative to earth is a function of time. The required hyperbola energy level also varies with time. In addition, the inertial orientation of the parking orbit is a function of time because of the perturbations caused by the Earth's oblateness. Thus, a coplanar injection onto the escape hyperbola can be made only at a point in time when the outgoing escape asymptote is contained by the plane of parking orbit. Even though this condition may be planned as a nominal situation, it will not generally represent the more probable injection geometry. The general case of an escape injection maneuver performed at a time other than the coplanar time will involve both a path angle and plane change and, therefore, a Delta(V) penalty. Usually, because of the Delta(V) penalty the actual departure injection window is smaller in duration than that determined by energy requirement alone. This report contains the formulation, characteristics, and test cases for five different launch window modes for Earth orbit. These modes are: (1) One impulsive maneuver from a Low Earth Orbit (LEO), (2) Two impulsive maneuvers from LEO, (3) Three impulsive maneuvers from LEO, (4) One impulsive maneuvers from a Highly Elliptical Orbit (HEO), (5) Two impulsive maneuvers from a Highly Elliptical Orbit (HEO) The formulation of these five different launch window modes provides a rapid means of generating realistic parametric data for space exploration studies. Also the formulation provides vector and geometrical data sufficient for use as a good starting point in detail trajectory analysis based on calculus of variations, steepest descent, or parameter optimization program techniques.

Young, Archie↗

Implementation of an Unsteady PSP System in the NASA TDT

An unsteady pressure-sensitive paint (uPSP) system has been developed to provide time-resolved pressure measurements in the NASA Langley Transonic Dynamics Tunnel (TDT). Obtaining these measurements necessitated the development of environmental enclosures to protect the high-speed camera and ultraviolet lights required for uPSP from the harsh environment present during tunnel operation. Since the facility main drive was non-functioning during the testing window, performance of the uPSP system was demonstrated using an impinging jet with a passive oscillator attachment to provide unsteady flow with a known frequency independent of amplitude. Measurements were obtained for tunnel pressures ranging from 565 to 2116 psf, and model angles of attack between -4 degrees and 4 degrees. Results indicate that the system is capable of measuring surface pressure differentials on the order of 0.01 psi at full scale with a camera frame rate of at least 10 kHz. Spectral analysis shows that the fundamental frequency of the oscillating jet is captured by the uPSP system, as are the second and third harmonics. Dynamic mode decomposition highlights the dominant coherent spatial structures of the surface pressure, along with the associated frequency and growth rate of each mode, allowing for a de-noised reconstruction of the uPSP measurements. The experimental campaign outlined within this report also confirmed compatibility of the uPSP system with the TDT facility data acquisition system, and verified the successful integration with existing processing capabilities within the NASA advanced supercomputing environment. Note: this presentation is an MP4 video with sound, color with a run time of 10 minutes 37 seconds.

unsteady pressure-sensitive paint↗

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗

An Investigation of the Elements which Contribute to Statical and Dynamical Stability, and of the Effects of Variation in Those Elements

Part 1 gives details of models tested and methods of testing of the Eiffel 36 wing alone and the JN2 aircraft. Characteristics and performance curves for standard JN are included. Part 2 presents a statistical analysis of the following: lift and drag contributed by body and chassis tested without wings; lift and drag contributed by tail, tested without wings; the effect on lift and drift of interference between the wings of a biplane combination; lift and drag contributed by the addition of body, chassis, and tail to a biplane combination; total parasite resistance; effect of varying size of tail, keeping angle of setting constant; effect of varying length of body and size of tail at the same time, keeping constant moment of tail surface about the center of gravity; forces on the tail and the effects of downwash; effect of size and setting of tail on statical longitudinal stability effects of length of body on stability; the effects of the various elements of an airplane on longitudinal stability and the placing of the force vectors. Part 3 presents the fundamental principals of dynamical stability; computations of resistance derivatives; solution of the stability equation; dynamical stability of the Curtiss JN2; tabulation of resistance derivatives; discussion of the resistance derivatives; formation and solution of stability equations; physical conceptions of the resistance derivatives; elements contributing to damping and an investigation of low speed conditions. Part 4 includes a summary of the results of the statistical investigation and a summary of the results for dynamic stability.

Klemin, Alexander↗

Topological analysis of computed three-dimensional viscous flow fields

Computed solutions of the time-dependent, Reynolds-averaged Navier-Stokes equations for three dimensional flows having thin shear layers are analyzed using topological concepts. Specific examples include the transonic flow over a body of revolution with conical afterbody at moderate angles of incidence to the free stream. Experimental flow-visualization techniques are simulated graphically to visualize the computed flow. Scalar and vector fluid dynamics properties such as pressure, shear stress, and vorticity on the body surface are presented as topological maps, and their relationship to one another in terms of orientation and singular points is discussed. The extrapolation from these surface topologies toward the understanding of external flow-field behavior is and demonstrated.

Deiwert, G. S.↗

Physical aging and its influence on the reliability of network epoxies and epoxy-matrix composites

The matrix-dominated physical and mechanical properties of a carbon fiber reinforced epoxy composite and a neat epoxy resin were found to be affected by sub-Tg annealing in nitrogen and dark atmosphere. Postcured specimens of Thornel 300 carbon-fiber/Fiberite 934 epoxy as well as Fiberite 934 epoxy resin were quenched from above Tg and given annealing at 140 C, 110 C, or 80 C, for time up to one-hundred thousand minutes. No weight loss was observed during annealing at these temperatures. Significant variations were found in density, modulus, hardness, damping, moisture absorption ability, thermal expansivity. Moisture-epoxy interactious were also studied. The kinetics of aging as well as the molecular aggregation during this densification process were monitored by differential scanning calorimetry, dynamic mechanical analysis, density gradient column, microhardness tester, Instron, and solid-state nuclear magnetic resonance spectroscopy.

Heinemann, K.↗

Transportation

A summary of tether transportation is given. Four steps were used over a period of time. First, theoretical engineering feasibility and technology requirements were determined. Then the survivors of that effort went into step two in the analysis of promising candidates. Those survivors went into the third phase which is engineering design and cost benefits. Survivors entered into the demonstration mission definition phase. Transportation studies have covered two kinds of deployments. First, steady state deployment was studied. Like the TSS, it's nearly vertical. It takes a long time to deploy and involves relatively high tether tension. Secondly, dynamic deployment was studied. Deployment started in an almost horizontal direction under a very shallow angle which allows a high deployment rate under very low tension. Momentum transfer here occurs by libration. Specific payloads were used to study tethered transportation benefits. Four transportation concepts were studied with regard to cost benefits. A tethered orbiter deboost from the space station, an OTV boost up from the Space Station, a science platform on a tether with a possible micro-g lab moving in between platform and station, and a tethered boost of payloads fromthe orbiter are the four concepts. These benefits are examined in detail.

Vontiesenhausen, G.↗

Displaying Data As Movies

NMSB Movie computer program displays large sets of data (more than million individual values). Presentation dynamic, rapidly displaying sequential image "frames" in main "movie" window. Any sequence of two-dimensional sets of data scaled between 0 and 255 (1-byte resolution) displayed as movie. Time- or slice-wise progression of data illustrated. Originally written to present data from three-dimensional ultrasonic scans of damaged aerospace composite materials, illustrates data acquired by thermal-analysis systems measuring rates of heating and cooling of various materials. Developed on Macintosh IIx computer with 8-bit color display adapter and 8 megabytes of memory using Symantec Corporation's Think C, version 4.0.

Moore, Judith G.↗

Dynamics of phase separation of binary fluids

The results of molecular-dynamics studies of surface-tension-dominated spinodal decomposition of initially well-mixed binary fluids in the absence and presence of gravity are presented. The growth exponent for the domain size and the decay exponent of the potential energy of interaction between the two species with time are found to be 0.6 +/- 0.1, inconsistent with scaling arguments based on dimensional analysis.

Ma, Wen-Jong↗

Acceleration of FDNS flow simulations using initial flowfields generated with a parabolized Navier-Stokes method

The capability of computational fluid dynamics (CFD) analysis to predict complex flowfields has been greatly advanced by the widespread use of flow simulation programs based on the Navier-Stokes (NS) equations. The flow physics are theoretically well represented, and with proper care the numerical solution should not introduce appreciable uncertainties. However, the computational cost of a typical simulation in terms of both memory and execution time are large by current standards. Therefore, the efficiency of the numerical algorithm in solving the set of model equations is one factor determining the usefulness of CFD tools. The objective of the present study is to reduce the execution time of the FDNS flow simulation code by using the parabolized Navier-Stokes (PNS) equations to provide a 'good' starting condition. The technique is not universal, however the PNS model can be applied to convection dominated flows with moderate deflection in geometries that are free of large obstructions.

Tenpas, Peter W.↗

Flight Investigation on a Fighter-type Airplane of Factors which Affect the Loads and Load Distributions on the Vertical Tail Surfaces During Rudder Kicks and Fishtails

Results are presented of a flight investigation conducted on a fighter-type airplane to determine the factors which affect the loads and load distributions on the vertical tail surfaces in maneuvers. An analysis is made of the data obtained in steady flight, rudder kicks, and fishtail maneuvers. For the rudder kicks, the significant loads were the "deflection load" resulting from an abrupt control deflection and the "dynamic load" consisting of a load corresponding to the new static equilibrium condition for the rudder deflected plus a load due to a transient overshoot. The minimum time to reach the maximum control deflection attainable by the pilot in any flight condition was found to be a constant. In the fishtail maneuvers, it was found that the pilot tends to deflect the rudder in phase with the natural frequency of the airplane. The maximum loads measured in fishtails were of the same order of magnitude as those from a rudder kick in which the rudder is returned to zero at the time of maximum sideslip.

Boshar, John↗

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗

Interface Finite Elements for the Analysis of Fracture Initiation and Progression

Progressive failure analyses (PFA) are important for the prediction of residual strength and damage tolerance of vehicle structures, and to predict the energy absorbing capability of vehicle structures under crash-type loads. Typically continuum damage mechanics (CDM) and fracture mechanics (FM) are the methods used for PFA. The method of interfacial damage mechanics (IDM) is used for PFA in this research. IDM has capabilities intermediate between CDM and FM, and is used to numerically model the initiation, growth, and arrest of cracks. IDM smooths the stress singularity at the crack tip, and is easily adaptable with other nonlinearties such as plasticity and material damage. IDM is implemented by user-defined interface elements in the ABAQUS/ Standard structural analysis software package. The structural components selected to demonstrate the effectiveness PFA using interface elements are, for the most part, those with published test data. These structural components were subjected to quasi-static loading in the tests. Thus, the ABAQUS analyses are used to predict geometrically and materially nonlinear equilibrium states. Impact loading, dynamic fracture, reflected stress waves, inertia, and time dependent material behavior are not considered.

Ambur, Damodar R.↗

Predicting survival in heart failure case and control subjects by use of fully automated methods for deriving nonlinear and conventional indices of heart rate dynamics

BACKGROUND: Despite much recent interest in quantification of heart rate variability (HRV), the prognostic value of conventional measures of HRV and of newer indices based on nonlinear dynamics is not universally accepted. METHODS AND RESULTS: We have designed algorithms for analyzing ambulatory ECG recordings and measuring HRV without human intervention, using robust methods for obtaining time-domain measures (mean and SD of heart rate), frequency-domain measures (power in the bands of 0.001 to 0.01 Hz [VLF], 0.01 to 0.15 Hz [LF], and 0.15 to 0.5 Hz [HF] and total spectral power [TP] over all three of these bands), and measures based on nonlinear dynamics (approximate entropy [ApEn], a measure of complexity, and detrended fluctuation analysis [DFA], a measure of long-term correlations). The study population consisted of chronic congestive heart failure (CHF) case patients and sex- and age-matched control subjects in the Framingham Heart Study. After exclusion of technically inadequate studies and those with atrial fibrillation, we used these algorithms to study HRV in 2-hour ambulatory ECG recordings of 69 participants (mean age, 71.7+/-8.1 years). By use of separate Cox proportional-hazards models, the conventional measures SD (P<.01), LF (P<.01), VLF (P<.05), and TP (P<.01) and the nonlinear measure DFA (P<.05) were predictors of survival over a mean follow-up period of 1.9 years; other measures, including ApEn (P>.3), were not. In multivariable models, DFA was of borderline predictive significance (P=.06) after adjustment for the diagnosis of CHF and SD. CONCLUSIONS: These results demonstrate that HRV analysis of ambulatory ECG recordings based on fully automated methods can have prognostic value in a population-based study and that nonlinear HRV indices may contribute prognostic value to complement traditional HRV measures.

Non-NASA Center↗

Bone mineral density, muscle strength, and recreational exercise in men

Muscle strength has been shown to predict bone mineral density (BMD) in women. We examined this relationship in 50 healthy men who ranged in age from 28 to 51 years (average 38.3 years). BMD of the lumbar spine, proximal femur, whole body, and tibia were measured by dual-energy x-ray absorptiometry (Hologic QDR 1000W). Dynamic strength using one repetition maximum was assessed for the biceps, quadriceps, and back extensors and for the hip abductors, adductors, and flexors. Isometric grip strength was measured by dynamometry. Daily walking mileage was assessed by 9 week stepmeter records and kinematic analysis of video filming. Subjects were designated as exercisers and nonexercisers. Exercisers participated in recreational exercise at least two times each week. The results demonstrated that BMD at all sites correlated with back and biceps strength (p < 0.01 to p = 0.0001). Body weight correlated with tibia and whole-body BMD (p < 0.001); age negatively correlated with Ward's triangle BMD (p < 0.01). In stepwise multiple regressions, back strength was the only independent predictor of spine and femoral neck density (R2 = 0.27). Further, back strength was the most robust predictor of BMD at the trochanter, Ward's triangle, whole body, and tibia, although biceps strength, age, body weight, and leg strength contributed significantly to BMD at these skeletal sites, accounting for 35-52% of the variance in BMD. Exercisers and nonexercisers were similar for walking (3.97 versus 3.94 miles/day), age (37.8 versus 38.5) years, and weight (80.0 versus 77.7 kg). However, BMD and muscle strength were significantly greater in exercises than in nonexercisers.(ABSTRACT TRUNCATED AT 250 WORDS).

NASA Center ARC↗

CFD Simulation on the J-2X Engine Exhaust in the Center-Body Diffuser and Spray Chamber at the B-2 Facility

A computational fluid dynamic (CFD) code is used to simulate the J-2X engine exhaust in the center-body diffuser and spray chamber at the Spacecraft Propulsion Facility (B-2). The CFD code is named as the space-time conservation element and solution element (CESE) Euler solver and is very robust at shock capturing. The CESE results are compared with independent analysis results obtained by using the National Combustion Code (NCC) and show excellent agreement.

Wang, Xiao-Yen↗