Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative systems”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

NASA Tech Briefs, March 2011

Topics covered include: Optimal Tuner Selection for Kalman-Filter-Based Aircraft Engine Performance Estimation; Airborne Radar Interferometric Repeat-Pass Processing; Plug-and-Play Environmental Monitoring Spacecraft Subsystem; Power-Combined GaN Amplifier with 2.28-W Output Power at 87 GHz; Wallops Ship Surveillance System; Source Lines Counter (SLiC) Version 4.0; Guidance, Navigation, and Control Program; Single-Frame Terrain Mapping Software for Robotic Vehicles; Auto Draw from Excel Input Files; Observation Scheduling System; CFDP for Interplanetary Overlay Network; X-Windows Widget for Image Display; Binary-Signal Recovery; Volumetric 3D Display System with Static Screen; MMIC Replacement for Gunn Diode Oscillators; Feature Acquisition with Imbalanced Training Data; Mount Protects Thin-Walled Glass or Ceramic Tubes from Large Thermal and Vibration Loads; Carbon Nanotube-Based Structural Health Monitoring Sensors; Wireless Inductive Power Device Suppresses Blade Vibrations; Safe, Advanced, Adaptable Isolation System Eliminates the Need for Critical Lifts; Anti-Rotation Device Releasable by Insertion of a Tool; A Magnetically Coupled Cryogenic Pump; Single Piezo-Actuator Rotary-Hammering Drill; Fire-Retardant Polymeric Additives; Catalytic Generation of Lift Gases for Balloons; Ionic Liquids to Replace Hydrazine; Variable Emittance Electrochromics Using Ionic Electrolytes and Low Solar Absorptance Coatings; Spacecraft Radiator Freeze Protection Using a Regenerative Heat Exchanger; Multi-Mission Power Analysis Tool; Correction for Self-Heating When Using Thermometers as Heaters in Precision Control Applications; Gravitational Wave Detection with Single-Laser Atom Interferometers; Titanium Alloy Strong Back for IXO Mirror Segments; Improved Ambient Pressure Pyroelectric Ion Source; Multi-Modal Image Registration and Matching for Localization of a Balloon on Titan; Entanglement in Quantum-Classical Hybrid; Algorithm for Autonomous Landing; Quantum-Classical Hybrid for Information Processing; Small-Scale Dissipation in Binary-Species Transitional Mixing Layers; Superpixel-Augmented Endmember Detection for Hyperspectral Images; Coding for Parallel Links to Maximize the Expected Value of Decodable Messages; and Microwave Tissue Soldering for Immediate Wound Closure.

Source record↗
NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗
NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗
NASA NTRS2005

Thermal Assessment of Swift BAT Instrument Thermal Control System in Flight

THE BAT is the primary instrument on the Swift spacecraft. The Swift mission is part of the National Aeronautics and Space Administration (NASA) Medium-Size Explorer (MIDEX) Program, and is managed by Goddard Space Flight Center (GSFC). It is designed to detect gamma ray burst over a broad region of the sky in a low Earth orbit of 600-km altitude and quickly align the telescopes on the spacecraft to the gamma ray source. It was successfully launched into orbit on November 20, 2004. The Swift mission is a first of its kind of multi-wavelength transient observatory for gamma ray burst astronomy. Its mission life is 2 years. The inclination is 22 deg maximum. The spacecraft bus voltage to the instruments is in the 24 V to 35 V range. The instruments will be turned off when the voltage is below 27 V. The BAT is mounted to the optical bench through five titanium flexures. The BAT has been developed at GSFC. Its telescope assembly consists of 256 Detector Modules (DMs) in the Detector Array. There are 16 Detector Array Blocks. Each Block holds 16 DMs, 3 Block Voltage Regulator (BVR) units and 3 Block Command & Data Handling (BCDH) units. The power dissipation of each Block has been measured to be 13 W. Therefore the total power dissipation of the 16 Blocks is 208 W. The DAP is 1.3 m (4.3 ft) x 1 m (3.3 ft), accommodates all the 16 Blocks. It also provides the mounting surface and the positional stability for the Blocks. The DMs are located at the top (+X side) of the DAP and is enclosed by graded-Z shields on the sides and a coded mask at the top. The BVRs and BCDHs are located at the bottom (-X side) of the DAP. Eight Blocks are located at the front (-Z side or radiator side) of the DAP, and eight are located at the rear (+Z side) of the DAP. The DMs and top of DAP are insulated with a 7-layer multi-layer insulation (MLI). There is a 5.08 cm (2 in) x 5.08 cm (2 in) MLI cutout over each Block heater controller so that heat radiates from the heater controller to the mask. The exterior of the mask, graded-Z shields and bottom of DAP is insulated with a 15-layer MLI.

