Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “JuMP model”

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 289 records · Page 16

Differential GPS and system integration of the Low Visibility Landing and Surface Operations (LVLASO) demonstration

The LVLASO Flight Demonstration of ASTA concepts (FDAC) integrates NASA-Langley's electronic moving map display and Transport Systems Research Vehicle (TSRV) (a modified Boeing 737 aircraft); ARINC's VHF data link, GPS ground station, and automated controller workstation; and Norden's surface radar/airport movement safety system. Aircraft location is shown on the electronic map display in the cockpit. An approved taxi route as well as other aircraft and surface traffic are also displayed. An Ashtech Z12 Global Positioning System (GPS) receiver on the TSRV estimates the aircraft's position. In Differential mode (DSPS), the Ashtech receiver accepts differential C/A code pseudorange corrections from a GPS ground station. The GPS ground station provides corrections up to ten satellites. The corrections are transmitted on a VHF data link at a 1 Hz. rate using the RTCM-104 format. DGPS position estimates will be within 5 meters of actual aircraft position. DGPS position estimates are blended with position, velocity, acceleration, and heading data from the TSRV Air Data/Inertial Reference System (ADIRS). The ADIRS data is accurate in the short-term, but drifts over time. The DGPS data is used to keep the ADIRS position accurate. Ownship position, velocity, heading, and turn rate are sent at a 20 Hz. rate to the electronic map display. Airport traffic is detected by the airport surface radar system. Aircraft and vehicles such as fuel trucks and baggage carts are detected. The traffic's location, velocity, and heading are sent to the TSRV. To prevent traffic symbology from jumping each second when a location update arrives, velocity and heading are used to predict a new traffic location for each display update. Possible runway incursions and collisions can be shown on the electronic map. Integrating the different systems used in the FDAC requires attention to the underlying coordinate systems. The airport diagram displayed on the electronic map is obtained from published navigational charts. The charts reference the North American Datum of 1927 (NAD27) or a local state-plane coordinate system. GPS uses the World Geodetic Standard of 1984 (WGS84). Both NAD27 and WGS84 model the Earth as an ellipsoid, however, they use a different origin and different size ellipsoids. Latitudes and longitudes given in these systems can be converted to a Cartesian system with the origin at the Earth's center. The surface radar detects traffic in a locally-level, rho-theta coordiante system. The electronic airport diagram is stored using a flat XY coordinate system. The map origin is at the tower and is referenced as True North up. All ownship and other traffic positions must be converted to the electronic map's frame of reference for display.

Rankin, James M.↗

Motion of the heliospheric termination shock. 2: Energy loss due to energetic particle acceleration

We present a generalization of earlier analysis of the motion of the heliospheric termination shock in response to heliospheric disturbances (Barnes, 1993) (paper 1), to allow jump conditions that include an energy sink at the shock front. The motivation for this study is that acceleration of the anomalous cosmic ray component may in fact represent such a sink. We have idealized the situation by assuming an infinitely thin shock parameterized by a quantity lambda(0 less than or equals lambda less than or equals 1), defined as the fraction of solar wind energy that is lost due to acceleration of the energetic particle component. If the sink is strong (consuming, say, 50% or more of the incident solar wind energy), the model leads to the following principal conclusions: (1) the shocked plasma would be much denser and cooler than in the standard gasdynamic case, thereby leading to more favorable conditions for direct observation of the shocked plasma; (2) the equilibrium shock position would be slightly farther (less than 10%) from the Sun than in the standard model; (3) as in the gasdynamic case, the shock would normally be in motion, so that the instantaneous position of the termination shock is not determined by interstellar conditions but by the recent history of the wind that has passed through the termination shock; and (4) the response of the shock to upstream disturbances would be similar to the response in the gasdynamic case, but the speed of the new termination shock would be somewhat smaller (probably by a f actor of 4 or less). We estimate that this speed is typically, approximately a few tens of kilometers per second corresponding to an inward or outward excursion of order of less than 1 to several AU, rather less than conventional estimates of several tens of AU.

Barnes, Aaron↗

