Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence rate”

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

Navier-Stokes simulation of the supersonic combustion flowfield in a ram accelerator

A computational study of the ram accelerator, a ramjet-in-tube device for accelerating projectiles to ultrahigh velocities, is presented. The analysis is performed using a fully implicit TVD scheme that efficiently solves the Reynolds-averaged Navier-Stokes equations and the species continuity equations associated with a finite rate combustion model. Previous analyses of this concept were based on inviscid assumptions. The present results indicate that viscous effects are of primary importance; in all the cases studied, shock-induced combustion always started in the boundary layer. The effects of Mach number, mixture composition, pressure, and turbulence are investigated for various configurations. Two types of combustion processes, one stable and the other unstable, were observed depending on the inflow conditions. In the unstable case, a detonation wave is formed, which propagates upstream and unstarts the ram accelerator. In the stable case, a solution that converges to steady-state is obtained, in which the combustion wave remains stationary with respect to the ram accelerator projectile. The possibility of stabilizing the detonation wave by means of a backward facing step is also investigated. In addition to these studies, two numerical techniques were tested. These two techniques are vector extrapolation to accelerate convergence, and a diagonal formulation that eliminates the expense of inverting large block matrices that arise in chemically reacting flows.

Yungster, Shaye↗

Control Demonstration of Multiple Doubly-Fed Induction Motors for Hybrid Electric Propulsion

The Convergent Aeronautics Solutions (CAS) High Voltage-Hybrid Electric Propulsion (HVHEP) task was formulated to support the move into future hybrid-electric aircraft. The goal of this project is to develop a new AC power architecture to support the needs of higher efficiency and lower emissions. This proposed architecture will adopt the use of the doubly-fed induction machine (DFIM) for propulsor drive motor application.The Convergent Aeronautics Solutions (CAS) High Voltage-Hybrid Electric Propulsion (HVHEP) task was formulated to support the move into future hybrid-electric aircraft. The goal of this project is to develop a new AC power architecture to support the needs of higher efficiency and lower emissions. This proposed architecture will adopt the use of the doubly-fed induction machine (DFIM) for propulsor drive motor application. DFIMs are attractive for several reasons, including but not limited to the ability to self-start, ability to operate sub- and super-synchronously, and requiring only fractionally rated power converters on a per-unit basis depending on the required range of operation. The focus of this paper is based specifically on the presentation and analysis of a novel strategy which allows for independent operation of each of the aforementioned doubly-fed induction motors. This strategy includes synchronization, soft-start, and closed loop speed control of each motor as a means of controlling output thrust; be it concurrently or differentially. The demonstration of this strategy has recently been proven out on a low power test bed using fractional horsepower machines. Simulation and hardware test results are presented in the paper.

CAS↗

A Hybrid-Learning Algorithm for Online Dynamic State Estimation in Multimachine Power Systems

With the increasing penetration of distributed generators in the smart grids, having knowledge of rapid real-time electromechanical dynamic states has become crucial to system stability control. Conventional Supervisory Control and Data Acquisition (SCADA)-based dynamic state estimation (DSE) techniques are limited by the slow sampling rates, while the emerging phasor measurement units (PMUs) technology enables rapid real-time measurements at network nodes. Using generator bus terminal voltages, we propose a hybrid-learning DSE (HL-DSE) algorithm to estimate the synchronous machine rotor angle and speed in real time. The HL-DSE takes the power system model into account and trains neuroestimators with real-time data in an online manner. Compared with traditional DSE methods, the HL-DSE overcomes limitations by using a data-driven approach in conjunction with the physical power system model. The time efficiency, accuracy, convergence, and robustness of the proposed algorithm are tested under noises and fault conditions in both small- and large-scale test systems. Simulation results show that the proposed HL-DSE is much more computationally efficient than widely used Kalman filter (KF)-based methods while maintaining comparable accuracy and robustness. In particular, HL-DSE is over 100 times faster than square-root unscented KF (SR-UKF) and 80 times faster than extended KF (EKF). The advantages and challenges of the HL-DSE are also discussed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