Choi, Michael K.↗
NASA NTRS2018

Orientation Effects in Two-Phase Microgap Flow

The high power density of emerging electronic devices is driving the transition from remote cooling, which relies onconduction and spreading, to embedded cooling, which extracts dissipated heat on-site. Two-phase microgap coolersemploy the forced flow of dielectric fluids undergoing phase change in a heated channel within or between devices. Such coolers must work reliably in all orientations for a variety of applications (e.g., vehicle-based equipment), as well as in microgravity and high-g for other applications (e.g., spacecraft and aircraft). The lack of acceptable models andcorrelations for orientation- and gravity-independent operation has limited the use of two-phase coolers in suchapplications. Previous research has revealed that gravitational acceleration plays a diminishing role in establishing flow regimes and transport rates as the channel size shrinks, but there is considerable variation among the proposed microscale criteria and limited research on two-phase flows in low aspect ratio microgap channels. Reliable criteria for achieving orientation- and gravity-independent flow boiling would enable emerging systems to exploit this thermal management technique and streamline the technology development process.

orientation↗
NASA NTRS1987

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.↗
NASA NTRS1990

Numerical solution of the incompressible Navier-Stokes equations

The current work is initiated in an effort to obtain an efficient, accurate, and robust algorithm for the numerical solution of the incompressible Navier-Stokes equations in two- and three-dimensional generalized curvilinear coordinates for both steady-state and time-dependent flow problems. This is accomplished with the use of the method of artificial compressibility and a high-order flux-difference splitting technique for the differencing of the convective terms. Time accuracy is obtained in the numerical solutions by subiterating the equations in psuedo-time for each physical time step. The system of equations is solved with a line-relaxation scheme which allows the use of very large pseudo-time steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. Numerous laminar test flow problems are computed and presented with a comparison against analytically known solutions or experimental results. These include the flow in a driven cavity, the flow over a backward-facing step, the steady and unsteady flow over a circular cylinder, flow over an oscillating plate, flow through a one-dimensional inviscid channel with oscillating back pressure, the steady-state flow through a square duct with a 90 degree bend, and the flow through an artificial heart configuration with moving boundaries. An adequate comparison with the analytical or experimental results is obtained in all cases. Numerical comparisons of the upwind differencing with central differencing plus artificial dissipation indicates that the upwind differencing provides a much more robust algorithm, which requires significantly less computing time. The time-dependent problems require on the order of 10 to 20 subiterations, indicating that the elliptical nature of the problem does require a substantial amount of computing effort.

Rogers, Stuart E.↗
NASA NTRS2010

Reply [to: Atlantic Tropical Cyclogenetic Processes during SOP-3 NAMMA in the GEOS-5 Global Data Assimilation and Forecast System]

