Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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,117 records · Page 62

Resistive-heating instability of a compressible, force-free plasma

The thermal stability of a resistive, force-free plasma in a two-dimensional, slab geometry is investigated numerically. Linearly unstable modes that arise as an energetic consequence of Ohmic heating are examined in relation to the dynamic magnetic-tearing instability. Plasma compression is found to give rise to a set of unstable thermal modes not seen in previous incompressible studies. Such modes can occur in parameter regimes where no purely growing tearing mode exists.

Sparks, L.↗

Progress in the development of advanced computational models for the analysis of generalized supersonic jet flowfields

This paper describes two three-dimensional computational codes (PARCH and CRAFT) for computing a generalized supersonic jet flowfields with subsonic external streams, which solve the Reynolds-averaged Navier-Stokes equations in a time asymptotic manner. The application of these codes to an imperfectly expanded jet flowfield is demonstrated. Special attention is given to numerical issues and to turbulence modeling. The application of the PARCH and CRAFT codes to other related flowfield problems, such as gas turbine and nozzle flows, and to generalized missile-plume flowfields is discussed.

Dash, S. M.↗

The effect of a turbulent wake on the stagnation point. I - Skin friction results

The response of a boundary layer in the stagnation region of a two-dimensional body to fluctuations in the freestream is examined. The analysis is restricted to laminar incompressible flow. The assumed form of the velocity distribution at the edge of the boundary layer represents both a pulsation of the incoming flow, and an oscillation of the stagnation point streamline. Both features are essential in accurately representing the effect which freestream spatial and temporal nonuniformities have upon the unsteady boundary layer. Finally, a simple model is proposed which relates the characteristic parameters in a turbulent wake to the unsteady boundary-layer edge velocity. Numerical results are presented for both an arbitrary two-dimensional geometry and a circular cylinder.

Wilson, Dennis E.↗

Closed-loop transfer recovery with observer-based controllers. II - Design

The paper focuses on the design problem related to the recovery of a target closed-loop transfer function using both full-order and reduced-order observer-based controllers. Design schemes for the reduced-order observer-based controllers are similar to those of full-order controllers; hence, for simplicity, only the case for full-order observer-based controller designs is given. Three types of design schemes are discussed in detail. The first one is an asymptotic time-scale eigenstructure assignment method, and the other two are optimization-based methods where a certain norm of the recovery matrix is minimized. Asymptotic behavior associated with each design method is highlighted. Relative advantages and disadvantages of both the ATEA and optimization-based schemes are discussed. A numerical example is used to illustrate the design possibilities from each of these methods, the available design freedom, and the significance of selecting the limit of the recovery matrix, especially in the case of a nonrecoverable design.

Chen, Ben M.↗

Supercritical blade design on stream surfaces of revolution with an inverse method

A method to solve the inverse problem of supercritical blade-to-blade flow on stream surfaces of revolution with variable radius and variable stream surface thickness in a relative system is described. Some aspects of shockless design and of leading edge resolution in the numerical procedure are depicted. Some supercritical compressor cascades were designed and their complete flow field results were compared with computations of two different analysis methods.

Schmidt, E.↗

Analysis of superconducting electromagnetic finite elements based on a magnetic vector potential variational principle

Electromagnetic finite elements are extended based on a variational principle that uses the electromagnetic four potential as primary variable. The variational principle is extended to include the ability to predict a nonlinear current distribution within a conductor. The extension of this theory is first done on a normal conductor and tested on two different problems. In both problems, the geometry remains the same, but the material properties are different. The geometry is that of a 1-D infinite wire. The first problem is merely a linear control case used to validate the new theory. The second problem is made up of linear conductors with varying conductivities. Both problems perform well and predict current densities that are accurate to within a few ten thousandths of a percent of the exact values. The fourth potential is then removed, leaving only the magnetic vector potential, and the variational principle is further extended to predict magnetic potentials, magnetic fields, the number of charge carriers, and the current densities within a superconductor. The new element produces good results for the mean magnetic field, the vector potential, and the number of superconducting charge carriers despite a relatively high system condition number. The element did not perform well in predicting the current density. Numerical problems inherent to this formulation are explored and possible remedies to produce better current predicting finite elements are presented.

