Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Saturation problem”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

A multidisciplinary approach to optimization of controlled space structures

A fundamental problem facing controls-structures analysts is a means of determining the trade-offs between structural design parameters and control design parameters in meeting some particular performance criteria. Developing a general optimization-based design methodology integrating the disciplines of structural dynamics and controls is a logical approach. The objective of this study is to develop such a method. Classical design methodology involves three phases. The first is structural optimization, wherein structural member sizes are varied to minimize structural mass, subject to open-loop frequency constraints. The next phase integrates control and structure design with control gains as additional design variables. The final phase is analysis of the 'optimal' integrated design phase considering 'real' actuators and 'standard' member sizes. The control gains could be further optimized for fixed structure, and actuator saturation constraints could be imposed. However, such an approach does not take full advantage of opportunities to tailor the structure and control system design as one system.

Woodard, Stanley E.↗
NASA NTRS2001

Preliminary Study of Electron Emission for Use in the PIC Portion of MAFIA

This memorandum summarizes a study undertaken to apply the program MAFIA to the modeling of an electron gun in a traveling wave tube (TWT). The basic problem is to emit particles from the cathode in the proper manner. The electrons are emitted with the classical Maxwell-Boltzmann (M-B) energy distribution; and for a small patch of emitting surface; the distribution with angle obeys Lambert's law. This states that the current density drops off as the cosine of the angle from the normal. The motivation for the work is to extend the analysis beyond that which has been done using older codes. Some existing programs use the Child-Langmuir, or 3/2 power law, for the description of the gun. This means the current varies as the 3/2 power of the anode voltage. The proportionality constant is termed the perveance of the gun. This is limited, however, since the 3/2 variation is only an approximation. Also, if the cathode is near saturation, the 3/2 law definitely will not hold. In most of the older codes, the electron beam is decomposed into current tubes, which imply laminar flow in the beam; even though experiments show the flow to be turbulent. Also, the proper inclusion of noise in the beam is not possible. These older methods of calculation do, however, give reasonable values for parameters of the electron beam and the overall gun, and these values will be used as the starting point for a more precise particle-in-cell (PIC) calculation. To minimize the time needed for a given computer run, all beams will use the same number of particles in a simulation. This is accomplished by varying the mass and charge of the emitted particles (macroparticles) in a certain manner, to be consistent with the desired beam current.

Freeman, Jon C.↗
NASA NTRS1989

Kinematics, controls, and path planning results for a redundant manipulator

The inverse kinematics solution, a modal position control algorithm, and path planning results for a 7 degree of freedom manipulator are presented. The redundant arm consists of two links with shoulder and elbow joints and a spherical wrist. The inverse kinematics problem for tip position is solved and the redundant joint is identified. It is also shown that a locus of tip positions exists in which there are kinematic limitations on self-motion. A computationally simple modal position control algorithm has been developed which guarantees a nearly constant closed-loop dynamic response throughout the workspace. If all closed-loop poles are assigned to the same location, the algorithm can be implemented with very little computation. To further reduce the required computation, the modal gains are updated only at discrete time intervals. Criteria are developed for the frequency of these updates. For commanding manipulator movements, a 5th-order spline which minimizes jerk provides a smooth tip-space path. Schemes for deriving a corresponding joint-space trajectory are discussed. Modifying the trajectory to avoid joint torque saturation when a tip payload is added is also considered. Simulation results are presented.

Gretz, Bruce↗
NASA NTRS1993

Robust nonlinear control of vectored thrust aircraft

An interdisciplinary program in robust control for nonlinear systems with applications to a variety of engineering problems is outlined. Major emphasis will be placed on flight control, with both experimental and analytical studies. This program builds on recent new results in control theory for stability, stabilization, robust stability, robust performance, synthesis, and model reduction in a unified framework using Linear Fractional Transformations (LFT's), Linear Matrix Inequalities (LMI's), and the structured singular value micron. Most of these new advances have been accomplished by the Caltech controls group independently or in collaboration with researchers in other institutions. These recent results offer a new and remarkably unified framework for all aspects of robust control, but what is particularly important for this program is that they also have important implications for system identification and control of nonlinear systems. This combines well with Caltech's expertise in nonlinear control theory, both in geometric methods and methods for systems with constraints and saturations.