This article is a Reply to a Comment by Scott Braun on a previously published article by O. Reale, K.-M. Lau, and E. Brin: "Atlantic tropical cyclogenetic processes during SOP-3 NAMMA in the GEOS-5 global data assimilation and forecast system", by Reale, Lau and Brin, hereafter referred to as RA09. RA09 investigated the role of the Saharan Air Layer (SAL) in tropical cyclogenetic processes associated with a non-developing easterly wave observed during the Special Observation Period (SOP-3) phase of the 2006 NASA African Monsoon Multidisciplinary Analyses (MAMMA). The wave was chosen because both interact heavily with Saharan air. Results showed: a) very steep moisture gradients are associated with the SAL in forecasts and analyses even at great distance from the Sahara; b) a thermal dipole (warm above, cool below) in the non-developing case. RA09A suggested that radiative effect of dust may play some role in producing a thermal structure less favorable to cyclogenesis, and also indicated that only global horizontal resolutions on the order of 20-30 kilometers can capture the large-scale transport and the fine thermal structure of the SAL Braun (2010) questions those results attributing the wave dissipation to midlatitude air. The core discussion is on a dry filament preceding the wave, on the presence of dust, and on the origin of the air contained in this dry filament. In the 'Reply', higher resolution analyses than the ones used by Braun, taken at almost coincident times with Aqua and Terra passes, are shown, to emphasize how the channel of dry air associated with W1 is indeed rich in dust. Backtrajectories on a higher resolution grid are also performed, leading to results drastically different from Braun (2010), and in particularly showing that there is a clear contribution of Saharan air. Finally, the 'Reply' presents evidence on that analyses at a horizontal resolution of one degree are inadequate to investigate such feature.

Reale, Oreste↗
NASA NTRS2003

Reliable Transport over SpaceWire for James Webb Space Telescope (JWST) Focal Plane Electronics (FPE) Network

NASA's James Webb Space Telescope (JWST) faces difficult technical and budgetary challenges to overcome before it is scheduled launch in 2010. The Integrated Science Instrument Module (ISIM), shares these challenges. The major challenge addressed in this paper is the data network used to collect, process, compresses and store Infrared data. A total of 114 Mbps of raw information must be collected from 19 sources and delivered to the two redundant data processing units across a twenty meter deployed thermally restricted interface. Further data must be transferred to the solid-state recorder and the spacecraft. The JWST detectors are kept at cryogenic temperatures to obtain the sensitivity necessary to measure faint energy sources. The Focal Plane Electronics (FPE) that sample the detector, generate packets from the samples, and transmit these packets to the processing electronics must dissipate little power in order to help keep the detectors at these cold temperatures. Separating the low powered front-end electronics from the higher-powered processing electronics, and using a simple high-speed protocol to transmit the detector data minimize the power dissipation near the detectors. Low Voltage Differential Signaling (LVDS) drivers were considered an obvious choice for physical layer because of their high speed and low power. The mechanical restriction on the number cables across the thermal interface force the Image packets to be concentrated upon two high-speed links. These links connect the many image packet sources, Focal Plane Electronics (FPE), located near the cryogenic detectors to the processing electronics on the spacecraft structure. From 12 to 10,000 seconds of raw data are processed to make up an image, various algorithms integrate the pixel data Loss of commands to configure the detectors as well as the loss of science data itself may cause inefficiency in the use of the telescope that are unacceptable given the high cost of the observatory. This combination of requirements necessitates a redundant, fault tolerant, high- speed, low mass, low power network with a low Bit error Rate(1E-9- 1E-12). The ISIM systems team performed many studies of the various network architectures that meeting these requirements. The architecture selected uses the Spacewire protocol, with the addition of a new transport and network layer added to implement end-to-end reliable transport. The network and reliable transport mechanism must be implemented in hardware because of the high average information rate and the restriction on the ability of the detectors to buffer data due to power and size restrictions. This network and transport mechanism was designed to be compatible with existing Spacewire links and routers so that existing equipment and designs may be leveraged upon. The transport layer specification is being coordinated with European Space Agency (ESA), Spacewire Working Group and the Consultative Committee for Space Data System (CCSDS) PlK Standard Onboard Interface (SOIF) panel, with the intent of developing a standard for reliable transport for Spacewire. Changes to the protocol presented are likely since negotiations are ongoing with these groups. A block of RTL VHDL that implements a multi-port Spacewire router with an external user interface will be developed and integrated with an existing Spacewire Link design. The external user interface will be the local interface that sources and sinks packets onto and off of the network (Figure 3). The external user interface implements the network and transport layer and handles acknowledgements and re-tries of packets for reliable transport over the network. Because the design is written in RTL, it may be ported to any technology but will initially be targeted to the new Actel Accelerator series (AX) part. Each link will run at 160 Mbps and the power will be about 0.165 Watt per link worst case in the Actel AX.

