Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “current conservation”

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 451 records · Page 25

Forecast of the Amplitude of Solar Cycle 24 Based on the Disturbed Days Precursor

R. J. Thompson (1993, Solar Physics 148, 383) exhibited a significant linear relationship between the number of geomagnetically disturbed days (those that exceed some threshold value of the Ap or nu index) in a solar cycle and the sum of the peak sunspot number in that cycle and the next cycle. Thus, the number of disturbed days during a full cycle, together with the peak sunspot number in that cycle, is a predictor of the amplitude of the next cycle. The work reported here applies Thompson's method to the current cycle. Linear relationships as described above are derived both for the pure aa record (1868-2006) and for a composite of aa (1868-1931) and Ap (1932-2006). For the composite record, the relationship between aa and Ap is determined cycle-by-cycle during the period of overlap. The method is tested for sensitivity to the adopted Ap (or equivalent aa) threshold. The highest smoothed monthly sunspot number for Cycle 24 is forecasted to be R(sub z) = 115 plus or minus 30, where the uncertainty is conservatively based on the full spread of the data around the fitted line in the sunspot number direction. In terms of smoothed monthly 10.7-cm radio flux, the forecast is F10.7 = 164 plus or minus 28.

Rabin, Douglas M.↗

A statistical model of carbon/carbon composite failure

A failure model which considers the stochastic nature of the damage accumulation process is essential to assess reliability and to accurately scale the results from standard test specimens to composite structures. A superior filamentary composite for high temperature applications is composed of carbon fibers in a carbon matrix. Carbon-carbon composites are the strongest known material at very high temperatures. Since there appears to be a significant randomness in C-C material strength which cannot be controlled or detected with current technology, a better model of the material failure based upon statistical principles should be used. Simple applications of the model based upon the limited data provide encouraging results that indicate that better design of test specimens would provide a substantially higher prediction for the design strength of C-C composites. An A-basis strength for the C-C tensile rings from a first stage D-5 billets was estimated. A statistical failure model was developed for these rings which indicates that this strength may be very conservative for larger C-C parts. The analysis may be improved by use of a heterogeneous/noncontinuum finite element approach on the minimechanical level.

Slattery, Kerry T.↗

Linearized form of implicit TVD schemes for the multidimensional Euler and Navier-Stokes equations

Linearized alternating direction implicit (ADI) forms of a class of total variation diminishing (TVD) schemes for the Euler and Navier-Stokes equations have been developed. These schemes are based on the second-order-accurate TVD schemes for hyperbolic conservation laws developed by Harten (1983, 1984). They have the property of not generating spurious oscillations across shocks and contact discontinuities. In general, shocks can be captured within 1-2 grid points. These schemes are relatively simple to understand and easy to implement into a new or existing computer code. One can modify a standard three-point central-difference code by simply changing the conventional numerical dissipation term into the one designed for the TVD scheme. For steady-state applications, the only difference in computation is that the current schemes require a more elaborate dissipation term for the explicit operator; no extra computation is required for the implicit operator. Numerical experiments with the proposed algorithms on a variety of steady-state airfoil problems illustrate the versatility of the schemes.

Yee, H. C.↗

A limitation of the generation of magnetic fields

Simple arguments are put forward to show that the currents produced in disks by radiative or other types of drag forces on electrons must give rise to extraordinarily small magnetic fields. The field strengths are consistent with earlier (more complex) treatments of this class of problem, but inconsistent with the claims of recent papers in the literature. The discrepancies involve the treatment of self-induction. Ion-electron inductive coupling limits the generated magnetic field to have an associated ion Larmor radius greater than or of order of the radius of the disk, a result which follows most directly from conservation of canonical ion momentum, p + eA/c. An explicit time-dependent model for the buildup of the field leads to the same conclusion. If internal velocity gradient scales are smaller than the ion Larmor radius, and the plasma is not collision dominated, standard dynamo amplification within the disk is hardly likely to be effective. But to explain the Galactic field, dynamo amplification in a collisional plasma is also likely to be problematic. The difficulties posed by the existence of ordered Galactic-scale magnetic fields are made all the more acute by the simplicity and the scope of the discussed field limitation.