Characterization of Dendritic Spatially Extended 3D Patterns in Directional Solidification: Microgravity Experiments in DECLIC-DSI onboard ISS and 3D Phase-field Simulations

To clarify and characterize the fundamental physical mechanisms active in the dynamical formation of three-dimensional (3D) arrays of dendrites under diffusive growth conditions, in situ monitoring of series of experiments on transparent model alloy succinonitrile – 0.46 wt% camphor was carried out under low gravity in the DECLIC Directional Solidification Insert onboard the International Space Station. These experiments offer the very unique opportunity to observe in situ and characterize the dynamics of the microstructure formation and evolution in extended 3D patterns under microgravity environment. The analyses of the dendritic patterns for a broad range of growth velocities displaying different levels of sidebranching will be presented. Especially, the time evolution of primary spacing, in case of solidifications at constant pulling rate as well as for experiments with pulling rate jump, will be compared to 3D phase-field simulations, and the results will be discussed in terms of stable spacing range.

Kaihua Ji↗

Characterization of Dendritic Spatially Extended 3D Patterns in Directional Solidification: Microgravity Experiments in DECLIC-DSI Onboard ISS and 3D Phase-field Simulations

To clarify and characterize the fundamental physical mechanisms active in the dynamical formation of three-dimensional (3D) arrays of dendrites under diffusive growth conditions, in situ monitoring of series of experiments on transparent model alloy succinonitrile – 0.46 wt% camphor was carried out under low gravity in the DECLIC Directional Solidification Insert onboard the International Space Station. These experiments offer the very unique opportunity to observe in situ and characterize the dynamics of the microstructure formation and evolution in extended 3D patterns under microgravity environment. The analyses of the dendritic patterns for a broad range of growth velocities displaying different levels of sidebranching will be presented. Especially, the time evolution of primary spacing, in case of solidifications at constant pulling rate as well as for experiments with pulling rate jump, will be compared to 3D phase-field simulations, and the results will be discussed in terms of stable spacing range.

Kaihua Ji↗

Wind Variability of B Supergiants: The Two-component Stellar Wind of gamma Arae - No. 2

The stellar wind of the rapidly rotating early-B supergiant, gamma Ara, is studied using time series, high-resolution IUE spectroscopy secured over approx. 6 days in 1993 March. Results are presented based on an analysis of several line species, including N(N), C(IV), Si(IV), Si(III), C(II), and Al(III). The wind of this star is grossly structured, with evidence for latitude-dependent mass loss which reflects the role of rapid rotation. Independent, co-existing time variable features are identified at low-velocity (redward of approx. -750 km/s) and at higher-speeds extending to approx. -1500 km/s. The interface between these structures is 'defined' by the appearance of a discrete absorption component which is extremely sharp (in velocity space). The central velocity of this 'Super DAC' changes only gradually, over several days, between approx. -400 and -750 km/s in most of the ions. However, its location is shifted redward by almost 400 km/s in Al(III) and C(II), indicating that the physical structure giving rise to this feature has a substantial velocity and ionization jump. Constraints on the relative ionization properties of the wind structures are discussed, together with results based on SEI line-profile-fitting methods. The overall wind activity in gamma Ara exhibits a clear ion dependence, such that low-speed features are promoted in low-ionization species, including Al(III), C(II), and Si(III). We also highlight that - in contrast to most OB stars - there are substantial differences in the epoch-to-epoch time-averaged wind profiles of gamma Ara. We interpret the results in terms of a two-component wind model for gamma Ara, with an equatorially compressed low ionization region, and a high speed, higher-ionization polar outflow. This picture is discussed in the context of the predicted bi-stability mechanism for line-driven winds in rapidly rotating early-B type stars, and the formation of compressed wind regions in rapidly rotating hot stars. The apparent absence of a substantial shift in the wind ionization mixture of gamma Ara, and the normal nature of its photospheric spectrum, suggests wind-compression as the likely dominant cause for the observed equatorial density enhancements.

Prinja, R. K.↗

Information access in a dual-task context: testing a model of optimal strategy selection