Rakow, Glenn↗
NASA NTRS2004

Making Precise Resonators for Mesoscale Vibratory Gyroscopes

An alternative approach to the design and fabrication of vibratory gyroscopes is founded on the use of fabrication techniques that yield best results in the mesoscopic size range, which is characterized by overall device dimensions of the order of a centimeter. This approach stands in contradistinction to prior approaches in (1) the macroscopic size range (the size range of conventional design and fabrication, characterized by overall device dimensions of many centimeters) and (2) the microscopic size range [the size range of microelectromechanical systems (MEMS), characterized by overall device dimensions of the order of a millimeter or less]. The mesoscale approach offers some of the advantage of the MEMS approach (sizes and power demands smaller than those of the macroscale approach) and some of the advantage of the macroscale approach (the possibility of achieving relative dimensional precision greater than that of the MEMS approach). Relative dimensional precision is a major issue in the operation of a vibratory gyroscope. The heart of a vibratory gyroscope is a mechanical resonator that is required to have a specified symmetry in a plane orthogonal to the axis about which rotation is to be measured. If the resonator could be perfectly symmetrical, then in the absence of rotation, a free vibration of the resonator could remain fixed along any orientation relative to its housing; that is, the gyroscope could exhibit zero drift. In practice, manufacturing imprecision gives rise to some asymmetry in mass, flexural stiffness or dissipation, resulting in a slight drift or beating motion of an initial vibration pattern that cannot be distinguished from rotation. In the mesoscale approach, one exploits the following concepts: For a given amount of dimensional error generated in manufacturing, the asymmetry and hence the rate-of-rotation drift of the gyroscope can be reduced by increasing the scale. The decrease in asymmetry also reduces coupling of vibrations to the external environment. Mechanical thermal noise and electronic measurement noise and drift can also be reduced by increasing the size of the resonator and its associated sensors.

Yang, Eui-Hyeok↗
NASA NTRS2012

Investigation of Transonic Wake Dynamics for Mechanically Deployable Entry Systems

A numerical investigation of transonic flow around a mechanically deployable entry system being considered for a robotic mission to Venus has been performed, and preliminary results are reported. The flow around a conceptual representation of the vehicle geometry was simulated at discrete points along a ballistic trajectory using Detached Eddy Simulation (DES). The trajectory points selected span the low supersonic to transonic regimes with freestream Mach numbers from 1:5 to 0:8, and freestream Reynolds numbers (based on diameter) between 2:09 x 10(exp 6) and 2:93 x 10(exp 6). Additionally, the Mach 0:8 case was simulated at angles of attack between 0 and 5 . Static aerodynamic coefficients obtained from the data show qualitative agreement with data from 70deg sphere-cone wind tunnel tests performed for the Viking program. Finally, the effect of choices of models and numerical algorithms is addressed by comparing the DES results to those using a Reynolds Averaged Navier-Stokes (RANS) model, as well as to results using a more dissipative numerical scheme.

Stern, Eric↗
NASA NTRS1999

Control of Meridional Flow in Circular Cylinders by a Travelling Axial Magnetic Field

