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 415 records · Page 23

Electromagnetic structure of the magnetopause and boundary layer

After a review of the properties and predictions of the closed and open models of the magnetopause, OGO-5 magnetometer data are used to illustrate various observed signatures of the magnetopause current layer and the adjacent plasma boundary layer. Among the topics touched upon are: fluctuations, diamagnetic effects, and field aligned currents in the boundary layer; one dimensionality of the magnetopause; presence and absence of a magnetic field component perpendicular to the magnetopause; finite ion gyroradius effects. A brief summary is given of existing Vlasov theory for the description of tangential, rotational, and contact discontinuities. Special attention is paid to the tangential momentum balance and the jump conditions at a rotational discontinuity. Low frequency fluctuations are discussed with emphasis on the signatures of the tearing mode.

Sonnerup, B. U. O.↗

Variability of the ocean-induced magnetic field predicted at sea surface and at satellite altitudes

Spatial and temporal variability of the magnetic field component induced by ocean circulation is investigated on the basis of a standard thin-shell approximation of electro- and magneto-static equations. Well-known difficulties of numerical solution of the governing equations are resolved by reducing the problem to an equation for the electric field potential,(ef) as opposed to a more conventional approach focused on the vertical jump, (psi), of the magnetic field potential across a combined ocean/ marine-sediment-layer spherical shell. The present formulation permits using more realistic input data on ocean currents and ultimately yields much greater (by at least an order of magnitude) values of the magnetic field at sea surface than predicted in earlier studies. Such large values are comparable to, and in some cases exceed, magnetic field variations caused by lithospheric and ionospheric sources on monthly to interannual timescales. At the 400-km altitude (of CHAMP satellite), the field attains 6 nT. The model predictions show favorable comparisons with some in situ measurements as well as with Challenging Minisatellite Payload (CHAMP) satellite magnetometer data.

ocean induced magnetic field↗

Rover Attitude and Pointing System Simulation Testbed

The MER (Mars Exploration Rover) Attitude and Pointing System Simulation Testbed Environment (RAPSSTER) provides a simulation platform used for the development and test of GNC (guidance, navigation, and control) flight algorithm designs for the Mars rovers, which was specifically tailored to the MERs, but has since been used in the development of rover algorithms for the Mars Science Laboratory (MSL) as well. The software provides an integrated simulation and software testbed environment for the development of Mars rover attitude and pointing flight software. It provides an environment that is able to run the MER GNC flight software directly (as opposed to running an algorithmic model of the MER GNC flight code). This improves simulation fidelity and confidence in the results. Further more, the simulation environment allows the user to single step through its execution, pausing, and restarting at will. The system also provides for the introduction of simulated faults specific to Mars rover environments that cannot be replicated in other testbed platforms, to stress test the GNC flight algorithms under examination. The software provides facilities to do these stress tests in ways that cannot be done in the real-time flight system testbeds, such as time-jumping (both forwards and backwards), and introduction of simulated actuator faults that would be difficult, expensive, and/or destructive to implement in the real-time testbeds. Actual flight-quality codes can be incorporated back into the development-test suite of GNC developers, closing the loop between the GNC developers and the flight software developers. The software provides fully automated scripting, allowing multiple tests to be run with varying parameters, without human supervision.

Vanelli, Charles A.↗

Nonlinear dynamics of global atmospheric and Earth-system processes

Researchers are continuing their studies of the nonlinear dynamics of global weather systems. Sensitivity analyses of large-scale dynamical models of the atmosphere (i.e., general circulation models i.e., GCM's) were performed to establish the role of satellite-signatures of soil moisture, sea surface temperature, snow cover, and sea ice as crucial boundary conditions determining global weather variability. To complete their study of the bimodality of the planetary wave states, they are using the dynamical systems approach to construct a low-order theoretical explanation of this phenomenon. This work should have important implications for extended range forecasting of low-frequency oscillations, elucidating the mechanisms for the transitions between the two wave modes. Researchers are using the methods of jump analysis and attractor dimension analysis to examine the long-term satellite records of significant variables (e.g., long wave radiation, and cloud amount), to explore the nature of mode transitions in the atmosphere, and to determine the minimum number of equations needed to describe the main weather variations with a low-order dynamical system. Where feasible they will continue to explore the applicability of the methods of complex dynamical systems analysis to the study of the global earth-system from an integrative viewpoint involving the roles of geochemical cycling and the interactive behavior of the atmosphere, hydrosphere, and biosphere.

Saltzman, Barry↗

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↗

The Role of Early Giant-planet Instability in Terrestrial Planet Formation

The terrestrial planets are believed to have formed by violent collisions of tens of lunar- to Mars-size protoplanets at time t < 200 Myr after the protoplanetary gas disk dispersal (t {sub 0}). The solar system giant planets rapidly formed during the protoplanetary disk stage and, after t {sub 0}, radially migrated by interacting with outer disk planetesimals. An early (t < 100 Myr) dynamical instability is thought to have occurred with Jupiter having gravitational encounters with a planetary-size body, jumping inward by ∼0.2–0.5 au, and landing on its current, mildly eccentric orbit. Here we investigate how the giant-planet instability affected the formation of the terrestrial planets. We study several instability cases that were previously shown to match many solar system constraints. We find that resonances with giant planets help to remove solids available for accretion near ∼1.5 au, thus stalling the growth of Mars. It does not matter, however, whether the giant planets are placed on their current orbits at t {sub 0} or whether they realistically evolve in one of our instability models; the results are practically the same. The tight orbital spacing of Venus and Earth is difficult to reproduce in our simulations, including cases where bodies grow from a narrow annulus at 0.7–1 au, because protoplanets tend to spread radially during accretion. The best results are obtained in the narrow-annulus model when protoplanets emerging from the dispersing gas nebula are assumed to have (at least) the mass of Mars. This suggests efficient accretion of the terrestrial protoplanets during the first ∼10 Myr of the solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

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

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗

Benchmark microgravity experiments and computations for 3D dendritic-array stability in directional solidification

In this study, we present a comprehensive quantitative analysis of stability bands for dendritic arrays during directional solidification of a transparent succinonitrile-0.46 wt % camphor alloy, spanning a broad range of pulling velocities. Taking advantage of the microgravity environment aboard the International Space Station where most convection effects are suppressed, we obtain unique measurements that quantify the stable primary spacing range of spatially extended three-dimensional dendritic array structures under purely diffusive growth conditions. Through carefully designed velocity jump experiments and detailed examination of sub-grain boundary dynamics, we characterize key instabilities, including elimination and tertiary branching, shedding new light on the mechanisms governing dynamic dendritic spacing selection in extended 3D arrays. Phase field simulations are performed to characterize the stability limits of dendritic array structures for quantitative comparison with the flight experiments. Although the simulations capture general trends, significant deviations are noted at the upper stability boundary, indicating the influence of additional, unexplored factors. These findings contribute to a deeper understanding of dendritic growth dynamics and offer valuable benchmark data that could aid in refining predictive models and improving control of dendritic microstructures in metallurgical applications.

36 MATERIALS SCIENCE↗

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↗