Balbus, Steven A.↗

Analytical Ion Thruster Discharge Performance Model

A particle and energy balance model of the plasma discharge in magnetic ring-cusp ion thrusters has been developed. The model follows the original work of Brophy in the development of global 0-D discharge models that utilize conservation of particles into and out of the thruster and conservation of energy into the discharge and out of the plasma in the form of charged particles to the walls and beam and plasma radiation. The present model is significantly expanded over Brophy's original work by including self-consistent calculations of the internal neutral pressure, electron temperature, primary electron density, electrostatic ion confinement (due to the ring-cusp fields), plasma potential, discharge stability, and time dependent behavior during recycling. The model only requires information on the thruster geometry, ion optics performance and electrical inputs such as discharge voltage and currents, etc. to produce accurate performance curves of discharge loss versus mass utilization efficiency. The model has been benchmarked against the NEXIS Laboratory Model (LM) and Development Model (DM) thrusters, and successfully predicts the thruster discharge loss as a function of mass utilization efficiency for a variety of thrusters. The discharge performance model will be presented and results showing ion thruster performance and stability given.

discharge chambers↗

The MELISSA pilot plant facility as as integration test-bed for advanced life support systems

The different advances in the Micro Ecological Life Support System Alternative project (MELISSA), fostered and coordinated by the European Space Agency, as well as in other associated technologies, are integrated and demonstrated in the MELISSA Pilot Plant laboratory. During the first period of operation, the definition of the different compartments at an individual basis has been achieved, and the complete facility is being re-designed to face a new period of integration of all these compartments. The final objective is to demonstrate the potentiality of biological systems such as MELISSA as life support systems. The facility will also serve as a test bed to study the robustness and stability of the continuous operation of a complex biological system. This includes testing of the associated instrumentation and control for a safe operation, characterization of the chemical and microbial safety of the system, as well as tracking the genetic stability of the microbial strains used. The new period is envisaged as a contribution to the further development of more complete biological life support systems for long-term manned missions, that should be better defined from the knowledge to be gained from this integration phase. This contribution summarizes the current status of the Pilot Plant and the planned steps for the new period. c2004 COSPAR. Published by Elsevier Ltd. All rights reserved.

Waste Management/instrumentation/methods↗

Extension of a semi-implicit shock-capturing algorithm for 3-D fully coupled, chemically reacting flows in generalized coordinates

A semiimplicit high-resolution shock-capturing method for multidimensional systems of hyperbolic conservation laws with stiff source terms has been developed by Yee and Shinn (1987). The goal of this work is to extend this method to solve the three-dimensional fully coupled Navier-Stokes equations for a hypersonic chemically reacting flow in generalized coordinates. In this formulation, the global continuity equation was replaced by all the species continuity equations. The shock-capturing technique is a second-order-accurate, symmetric total-variation-diminishing method which accounts fully and directly for the coupling among the fluid and all the species. To verify the current approach, it was implemented into an existing computer code which contained the MacCormack method. Test results for a five-species reacting flow are shown to be oscillation-free around the shock, and the time spent per iteration only doubles when compared to the result using classical way of supplying numerical dissipation. The extra computation is more than justified by the elimination of spurious oscillation and nonlinear instability associated with the classical shock-capturing schemes in computing hypersonic reacting flows.

Shinn, J. L.↗

On Formulations of Discontinuous Galerkin and Related Methods for Conservation Laws