Convective flow in a Bridgman or float zone configuration significantly affects the interface shape and segregation phenomena. While the primary causative factor for this flow is buoyancy induced convection in an enclosed Bridgman melt, the presence of a free surface gives rise to surface tension driven flows in the floating zone processing of melts. It is of interest to curtail these flows in order to realize near quiescent growth conditions that have shown to result in crystals with good longitudinal and radial homogeneity and thereby of better overall quality. While buoyancy effects can be reduced by careful processing in a low gravity (space) environment, the reduction of Marangoni flows due to surface tension variations is not that straight forward. Attempts have been made with some limited success with the use of external fields to affect the melt thermo-fluid behavior. The use of a static magnetic field that reduces convective contamination through the effects of a non-intrusively induced, dissipative Lorentz force in an electrically conducting melt is one such approach. Experiments have shown that axial fields of the order of 5 Tesla can significantly eliminate convection and yield close to diffusion limited crystal growth conditions. The generation and use of such high magnetic fields require substantial hardware and incur significant costs for its operation. Lately, the use of rotating magnetic fields has been tested in semiconductor crystal growth. The method is fairly well known and commonly used in metal processing but its adaptation to crystal growth of semiconductors is fairly recent. The elegance of the technique rests in its low power requirement (typically 10-20 milli-Tesla at 50-400 Hz) and its efficacy in curtailing deleterious temperature fluctuations in the melt. A rotating magnetic field imposes a rotational force and thereby induces a circulation within the melt that tends to dominate other sporadic convective effects. Thus a known low level of convective flow is introduced into the system. A new novel variation of the Lorentz force mechanism is proposed and investigated in this study. Since one of the desired process conditions in melt crystal growth is the minimization of convective effects, this investigation examines the use of an external field of magnetic origin to counteract existing convective flow within the melt. This is accomplished by utilizing a running or traveling axial magnetic wave in the system. The concept is similar to the use of vibrational means in order to induce streaming flows that oppose buoyant or surface tension driven convection in the system. The rotation direction as well as the magnitude (strength) of this circulation can be easily controlled by external inputs thus affording a direct means of controlling the developing shape of the crystallizing front (interface). The theoretical model of this technique is fully developed and presented in this paper. Results from the solution of the developed governing equations and boundary conditions are also presented. An experimental demonstration of the concept is presented through the suppression of natural convective flow in a mercury column. Implications to crystal growth systems will be fully explored in the final manuscript.

Mazuruk, K.↗
NASA NTRS2022

Space Induced Discharge Model to Be Used on Interface Circuit Vulnerability Assessment

Charging can occur in space due to plasma and high energy particles. Charging is especially severe in certain earth orbits such as geostationary orbit, medium earth orbit and highly inclined orbit. A spacecraft that travels to other planets such as Jupiter also encounters high level of charging due to Jupiter’s radiation belt. The major concern for spacecraft operation due to charging is the survival of sensitive electrical circuits. Dielectrics and floating metal can charge up and discharge onto interface circuits. The electrical circuits need to be evaluated for survival from the electrostatic-discharge (ESD) threat. Both ESD induced upset and damage need to be assessed. The interface circuits are often located inside a shielded unit in the spacecraft structure or a shielded vault with sufficient radiation shielding to avoid surface and internal charging. The materials external to the highly shielded area can charge up and discharge onto wires that can directly dissipate the ESD onto the interface circuits. Figure 1 shows a possible discharge mechanism onto the interfacing circuits. If the wire has an external shield or a coaxial line, some or most of the discharge can flow on the shield to chassis ground. There is always coupling from external discharge. If the wire has no shield, the discharge flows onto the interface circuit. Fig. 1 Possible discharge paths onto the interface circuits. The purpose of this article is to address how to model the ESD threat to be used as the discharge source onto the interface circuit based on the electron beam test data. Different models are needed for different impedances associated with the interface circuits. Once those models are derived, one needs to analyze the interface circuit with the corresponding model using an electrical simulation tool. Then, the circuit’s front-end components can be assessed for damage and upset. Upset is highly dependent on the circuit operation and system requirements, and thus it is not addressed any further. Damage of electrical parts is often rated under Human-Body-Model (HBM) ESD. One can utilize this industry-wide available information to assess the damage vulnerability due to space discharges. However, the HBM waveform does not always match the discharge waveform. To estimate the damage vulnerability of the circuit and set margin, one needs to compare the equivalent energy delivered to the victim between a HBM discharge and the results from an ESD source model based on actual discharges. This paper presents such a comparison for a few select cases.