Doyle, John C.↗
NASA NTRS1999

Computer Modeling of Sand Transport on Mars Using a Compart-Mentalized Fluids Algorithm (CFA)

It has been postulated that aeolian transport on Mars may be significantly different from that on Earth. From laboratory experiments simulating martian grain transport [2], it has been observed that (saltating) grains striking the bed can cause hundreds of secondary reptation trajectories when impact occurs at speeds postulated for Mars. Some of the ballistically induced trajectories "die ouf' and effectively join the ranks on the creep population that is merely nudged along by impact. Many of the induced reptation trajectories, however, are sufficiently high for the grains to become part of the saltation load (it is irrelevant to the boundary layer how a grain attained its initial lift force). When these grains, in turn, strike the surface, they too are capable of inducing more reptating grains. This cascading effect has been discussed in connection with terrestrial aeolian transport in an attempt to dispel the notion that sand motion is divisible only into creep and saltation loads. On Earth, only a few grains are splashed by impact. On Mars, it may be hundreds. We developed a computer model to address this phenomenon because there are some important ramifications: First, this ratio may mean that martian aeolian transport is dominated by reptation flux rather than saltation. On Earth, the flux would be a roughly balanced mixture between reptation/creep and saltation. On Venus, there would be no transport other than by saltation. In other words, an understanding of planetary aeolian processes may not be necessarily understood by extrapolating from the "Earth case", with only gravity and atmospheric density/viscosity being considered as variables. Second, the reptation flux on Mars may be self sustaining, so that little input is required by the wind once transport has been initiated. The number of grains saturating the boundary layer near the bed may mean that average grain speed on Mars might conceivably be less than that on Earth. This would say much for models of sand comminution on Mars. A multiple-grain transport model using just the equations of grain motion describing lift and drag is impossible to develop owing to stochastic effects --the very effects we wish to model. Also, unless we were to employ supercomputing techniques and extremely complex computer codes that could deal with millions of grains simultaneously, it would also be difficult to model grain transport if we attempted to consider every grain in motion. No existing computer models were found that satisfactorily used the equations of motion to arrive at transport flux numbers for the different populations of saltation and reptation. Modeling all the grains in a transport system was an intractable problem within our resources, and thus we developed what we believe to be a new modeling approach to simulating grain transport. The CFA deals with grain populations, but considers them to belong to various compartmentalized fluid units in the boundary layer. In this way, the model circumvents the multigrain problem by dealing primarily with the consequences of grain transport --momentum transfer between air and grains, which is the physical essence of a dynamic grain-fluid mixture. We thus chose to model the aeolian transport process as a superposition of fluids. These fluids include the air as well as particle populations of various properties. The prime property distinguishing these fluids is upward and downward grain motion. In a normal saltation trajectory, a grain's downwind velocity increases with time, so a rising grain will have a smaller downwind velocity than a failing grain. Because of this disparity in rising and falling grain proper-ties, it seemed appropriate to track these as two separate grain populations within the same physical space. The air itself can be considered a separate fluid superimposed within and interacting with the various grain-cloud "fluids". Additional informaiton is contained in the original.

Marshall, J.↗
NASA NTRS1991

The Generic Spacecraft Analyst Assistant (GenSAA): A tool for automating spacecraft monitoring with expert systems

Flight Operations Analysts (FOAs) in the Payload Operations Control Center (POCC) are responsible for monitoring a satellite's health and safety. As satellites become more complex and data rates increase, FOAs are quickly approaching a level of information saturation. The FOAs in the spacecraft control center for the COBE (Cosmic Background Explorer) satellite are currently using a fault isolation expert system named the Communications Link Expert Assistance Resource (CLEAR), to assist in isolating and correcting communications link faults. Due to the success of CLEAR and several other systems in the control center domain, many other monitoring and fault isolation expert systems will likely be developed to support control center operations during the early 1990s. To facilitate the development of these systems, a project was initiated to develop a domain specific tool, named the Generic Spacecraft Analyst Assistant (GenSAA). GenSAA will enable spacecraft analysts to easily build simple real-time expert systems that perform spacecraft monitoring and fault isolation functions. Lessons learned during the development of several expert systems at Goddard, thereby establishing the foundation of GenSAA's objectives and offering insights in how problems may be avoided in future project, are described. This is followed by a description of the capabilities, architecture, and usage of GenSAA along with a discussion of its application to future NASA missions.