A formulation for the discontinuous Galerkin (DG) method that leads to solutions using the differential form of the equation (as opposed to the standard integral form) is presented. The formulation includes (a) a derivative calculation that involves only data within each cell with no data interaction among cells, and (b) for each cell, corrections to this derivative that deal with the jumps in fluxes at the cell boundaries and allow data across cells to interact. The derivative with no interaction is obtained by a projection, but for nodal-type methods, evaluating this derivative by interpolation at the nodal points is more economical. The corrections are derived using the approximate (Dirac) delta functions. The formulation results in a family of schemes: different approximate delta functions give rise to different methods. It is shown that the current formulation is essentially equivalent to the flux reconstruction (FR) formulation. Due to the use of approximate delta functions, an energy stability proof simpler than that of Vincent, Castonguay, and Jameson (2011) for a family of schemes is derived. Accuracy and stability of resulting schemes are discussed via Fourier analyses. Similar to FR, the current formulation provides a unifying framework for high-order methods by recovering the DG, spectral difference (SD), and spectral volume (SV) schemes. It also yields stable, accurate, and economical methods.

Huynh, H. T.↗

A CFD Approach to Modeling Spacecraft Fuel Slosh

Energy dissipation and resonant coupling from sloshing fuel in spacecraft fuel tanks is a problem that occurs in the design of many spacecraft. In the case of a spin stabilized spacecraft, this energy dissipation can cause a growth in the spacecrafts' nutation (wobble) that may lead to disastrous consequences for the mission. Even in non-spinning spacecraft, coupling between the spacecraft or upper stage flight control system and an unanticipated slosh resonance can result in catastrophe. By using a Computational Fluid Dynamics (CFD) solver such as Fluent, a model for this fuel slosh can be created. The accuracy of the model must be tested by comparing its results to an experimental test case. Such a model will allow for the variation of many different parameters such as fluid viscosity and gravitational field, yielding a deeper understanding of spacecraft slosh dynamics. In order to gain a better understanding of the dynamics behind sloshing fluids, the Launch Services Program (LSP) at the NASA Kennedy Space Center (KSC) is interested in finding ways to better model this behavior. Thanks to past research, a state-of-the-art fuel slosh research facility was designed and fabricated at Embry Riddle Aeronautical University (ERAU). This test facility has produced interesting results and a fairly reliable parameter estimation process to predict the necessary values that accurately characterize a mechanical pendulum analog model. The current study at ERAU uses a different approach to model the free surface sloshing of liquid in a spherical tank using Computational Fluid Dynamics (CFD) methods. Using a software package called Fluent, a model was created to simulate the sloshing motion of the propellant. This finite volume program uses a technique called the Volume of Fluid (VOF) method to model the interaction between two fluids [4]. For the case of free surface slosh, the two fluids are the propellant and air. As the fuel sloshes around in the tank, it naturally displaces the air. Using the conservation of mass, momentum, and energy equations, as well as the VOF equations, one can predict the behavior of the sloshing fluid and calculate the forces, pressure gradients, and velocity field for the entire liquid as a function of time.

Marsell, Brandon↗

Active Control Analysis for Aeroelastic Instabilities in Turbomachines

Turbomachines onboard aircraft operate in a highly complex and harsh environment. The unsteady flowfield inherent to turbomachines leads to several problems associated with safety, stability, performance and noise. In-flight surge or flutter incidents could be catastrophic and impact the safety and reliability of the aircraft. High-Cycle-Fatigue (HCF), on the other hand, can significantly impact safety, readiness and maintenance costs. To avoid or minimize these problems generally a more conservative design method must be initiated which results in thicker blades and a loss of performance. Actively controlled turbomachines have the potential to reduce or even eliminate the instabilities by impacting the unsteady aerodynamic characteristics. By modifying the unsteady aerodynamics, active control may significantly improve the safety and performance especially at off-design conditions, reduce noise, and increase the range of operation of the turbomachine. Active control can also help improve reliability for mission critical applications such as the Mars Flyer. In recent years, HCF has become one of the major issues concerning the cost of operation for current turbomachines. HCF alone accounts for roughly 30% of maintenance cost for the United States Air-Force. Other instabilities (flutter, surge, rotating-stall, etc.) are generally identified during the design and testing phase. Usually a redesign overcomes these problems, often reducing performance and range of operation, and resulting in an increase in the development cost and time. Despite a redesign, the engines do not have the capabilities or means to cope with in-flight unforeseen vibration, stall, flutter or surge related instabilities. This could require the entire fleet worldwide to be stood down for expensive modifications. These problems can be largely overcome by incorporating active control within the turbomachine and its design. Active control can help in maintaining the integrity of the system in unforeseen events and provide for more aggressive designs to reduce the weight and improve efficiency of the turbomachine. Another area where active control can be useful is in controlling and suppressing rotating stall and surge in compressors, thereby increasing its operating range. Although some of these benefits will be offset by the added cost and weight penalty of the control system, the potential benefits in safety, reliability, performance, and noise characteristics are significant enough to warrant research in the area of active control of turbomachines. There is renewed interest within industry to understand unsteady aerodynamic behavior. This improved understanding not only leads to better design of turbomachines, which avoid instabilities but also which helps in understanding the controllability of the instabilities. The proliferation of micro-electro-mechanical-system (MEMS) devices has made available new tools to designers for employing feedback controls at reasonable costs. MEMS have also made the control devices small and unobtrusive enough to be implemented within the turbomachine without significant obstruction to the flow path. This has made active-control very attractive especially for systems requiring extreme confidence.