An extended Kalman filter for spinning spacecraft attitude estimation

An extended Kalman filter for real-time ground attitude estimation of a gyro-less spinning spacecraft was developed and tested. The filter state vector includes the angular momentum direction, phase angle, inertial nutation angle, and inertial and body nutation rates. The filter solves for the nutating three-axis attitude and accounts for effects due to principle axes offset from the body axes. The attitude is propagated using the kinematics of a rigid body symmetric about the principle spin axis; disturbance torques are assumed to be small. Filter updates consist only of the measured angles between celestial objects (Sun, Earth, etc.) and the nominal spin axis, and the times these angles were measured. Both simulated data and real data from the Dynamics Explorer -A (DE-A) spacecraft were used to test the filter; the results are presented. Convergence was achieved rapidly from a wide range of a priori state estimates, and sub-degree accuracy was attained. Systematic errors affecting the solution accuracy are discussed, as are the results of an attempt to solve for sensor measurement angle biases in the state vector.

Baker, David F.↗

Simulation of preburner sprays, volumes 1 and 2

The present study considered characteristics of sprays under a variety of conditions. Control of these sprays is important as the spray details can control both rocket combustion stability and efficiency. Under the present study Imperial College considered the following: (1) Measurement of the size and rate of spread of the sprays produced by single coaxial airblast nozzles with axial gaseous stream. The local size, velocity, and flux characteristics for a wide range of gas and liquid flowrates were measured, and the results were correlated with the conditions of the spray at the nozzle exit. (2) Examination of the effect of the geometry of single coaxial airblast atomizers on spray characteristics. The gas and liquid tube diameters were varied over a range of values, the liquid tube recess was varied, and the shape of the exit of the gaseous jet was varied from straight to converging. (3) Quantification of the effect of swirl in the gaseous stream on the spray characteristics produced by single coaxial airblast nozzles. (4) Quantification of the effect of reatomization by impingement of the spray on a flat disc positioned around 200 mm from the nozzle exit. This models spray impingement on the turbopump dome during the startup process of the preburner of the SSME. (5) Study of the interaction between multiple sprays without and with swirl in their gaseous stream. The spray characteristics of single nozzles were compared with that of three identical nozzles with their axis at a small distance from each other. This study simulates the sprays in the preburner of the SSME, where there are around 260 elements on the faceplate of the combustion chamber. (6) Design an experimental facility to study the characteristics of sprays at high pressure conditions and at supercritical pressure and temperature for the gas but supercritical pressure and subcritical temperature for the liquid.

Hardalupas, Y.↗

Simulation of preburner sprays, volume 1

Under the present contract a study was made on the characteristics of the sprays of coaxial injectors with particular emphasis on those aspects relevant to the performance of rocket engines. The technical effort was performed primarily by Imperial College of Science and Technology and the details of this effort are contained in the accompanying technical report by Hardalupas and Whitelaw, which forms Vol. 2 of the report. The role of Scientific Research Associates was review of material, initial direction of the technical effort and review of all reports. The study considered characteristics of sprays under a variety of conditions. Control of these sprays is important as the spray details can control both rocket combustion stability and efficiency. Imperial College considered the following: first, measurement of the size and rate of spread of the sprays produced by single coaxial airblast nozzles with axial gaseous stream. The local size, velocity and flux characteristics for a wide range of gas and liquid flowrates were measured and the results were correlated with the conditions of the spray at the nozzle exit. Second, examination of the effect of the geometry of single coaxial airblast atomizers on spray characteristics. The gas and liquid tube diameters were varied over a range of values, the liquid tube recess was varied and the shape of the exit of the gaseous jet was varied from straight to converging. Third, quantification of the effect of swirl in the gaseous stream on the spray characteristics produced by single coaxial airblast nozzles. Fourth, quantification of the effect of reatomization by impingement of the spray on a flat disc positioned around 200 mm from the nozzle exit. This models spray impingement on the turbopump dome during the startup process of the preburner of the Space Shuttle Main Engine (SSME). Fifth, study of the interaction between multiple sprays without and with swirl in their gaseous stream. The spray characteristics of single nozzles were compared with that of three identical nozzles with their axis at a small distance from each other. This study simulates the sprays in the preburner of the SSME, where there are around 260 elements on the faceplate of the combustion chamber. And lastly, design an experimental facility to study the characteristics of sprays at high pressure conditions and at supercritical pressure and temperature for the gas but supercritical pressure and subcritical temperature for the liquid.