Hughes, Peter M.↗
NASA NTRS1992

Evolution of the shock front and turbulence structures in the shock/turbulence interaction

The interaction of a weak shock front with isotropic turbulence has been investigated using Direct Numerical Simulation (DNS). Two problems were considered: the ability of the field equation (the equation for a propagating surface) to model the shock; and a quantitative study of the evolution of turbulence structure using the database generated by Lee et al. Field equation model predictions for front shape have been compared with DNS results; good agreement is found for shock wave interaction with 2D turbulence and for a single steady vorticity wave. In the interaction of 3D isotropic turbulence with a normal shock, strong alignment of vorticity with the intermediate eigenvector of the rate of strain tensor (S(sup *)(sub ij) = S(sub ij) - (1/3)(delta(sub ij))(S(sub kk))) is seen to develop upstream of the shock and to be further amplified on passage through the shock. Vorticity tends to align at 90 deg to the largest eigenvector, but there is no preferred alignment with the smallest eigenvector. Upstream of the shock, the alignments continue to develop even after the velocity derivative skewness saturates. There is a significant tendency, which increases with time throughout the computational domain, for velocity to align with vorticity. The alignment between velocity and vorticity is strongest in eddy regions and weakest in convergence regions.

Kevlahan, N.↗
NASA NTRS2010

Comparing Spatial Distributions of Solar Prominence Mass Derived from Coronal Absorption

In the present work we extend the use of this mass-inference technique to a sample of prominences observed in at least two coronal lines. This approach, in theory, allows a direct calculation of prominence mass and helium abundance and how these properties vary spatially and temporally. Our motivation is two-fold: to obtain a He(exp 0)/H(exp 0) abundance ratio, and to determine how the relative spatial distribution of the two species varies in prominences. The first of these relies on the theoretical expectation that the amount of absorption at each EUV wavelength is well-characterized. However, in this work we show that due to a saturation of the continuum absorption in the 625 A and 368 A lines (which have much higher opacity compared to 195 A-) the uncertainties in obtaining the relative abundances are too high to give meaningful estimates. This is an important finding because of its impact on future studies in this area. The comparison of the spatial distribution of helium and hydrogen presented here augments previous observational work indicating that cross-field diffusion of neutrals is an important mechanism for mass loss. Significantly different loss timescales for neutral He and H (helium drains much more rapidly than hydrogen) can impact prominence structure, and both the present and past studies suggest this mechanism is playing a role in structure and possibly dynamics. Section 2 of this paper contains a description of the observations and Section 3 summarizes the method used to infer mass along with the criteria imposed in choosing prominences appropriate for this study. Section 3 also contains a discussion of the problems due to limitations of the available data and the implications for determining relative abundances. We present our results in Section 4, including plots of radial-like scans of prominence mass in different lines to show the spatial distribution of the different species. The last section contains a discussion summarizing the importance of the qualitative results found in this work. The Appendix provides a detailed derivation of how to obtain prominence mass and helium abundance (A 1) and includes the data for all prominences studied (A2).

Gilbert, Holly↗
NASA NTRS2013

NASA Technology Evaluation for Environmental Risk Mitigation Remediation Technology Collaboration Development