Pilots were required to access information from a hierarchical aviation database by navigating under single-task conditions (Experiment 1) and when this task was time-shared with an altitude-monitoring task of varying bandwidth and priority (Experiment 2). In dual-task conditions, pilots had 2 viewports available, 1 always used for the information task and the other to be allocated to either task. Dual-task strategy, inferred from the decision of which task to allocate to the 2nd viewport, revealed that allocation was generally biased in favor of the monitoring task and was only partly sensitive to the difficulty of the 2 tasks and their relative priorities. Some dominant sources of navigational difficulties failed to adaptively influence selection strategy. The implications of the results are to provide tools for jumping to the top of the database, to provide 2 viewports into the common database, and to provide training as to the optimum viewport management strategy in a multitask environment.

User-Computer Interface↗

Low density gas dynamic wall boundary conditions

Low density nozzles or large expansion ratio nozzles used in space experience rarefaction effects near their exit in the form of velocity slip and temperature jump at the walls. In addition, the boundary layers become very thick and there is a very strong viscous/inviscid interaction. For these reasons no existing design technique has been found to accurately predict the nozzle flow properties up to the nozzle exit. The objective of this investigation was to examine the slip boundary conditions and formulate them in a form appropriate for use with a full Navier-Stokes numerical code. The viscous/inviscid interaction would automatically be accounted for by using a compressible Navier-Stokes code. Through examination of the interaction of molecules with solid surfaces, a model for the distribution function of the reflected molecules has been determined and this distribution function has been used to develop a new slip boundary condition that can be shown to yield more realistic surface boundary conditions.

Collins, F. G.↗

Time-Variable Gravity Signal due to Extratropic Pacific Water Mass Redistribution

Using the satellite-laser-ranging (SLR) data, Cox and Chao [2002] reported the detection of a large post-1998 anomaly (in the form of a positive jump) in the time series of Earth s lowest-degree gravity harmonic 52, or the dynamic oblateness. Among several groups now examining the mass redistribution in the global geophysical fluids in search of the cause(s), we report here a temporally coinciding anomalies found in the extratropic north + south Pacific basins. Clearly seen in the leading EOFPC mode for extratropic Pacific, these anomalies occurred in sea-surface height, sea-surface temperature, and temperature- and salinity-depth profiles. We based our analysis on two different data sources: TOPEX/Poseidon altimetry, and the ECCO ocean general circulation model output assimilating T/P data. The magnitude of these changes, when converted to equivalent J2 change, appears to be a few times too small to explain the observed J2 directly. These findings, and the fact that the anomalies occurred following the strong 1997-98 El Nino, suggest strong geophysical connection of the interannual-to-decadal variation of 52 with the Pacific Decadal Oscillation (PDO) and the ultimate global-change processes that cause PDO. More work is underway, and additional independent data sources are examined, paying close attention to the fact that the J2 anomaly has been reversing back to normal since 2001. These include: (1) cryospheric contributions (melting of glaciers and ice sheets); (2) land hydrological contributions; (3) polar sea influences ( e g , via deep flow); (4) fluid flow in Earth's core; (5) time-variable gravity signals from SLR in higher harmonic degree/order, including J3,J4, (2,1), and (2,2) coefficients, considering their lower signal-to-noise ratios; (6) Earth rotation data in terms of length-of-day and polar motion.

Chao, B. F.↗

Runway Independent Aircraft Extremely Short Takeoff and Landing

This report Summarizes the work done is support of NASA/Ames Runway Independent Aircraft Research during the summer of 2003. This work centered on the tasks laid out by the Statement of Work, which was to: Identify and assess operational scenarios including airport air and ground operations and how RIA operations would interface; 2) Identify critical technologies and create a list of technologies that might be pushed to provide a quantum jump in operating economy, reliability, and safety should sufficient finding be available; 3) Create public domain powered high lift methodologies; and 4) Identify and assess vehicle concepts that provide innovative approaches to RIA operations. All these tasks were accomplished, with certain areas needing additional exploration in future grant work. Three designs were analyzed to provide strawman configurations for the RIA operations. All three aircraft carried 60 passengers, with a stage length of 1,000 nautical miles. They were capable of operating with a balanced field length of 2000 feet or less. Three different technology approaches were explored. The first, the Model 115, was a mid-wing USB design, developed as a near-term, low risk concept. The second aircraft, the EMAX, used a directed thrust system, was a far-term, high-risk approach. The third configuration was the Model 114, whose development began in summer 2002. In addition, further research was conducted on issues related to STOL operations, such as noise concerns, SNI operations, and other areas of interest.