Source record↗

Dirty dishes or dirty laundry? Comparing two methods for quantifying American consumers' preferences for load management in a smart home

One challenge of transitioning to renewable energy is that household electricity use and renewable generation are often misaligned. Smart home energy management systems hold promise for shifting usage to match generation, but these systems need to be designed with the occupants’ preferences in mind. The purpose of the present research is to compare two approaches for collecting and modeling consumers’ load management preferences, both of which are amenable to use in a home energy management system. Specifically, we examine the performance of Simple Multi-Attribute Rating Technique Exploiting Ranks (SMARTER) and Analytic Hierarchy Process (AHP) in quantifying consumers’ preferences regarding air temperature (air conditioning and heating), water heating, dishwashing, clothes washing and drying, monetary costs, environmental impacts, and comfort/convenience. Two studies are presented: Study 1 examines the SMARTER approach, and Study 2 focuses on the AHP approach. In both studies, online surveys (N SMARTER = 956 and N AHP = 1023) were conducted to elicit preferences from participants across the United States. The preferences modeled by both approaches were validated based on (a) their ability to predict participants’ choices in a Discrete Choice Experiment and (b) their convergence with previous research on load-shifting behavior. The validation procedure suggests that the SMARTER approach is superior in modeling consumers’ preferences for load management. Overall, this research lays the groundwork for designing a smart home interface capable of collecting occupants’ preferences and using those preferences to deliver improved occupant comfort, lower operating costs, reduced environmental impact, and more significant demand response than exists today.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Venus tectonic styles and crustal differentiation

Two of the most important constraints are known from Pioneer Venus data: the lack of a system of spreading rises, indicating distributed deformation rather than plate tectonics; and the high gravity/topography ratio, indicating the absence of an asthenosphere. In addition, the high depth/diameter ratios of craters on Venus indicate that Venus probably has no more crust than Earth. The problems of the character of tectonics and crustal formation and recycling are closely coupled. Venus appears to lack a recycling mechanism as effective as subduction, but may also have a low rate of crustal differentiation because of a mantle convection pattern that is more distributed, less concentrated, than Earth's. Distributed convection, coupled with the nonlinear dependence of volcanism on heat flow, would lead to much less magmatism, despite only moderately less heat flow, compared to Earth. The plausible reason for this difference in convective style is the absence of water in the upper mantle of Venus. We have applied finite element modeling to problems of the interaction of mantle convection and crust on Venus. The main emphasis has been on the tectonic evolution of Ishtar Terra, as the consequence of convergent mantle flow. The early stage evolution is primarily mechanical, with crust being piled up on the down-stream side. Then the downflow migrates away from the center. In the later stages, after more than 100 m.y., thermal effects develop due to the insulating influence of the thickened crust. An important feature of this modeling is the entrainment of some crustal material in downflows. An important general theme in both convergent and divergent flows is that of mixing vs. stratification. Models of multicomponent solid-state flow obtain that lower-density crustal material can be entrained and recycled, provided that the ration of low-density to high-density material is small enough (as in subducted slabs on Earth). The same considerations should apply in upflows; a small percent of partial melt may be carried along with its matrix and never escape to the surface. Models that assume melt automatically rising to the crust and no entrainment or other mechanism of recycling lower-density material obtain oscillatory behavior, because it takes a long time for heat to build up enough to overcome a Mg-rich low-density residuum. However, these models develop much thicker crust than consistent with estimates from crater depth/diameter ratios.