NASA is committed to finding solutions to agency cleanup problems that are better, cheaper, and more effective than the status quo. Unfortunately, some potential solutions involve innovative technologies for which NASA remediation managers may not have a high level of understanding or confidence. Since 2004, NASA's Stennis Space Center (SSC) in Mississippi has been pumping groundwater contaminated with trichloroethylene (TCE) and other halogenated volatile organic compounds (HVOC) from their cleanup location designated "Area G" through extraction wells to an aboveground treatment system. Over time, however, the effectiveness of this treatment strategy has diminished and an alternative approach is needed. In 2012, professionals from NASA's Principal Center for Technology Evaluation for Environmental Risk Mitigation (TEERM) introduced SSC managers to an innovative technology for enhancing the performance of SSC's existing pump and treat system. The technology, generally referred to as in situ chemical oxidation (ISCO), involves slowly and continuously injecting a strong but safe chemical oxidant into the groundwater. Treatment is enhanced by a "surfactant-type effect" which causes residual contamination from saturated soil to be released into the dissolved-phase where it can be readily oxidized. Any dissolved-phase contamination that was not oxidized can be collected by the extraction well network and treated aboveground. SSC was not familiar with the technology so to increase their confidence, TEERM identified a contractor who was willing to demonstrate their product and process at a significantly reduced price. An initial, small-scale demonstration of ISCO began at sse in March 2012 and completed in August 2012. This successful demonstration was followed by three larger-scale ISCO demonstrations between August and December 2012. The contractor's innovative Continuous Injection System (CIS) incorporated "green" and sustainable technologies and practices. A slow injection rate was maintained autonomously by the CIS, eliminating the need for multiple mobilizations of personnel and powered equipment. The CIS was calibrated to deliver only as much reagent as the formation would accept without "short circuiting", minimizing material waste. Public utility water pressure was used to mix, dilute, and inject the reagent. NASA personnel were trained to operate and maintain the system and remote monitoring and injection control capabilities were developed, further reducing the need for contractor mobilizations to the site. An integrated solar photovoltaic panel was used to power the control valves and web monitoring telemetry. SSC provided hydrogen peroxide for the process using their existing supplier for the aboveground pump and treat systems. The only consumable required from the contractor was their proprietary activator to catalyze SSC's peroxide.

Romeo, James↗
NASA NTRS2011

Validating Phasing and Geometry of Large Focal Plane Arrays

The Kepler Mission is designed to survey our region of the Milky Way galaxy to discover hundreds of Earth-sized and smaller planets in or near the habitable zone. The Kepler photometer is an array of 42 CCDs (charge-coupled devices) in the focal plane of a 95-cm Schmidt camera onboard the Kepler spacecraft. Each 50x25-mm CCD has 2,200 x 1,024 pixels. The CCDs accumulate photons and are read out every six seconds to prevent saturation. The data is integrated for 30 minutes, and then the pixel data is transferred to onboard storage. The data is subsequently encoded and transmitted to the ground. During End-to-End Information System (EEIS) testing of the Kepler Mission System (KMS), there was a need to verify that the pixels requested by the science team operationally were correctly collected, encoded, compressed, stored, and transmitted by the FS, and subsequently received, decoded, uncompressed, and displayed by the Ground Segment (GS) without the outputs of any CCD modules being flipped, mirrored, or otherwise corrupted during the extensive FS and GS processing. This would normally be done by projecting an image on the focal plane array (FPA), collecting the data in a flight-like way, and making a comparison between the original data and the data reconstructed by the science data system. Projecting a focused image onto the FPA through the telescope would normally involve using a collimator suspended over the telescope opening. There were several problems with this approach: the collimation equipment is elaborate and expensive; as conceived, it could only illuminate a limited section of the FPA (.25 percent) during a given test; the telescope cover would have to be deployed during testing to allow the image to be projected into the telescope; the equipment was bulky and difficult to situate in temperature-controlled environments; and given all the above, test setup, execution, and repeatability were significant concerns. Instead of using this complicated approach of projecting an optical image on the FPA, the Kepler project developed a method using known defect features in the CCDs to verify proper collection and reassembly of the pixels, thereby avoiding the costs and risks of the optical projection approach. The CCDs composing the Kepler FPA, as all CCDs, had minor defects. At ambient temperature, some pixels look far brighter than they should. These ghot h pixels have a higher rate of charge leakage than the others due to manufacturing variations. They are usually stable over time, and appear at temperatures above 5 oC. The hot pixels on the Kepler FPA were mapped before photometer assembly during module testing. Selected hot pixels were used as target gstars h for the purposes of EEIS testing. gDead h pixels are permanently off, producing a permanently black pixel. These can also be used if there is some illumination of the FPA. During EEIS testing, Dark Current Full Frame Images (FFIs) taken at room temperature were used to create the hot pixel maps for all 84 Kepler photometer CCD channels. Data from two separate nights were used to create two hot pixel maps per channel, which were cross-correlated to remove cosmic ray events which appear to be hot pixels. These hot pixel maps obtained during EEIS testing were compared to the maps made during module testing to verify that the end-to-end data flow was correct.