Hall, David W.↗

Particle drift, diffusion, and acceleration at shocks

The gradient and curvature drifts implicit in change of the ambient magnetic field at a hydromagnetic shock wave are incorporated into the diffusive theory of shock acceleration of charged particles. The conventional jump condition at the shock is modified by a term incorporating the large drift along the shock plane. This term vanished identically for one-dimensional systems, but must be included in general for shocks which are finite in transverse extent or which have transverse structure. It is found that the effect of the drift is such that the transverse drift rate is proportional to the acceleration rate, and for perpendicular shocks is exactly equal to the rate of change of energy in the V x B electric field observed in the shock frame. This establishes a connection with the 'shock drift' models which neglect diffusion.

Jokipii, J. R.↗

Modeling of vortex-induced oscillations based on indicial response approach

The indicial response approach, a modeling approach originally used for studying nonlinear problems in flight dynamics is applied to the study of vortex-induced oscillations phenomena. The indicial response of the velocity field is derived for the problem studied with emphasis on physical postulates involved. A full account of fluid dynamics effects is taken by considering the incompressible Navier-Stokes equations. The theory is applied to the particular case of flow past a cylinder in periodically forced motion to derive some salient features. The indicial response approach is shown to be equivalent to a currently popular approach based on the use of the amplitude equation. Jump and hysteresis phenomena that experiments indicate occur within the lock-in regime (where the velocity field oscillates at the forcing frequency) are captured by the approach.

Truong, Khiem Van↗

The Radiated Field Generated by a Monopole Source in a Short, Rigid, Rectangular Duct

This thesis develops a method to model the acoustic field generated by a monopole source placed in a moving rectangular duct. The walls of the duct are assumed to be infinitesimally thin and the source is placed at the center of the duct. The total acoustic pressure is written in terms of the free-space pressure, or incident pressure, and the scattered pressure. The scattered pressure is the augmentation to the incident pressure due to the presence of the duct. It satisfies a homogeneous wave equation and is discontinuous across the duct walls. Utilizing an integral representation of the scattered pressure, a set of singular boundary integral equations governing the unknown jump in scattered pressure is derived. This equation is solved by the method of collocation after representing the jump in pressure as a double series of shape functions. The solution obtained is then substituted back into the integral representation to determine the scattered pressure, and the total acoustic pressure at any point in the field. A few examples are included to illustrate the influence of various geometric and kinematic parameters on the radiated sound field.

Lakota, Barbara Anne↗

Locomotion training of legged robots using hybrid machine learning techniques

In this study artificial neural networks and fuzzy logic are used to control the jumping behavior of a three-link uniped robot. The biped locomotion control problem is an increment of the uniped locomotion control. Study of legged locomotion dynamics indicates that a hierarchical controller is required to control the behavior of a legged robot. A structured control strategy is suggested which includes navigator, motion planner, biped coordinator and uniped controllers. A three-link uniped robot simulation is developed to be used as the plant. Neurocontrollers were trained both online and offline. In the case of on-line training, a reinforcement learning technique was used to train the neurocontroller to make the robot jump to a specified height. After several hundred iterations of training, the plant output achieved an accuracy of 7.4%. However, when jump distance and body angular momentum were also included in the control objectives, training time became impractically long. In the case of off-line training, a three-layered backpropagation (BP) network was first used with three inputs, three outputs and 15 to 40 hidden nodes. Pre-generated data were presented to the network with a learning rate as low as 0.003 in order to reach convergence. The low learning rate required for convergence resulted in a very slow training process which took weeks to learn 460 examples. After training, performance of the neurocontroller was rather poor. Consequently, the BP network was replaced by a Cerebeller Model Articulation Controller (CMAC) network. Subsequent experiments described in this document show that the CMAC network is more suitable to the solution of uniped locomotion control problems in terms of both learning efficiency and performance. A new approach is introduced in this report, viz., a self-organizing multiagent cerebeller model for fuzzy-neural control of uniped locomotion is suggested to improve training efficiency. This is currently being evaluated for a possible patent by NASA, Johnson Space Center. An alternative modular approach is also developed which uses separate controllers for each stage of the running stride. A self-organizing fuzzy-neural controller controls the height, distance and angular momentum of the stride. A CMAC-based controller controls the movement of the leg from the time the foot leaves the ground to the time of landing. Because the leg joints are controlled at each time step during flight, movement is smooth and obstacles can be avoided. Initial results indicate that this approach can yield fast, accurate results.