Kaula, W. M.↗

Global Oceanic Precipitation: A Joint View by TOPEX and the TOPEX Microwave Radiometer

The TOPEX/POSEIDON mission offers the first opportunity to observe rain cells over the ocean by a dual-frequency radar altimeter (TOPEX) and simultaneously observe their natural radiative properties by a three-frequency radiometer (TOPEX microwave radiometer (TMR)). This work is a feasibility study aimed at understanding the capability and potential of the active/passive TOPEX/TMR system for oceanic rainfall detection. On the basis of past experiences in rain flagging, a joint TOPEX/TMR rain probability index is proposed. This index integrates several advantages of the two sensors and provides a more reliable rain estimate than the radiometer alone. One year's TOPEX/TMR data are used to test the performance of the index. The resulting rain frequency statistics show quantitative agreement with those obtained from the Comprehensive Ocean-Atmosphere Data Set (COADS) in the Intertropical Convergence Zone (ITCZ), while qualitative agreement is found for other regions of the world ocean. A recent finding that the latitudinal frequency of precipitation over the Southern Ocean increases steadily toward the Antarctic continent is confirmed by our result. Annual and seasonal precipitation maps are derived from the index. Notable features revealed include an overall similarity in rainfall pattern from the Pacific, the Atlantic, and the Indian Oceans and a general phase reversal between the two hemispheres, as well as a number of regional anomalies in terms of rain intensity. Comparisons with simultaneous Global Precipitation Climatology Project (GPCP) multisatellite precipitation rate and COADS rain climatology suggest that systematic differences also exist. One example is that the maximum rainfall in the ITCZ of the Indian Ocean appears to be more intensive and concentrated in our result compared to that of the GPCP. Another example is that the annual precipitation produced by TOPEX/TMR is constantly higher than those from GPCP and COADS in the extratropical regions of the northern hemisphere, especially in the northwest Pacific Ocean. Analyses of the seasonal variations of prominent rainy and dry zones in the tropics and subtropics show various behaviors such as systematic migration, expansion and contraction, merging and breakup, and pure intensity variations. The seasonality of regional features is largely influenced by local atmospheric events such as monsoon, storm, or snow activities. The results of this study suggest that TOPEX and its follow-on may serve as a complementary sensor to the special sensor microwave/imager in observing global oceanic precipitation.

Chen, Ge↗

Vortex simulation of an inviscid shear layer

The accuracy of the vortex-blob method was tested by simulating a free-shear-layer instability, Kirchhoff's elliptical vortex, and a circular vortex. The main numerical parameters in the vortex-blob method are the density of the vortices, and the distribution of vorticity within each vortex core. The growth rate of a periodic unstable mode of the shear layer was calculated numerically and compared with the exact result. The error is only a few percent for about 10 rows of vortex blobs. The error is reduced by decreasing the spacing between vortices and, correspondingly, the core size. In the simulation of the motion of the elliptical vortex, the rotation of the boundary, without change of shape, and the circular particle paths of the vortical fluid were well simulated. For the circular vortex, optimum sets of parameters were obtained by comparing them with the exact velocity. The results are consistent with convergence theories of the vortex-blob method. In particular, second-order convergence is observed with a Gaussian core from velocity calculation.

Nakamura, Y.↗

An assessment of the utility of multirate time integration for the modeling of laser powder bed fusion