Standley, Shaun P.↗
NASA NTRS1993

A Report on Stochastic Fairness Queueing (SFQ) Experiments

SRI International (SRI) has developed an improved queueing algorithm, known as Stochastic Fairness Queueing (SFQ), for best-effort traffic (i.e., traffic that does not require any guaranteed service). SFQ is a probablistic variant of strict fair queueing where instead of a single queue being allocated per flow, a fixed number of queues are used and a hash function maps the IP source and destination to a particular queue. A seed to the hash function is also perturbed occasionally to help distribute the flows amongst different queues when more than one flow maps to the same queue during the lifetime of the flow. SFQ provides 'fair' access by trying to ensure that each flow from source to destination host obtains equal access to the available bandwidth. This report covers a series of experiments performed on DARTnet evaluating the behavior and performance of SFQ against a FIFO queueing discipline. These experiments were designed to show SFQ's advantages and performance, and include tests demonstrating: Fair utilization of available resources; Starvation prevention; Graceful degradation under overload conditions; and Resource usage. In general, the experiments do show that SFQ is better than FIFO queueing at allocating bandwidth equally among a set of flows. SFQ also prevents a stream from dominating the available bandwidth, which seems to be a tendency with FIFO queueing (i.e., if a flow demands more than its share of the available bandwidth, with FIFO queueing that stream receives a disproportionate amount when compared to flows demanding less than their share). Furthermore, SFQ seems to reward 'nice' users of the network by providing a lower variance in delay and more throughput when their resource demand is less than their available share. Both SFQ and FIFO queueing seem to degrade fairly well as the network becomes saturated and to recover well as the network becomes less congested. Not unexpectedly, FIFO queueing is a little more efficient than SFQ-the delays are less and the throughput slightly higher because SFQ requires more processing. However, the performance difference between the two queueing disciplines is relatively small. However, the experiments do point out some interesting behavior. FIFO queueing can behave better than SFQ with seed perturbation. We recommend further evaluation of the hash function and the seed perturbation technique. There are probably weaknesses in their current selection that cause this unexpected behavior. SFQ also seems to possess good scaling properties. To verify this, more experiments with a larger number of streams from more hosts need to be executed and examined, including the staggered introduction of streams. Staggering the streams may prove important, because graphs in the degradation experiment revealed some unexpected increases and decreases in throughput, which should be examined. This may again be due to the interaction of the hash function with the seed perturbation but it may also be related to some other unknown problem.

Denny, Barbara A.↗
NASA NTRS2022

Coronal Polarimetry: Determining the Magnetic Origins of Coronal Mass Ejections