Xie, Julie↗
NASA NTRS1997

Aerodynamics of a Simple Three-Dimensional High-Lift System

Experimental results are presented for a three-element high-lift system with a half-span flap and both a full- and three-quarter span slat. The simplified wing spanned the wind-tunnel test section with no wing sweep or wing tip. Full- and part-span slat configurations were tested to investigate the three-dimensional effects of a finite-span slat in a landing configuration. With a part-span slat, flow visualization and surface pressures revealed two vortices at the flap tip and one at the slat tip. Compared to the full-span slat case, the part-span slat configuration stalled at a lower angle of attack and yielded reduced lift at all angles of attack. At a representative approach attitude (10 deg), wake surveys indicated comparable induced drag for the full- and part-span slat configurations while the profile drag was more than doubled with a part-span slat. The vorticity contribution of the slat-tip vortex was minimal, suggesting that the vortex dissipates through viscous interaction with the main-element boundary layer. The addition of a lower-surface fence at the slat tip yielded two co-rotating vortices and partially alleviated the profile drag increase of the part-span slat.

Storms, Bruce↗
NASA NTRS2013

Spurious Behavior of Shock-Capturing Methods: Problems Containing Stiff Source Terms and Discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The smearing introduces a nonequilibrium state into the calculation. Thus as soon as a nonequilibrium value is introduced in this manner, the source term turns on and immediately restores equilibrium, while at the same time shifting the discontinuity to a cell boundary. The present study is to show that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Moreover, employing finite time steps and grid spacings that are below the standard Courant-Friedrich-Levy (CFL) limit on shockcapturing methods for compressible Euler and Navier-Stokes equations containing stiff reacting source terms and discontinuities reveals surprising counter-intuitive results. Unlike non-reacting flows, for stiff reactions with discontinuities, employing a time step and grid spacing that are below the CFL limit (based on the homogeneous part or non-reacting part of the governing equations) does not guarantee a correct solution of the chosen governing equations. Instead, depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The present investigation for three very different stiff system cases confirms some of the findings of Lafon & Yee (1996) and LeVeque & Yee (1990) for a model scalar PDE. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

Yee, Helen M. C.↗
NASA NTRS1996

Parallel Simulation of Unsteady Turbulent Flames

Time-accurate simulation of turbulent flames in high Reynolds number flows is a challenging task since both fluid dynamics and combustion must be modeled accurately. To numerically simulate this phenomenon, very large computer resources (both time and memory) are required. Although current vector supercomputers are capable of providing adequate resources for simulations of this nature, the high cost and their limited availability, makes practical use of such machines less than satisfactory. At the same time, the explicit time integration algorithms used in unsteady flow simulations often possess a very high degree of parallelism, making them very amenable to efficient implementation on large-scale parallel computers. Under these circumstances, distributed memory parallel computers offer an excellent near-term solution for greatly increased computational speed and memory, at a cost that may render the unsteady simulations of the type discussed above more feasible and affordable.This paper discusses the study of unsteady turbulent flames using a simulation algorithm that is capable of retaining high parallel efficiency on distributed memory parallel architectures. Numerical studies are carried out using large-eddy simulation (LES). In LES, the scales larger than the grid are computed using a time- and space-accurate scheme, while the unresolved small scales are modeled using eddy viscosity based subgrid models. This is acceptable for the moment/energy closure since the small scales primarily provide a dissipative mechanism for the energy transferred from the large scales. However, for combustion to occur, the species must first undergo mixing at the small scales and then come into molecular contact. Therefore, global models cannot be used. Recently, a new model for turbulent combustion was developed, in which the combustion is modeled, within the subgrid (small-scales) using a methodology that simulates the mixing and the molecular transport and the chemical kinetics within each LES grid cell. Finite-rate kinetics can be included without any closure and this approach actually provides a means to predict the turbulent rates and the turbulent flame speed. The subgrid combustion model requires resolution of the local time scales associated with small-scale mixing, molecular diffusion and chemical kinetics and, therefore, within each grid cell, a significant amount of computations must be carried out before the large-scale (LES resolved) effects are incorporated. Therefore, this approach is uniquely suited for parallel processing and has been implemented on various systems such as: Intel Paragon, IBM SP-2, Cray T3D and SGI Power Challenge (PC) using the system independent Message Passing Interface (MPI) compiler. In this paper, timing data on these machines is reported along with some characteristic results.