Keith, Theo G., Jr.↗

Current concepts of space flight induced changes in hormonal control of fluid and electrolyte metabolism

A systematic analysis of body fluid and renal dynamics during simulated space flight (head-down bedrest) was undertaken to increase understanding of the physiologic effects of acute cephalad fluid shifts. The earliest effects were increases in central venous pressure and decreases in plasma aldosterone, epinephrine and norepinephrine and glomerular filtration rate, 2 h after the beginning of bedrest. Decreases in plasma angiotensin I at 6 h may have resulted from the increased effective pressure and decreased sympathetic activity seen earlier in bedrest. The early decrease in aldosterone and ADH is thought to contribute to an increase, by 6 h, in urinary excretion of salt and water. Fluid and electrolyte losses occur during space flight, and analysis of body fluids from Space Shuttle crewmembers has indicated that conservation of these substances is begun almost immediately upon cessation of weightlessness. Operational medicine measures to counteract dehydration and electrolyte loss resulted in a less extreme physiologic response to the flight.

Leach, C. S.↗

The method of space-time conservation element and solution element-applications to one-dimensional and two-dimensional time-marching flow problems

A nontraditional numerical method for solving conservation laws is being developed. The new method is designed from a physicist's perspective, i.e., its development is based more on physics than numerics. Even though it uses only the simplest approximation techniques, a 2D time-marching Euler solver developed recently using the new method is capable of generating nearly perfect solutions for a 2D shock reflection problem used by Helen Yee and others. Moreover, a recent application of this solver to computational aeroacoustics (CAA) problems reveals that: (1) accuracy of its results is comparable to that of a 6th order compact difference scheme even though nominally the current solver is only of 2nd-order accuracy; (2) generally, the non-reflecting boundary condition can be implemented in a simple way without involving characteristic variables; and (3) most importantly, the current solver is capable of handling both continuous and discontinuous flows very well and thus provides a unique numerical tool for solving those flow problems where the interactions between sound waves and shocks are important, such as the noise field around a supersonic over- or under-expansion jet.

Chang, Sin-Chung↗

Origin of the Greenland Ice Sheet

A climate model resolving the seasonal cycle and the two horizontal dimensions was developed and applied to several problems of current interest. Models of this type are useful when for various reasons a general circulation model experiment is not warranted or not feasible. For example, in cases where the signal to natural variability is small it may be advantageous to first consider such a statistical dynamical model because extremely long runs may be necessary in the application. In this case the simpler statistical dynamical model serves as a pilot study device. The model developed is a thermodynamic model whose solution yields the equilibrium seasonal cycle for the surface temperature field over the globe. The model is essentially a statement of the conservation of heat energy for individual columns of the earth atmosphere system. Various terms such as the infrared radiation flux to space are parameterized with earth radiation budget data from satellites such as Nimbus 6. The primary agent modulating the seasonal cycle amplitude is the heat capacity per unit area which is a strong function of surface type -- ocean surface can store 60 times more heat per unit time than land. By adjusting its few empirical parameters the model can be brought into remarkable agreement with the observed seasonal cycle. The model is then very useful for looking at the dependence of the seasonal cycle of the temperature on such externally defined variables as the Earth's orbital elements (eccentricity, tilt, precession of equinoxes) or the configuration of land-sea geography which can be changed by continental drift.