Simon, William E.↗

Unexpected Nonlinear Effects in Superconducting Transition-Edge Sensors

When a normal metal transitions into the superconducting state the DC resistance drops from a finite value to zero over some finite transition width in temperature, current, and magnetic field. Superconducting transition-edge sensors (TESs) operate within this transition region and uses resistive changes to measure deposited thermal energy. This resistive transition is not perfectly smooth and a wide range of TES designs and materials show sub-structure in the resistive transition (as seen in smooth nonmonotonic behavior, jump discontinuities, and hysteresis in the devices current-voltage relation and derivatives of the resistance with respect to temperature, bias current, and magnetic field). TES technology has advanced to the point where for many applications this structure is the limiting factor in performance and optimization consists of finding operating points away from these structures. For example, operating at or near this structure can lead to nonlinearity in the detectors response and gain scale, limit the spectral range of the detector by limiting the usable resistive range, and degrade energy resolution. The origin of much of this substructure is unknown. This presentation investigates a number of possible sources in turn. First we model the TES as a superconducting weak-link and solve for the characteristic differential equations current and voltage time dependence. We find:(1) measured DC biased current-voltage relationship is the time-average of a much higher frequency limit cycle solution.(2) We calculate the fundamental frequency and estimate the power radiated from the TES treating the bias leads as an antennae.(3) The solution for a set of circuit parameters becomes multivalued leading to current transitions between levels.(4)The circuit parameters can change the measure resistance and mask the true critical current. As a consequence the TES resistance surface is not just a function of temperature, current, and magnetic field but is also a function of the circuit elements (such as shunt resistor, SQUID inductance, and capacitor values). In other words, same device measured in different electrical circuits will have a different resistive surface in temperature, current, and magnetic field. Next we consider that at the transition temperature of a superconductor both the magnetic penetration depth and coherence length are divergent. As a consequence these important characteristic length scales are changing with operating point. We present measurements on devices showing commensurate behavior between these characteristic lengths and the length scale of added normal metal structures. Reordering of proximity vortices leads to discontinuities and irreversibility of the current-voltage curves. Last we consider a weak-link TES including both thermal activated resistance effects and the effect of the magnetic penetration depth being a function of temperature and magnetic field. We derive its impact on the resistive transition surface and the important device parameters a and b.

Superconducting Detectors↗

Marangoni Effects of a Drop in an Extensional Flow: The Role of Surfactant Physical Chemistry

While the changes in stresses caused by surfactant adsorption on non-deforming interfaces have been fairly well established, prior to this work, there were few studies addressing how surfactants alter stresses on strongly deforming interfaces. We chose the model problem of a drop in a uniaxial extensional flow to study these stress conditions To model surfactant effects at fluid interfaces, a proper description of the dependence of the surface tension on surface concentration, the surface equation of state, is required. We have adopted a surface equation of state that accounts for the maximum coverage limit; that is, because surfactants have a finite cross sectional area, there is an upper bound to the amount of surfactant that can adsorb in a monolayer. The surface tension reduces strongly only when this maximum coverage is approached. Since the Marangoni stresses go as the derivative of the surface equation of state times the surface concentration gradient, the non-linear equation of state determines both the effect of surfactants in the normal stress jump, (which is balanced by the product of the mean curvature of the interface times the surface tension), and the tangential stress jump, which is balanced by Marangoni stresses. First, the effects of surface coverage and intermolecular interactions among surfactants which drive aggregation of surfactants in the interface were studied. (see Pawar and Stebe, Physics of Fluids).