Synopsis: The mechanism for the release of stored magnetic energy in solar eruptions remains a major unsolved problem of Heliophysics. Choosing between triggers requires knowledge of pre-eruptive magnetic fields (B). Although linear polarization in visible/infrared (VIR) coronal emission lines reveals intriguing clues about coronal mass ejection (CME) precursor topology, small telescope apertures limit current capability for measuring vector field. Current coronal observations are insufficient to diagnose 3D coronal magnetic fields in CME precursors. Coronal cavities are the ideal candidates for CME precursor studies. BLOS in cavities is a direct measure of stored magnetic free energy, and the presence and location of topological X-points (reconnection locations) and O-points (circulation of B about axis) distinguish predictions of the flux rope-torus-instability and sheared-arcade-breakout CME models. - Finding: Coronal cavities are ubiquitous throughout the solar cycle, erupt as CMEs, and have coronal polarimetric signatures that distinguish between model predictions. - Recommendation: Coronal cavities should be targeted in a comprehensive multiwavelength study of CME precursors & associated eruption-driving mechanisms. Large ground-based VIR telescopes can measure both coronal magnetic field topology and strength through a combination of the saturated Hanle and BLOS -sensitive Zeeman effects. - Finding: The 4m Daniel K. Inouye Solar Telescope (DKIST) and proposed 1.5m Coronal Solar Magnetism Observatory Large Coronagraph (COSMO-LC) make a major leap forward in VIR coronal sensitivity, enabling measurement of vector B. - Recommendation: A dedicated coronal synoptic telescope with a large field-of-view (i.e., COSMO-LC) is needed to identify the dominant CME driving mechanisms, by measuring coronal magnetic fields from precursor state into eruption. In the UV, a new opportunity has arisen to make use of the strong H I Lyman-α coronal line, which obtains a measurement of B that is independent of the Zeeman effect. - Finding: Small-telescope spectropolarimetric capability in the unsaturated Hanle regime provides an independent and complementary coronal magnetic diagnostic to large ground-based telescope measurements. - Recommendation: The proposed 12-cm Coronal Lyman-α Resonance Observatory (CLARO) spectropolarimetric coronagraph demonstrates a path forward for space-based observations of the coronal magnetic field that should be incorporated into future missions away from the Sun-Earth line (e.g., COMPLETE).

magnetic energy↗
NASA NTRS2015

Use of Generalized Fluid System Simulation Program (GFSSP) for Teaching and Performing Senior Design Projects at the Educational Institutions

This paper describes the experience of the authors in using the Generalized Fluid System Simulation Program (GFSSP) in teaching Design of Thermal Systems class at University of Alabama in Huntsville. GFSSP is a finite volume based thermo-fluid system network analysis code, developed at NASA/Marshall Space Flight Center, and is extensively used in NASA, Department of Defense, and aerospace industries for propulsion system design, analysis, and performance evaluation. The educational version of GFSSP is freely available to all US higher education institutions. The main purpose of the paper is to illustrate the utilization of this user-friendly code for the thermal systems design and fluid engineering courses and to encourage the instructors to utilize the code for the class assignments as well as senior design projects. The need for a generalized computer program for thermofluid analysis in a flow network has been felt for a long time in aerospace industries. Designers of thermofluid systems often need to know pressures, temperatures, flow rates, concentrations, and heat transfer rates at different parts of a flow circuit for steady state or transient conditions. Such applications occur in propulsion systems for tank pressurization, internal flow analysis of rocket engine turbopumps, chilldown of cryogenic tanks and transfer lines, and many other applications of gas-liquid systems involving fluid transients and conjugate heat and mass transfer. Computer resource requirements to perform time-dependent, three-dimensional Navier-Stokes computational fluid dynamic (CFD) analysis of such systems are prohibitive and therefore are not practical. Available commercial codes are generally suitable for steady state, single-phase incompressible flow. Because of the proprietary nature of such codes, it is not possible to extend their capability to satisfy the above-mentioned needs. Therefore, the Generalized Fluid System Simulation Program (GFSSP1) has been developed at NASA Marshall Space Flight Center (MSFC) as a general fluid flow system solver capable of handling phase changes, compressibility, mixture thermodynamics and transient operations. It also includes the capability to model external body forces such as gravity and centrifugal effects in a complex flow network. The objectives of GFSSP development are: a) to develop a robust and efficient numerical algorithm to solve a system of equations describing a flow network containing phase changes, mixing, and rotation; and b) to implement the algorithm in a structured, easy-to-use computer program. The analysis of thermofluid dynamics in a complex network requires resolution of the system into fluid nodes and branches, and solid nodes and conductors as shown in Figure 1. Figure 1 shows a schematic and GFSSP flow circuit of a counter-flow heat exchanger. Hot nitrogen gas is flowing through a pipe, colder nitrogen is flowing counter to the hot stream in the annulus pipe and heat transfer occurs through metal tubes. The problem considered is to calculate flowrates and temperature distributions in both streams. GFSSP has a unique data structure, as shown in Figure 2, that allows constructing all possible arrangements of a flow network with no limit on the number of elements. The elements of a flow network are boundary nodes where pressure and temperature are specified, internal nodes where pressure and temperature are calculated, and branches where flowrates are calculated. For conjugate heat transfer problems, there are three additional elements: solid node, ambient node, and conductor. The solid and fluid nodes are connected with solid-fluid conductors. GFSSP solves the conservation equations of mass and energy, and equation of state in internal nodes to calculate pressure, temperature and resident mass. The momentum conservation equation is solved in branches to calculate flowrate. It also solves for energy conservation equations to calculate temperatures of solid nodes. The equations are coupled and nonlinear; therefore, they are solved by an iterative numerical scheme. GFSSP employs a unique numerical scheme known as simultaneous adjustment with successive substitution (SASS), which is a combination of Newton-Raphson and successive substitution methods. The mass and momentum conservation equations and the equation of state are solved by the Newton-Raphson method while the conservation of energy and species are solved by the successive substitution method. GFSSP is linked with two thermodynamic property programs, GASP2 and WASP3 and GASPAK4, that provide thermodynamic and thermophysical properties of selected fluids. Both programs cover a range of pressure and temperature that allows fluid properties to be evaluated for liquid, liquid-vapor (saturation), and vapor region. GASP and WASP provide properties of 12 fluids. GASPAK includes a library of 36 fluids. GFSSP has three major parts. The first part is the graphical user interface (GUI), visual thermofluid analyzer of systems and components (VTASC). VTASC allows users to create a flow circuit by a 'point and click' paradigm. It creates the GFSSP input file after the completion of the model building process. GFSSP's GUI provides the users a platform to build and run their models. It also allows post-processing of results. The network flow circuit is first built using three basic elements: boundary node, internal node, and branch.