North, G. R.↗

A survey on the structured singular value

The structured singular value, U, is an important linear algebra tool to study a class of matrix perturbation problems. It is useful for analyzing the robustness of stability and performance of uncertain, (nominally) linear systems. Computation of (M) is difficult, and usually, upper and lower bounds are all that can be reliably computed. Upper bounds give conservative estimates of the sizes of allowable perturbations. The maximum singular value of a matrix M is an upper bound for (M). As an upper bound, it can be improved by finding a transformations to the data (i.e. M) which do not change the structured singular value, but do reduce the maximum singular value. Typically, upper bound algorithms involve searches over sets of transformations to yield the tightest bound. Lower bound algorithms are intelligent searches for minimum-norm solutions to multivariable polynomial equations, and are based on various optimality conditions that hold at the global (and, unfortunately, some local) minima. The current methods to compute both of these types of bounds are reviewed. Theoretical justification and extensive numerical experience with the various algorithms are covered.

Packard, Andy↗

Double layer propagation in experiments with electron beam injection

Electron beam injection into a plasma is investigated using the analytical inverted Bernstein-Green-Kruskal method. Particle number and momentum conservation laws are applied to evaluate the propagation velocity and potential drop on the leading edge of the beam. Electric potential is supposed to be monotonic, thus the leading front has a double-layer-like structure. For the case of cold particles, analytical expressions for the double layer velocity and potential drop are obtained. It is pointed out that double layer velocity differs from the initial electron speed: even for weak beams a noticeable deceleration takes place. Strong beams are found incapable of penetrating into plasma - their propagation velocity is very small. Ambient electrons undergo a considerable acceleration forming a return current which neutralizes the injector. Possible instability of the distribution functions is discussed.

Bruskin, L. G.↗

The Finer Details: Climate Modeling

If you want to know whether you will need sunscreen or an umbrella for tomorrow's picnic, you can simply read the local weather report. However, if you are calculating the impact of gas combustion on global temperatures, or anticipating next year's rainfall levels to set water conservation policy, you must conduct a more comprehensive investigation. Such complex matters require long-range modeling techniques that predict broad trends in climate development rather than day-to-day details. Climate models are built from equations that calculate the progression of weather-related conditions over time. Based on the laws of physics, climate model equations have been developed to predict a number of environmental factors, for example: 1. Amount of solar radiation that hits the Earth. 2. Varying proportions of gases that make up the air. 3. Temperature at the Earth's surface. 4. Circulation of ocean and wind currents. 5. Development of cloud cover. Numerical modeling of the climate can improve our understanding of both the past and, the future. A model can confirm the accuracy of environmental measurements taken. in, the past and can even fill in gaps in those records. In addition, by quantifying the relationship between different aspects of climate, scientists can estimate how a future change in one aspect may alter the rest of the world. For example, could an increase in the temperature of the Pacific Ocean somehow set off a drought on the other side of the world? A computer simulation could lead to an answer for this and other questions. Quantifying the chaotic, nonlinear activities that shape our climate is no easy matter. You cannot run these simulations on your desktop computer and expect results by the time you have finished checking your morning e-mail. Efficient and accurate climate modeling requires powerful computers that can process billions of mathematical calculations in a single second. The NCCS exists to provide this degree of vast computing capability.

Source record↗

Impeller leakage flow modeling for mechanical vibration control