Finite element simulation of the laser powder bed fusion process is made challenging by the inherently multiscale nature of the process. When using the typical global time stepping techniques, slowly-evolving regions of the domain receive the same numerical treatment as the regions with the highest temperature rates. The current work details the implementation and evaluation of an implicit multirate method which is able to advance different regions of the domain with distinct time steps, depending on their current solution rate. Previous work indicates that this representation of the temporal scales of the problem can lead to significant reductions in wall clock run time, and it is shown herein that multirate time integration, when used with a uniform process-scale mesh, can result in speedups of approximately 19, 42, and 87 times, for domains with edge lengths of approximately 1 mm, 2 mm, and 5 mm, respectively. When used in conjunction with $\ h$-refinement (limited to two levels of refinement), the resulting speedups (taken relative to the uniform mesh) are around 28, 85, and 665 times. The method is demonstrated to converge as indicated by the literature, and its use with an AM-Bench domain is demonstrated.

42 ENGINEERING↗

Off-design operation and performance of pumped thermal energy storage

In this article, we describe off-design models and control strategies for a Pumped Thermal Energy Storage (PTES) system that uses liquid thermal energy storage: specifically molten salt for hot storage and methanol for cold storage. Off-design conditions arise when load-following, or due to variations in storage tank temperatures or ambient temperatures. We propose a control strategy that uses inventory control to manage the mass flow rate in the thermodynamic cycles, which facilitates load following. We also propose a control strategy for the storage fluid mass flow rates, which are varied to ensure the molten salt is maintained at its design temperature. This maximizes efficiency and minimizes problems with salt freezing or degradation. The cold storage fluid mass flow rate is varied so that the cold tanks have the same state-of-charge as the hot tanks. This leads to variations in cold fluid temperature, but these variations are shown to be acceptably small (e.g. 7.5% increase), and this control method is shown to be simpler and more efficient than an alternative strategy where tanks become unbalanced. The ambient temperature and storage tank temperatures are moved ±50 °C from the design values and the impact on power, duration, and tank temperatures is quantified. Results demonstrate that the proposed control strategy is stable and self-correcting - that is, storage temperatures converge on stable values after two-to-three charge-discharge cycles. When inputs return to design values, the system returns to its design point after two charge-discharge cycles. We also demonstrate that inventory control enables delivery of the target power output even when off-design conditions exist that would normally reduce the power output.

25 ENERGY STORAGE↗

Automating Bayesian inference and design to quantify acoustic particle levitation

Self-propulsion of micro- and nanoparticles powered by ultrasound provides an attractive strategy for the remote manipulation of colloidal matter using biocompatible energy inputs. Quantitative understanding of particle motion and its dependence on size, shape, and composition requires accurate characterization of the acoustic field, which depends sensitively on the experimental setup. Here, we show how automated experiments based on Bayesian inference and design can accurately and efficiently characterize the acoustic field within resonant chambers used to propel acoustic nanomotors. Repeated cycles of observation, inference, and design (OID) are guided by a physical model that describes the rate at which levitating particles approach the nodal plane. Using video microscopy, we observe the relaxation of tracer particles to this plane following the application of the acoustic field. We use sequential Monte Carlo methods to infer model parameters such as the amplitude and frequency of the resonant chamber while accounting for particle-level measurement noise and population-level heterogeneity in the field. Guided by simulated outcomes, we select the optimal design for the next experiment as to maximize the information gain in the relevant parameters. We show how this iterative process serves to discriminate between competing hypotheses and efficiently converges to accurate parameter estimates using only few automated experiments. We discuss the need for model criticism to ensure the validity of the guiding model throughout automated cycles of observation, inference, and design. Furthermore, this work demonstrates how Bayesian methods can learn the parameters of nonlinear, hierarchical models used to describe video microscopy data of active colloids.

42 ENGINEERING↗

Control Demonstration of Multiple Doubly Fed Induction Motors for Hybrid Electric Propulsion