Schuler, James J.↗

Modeling of compressible turbulent shear flows

Despite all the recent developments in computer technologies and numerical algorithms, full numerical simulations of turbulent flows are feasible only at moderate Reynolds numbers and for flows with relatively simple geometries. The main goal of this research is to develop new second order moment closures for compressible turbulence. It has been shown that the models based on the extension of those developed originally for incompressible flows fail to adequately predict turbulent flows at high Mach numbers. In this attempt, the compressibility effects are explicitly considered. A successful development of these models that directly takes into account the compressibility effects may have a range of technological implications in the design of supersonic and hypersonic vehicles.

Liou, William W.↗

Scaling of fiber composite structures

Structural scaling is being investigated using graphite/epoxy Z-section stiffeners loaded in compression. Sections of various web, flange and web-to-flange corner radii dimensions are being considered. Initial local buckling loads and modes have been determined analytically, numerically and by experiment. A normalizing parameter has been developed for initial buckling data. Postbuckling behavior has been evaluated numerically and experimentally. The nature of the load redistribution following buckling is a key factor in relating the response of the various section sizes. Additional experiments and the development of scaling parameters for postbuckling and crippling are underway.

Wieland, Todd M.↗

Recent advances in MPD thruster research at Princeton

A summary of last years anode, plasma, and cathode findings is presented. A summary of this years activities and findings is also presented. A brief discussion of previous and current understanding is given and covers the following topics: existence of microinstabilities; the scaling of Va with the Hall parameter; the scaling of anomalous resistivity with the Hall parameter; the relation between anomalous resistivity and the anode drop; the presence of turbulence in the anode region; numerical simulation with anomalous transport; the use of magnets to decrease dissipation; performance testing with the new anode; the mechanisms behind the ionization sink; and lithiated cathode research.

Jahn, Robert G.↗

Gravity and count probabilities in an expanding universe

The time evolution of nonlinear clustering on large scales in cold dark matter, hot dark matter, and white noise models of the universe is investigated using N-body simulations performed with a tree code. Count probabilities in cubic cells are determined as functions of the cell size and the clustering state (redshift), and comparisons are made with various theoretical models. We isolate the features that appear to be the result of gravitational instability, those that depend on the initial conditions, and those that are likely a consequence of numerical limitations. More specifically, we study the development of skewness, kurtosis, and the fifth moment in relation to variance, the dependence of the void probability on time as well as on sparseness of sampling, and the overall shape of the count probability distribution. Implications of our results for theoretical and observational studies are discussed.

Bouchet, Francois R.↗

On the stability of nongyrotropic ion populations - A first (analytic and simulation) assessment

The wave and dispersion equations for perturbations propagating parallel to an ambient magnetic field in magnetoplasmas with nongyrotropic ion populations show, in general, the occurrence of coupling between the parallel (left- and right-hand circularly polarized electromagnetic and longitudinal electrostatic) eigenmodes of the associated gyrotropic medium. These interactions provide a means to driving linearly one mode with free-energy sources of other modes in homogeneous media. Different types of nongyrotropy bring about distinct classes of coupling. The stability of a hydrogen magnetoplasma with anisotropic, nongyrotropic protons that only couple the electromagnetic modes to each other is investigated analytically (via solution of the derived dispersion equation) and numerically (via simulation with a hybrid code). Nongyrotropy enhances growth and enlarges the unstable spectral range relative to the corresponding gyrotropic situation. The relevance of the properties of nongyrotropic populations to space plasma environments is also discussed.

Brinca, A. L.↗

Time-accurate Navier-Stokes calculations with multigrid acceleration