Menon, Suresh↗
NASA NTRS1993

On the dynamics of approximating schemes for dissipative nonlinear equations

Since one can rarely write down the analytical solutions to nonlinear dissipative partial differential equations (PDE's), it is important to understand whether, and in what sense, the behavior of approximating schemes to these equations reflects the true dynamics of the original equations. Further, because standard error estimates between approximations of the true solutions coming from spectral methods - finite difference or finite element schemes, for example - and the exact solutions grow exponentially in time, this analysis provides little value in understanding the infinite time behavior of a given approximating scheme. The notion of the global attractor has been useful in quantifying the infinite time behavior of dissipative PDEs, such as the Navier-Stokes equations. Loosely speaking, the global attractor is all that remains of a sufficiently large bounded set in phase space mapped infinitely forward in time under the evolution of the PDE. Though the attractor has been shown to have some nice properties - it is compact, connected, and finite dimensional, for example - it is in general quite complicated. Nevertheless, the global attractor gives a way to understand how the infinite time behavior of approximating schemes such as the ones coming from a finite difference, finite element, or spectral method relates to that of the original PDE. Indeed, one can often show that such approximations also have a global attractor. We therefore only need to understand how the structure of the attractor for the PDE behaves under approximation. This is by no means a trivial task. Several interesting results have been obtained in this direction. However, we will not go into the details. We mention here that approximations generally lose information about the system no matter how accurate they are. There are examples that show certain parts of the attractor may be lost by arbitrary small perturbations of the original equations.

Jones, Donald A.↗
NASA NTRS2012

Dielectric Performance of High Purity HTCC Alumina at High Temperatures - A Comparison Study with Other Polycrystalline Alumina

A very high purity (99.99+) high temperature co-fired ceramic (HTCC) alumina has recently become commercially available. The raw material of this HTCC alumina is very different from conventional HTCC alumina, and more importantly there is no glass additive in this co-fired material. Previously, selected HTCC and LTCC (low temperature co-fired ceramic) alumina materials were evaluated at high temperatures as dielectric and compared to a regularly sintered 96 polycrystalline alumina (96 Al2O3), where 96 alumina was used as the benchmark. A prototype packaging system based on regular 96 alumina with Au thick-film metallization successfully facilitated long term testing of high temperature silicon carbide (SiC) electronic devices for over 10,000 hours at 500C. In order to evaluate this new HTCC alumina for possible high temperature packaging applications, the dielectric properties of this HTCC alumina substrate were measured and compared with those of 96 alumina and a LTCC alumina from room temperature to 550C at frequencies of 120 Hz, 1 KHz, 10 KHz, 100 KHz, and 1 MHz. A parallel-plate capacitive device with dielectric of the HTCC alumina and precious metal electrodes were used for measurements of the dielectric constant and dielectric loss of the co-fired alumina material in the temperature and frequency ranges. The capacitance and AC parallel conductance of the capacitive device were directly measured by an AC impedance meter, and the dielectric constant and parallel AC conductivity of the dielectric were calculated from the capacitance and conductance measurement results. The temperature and frequency dependent dielectric constant, AC conductivity, and dissipation factor of the HTCC alumina substrate are presented and compared to those of 96 alumina. Other technical advantages of this new co-fired material for possible high packaging applications are also discussed.

High temperature↗