Majumdar, A. K.↗
NASA NTRS1996

The infrared spectrum of the Galactic center and the composition of interstellar dust

We have obtained 5-8 micrometers spectra of the Galactic center from the Kuiper Airborne Observatory at resolving powers of approximately 50, approximately 150, and approximately 300. These spectra show absorption features at 5.5, 5.8, 6.1, and 6.8 micrometers. Together with previously observed features in the 3 micrometers region, these features are compared with laboratory spectra of candidate materials. The 3.0 and 6.1 micrometers features are due to the OH stretching and bending variations of H2O and are well fitted by water of hydration in silicates (e.g., talc). The 3.0 micrometer band is equally well fitted by ice mixtures containing 30% H2O, but such mixtures do not provide a good fit to the observed 6.1 micrometer band. The 3.4 and 6.8 micrometers features are identified with the CH stretching and deformation modes in CH2 and CH3 groups in saturated aliphatic hydrocarbons. The 6.1 micrometer band shows a short wavelength shoulder centered on 5.8 micrometer, attributed to carbonyl (C double bond O) groups in this interstellar hydrocarbon dust component. Finally, the narrow 5.5 micrometer feature is also attributed to carbonyl groups, but in the form of metal carbonyls [e.g., Fe(CO)4]. We have derived column densities and abundances along the line of sight toward the Galactic center for the various identified dust components. This analysis shows that hydrocarbon grains contain only 0.08 of the elemental abundance of C and contribute only a relatively minor fraction (0.1) of the total dust volume. Most of the interstellar dust volume is made up of silicates (approximately 0.6). Small graphite grains, responsible for the 2200 angstroms bump, account for 0.07 of the total dust volume. The remaining one-quarter of the interstellar dust volume consists of a material(s) without strong IR absorption features. Likely candidates include large graphite grains, diamonds, or amorphous carbon grains, which all have weak or no IR active modes. Finally, various models for the origin of the hydrocarbon dust component of the interstellar dust are discussed. All of them face some problems in explaining the observations, in particular, the absence of the spectroscopic signature of hydrocrbon grains in sources associated with molecular clouds.

NASA Discipline Exobiology↗
NASA NTRS1996

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗
NASA NTRS2016

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗
NASA NTRS1995

Nonequilibrium, large-amplitude MHD fluctuations in the solar wind