Stebe, Kathleen J.↗

Approximate scaling properties of RNA free energy landscapes

RNA free energy landscapes are analysed by means of "time-series" that are obtained from random walks restricted to excursion sets. The power spectra, the scaling of the jump size distribution, and the scaling of the curve length measured with different yard stick lengths are used to describe the structure of these "time series". Although they are stationary by construction, we find that their local behavior is consistent with both AR(1) and self-affine processes. Random walks confined to excursion sets (i.e., with the restriction that the fitness value exceeds a certain threshold at each step) exhibit essentially the same statistics as free random walks. We find that an AR(1) time series is in general approximately self-affine on timescales up to approximately the correlation length. We present an empirical relation between the correlation parameter rho of the AR(1) model and the exponents characterizing self-affinity.

Signal Processing, Computer-Assisted↗

Monte Carlo Bayesian Inference on a Statistical Model of Sub-gridcolumn Moisture Variability Using High-resolution Cloud Observations : Sensitivity Tests and Results - Part II

Part I presented a Monte Carlo Bayesian method for constraining a complex statistical model of GCM sub-gridcolumn moisture variability using high-resolution MODIS cloud data, thereby permitting large-scale model parameter estimation and cloud data assimilation. This part performs some basic testing of this new approach, verifying that it does indeed significantly reduce mean and standard deviation biases with respect to the assimilated MODIS cloud optical depth, brightness temperature and cloud top pressure, and that it also improves the simulated rotational-Ramman scattering cloud optical centroid pressure (OCP) against independent (non-assimilated) retrievals from the OMI instrument. Of particular interest, the Monte Carlo method does show skill in the especially difficult case where the background state is clear but cloudy observations exist. In traditional linearized data assimilation methods, a subsaturated background cannot produce clouds via any infinitesimal equilibrium perturbation, but the Monte Carlo approach allows finite jumps into regions of non-zero cloud probability. In the example provided, the method is able to restore marine stratocumulus near the Californian coast where the background state has a clear swath. This paper also examines a number of algorithmic and physical sensitivities of the new method and provides guidance for its cost-effective implementation. One obvious difficulty for the method, and other cloud data assimilation methods as well, is the lack of information content in the cloud observables on cloud vertical structure, beyond cloud top pressure and optical thickness, thus necessitating strong dependence on the background vertical moisture structure. It is found that a simple flow-dependent correlation modification due to Riishojgaard (1998) provides some help in this respect, by better honoring inversion structures in the background state.

Monte Carlo Bayesian inference↗

Fast-shock formation in line-tied magnetic reconnection models of solar flares

In a previous study by the author, an approximately stationary fast shock was tentatively identified in a numerical experiment designed to study line-tied magnetic reconnection. Here the evidence for the occurrence of a stationary fast shock is reexamined, and the previous identification is confirmed. In the numerical experiment, line-tied reconnection is modeled by a configuration which produces two supermagnetosonic outflow jets - one directed upward, away from the photosphere, and one directed downward, toward an arcade of closed magnetic loops tied to the photosphere. The fast shock occurs when the downward-directed jet encounters the obstacle formed by the closed loops. Although the existence of a stationary, or nearly stationary, fast shock is confirmed, the transition from the supermagnetosonic flow region upstream of the shock to the nearly static region downstream of the shock is more complicated than was previously thought. Immediately downstream of the shock, there exists a deflection sheath in which the submagnetosonic flow coming out of the shock is diverted around the region of static closed loops. The MHD jump conditions are used to investigate the characteristics of the fast shock and to show that a stationary shock cannot exist unless accompanied by a deflection sheath. Analysis of the shock's location and dimensions suggests that such fast shocks may contribute to particle acceleration and to thermal condensation in flares.

Forbes, T. G.↗