Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time step”

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,207 records · Page 67

Multifunction Vacuum Chamber for IC Metallization

Vacuum system chamber processing multilayer metallization on integrated circuits (IC's) performs four operations ordinarily done in separated equipment. Chamber etches holes, removes photoresist, cleans by sputter etching, and deposits final layer of metal. Combinedfunction chamber costs less than separate equipment. Chamber avoids exposing integrated circuits to room air and, to oxidation and dust, between steps. Eliminates time spent in transferring circuits from one apparatus to next.

Routh, D. E.↗

A numerical study of ENO and TVD schemes for shock capturing

The numerical performance of a second-order upwind-based total variation diminishing (TVD) scheme and that of a uniform second-order essentially non-oscillatory (ENO) scheme for shock capturing are compared. The TVD scheme used is a modified version of Liou, using the flux-difference splitting (FDS) of Roe and his superbee function as the limiter. The construction of the basic ENO scheme is based on Harten, Engquist, Osher, and Chakravarthy, and the 2-D extensions are obtained by using a Strang-type of fractional-step time-splitting method. Numerical results presented include both steady and unsteady, 1-D and 2-D calculations. All the chosen test problems have exact solutions so that numerical performance can be measured by comparing the computer results to them. For 1-D calculations, the standard shock-tube problems of Sod and Lax are chosen. A very strong shock-tube problem, with the initial density ratio of 400 to 1 and pressure ratio of 500 to 1, is also used to study the behavior of the two schemes. For 2-D calculations, the shock wave reflection problems are adopted for testing. The cases presented in this report include flows with Mach numbers of 2.9, 5.0, and 10.0.

Chang, Shih-Hung↗

Unsteady flow computation of oscillating flexible wings

The problem of unsteady flow around flexible wings is solved using the unsteady, compressible, thin-layer Navier-Stokes equations in conjunction with the unsteady, linearized, Navier-displacement equations. Starting with the initial shape of the wing, the Navier-Stokes equations are solved on an initial structured grid to obtain the steady-flow solution which is used for the initial conditions. The forced deformation motion of the wing boundaries is then applied, and the problem is solved accurately in time. During the time-accurate stepping, the Navier-displacement equations are used to solve for the grid deformation and sequently, the Navier-Stokes equations are used to solve for the flowfield. Two applications are presented; the first is for a pulsating oscillation of a bending-mode airfoil in transonic flow, and the second is for a bending-mode oscillation of a sharp-edged delta wing in supersonic flow.

Kandil, Osama A.↗

Dynamics of surface tension in microgravity environment

Time-dependent evolutions of the profile of free surface (bubble shapes) for a cylindrical container partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry, in low and microgravity environments, have been studied. Numerical computations of the dynamics of bubble shapes have been carried out (1) linear time-dependent functions of spin-up and spin-down in low and microgravity environments, (2) linear time-dependent functions of increasing and decreasing gravity environment at high and low rotating cylinder speeds, (3) time-dependent step functions of spin-up and spin-down in a low-gravity environment, and (4) sinusoidal function oscillation of the gravity environment at high and low rotating cylinder speeds.

Hung, R. J.↗

Evaluation of electrolytic tilt sensors for measuring model angle of attack in wind tunnel tests

The results of a laboratory evaluation of electrolytic tilt sensors as potential candidates for measuring model attitude or angle of attack in wind tunnel tests are presented. The performance of eight electrolytic tilt sensors was compared with that of typical servo accelerometers used for angle-of-attack measurements. The areas evaluated included linearity, hysteresis, repeatability, temperature characteristics, roll-on-pitch interaction, sensitivity to lead-wire resistance, step response time, and rectification. Among the sensors being evaluated, the Spectron model RG-37 electrolytic tilt sensors have the highest overall accuracy in terms of linearity, hysteresis, repeatability, temperature sensitivity, and roll sensitivity. A comparison of the sensors with the servo accelerometers revealed that the accuracy of the RG-37 sensors was on the average about one order of magnitude worse. Even though a comparison indicates that the cost of each tilt sensor is about one-third the cost of each servo accelerometer, the sensors are considered unsuitable for angle-of-attack measurements. However, the potential exists for other applications such as wind tunnel wall-attitude measurements where the errors resulting from roll interaction, vibration, and response time are less and sensor temperature can be controlled.

Wong, Douglas T.↗

Temporal averaging in a turbulent environment - Compensation for phase drifts in phase shifting interferometry

In an effort to improve the signal to noise in an interference experiment, we have developed a method to remove systematic phase drift between data sets acquired over long time intervals. Using this technique, it is possible to average repeatedly acquired phase measurements and improve the phase estimate without sacrificing spatial resolution. Results from tests using real-time phase stepping holographic interferometry applied to cantilever bending of a piezoelectric bimorph indicate that white noise has been reduced from 3 to less than 1 deg (lambda/360) by averaging 36 phase compensated data sets before object bending and 36 data sets after bending.

Ovryn, B.↗

Generic mission planning concepts for space astronomy missions

The past two decades have seen the rapid development of space astronomy, both manned and unmanned, and the concurrent proliferation of the operational concepts and software that have been produced to support each individual project. Having been involved in four of these missions since the '70's and three yet to fly in the present decade, the authors believe it is time to step back and evaluate this body of experience from a macro-systems point of view to determine the potential for generic mission planning concepts that could be applied to future missions. This paper presents an organized evaluation of astronomy mission planning functions, functional flows, iteration cycles, replanning activities, and the requirements that drive individual concepts to specific solutions. The conclusions drawn from this exercise are then used to propose a generic concept that could support multiple missions.

Guffin, O. T.↗

Planning applications in image analysis

We describe two interim results from an ongoing effort to automate the acquisition, analysis, archiving, and distribution of satellite earth science data. Both results are applications of Artificial Intelligence planning research to the automatic generation of processing steps for image analysis tasks. First, we have constructed a linear conditional planner (CPed), used to generate conditional processing plans. Second, we have extended an existing hierarchical planning system to make use of durations, resources, and deadlines, thus supporting the automatic generation of processing steps in time and resource-constrained environments.

Boddy, Mark↗

A Hybrid Boundary Element-Finite Volume Method for Unsteady Transonic Airfoil Flows

A hybrid boundary element finite volume method for unsteady transonic flow computation has been developed. In this method, the unsteady Euler equations in a moving frame of reference are solved in a small embedded domain (inner domain) around the airfoil using an implicit finite volume scheme. The unsteady full-potential equation, written in the same frame of reference and in the form of the Poisson equation. is solved in the outer domain using the integral equation boundary element method to provide the boundary conditions for the inner Euler domain. The solution procedure is a time-accurate stepping procedure, where the outer boundary conditions for the inner domain are updated using the integral equation -- boundary element solution over the outer domain. The method is applied to unsteady transonic flows around the NACA0012 airfoil undergoing pitching oscillation and ramp motion. The results are compared with those of an implicit Euler equation solver, which is used throughout a large computational domain, and experimental data.

Hu, Hong↗

Adjoint Techniques for Topology Optimization of Structures Under Damage Conditions

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation (Haftka and Gurdal, 1992) in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers (Akgun et al., 1998a and 1999). It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages (Haftka et al., 1983). A common method for topology optimization is that of compliance minimization (Bendsoe, 1995) which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system.. Shennan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this (Akgun et al., 1998b). SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

Optimization of Aerospace Structure Subject to Damage Tolerance Criteria

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers. It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages. A common method for topology optimization is that of compliance minimization which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system. Sherrnan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this. SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Wang, Chi R.↗

Convective Weather Forecast Quality Metrics for Air Traffic Management Decision-Making

Since numerical weather prediction models are unable to accurately forecast the severity and the location of the storm cells several hours into the future when compared with observation data, there has been a growing interest in probabilistic description of convective weather. The classical approach for generating uncertainty bounds consists of integrating the state equations and covariance propagation equations forward in time. This step is readily recognized as the process update step of the Kalman Filter algorithm. The second well known method, known as the Monte Carlo method, consists of generating output samples by driving the forecast algorithm with input samples selected from distributions. The statistical properties of the distributions of the output samples are then used for defining the uncertainty bounds of the output variables. This method is computationally expensive for a complex model compared to the covariance propagation method. The main advantage of the Monte Carlo method is that a complex non-linear model can be easily handled. Recently, a few different methods for probabilistic forecasting have appeared in the literature. A method for computing probability of convection in a region using forecast data is described in Ref. 5. Probability at a grid location is computed as the fraction of grid points, within a box of specified dimensions around the grid location, with forecast convection precipitation exceeding a specified threshold. The main limitation of this method is that the results are dependent on the chosen dimensions of the box. The examples presented Ref. 5 show that this process is equivalent to low-pass filtering of the forecast data with a finite support spatial filter. References 6 and 7 describe the technique for computing percentage coverage within a 92 x 92 square-kilometer box and assigning the value to the center 4 x 4 square-kilometer box. This technique is same as that described in Ref. 5. Characterizing the forecast, following the process described in Refs. 5 through 7, in terms of percentage coverage or confidence level is notionally sound compared to characterizing in terms of probabilities because the probability of the forecast being correct can only be determined using actual observations. References 5 through 7 only use the forecast data and not the observations. The method for computing the probability of detection, false alarm ratio and several forecast quality metrics (Skill Scores) using both the forecast and observation data are given in Ref. 2. This paper extends the statistical verification method in Ref. 2 to determine co-occurrence probabilities. The method consists of computing the probability that a severe weather cell (grid location) is detected in the observation data in the neighborhood of the severe weather cell in the forecast data. Probabilities of occurrence at the grid location and in its neighborhood with higher severity, and with lower severity in the observation data compared to that in the forecast data are examined. The method proposed in Refs. 5 through 7 is used for computing the probability that a certain number of cells in the neighborhood of severe weather cells in the forecast data are seen as severe weather cells in the observation data. Finally, the probability of existence of gaps in the observation data in the neighborhood of severe weather cells in forecast data is computed. Gaps are defined as openings between severe weather cells through which an aircraft can safely fly to its intended destination. The rest of the paper is organized as follows. Section II summarizes the statistical verification method described in Ref. 2. The extension of this method for computing the co-occurrence probabilities in discussed in Section HI. Numerical examples using NCWF forecast data and NCWD observation data are presented in Section III to elucidate the characteristics of the co-occurrence probabilities. This section also discusses the procedure for computing throbabilities that the severity of convection in the observation data will be higher or lower in the neighborhood of grid locations compared to that indicated at the grid locations in the forecast data. The probability of coverage of neighborhood grid cells is also described via examples in this section. Section IV discusses the gap detection algorithm and presents a numerical example to illustrate the method. The locations of the detected gaps in the observation data are used along with the locations of convective weather cells in the forecast data to determine the probability of existence of gaps in the neighborhood of these cells. Finally, the paper is concluded in Section V.

Chatterji, Gano B.↗

"Far" and "Near" Visual Acuity While Walking and the Collective Contributions of Non-Ocular Mechanisms to Gaze Stabilization

Gaze stabilization was quantified in subjects (n=11) as they walked on a motorized treadmill (1.8 m/s) and viewed visual targets at two viewing distances. A "far" target was positioned at 4 m (FAR) in front of the subject and the "near" target was placed at a distance of 0.5 m (NEAR). A direct measure of visual acuity was used to assess the overall effectiveness of the gaze stabilization system. The contributions of nonocular mechanisms to the gaze goal were also quantified using a measure of the distance between the subject and point in space where fixation of the visual target would require the least eye movement amplitude (i.e. the head fixation distance (HFD)). Kinematic variables mirrored those of previous investigations with the vertical trunk translation and head pitch signals, and the lateral translation and head yaw signals maintaining what appear as antiphase relationships. However, an investigation of the temporal relationships between the maxima and minima of the vertical translation and head pitch signals show that while the maximum in vertical translation occurs at the point of the minimum head pitch signal, the inverse is not true. The maximum in the head pitch signal lags the vertical translation minimum by an average of greater than 12 percent of the step cycle time. Three HFD measures, one each for data in the sagittal and transverse planes, and one that combined the movements from both planes, all revealed changes between the FAR and NEAR target viewing conditions. This reorganization of the nonocular degrees of freedom while walking was consistent with a strategy to reduce the magnitude of the eye movements required when viewing the NEAR target. Despite this reorganization, acuity measures show that image stabilization is not occurring while walking and viewing the NEAR target. Group means indicate that visual acuity is not affected while walking in the FAR condition, but a decrement of 0.15 logMAR (i.e. 1.5 eye chart lines) exists between the standing and walking acuity measures when viewing the NEAR target.

Peters, Brian T.↗

Fabrication and Modification of Nanoporous Silicon Particles

Silicon-based nanoporous particles as biodegradable drug carriers are advantageous in permeation, controlled release, and targeting. The use of biodegradable nanoporous silicon and silicon dioxide, with proper surface treatments, allows sustained drug release within the target site over a period of days, or even weeks, due to selective surface coating. A variety of surface treatment protocols are available for silicon-based particles to be stabilized, functionalized, or modified as required. Coated polyethylene glycol (PEG) chains showed the effective depression of both plasma protein adsorption and cell attachment to the modified surfaces, as well as the advantage of long circulating. Porous silicon particles are micromachined by lithography. Compared to the synthesis route of the nanomaterials, the advantages include: (1) the capability to make different shapes, not only spherical particles but also square, rectangular, or ellipse cross sections, etc.; (2) the capability for very precise dimension control; (3) the capacity for porosity and pore profile control; and (4) allowance of complex surface modification. The particle patterns as small as 60 nm can be fabricated using the state-of-the-art photolithography. The pores in silicon can be fabricated by exposing the silicon in an HF/ethanol solution and then subjecting the pores to an electrical current. The size and shape of the pores inside silicon can be adjusted by the doping of the silicon, electrical current application, the composition of the electrolyte solution, and etching time. The surface of the silicon particles can be modified by many means to provide targeted delivery and on-site permanence for extended release. Multiple active agents can be co-loaded into the particles. Because the surface modification of particles can be done on wafers before the mechanical release, asymmetrical surface modification is feasible. Starting from silicon wafers, a treatment, such as KOH dipping or reactive ion etching (RIE), may be applied to make the surface rough. This helps remove the nucleation layer. A protective layer is then deposited on the wafer. The protective layer, such as silicon nitride film or photoresist film, protects the wafer from electrochemical etching in an HF-based solution. A lithography technique is applied to pattern the particles onto the protective film. The undesired area of the protective film is removed, and the protective film on the back side of the wafer is also removed. Then the pattern is exposed to HF/surfactant solution, and a larger DC electrical current is applied to the wafers for a selected time. This step removes the nucleation layer. Then a DC current is applied to generate the nanopores. Next, a large electrical current is applied to generate a release layer. The particles are mechanically suspended in the solvent and collected by filtration or centrifuge.

Ferrari, Mauro↗

Telescope Innovations Improve Speed, Accuracy of Eye Surgery

One of the main components of NASA's vision for the future of space exploration will actually have a keen eye for the past. The James Webb Space Telescope (JWST), scheduled to launch in 2018, will have spectacular sight, after it reaches orbit, one of its main goals is to observe the first galaxies that formed in the early universe. "JWST offers new capabilities in the infrared well beyond what we can see from current telescopes, either on the ground or in space. It will let us explore the early universe, extrasolar planets, and really, all branches of astrophysics," says Lee Feinberg, optical telescope element manager for the JWST at Goddard Space Flight Center. Building such a keen space telescope is an astronomic task. Because JWST will gaze over such incredible distances, it requires very large mirrors. In fact, the primary mirror will be more than two stories in diameter and consists of 18 separate segments. Each segment must be perfectly smooth, flat, and scratch-free in order to deliver a view 13 billion light years away. Construction of the 18 mirror segments involved measuring, grinding, polishing, and testing - and more measuring, grinding, polishing, and testing - and more measuring, grinding, polishing, and testing (you get the idea). One of the most time consuming steps of the mirror development process, the grinding phase, can take years.

Source record↗

LiveCheckHSI: a Hardware/Software Co-verification Tool for Hyperspectral Imaging Systems with Embedded System-on-Chip Instrument Avionics

The emergent technology of system-on-chip (SoC) devices promises lighter, smaller, cheaper, and more capable and reliable space electronic systems that could help to unveil some of the most treasured secrets in our universe. This technology is an improvement over the technology that is currently used in space applications, which lags behind state-of-the-art commercial-off-the-shelf (COTS) equipment by several generations. SoC technology integrates all computational power required by next-generation space exploration science instruments onto a single chip. This paper describes hardware/software co-verification tools for the Xilinx Zynq-based control and data handling system that have been developed at the Jet Propulsion Laboratory (JPL) for visible-infrared imaging spectrometers. The system acquires and compresses images in real-time, in addition to programming the spectrometer (frame rate, exposure time), focus step motor, and heaters and reporting telemetry.

Dolman, David↗