The Convergent Aeronautics Solutions (CAS) High Voltage Hybrid Electric Propulsion (HVHEP) activity was formulated to support the move into future hybrid-electric aircraft. The goal of this high-voltage activity is to develop a new alternating current (AC) power architecture to support the needs of higher efficiency and lower emissions aircraft. This proposed architecture will adopt the use of the doubly fed induction machine (DFIM) for propulsor drive motor application. The DFIMs are attractive for several reasons, including but not limited to, the ability to self-start, the ability to operate subsynchronously and supersynchronously, and the requirement for power converters rated at a fraction of those required in a corresponding direct current system, dependent on the range of operation. The focus of this report is based specifically on the presentation and analysis of a novel strategy, which allows for independent operation of multiple doubly fed induction motors from a common AC bus

David J Sadey↗

Ma'adim Vallis Estuarine Delta in Elysium Basin and Its Relevance as a Landing Site for Exobiology Exploration on Mars

The debouche of Ma'adim Vallis in the Elysium Basin generated a transitional transported sediment structure, which planimetric shape is controlled by the enclosing topography of a deep reentrant gulf of the Basin into the highland. We defined it as an estuarine delta. The location and the importance of this estuarine delta is supported by the theoretical model of graded profile constructed for Ma'adim Vallis, and by two approaches: (i) the reconstruction of Ma'adim Vallis downstream course from Gusev to Elysium Basin, and (ii) the survey of the sediment deposit in the alleged estuary. The longitudinal graded profile of Ma'adim Vallis finds its base-level in the Elysium Basin, at a about 1000 m elevation, which is in agreement with the observed Basin shoreline. This model is supported by observational evidence of flow between the northern rim of Gusev crater, and the Elysium Basin shoreline. This downstream course of Ma'adim Vallis can be divided into three hydrogeologic regions. into three hydrogeologic regions. (a) The first region is a flooded plain (Zephiria Mensae), consisting in chaotic terrain formed by highland rocks, and disintegrated lava of the western flank of Apollinaris. Morphologic indicators of the flood process are: (1) the sediment deposit over the Gusev crater northern rim that reflects the overspilling of the crater-lake water through a 40-km wide gap provided by an ancient impact crater, (2) the tear-drop shaped feature on the northeastern flank of Apollinaris Patera, and (3) the chaotic terrain that suggest the emergence of ground water generated by the seepage of the crater lake through high-permeable broken rampart material. This underground water circulation sustained by the hydrostatic pressure of the crater-lake has likely generated a hydrothermal system in the volcanic environment of Apollinaris Patera. The stratigraphy of the flooded area is identified as Hesperian age, with occurrences of Noachian hilly individual features, and as Amazonian flooded plain and chaotic material, (b) The second region is located on the western flank of Apollinaris Patera. It is surrounded by relics of deep valleys that suggest a former downstream course of Ma'adim Vallis. The geologic setting of this region (Lucus Planum) is interpreted to be an Amazonian formation composed by the middle and lower members of the Medusae Formation., c) The third region corresponds to the convergence of the west and east branches of Ma'adirn Vallis into a deep re-entrant wide gulf that penetrates about 100 km into the highland. This topographic depression is delineated by the 1000 in elevation contour. This gulf has formed an estuarine configuration centered at 3S/190W within the Elysium Basin. This configuration has favored the formation of a estuarine sedimentary delta, because of topographically controlled lateral migration. This estuarine structure is strongly dominated by the incoming supply of Ma'adim Vallis fluvial sediment extracted from Zephiria Mensae and Lucus Planum. The obtuse-angle geometry of the estuary increases the sedimentation rate, which is higher than in the course of the channel. The sediment deposition process is governed by the estuarine water circulation. The inflowing loaded fluvial water enters the estuary as a bottom current, and mixes with the relatively less-loaded water of the receiving basin. When they mixed. the inflowing fluvial material, and the landward basin circulating water generate an accumulation of highly-diversified estuarine deposit stratification. This accumulation of material is mostly centered in the transitional zone of the delta. The sediment trapping efficiency of the estuary is function of the energy balance between the inflowing fluvial water, and the ingoing basin current. The submergence of the delta by the rising of the water-level increases the estuary water-depth, and consequently the sediment entrapment is favored. The locus of sediment accumulation moves landward in the zone of inflowing fluvial water. This results in the rising of the channel base-level, thus in the increase of the length of the longitudinal graded-profile. The sediment deposit facies of the zone A shows a generally smooth surface. The longitudinal deposit is bordered by alluvial terraces that reflect the variations of the channel level. The waning of the Elysium Basin caused the erosion of the Basin estuarine zone by small channels, this episode being characterized by dissected tear-drop shaped mesa-like morphologies in the delta. Our estuarine delta model predicts a lithostratigraphic depositional sequence associated with the water submergence and the transgression of Elysium Basin. The thickness of the estuarine sediment corresponds to the Elysium Basin levels changes relatively to the bed floor of the estuary, The depositional sequence of Ma'adim Vallis are described: (1) a pro-current filled region (A), where fluvial are longitudinally accumulated by the inflowing water, (2) inverse current from Elysium Basin (B), where fluvial and lacustrine sediments are accumulated, and (3) zone of current equilibrium (C), where the sediments are distributed as a shoreline at the boundary of the estuarine delta. The estuary sedimentology dynamics collects and keeps the record of the geologic unit material crossed by Ma'adim Vallis, and those of the lakebed deposit of Elysium Basin. The predicted mixed stratigraphic sequence from fluvial and lacustrine sediment makes this site an exceptional environment to concentrate potential multi-origin biologic records. We envision four possible strategies to explore this sedimentologic record: (1) longitudinal surface and subsurface traverses in region A to investigate outcrop levees, (2) exploration of the mesa walls in region B, (3) deep drilling hole lodging of the sequential deposits in the zones A and B, and (4) surface and subsurface exploration of the shoreline delta. The expected results for each of these strategies are: (1) in the deepest layers of region A are predicted frequent and abundant coarse material, sandy lenses lamination grading downward from sand to cobbles. Volcanic debris from the Noachian crustal Plateau unit material, hydrothermal altered rocks, carbonates, Hesperian and possibly Amazonian volcanic material, from Apollinaris Patera, altered rocks and carbonates from Zephiria Mensae are expected. As a favorable environment for inception of life, possible biological records are expected in transported rock, (2) At the surface, and subsurface (<=100 m), large deposits sandy to silted material from Elysium paleolake basin mixed with fine-grained sediments from Ma'adim Vallis are expected mostly in the upstream part of region B, (3) on the shoreline of the estuarine delta, abundant fine material from Elysium paleolake basin (evaporites, carbonates), mostly Amazonian in age are expected. The Ma'adini estuary is a favorable landing site for all the above mentioned science aspects, and .for its location. The site lies near the equator, which is favorable for the rover solar power supply, and at 1000m elevation, which is a favorable configuration for the descent system braking. Another advantage is the extent of the area of high scientific interest (33,000 sq km), which provides a good ellipse, and potential long study traverses.