HPOTP and HPFTP vibration test results have exhibited transient and steady characteristics which may be due to impeller leakage path (ILP) related forces. For example, an axial shift in the rotor could suddenly change the ILP clearances and lengths yielding dynamic coefficient and subsequent vibration changes. ILP models are more complicated than conventional-single component-annular seal models due to their radial flow component (coriolis and centrifugal acceleration), complex geometry (axial/radial clearance coupling), internal boundary (transition) flow conditions between mechanical components along the ILP and longer length, requiring moment as well as force coefficients. Flow coupling between mechanical components results from mass and energy conservation applied at their interfaces. Typical components along the ILP include an inlet seal, curved shroud, and an exit seal, which may be a stepped labyrinth type. Von Pragenau (MSFC) has modeled labyrinth seals as a series of plain annular seals for leakage and dynamic coefficient prediction. These multi-tooth components increase the total number of 'flow coupled' components in the ILP. Childs developed an analysis for an ILP consisting of a single, constant clearance shroud with an exit seal represented by a lumped flow-loss coefficient. This same geometry was later extended to include compressible flow. The objective of the current work is to: supply ILP leakage-force impedance-dynamic coefficient modeling software to MSFC engineers, base on incompressible/compressible bulk flow theory; design the software to model a generic geometry ILP described by a series of components lying along an arbitrarily directed path; validate the software by comparison to available test data, CFD and bulk models; and develop a hybrid CFD-bulk flow model of an ILP to improve modeling accuracy within practical run time constraints.

Palazzolo, Alan B.↗

Resolution of the buoyancy in the 8-foot high temperature tunnel combustor

Currently, the 8-Foot High Temperature Tunnel (8-Ft. HTT) combustor produces a good profile at only one point (2000 psia and 3650 R with oxygen enrichment). Air is enriched with oxygen (liquid) so that the combustor product gas will contain the volumetric amount of oxygen normally found in air. The oxygen enriched air has a large fraction that is not reacted and flows through the outer periphery of the fuel injector. This ring of cold air in addition to the relatively cold walls of the combustor set up buoyancy forces that produce a segregation of relatively cool gases at the bottom of the combustor exit. The basic problem is to produce a test gas that has uniform properties at all combustor conditions. The combustor temperature may be as high as 3700 R or as low as 2000 R. Combustor pressures can be as high as 3500 psia (no oxygen enrichment) and as low as 600 psia. The segregation is most severe with oxygen enriched air, since its temperature is lower and its density is high. The combustor is lined with nickel 201 and can be operated at about 1600 R maximum. A global mixing process is desired that produces an acceptable profile of temperature, species, and velocity at the exit of the combustor. The ultimate goal is a temperature profile with about 100 R variance and about 2 percent variance in oxygen. The exit total temperature must not be lowered significantly by the mixing apparatus or mechanisms employed. If immersed bodies are used, they must also be kept very hot. All combustor wall modifications must be able to survive the heat and structural conditions of the varied operating conditions. Our approach to resolving this issue is being conducted in three stages: (1) Consider mixing exclusively, (2) Resolve the heat transfer concerns resulting from the chosen mixing strategy, and (3) Solve the material and structural problems resulting from stages (1) and (2). Since the 8-Ft. HTT is unavailable for experimentation, the study is conducted exclusively with computational fluid dynamic (CPD) codes (Fluent/Uns and Rampant 3.1) using unstructured grid through body fitted coordinates. Both CFD codes are general purpose Navier-Stokes solution packages that can solve integral conservation equations for conservation of mass, momentum and energy. The governing equations are discretized using a control-volume finite-element method on unstructured triangular 2-D grids. In the interest of time, a 3-D tetrahedral grid was used to check the 2-D results on one mixing strategy and the 2-D results were confirmed. Preliminary results indicate that excellent mixing can be achieved with a body placed in the center of the flowing hot fluid with a minor modification to the combustor wall similar to a model positioned in the test section of a wind tunnel. The concept here, is to create longitudinal vortices strong enough to bring both fluids into intimate contact with each other near the wall where the cold fluid resides. However, there may be a trade-off in the choice of mixing strategy, heat transfer and structural requirements. Currently, we are examining the more promising geometries for heat transfer concerns and developing strategies for the material of construction of the center body.

Loney, Norman W.↗