A numerical scheme to solve the unsteady Navier-Stokes equations is described. The scheme is implemented by modifying the multigrid-multiblock version of the steady Navier-Stokes equations solver, TLNS3D. The scheme is fully implicit in time and uses TLNS3D to iteratively invert the equations at each physical time step. The design objective of the scheme is unconditional stability (at least for first- and second-order discretizations of the physical time derivatives). With unconditional stability, the choice of the time step is based on the physical phenomena to be resolved rather than limited by numerical stability which is especially important for high Reynolds number viscous flows, where the spatial variation of grid cell size can be as much as six orders of magnitude. An analysis of the iterative procedure and the implementation of this procedure in TLNS3D are discussed. Numerical results are presented to show both the capabilities of the scheme and its speed up relative to the use of global minimum time stepping. Reductions in computational times of an order of magnitude are demonstrated.

Melson, N. Duane↗

Time-dependent grid adaptation for meshes of triangles and tetrahedra

This paper presents in viewgraph form a method of optimizing grid generation for unsteady CFD flow calculations that distributes the numerical error evenly throughout the mesh. Adaptive meshing is used to locally enrich in regions of relatively large errors and to locally coarsen in regions of relatively small errors. The enrichment/coarsening procedures are robust for isotropic cells; however, enrichment of high aspect ratio cells may fail near boundary surfaces with relatively large curvature. The enrichment indicator worked well for the cases shown, but in general requires user supervision for a more efficient solution.

Russ D Rausch↗

Flux pinning characteristics and irreversibility line in high temperature superconductors

The flux pinning properties in high temperature superconductors are strongly influenced by thermally activated flux motion. The scaling relation of the pinning force density and the irreversibility line in various high temperature superconductors are numerically analyzed in terms of the flux creep model. The effect of two factors, i.e., the flux pinning strength and the dimensionality of the material, on these properties are investigated. It is speculated that the irreversibility line in Bi-2212 superconductors is one order of magnitude smaller than that in Y-123, even if the flux pinning strength in Bi-2212 is improved up to the level of Y-123. It is concluded that these two factors are equally important in determination of the flux pinning characteristics at high temperatures.

Matsushita, T.↗

Self-Disclosure as a Predictor of EAP Supervisory Utilization

The value of Employee Assistance Programs (EAPS) has been cited in a variety of published papers and articles. An important managerial element relative to the assessment and referral of troubled employees has been supervisory training. There has been numerous studies highlighting the various factors and circumstances associated with supervisory behavior and EAP referrals. The inclusion of emotional awareness factors in EAP supervisory utilization has not been thoroughly investigated, although frequently found in the literature as a training and development objective for managers in business and education. The present study sought to determine what role supervisory denial and anxiety avoidance plays in confrontation of troubled employees and if admission of specific, internal emotional events is a characteristic among EAP utilizing supervisors.

Donohoe, Timothy L.↗

The Generation of Smooth High Speed Solar Wind From Plume-Interplume Mixing

Ulysses has shown that fast solar wind is extremely smooth, with a variance of less than 5%, in contrast to slow wind with a variance of approximately 30%. Now UVCS has produced the surprising result that the flow speed within coronal holes, the source of fast wind, is not at all smooth. Specifically, Giordano et al. (ApJ, v531, L79-L82, 2000) report that at 1.7 R(sub SUN) the interplume flow speed is typically more than twice the plume flow speed. Other less direct evidence supports this same result, with speeds from less than 300 to over 1000 km/s reported at approximately 5 R(sub SUN). This presents the paradox of how strongly differing plume and interplume flow speeds can exist near the Sun and be absent far from the Sun. The only answer is that plume and interplume material or momentum must be strongly mixed and that the mixing must occur mainly inside 0.3 AU to be consistent with Helios observations of smooth fast wind. Pressure balance structures (PBSs) and He abundance anomalies (Reisenfeld, et al., GRL, v26(13), 1805-1808, 1999) have been identified as interplanetary remnants of plumes, implying momentum mixing is the dominant coronal process. One possible source for plume/interplume momentum mixing is MHD Kelvin-Helmholtz (KH) shear instabilities occurring on the velocity shear interfaces. The velocity shear is a source of free energy and KH fluctuations could, through nonlinear cascade, provide the forcing required for the plasma oscillations (Cranmer, ApJ, v532, 1197-1208, 2000) reported to exist in coronal holes. The physical properties in coronal holes are now sufficiently well known that we can show plume/interplume shear interfaces become unstable to the KH instability at 5 - 10 R(sub SUN). The KH dispersion relation can be used to analyze marginal stability, the most unstable wavelengths, and linear growth rates. Numerical simulations can be used to verify results from the linear analysis and study the nonlinear development of KH modes. Here we will describe the marginal stability criterion for the KH instability, how and where this condition is met in the corona, and the possible character of the resultant fluctuations.