Grin, E. A.↗

NUMERICAL MODELING OF A SOLID OXIDE FUEL CELL FOR USE IN REAL-TIME SIMULATION AND CYBER-PHYSICAL SYSTEMS

Cyber-physical systems provide a mechanism with which to investigate the physical phenomena and behavior of traditionally cost-prohibitive or otherwise fragile equipment. For the National Energy Technology Laboratory (NETL), this approach resulted in the Hybrid Performance (Hyper) facility which features a gas turbine-SOFC hybrid cycle utilizing real turbomachinery and a simulated SOFC stack. This allows for the investigation of combined cycle performance and control strategies, in an exhaustive manner, both without fear of destroying delicate state-of-the-art fuel cells, and with the full accuracy of real-world turbomachinery. Issues arose between the transient response of the SOFC model being limited to a sample time of 80 milliseconds, due to the calculation time of the SOFC model taking on average 40 milliseconds to calculate for a given timestep with spikes in calculation time reaching the 80 millisecond threshold. In order to be able to match the speed of transients from the turbomachinery and likewise better discern transient behavior, it was determined that the SOFC model must be optimized to operate at a sample time of 5 milliseconds. Therefore, it is necessary to optimize the SOFC model in order to decrease the calculation time from around 40 milliseconds, down to at the most 5 milliseconds. To do this, both the electrochemical algorithm and the thermal algorithm used to simulate the physical behavior of the SOFC are investigated to determine where improvements can be made. To this end the rootfinding numerical recipes of the electrochemical algorithm are investigated as the complex electrochemistry requires a highly iterative nested dual convergence loop to resolve the voltage-current relationship, and likewise the temporal discretization of the thermal algorithm is modified for the sake of higher accuracy and stability. Ultimately the new electrochemical algorithm featuring higher order rootfinding schemes proves to be efficient enough to reach the sub 5 millisecond target, signifying an order of magnitude reduction in calculation time, and when coupled with the new temporal discretization similar calculation time characteristics show that a fully implicit, higher order temporal discretization can also successfully be used if desired. Ultimately this result means that the cyber-physical simulation system can operate at higher sample rates, and resolve transient events at significantly higher resolution and fidelity.