Compressible MHD simulations in one dimension with three-dimensional vectors are used to investigate a number of processes relevant to problems in interplanetary physics. The simulations indicate that a large-amplitude nonequilibrium (e.g., linearly polarized) Alfvenic wave, which always starts with small relative fluctuations in the magnitude B of the magnetic field, typically evolves to flatten the magnetic profile in most regions. Under a wide variety of conditions B and the density rho become anticorrelated on average. If the mean magnetic field is allowed to decrease in time, the point where the transverse magnetic fluctuation amplitude delta B(sub T) is greater than the mean field B(sub 0) is not special, and large values of delta B(sub T)/B(sub 0) do not cause the compressive thermal energy to increase remarkably or the wave energy to dissipate at an unusually high rate. Nor does the 'backscatter' of the waves that occurs when the sound speed is less than the Alfven speed result, in itself, in substantial energy dissipation, but rather primarily in a phase change between the magnetic and velocity fields. For isolated wave packets the backscatter does not occur for any of the parameters examined; an initial radiation of acoustic waves away from the packet establishes a stable traveling structure. Thus these simulations, although greatly idealized compared to reality, suggest a picture in which the interplanetary fluctuations should have small deltaB and increasingly quasi-pressure balanced compressive fluctuations, as observed, and in which the dissipation and 'saturation' at delta B(sub T)/B(sub 0) approximately = 1 required by some theories of wave acceleration of the solar wind do not occur. The simulations also provide simple ways to understand the processes of nonlinear steepening and backscattering of Alfven waves and demonstrate the existence of previously unreported types of quasi-steady MHD states.

Roberts, D. Aaron↗
NASA NTRS2017

Integrated Demand Management: Minimizing Unanticipated Excessive Departure Delay while Ensuring Fairness from a Traffic Management Initiative

This paper introduces NASA's Integrated Demand Management (IDM) concept and presents the results from an early proof-of-concept evaluation and an exploratory experiment. An initial development of the concept was focused on integrating two systems - i.e. the FAA's newly deployed Traffic Flow Management System (TFMS) tool called the Collaborative Trajectory Options Program (CTOP) and the Time-Based Flow Management (TBFM) system with Extended Metering (XM) capabilities to manage projected heavy traffic demand into a capacity-constrained airport. A human-in-the-loop (HITL) simulation experiment was conducted to demonstrate the feasibility of the initial development of the concept by adapting it to an arrival traffic problem at Newark Liberty International Airport (EWR) during clear weather conditions. In this study, the CTOP was utilized to strategically plan the arrival traffic demand by controlling take-off times of both short- and long-haul flights (long-hauls specify aircraft outside TBFM regions and short-hauls specify aircraft within TBFM regions) in a way that results in equitable delays among the groups. Such strategic planning allows less airborne delay to occur within TBFM by feeding manageable long-haul traffic demand while reserving sufficient slots in the overhead streams for the short-haul departures. The manageable traffic demand indicates the TBFM scheduler assigns no more airborne delay than its assigned airspace is capable of absorbing. TBFM then uses its time-based metering capabilities to deliver the desirable throughput by tactically rescheduling the TBFM entered long-haul flights and short-haul departures. Additional research was also performed to explore use of Required Time of Arrival (RTA) capabilities as a potential control mechanism for the airborne flights to improve arrival traffic delivery accuracy of scheduled long-haul traffic demand. The study results show that both short- and long-haul flights received similar ground delays. In addition, there was a noticeable reduction in the total amount of excessive unanticipated last-minute ground delays, i.e. delays that are frequently imposed on the short-haul flight in current day operations due to saturation in the overhead stream, commonly referred to as 'double penalty'. Furthermore, the concept achieved the target throughput while minimizing the expected cost associated with overall delays in arrival traffic. Assessment of the RTA capabilities showed that there was indeed improvement of the scheduled entry times into TBFM regions by using RTA capabilities. However, with respect to reduction in delays incurred within TBFM, there was no observable benefit of improving the precision of long-haul flights entry times.

Collaborative Trajectory Options Program (CTOP)↗