Suess, Steve↗

Recommendations for Achieving Accurate Numerical Simulation of Tip Clearance Flows in Transonic Compressor Rotors

The tip clearance flows of transonic compressor rotors are important because they have a significant impact on rotor and stage performance. While numerical simulations of these flows are quite sophisticated. they are seldom verified through rigorous comparisons of numerical and measured data because these kinds of measurements are rare in the detail necessary to be useful in high-speed machines. In this paper we compare measured tip clearance flow details (e.g. trajectory and radial extent) with corresponding data obtained from a numerical simulation. Recommendations for achieving accurate numerical simulation of tip clearance flows are presented based on this comparison. Laser Doppler Velocimeter (LDV) measurements acquired in a transonic compressor rotor, NASA Rotor 35, are used. The tip clearance flow field of this transonic rotor was simulated using a Navier-Stokes turbomachinery solver that incorporates an advanced k-epsilon turbulence model derived for flows that are not in local equilibrium. Comparison between measured and simulated results indicates that simulation accuracy is primarily dependent upon the ability of the numerical code to resolve important details of a wall-bounded shear layer formed by the relative motion between the over-tip leakage flow and the shroud wall. A simple method is presented for determining the strength of this shear layer.

VanZante, Dale E.↗

Effects of training on muscle O2 transport at VO2max

To quantify the relative contributions of convective and peripheral diffusive components of O2 transport to the increase in leg O2 uptake (VO2leg) at maximum O2 uptake (VO2max) after 9 wk of endurance training, 12 sedentary subjects (age 21.8 +/- 3.4 yr, VO2max 36.9 +/- 5.9 ml.min-1.kg-1) were studied. VO2max, leg blood flow (Qleg), and arterial and femoral venous PO2, and thus VO2leg, were measured while the subjects breathed room air, 15% O2, and 12% O2. The sequence of the three inspirates was balanced. After training, VO2max and VO2leg increased at each inspired O2 concentration [FIO2; mean over the 3 FIO2 values 25.2 +/- 17.8 and 36.5 +/- 33% (SD), respectively]. Before training, VO2leg and mean capillary PO2 were linearly related through the origin during hypoxia but not during room air breathing, suggesting that, at 21% O2, VO2max was not limited by O2 supply. After training, VO2leg and mean capillary PO2 at each FIO2 fell along a straight line with zero intercept, just as in athletes (Roca et al. J. Appl. Physiol. 67: 291-299, 1989). Calculated muscle O2 diffusing capacity (DO2) rose 34% while Qleg increased 19%. The relatively greater rise in DO2 increased the DO2/Qleg, which led to 9.9% greater O2 extraction. By numerical analysis, the increase in Qleg alone (constant DO2) would have raised VO2leg by 35 ml/min (mean), but that of DO2 (constant Qleg) would have increased VO2leg by 85 ml/min, more than twice as much. The sum of these individual effects (120 ml/min) was less (P = 0.013) than the observed rise of 164 ml/min (mean). This synergism (explained by the increase in DO2/Qleg) seems to be an important contribution to increases in VO2max with training.

NASA Discipline Musculoskeletal↗