Arias, Jesus↗

The sensory transduction pathways in bacterial chemotaxis

Bacterial chemotaxis is a useful model for investigating in molecular detail the behavioral response of cells to changes in their environment. Peritrichously flagellated bacteria such as coli and typhimurium swim by rotating helical flagella in a counterclockwise direction. If flagellar rotation is briefly reversed, the bacteria tumble and change the direction of swimming. The bacteria continuously sample the environment and use a temporal sensing mechanism to compare the present and immediate past environments. Bacteria respond to a broad range of stimuli including changes in temperature, oxygen concentration, pH and osmotic strength. Bacteria are attracted to potential sources of nutrition such as sugars and amino acids and are repelled by other chemicals. In the methylation-dependent pathways for sensory transduction and adaptation in E. coli and S. typhimurium, chemoeffectors bind to transducing proteins that span the plasma membrane. The transducing proteins are postulated to control the rate of autophosphorylation of the CheA protein, which in turn phosphorylates the CheY protein. The phospho-CheY protein binds to the switch on the flagellar motor and is the signal for clockwise rotation of the motor. Adaptation to an attractant is achieved by increasing methylation of the transducing protein until the attractant stimulus is cancelled. Responses to oxygen and certain sugars involve methylation-independent pathways in which adaption occurs without methylation of a transducing protein. Taxis toward oxygen is mediated by the electron transport system and changes in the proton motive force. Recent studies have shown that the methylation-independent pathway converges with the methylation-dependent pathway at or before the CheA protein.

Taylor, Barry L.↗

Characteristic momentum of Hydro + and a bound on the enhancement of the speed of sound near the QCD critical point

Near the critical point in the QCD phase diagram, hydrodynamics breaks down at a momentum where the frequency of the fastest hydrodynamic mode becomes comparable with the decay rate of the slowest nonhydrodynamic mode. Hydro+ was developed as a framework which extends the range of validity of hydrodynamics beyond that momentum value. This was achieved through coupling the hydrodynamic modes to the slowest nonhydrodynamic mode. In this work, analyzing the spectrum of linear perturbations in Hydro+, we find that a slow mode falls out of equilibrium if its momentum is greater than a characteristic momentum value. That characteristic momentum turns out to be set by the branch points of the dispersion relations. These branch points occur at the critical momenta of so-called spectral curves and are related to the radius of convergence of the derivative expansion. The existence of such a characteristic momentum scale suggests that a particular class of slow modes has no remarkable effect on the flow of the plasma. Based on these results, we show that there is an enhancement of the speed of sound due to the Hydro+ corrections compared to the known thermodynamic corrections near the critical point. This enhancement can not exceed a temperature-dependent upper bound which we derive near the critical point in the QCD phase diagram.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