Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BOUND STATE”

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 829 records · Page 46

Design of Energy Storage Reactors for Dc-To-Dc Converters

Two methodical approaches to the design of energy-storage reactors for a group of widely used dc-to-dc converters are presented. One of these approaches is based on a steady-state time-domain analysis of piecewise-linearized circuit models of the converters, while the other approach is based on an analysis of the same circuit models, but from an energy point of view. The design procedure developed from the first approach includes a search through a stored data file of magnetic core characteristics and results in a list of usable reactor designs which meet a particular converter's requirements. Because of the complexity of this procedure, a digital computer usually is used to implement the design algorithm. The second approach, based on a study of the storage and transfer of energy in the magnetic reactors, leads to a straightforward design procedure which can be implemented with hand calculations. An equation to determine the lower-bound volume of workable cores for given converter design specifications is derived. Using this computer lower-bound volume, a comparative evaluation of various converter configurations is presented.

Chen, D. Y.↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

Far field conditions for compressible flows

The numerical treatment of exterior regions requires a method to convert the problem to one in a bounded region. One method is to truncate the unbounded region at some finite, artificial surface. This creates a finite computational region at the expense of imposing boundary conditions at the artificial boundary. If the problem has wave-like solutions near infinity, then these boundary conditions must simulate the radiation of energy out of the computational domain and towards infinity. Incorrect specification of these radiation boundary conditions can cause spurious reflected waves to be generated at the artificial boundary. Radiation boundary conditions which have the following properties were constructed: to accurately simulate the radiation of energy out of the computational domain; to accomplish the simulation with an accuracy which improves as fast as feasible as the artificial surface is moved outward; and to accelerate convergence to the steady state (by minimizing spurious reflections).

Bayliss, A.↗

Magnetic reconnection in the presence of sheared flow and density asymmetry: Applications to the Earth's magnetopause

Classical models of magnetic reconnection consist of a small diffusion region bounded by two symmetric slow shocks, across which the plasma is accelerated. Asymmetries often present in space plasmas are sheared plasma flow and dissimilar plasma densities on the two sides of current sheets. In this paper, we investigate magnetic reconnection in the presence of a shear flow and an asymmetric density across the current sheet using two-dimensional magnetohydrodynamic (MHD) simulations. The results demonstrate that magnetic reconnection can occur only for a plasma flow velocity (in the frame of the X line) which is below the Alfven speed in each inflow region. This limits the velocity of the X line to a certain range for a given flow shear and provides an upper limit to the total velocity shear at which reconnection ceases to operate. Depending on the direction of the flow in the adjacent inflow region, the effects from the sheared flow and from the density asymmetry will compete with or enhance each other in respect to the magnitude and location of the currents which bound the outflow regions. The results are applied to the dayside and flank regions of the magnetosphere. For the dayside region where the magnetosheath flow is slow, the magnetic field transition region is thin and the accelerated flow is earthward of the sharp current layer (magnetopause). At the flanks tailward of the X line, shear flow and density asymmetry effects compete making the magnetic field transition layer broad with the high-speed flow contained within the transition region which explains corresponding observations. At the flanks sunward of the X line, shear flow and density asymmetry effects enhance each other and lead to a strong current sheet on the magnetosheath side of the accelerated flow. The total volume affected by magnetic reconnection is much larger than the steady state region. A large bulge region precedes the steady state region. Qualitatively, the bulge and the steady state region have similar signatures and both can explain observations. We provide criteria in order to distinguish between the bulge and the steady state region in observations.

La Belle-Hamer, A. L.↗

Forced baroclinic ocean motions. II - The linear equatorial bounded case

Results obtained in the first part of this paper (1976), which studied forced baroclinic ocean motions in the linear equatorial unbounded case, are extended to the case of an ocean bounded by two meridians. A complete solution for the linear spin-up in response to a switched-on, x-dependent wind stress is obtained in terms of five components: the unbounded nonoscillatory response to the wind stress forcing; the inertia-gravity waves generated when the forcing switches on, together with their reflections; the quasi-geostrophic Rossby modes forming the eastern boundary response and a western boundary layer; and the equatorial Kelvin wave component of the western boundary response. Kelvin waves are generated if the forcings include a zonal wind component that is symmetric about the equator or an antisymmetric meridional wind stress. For non-Kelvin symmetries, spin-up occurs entirely by the effects of Rossby waves emanating from the eastern boundary. The state approaches the steady solution as increasing numbers of these reach an interior point. In Kelvin symmetries, the Rossby waves act to bring the sea surface tilt to the steady value and the spin-up of the tilt proceeds as in non-Kelvin cases.

Cane, M. A.↗

Information Science Panel joint meeting with Imaging Science Panel

Specific activity in information extraction science (taken to include data handling) is needed to: help identify the bounds of practical missions; identify potential data handling and analysis scenarios; identify the required enabling technology; and identify the requirements for a design data base to be used by the disciplines in determining potential parameters for future missions. It was defined that specific analysis topics were a function of the discipline involved, and therefore no attempt was made to define any specific analysis developments required. Rather, it was recognized that a number of generic data handling requirements exist whose solutions cannot be typically supported by the disciplines. The areas of concern were therefore defined as: data handling aspects of system design considerations; enabling technology for data handling, with specific attention to rectification and registration; and enabling technology for analysis. Within each of these areas, the following topics were addressed: state of the art (current status and contributing factors); critical issues; and recommendations for research and/or development.

Source record↗

Oxidation of the Mantle Wedge by H+ in Aqueous Fluids: A New Interdisciplinary Approach

The observation that arc magmas are the most oxidized on Earth have led petrologists to question whether the subduction process might cause oxidation of the sub-arc mantle source. A strong correlation between the input of slab-derived aqueous fluid and Fe3(+)/ΣFe in arc magmas has led to the hypothesis that slab fluids may facilitiate the transfer of redox potential from oxidized slab material to mantle wedge and subsequently to primary arc melts. Despite this intuitive link, identifying an efficient and ubiquitous chemical process to transfer oxidation state in slab fluids has been challenging. Pure H2O alone is an inefficient oxidizer in the mantle, necessitating another oxidant within the fluid. Commonly invoked components include S or C, which are likely heterogeneously distributed within subducting slabs, and direct transfer of Fe3(+), which is only possible in very solute-rich fluids. Here we present a new mechanism to explain the oxidation of the sub-arc mantle by slab-derived aqueous fluid, wherein dissolved H(+)(aq) is reduced to H2(aq). Redox equilibrium is satisfied by electron transfer from Fe2(+)(s), which is oxidized in the rock to Fe3(+)(s) in the reaction: 2Fe2(+)(s) + 2H(+)(aq) = 2Fe3(+)(s) + H2(aq). This simple reaction was previously overlooked due to the assumption that all H must be bound by H2O molecules (i.e., no dissolved H), applicable to a pure H2O fluid, but not one with a solute load as addressed here. The application of thermodynamic modeling tools from aqueous geochemistry (DEW, EQ3/6) to mantle petrology predicts this reaction between mantle rock and slab fluid with our own experimentally constrained chemistry and fO2 (QFM+2). Mass transfer model results show an increase in Fe3(+)/ΣFe in mantle rock from MORB- like (0.15) to arc-like (0.2–0.3) values at reasonable fluid fluxes, co-incident with increasing H2 activity in the fluid. This can explain the stark correlation between slab fluid input into the mantle wedge and the observed fO2 of arc magmas. Moreover, it does not require diffusion of hydrogen out of the system. Instead, H2 remains dissolved in fluid in equilibrium with oxidized mantle rock.

Iacovino, Kayla↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand while simultaneously maintaining air travel as one of the safest forms of transportation. One of the reasons for this success is the ability of the air traffic control system and the operators to adapt and accommodate to situations that routinely disrupt normal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators’ ability to control. These factors can lead to states where automation is unable to properly handle these issues, and therefore air traffic controllers and pilots have to intervene — ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions, complexity increases. This is because, under these conditions humans are required to make tactical decisions in response to external factors. This results in a departure from the original strategic plan where operations would be more efficiently managed. Human operators manage airspace complexity under rigid regulations but in a constantly changing environment. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. Some prior studies devised airspace complexity metrics in commercial aviation and related these metrics to controller workload. The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic — including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to ours that identifies such contributing factors or precursor patterns.

Precursor↗

Direct simulation of isothermal-wall supersonic channel flow

The motivation for this work is the fact that in turbulent flows where compressibility effects are important, they are often poorly understood. A few examples of such flows are those associated with astrophysical phenomena and those found in combustion chambers, supersonic diffusers and nozzles, and over high-speed airfoils. For this project, we are primarily interested in compressibility effects near solid surfaces. Our main objective is an improved understanding of the fundamentals of compressible wall-bounded turbulence, which can in turn be used to cast light upon modeling concepts such as the Morkovin hypothesis and the Van Driest transformation. To this end, we have performed a direct numerical simulation (DNS) study of supersonic turbulent flow in a plane channel with constant-temperature walls. All of the relevant spatial and temporal scales are resolved so that no sub grid scale or turbulence model is necessary. The channel geometry was chosen so that finite Mach number effects can be isolated by comparing the present results to well established incompressible channel data. Here the fluid is assumed to be an ideal gas with constant specific heats, constant Prandtl number, and power-law temperature-dependent viscosity. Isothermal-wall boundary conditions are imposed so that a statistically stationary state may be obtained. The flow is driven by a uniform (in space) body force (rather than a mean pressure gradient) to preserve stream wise homogeneity, with the body force defined so that the total mass flux is constant.

Coleman, Gary N.↗

Development Status of the Carbon Dioxide and Moisture Removal Amine Swing-Bed System (CAMRAS)

Under a cooperative agreement with NASA, Hamilton Sundstrand has successfully designed, fabricated, tested and delivered three, state-of-the-art, solid amine prototype systems capable of continuous CO2 and humidity removal from a closed, habitable atmosphere. Two prototype systems (CAMRAS #1 and #2) incorporated a linear spool valve design for process flow control through the sorbent beds, with the third system (CAMRAS #3) employing a rotary valve assembly that improves system fluid interfaces and regeneration capabilities. The operational performance of CAMRAS #1 and #2 has been validated in a relevant environment, through both simulated human metabolic loads in a closed chamber and through human subject testing in a closed environment. Performance testing at Hamilton Sundstrand on CAMRAS #3, which incorporates a new valve and modified canister design, showed similar CO2 and humidity removal performance as CAMRAS #1 and #2, demonstrating that the system form can be modified within certain bounds with little to no effect in system function or performance. Demonstration of solid amine based CO2 and humidity control is an important milestone in developing this technology for human spaceflight. The systems have low power requirements; with power for air flow and periodic valve actuation and indication the sole requirements. Each system occupies the same space as roughly four shuttle non-regenerative LiOH canisters, but have essentially indefinite CO2 removal endurance provided a regeneration pathway is available. Using the solid amine based systems to control cabin humidity also eliminates the latent heat burden on cabin thermal control systems and the need for gas/liquid phase separation in a low gravity environment, resulting in additional simplification of vehicle environmental control and life support system process requirements.

Papale, William↗

Algorithmic information theory and the hidden variable question

The admissibility of certain nonlocal hidden-variable theories are explained via information theory. Consider a pair of Stern-Gerlach devices with fixed nonparallel orientations that periodically perform spin measurements on identically prepared pairs of electrons in the singlet spin state. Suppose the outcomes are recorded as binary strings l and r (with l sub n and r sub n denoting their n-length prefixes). The hidden-variable theories considered here require that there exists a recursive function which may be used to transform l sub n into r sub n for any n. This note demonstrates that such a theory cannot reproduce all the statistical predictions of quantum mechanics. Specifically, consider an ensemble of outcome pairs (l,r). From the associated probability measure, the Shannon entropies H sub n and H bar sub n for strings l sub n and pairs (l sub n, r sub n) may be formed. It is shown that such a theory requires that the absolute value of H bar sub n - H sub n be bounded - contrasting the quantum mechanical prediction that it grow with n.

Fuchs, Christopher↗

A Redetermination of the Dissociation Energy of MgO(+)

In 1986, we reported a dissociation energy (D(sub 0) of 2.31 eV for the X(sup 2)Pi ground state of MgO(+). This value was determined by computing the dissociation energy to the Mg(2+) + O(-) limit and adjusting the value to the Mg(+) + O limit using the experimental Ionization Potential (IP) of Mg(+) and the Electron Affinity (EA) of O. The success of this method relies on the assumption that there is little covalent contribution to the bonding. The very small (0.04 eV) correlation contribution to the binding energy was taken as corroboration for the validity of this approach. Our earlier theoretical value was estimated to be accurate to at least 0.2 eV. It is in excellent agreement with the subsequent value of 2.30 +/- 0.13 eV determined by Freiser and co-workers from photodissociation experiments. It is also consistent with the upper (less than 3.1 eV) and lower (greater than 1.1 eV) bounds determined by Rowe obtained by studying the reactions of Mg(+) with 03 and NO2. However, it is inconsistent with an upper bound of 1.7 eV reported by Kappes and Staley based on their failure to observe MgO(+) in the reaction of Mg(+) with N2O. The picture became somewhat clouded, however, by the recent guided-ion beam mass spectrometric studies of Dalleska and Armentrout. Their initial analysis of the reaction data for Mg(+) + O2 lead to a bond dissociation energy of 2.92 +/- 0.25 eV, which is considerably larger than the value of 2.47 +/- 0.06 eV deduced from their studies of the Mg(+)+NO2 reaction.

Bauschlicher, Charles W., Jr.↗

Thermal State-of-Charge in Solar Heat Receivers

A theoretical framework is developed to determine the so-called thermal state-of-charge (SOC) in solar heat receivers employing encapsulated phase change materials (PCMS) that undergo cyclic melting and freezing. The present problem is relevant to space solar dynamic power systems that would typically operate in low-Earth-orbit (LEO). The solar heat receiver is integrated into a closed-cycle Brayton engine that produces electric power during sunlight and eclipse periods of the orbit cycle. The concepts of available power and virtual source temperature, both on a finite-time basis, are used as the basis for determining the SOC. Analytic expressions for the available power crossing the aperture plane of the receiver, available power stored in the receiver, and available power delivered to the working fluid are derived, all of which are related to the SOC through measurable parameters. Lower and upper bounds on the SOC are proposed in order to delineate absolute limiting cases for a range of input parameters (orbital, geometric, etc.). SOC characterization is also performed in the subcooled, two-phase, and superheat regimes. Finally, a previously-developed physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) system is used in order to predict the SOC as a function of measurable parameters.

Hall, Carsie A., Jr.↗

Thermal State-of-Charge in Solar Heat Receivers

A theoretical framework is developed to determine the so-called thermal state-of-charge (SOC) in solar heat receivers employing encapsulated phase change materials (PCMS) that undergo cyclic melting and freezing. The present problem is relevant to space solar dynamic power systems that would typically operate in low-Earth-orbit (LEO). The solar heat receiver is integrated into a closed-cycle Brayton engine that produces electric power during sunlight and eclipse periods of the orbit cycle. The concepts of available power and virtual source temperature, both on a finite-time basis, are used as the basis for determining the SOC. Analytic expressions for the available power crossing the aperture plane of the receiver, available power stored in the receiver, and available power delivered to the working fluid are derived, all of which are related to the SOC through measurable parameters. Lower and upper bounds on the SOC are proposed in order to delineate absolute limiting cases for a range of input parameters (orbital, geometric, etc.). SOC characterization is also performed in the subcooled, two-phase, and superheat regimes. Finally, a previously-developed physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) system is used in order to predict the SOC as a function of measurable parameters.

Hall, Carsie, A., III↗

Uncertainty in Soil Moisture Retrievals: an Ensemble Approach Using SMOS L-Band Microwave Data

The uncertainty of soil moisture (SM) retrievals from satellite brightness temperature (TB) observations depends primarily on the choice of radiative transfer model (RTM) parameters, prior SM information and TB inputs. This paper studies the sensitivity of several (quasi-)operational and experimental SM retrieval products from the Soil Moisture Ocean Salinity (SMOS) mission to these choices at 11 reference sites, located in 7 watersheds across the United States (US). Different literature-based RTM parameter sets cause large biases between retrievals. Whereas typical RTM parameter sets are calibrated for SM retrievals, it is shown that a parameter set carefully optimized for TB forward modeling can also be used for retrieving SM. It is also shown that the inclusion of dynamic prior SM estimates in a Bayesian retrieval scheme can strongly improve SM retrievals, regardless of the choice of RTM parameters, and that the use of multi-angular and multi-polarization TB does not necessarily lead to superior retrievals compared to retrievals based on TB data at a single incidence angle and polarization. The second part of this paper evaluates ensemble uncertainty metrics for SM retrievals obtained by propagating a wide range of RTM parameters through the RTM. As expected for bounded variables, the spread in the ensemble SM retrievals is smallest for wet and dry SM values and highest for intermediate SM values. After removal of the strong long-term mean bias associated with the RTM parameter values for individual ensemble members, the remaining anomaly ensemble SM spread of 0.037 cu m/cu m approximates the actual time series unbiased root-mean-square-difference of 0.042 cu m/cu m between ensemble mean retrievals and in situ reference data across the reference sites. However, the temporal variability in the anomaly ensemble spread reveals higher-order biases in the retrieval error, which should be accounted for when characterizing retrieval error.

Jan Quets↗

Point to Point Control of the Hydrogen Mixer

We continue to build on the theoretical modelling of a liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem. The mixer described in this work is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. A small signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasureable but bounded additive step perturbations in the valve positions.

Barbieri, Enrique↗

Stellar Laboratories II. New Zn Iv and Zn v Oscillator Strengths and Their Validation in the Hot White Dwarfs G191-B2B and RE0503-289

Context. For the spectral analysis of high-resolution and high-signal-to-noise (SN) spectra of hot stars, state-of-the-art non-local thermodynamic equilibrium (NLTE) model atmospheres are mandatory. These are strongly dependent on the reliability of the atomic data that is used for their calculation. In a recent analysis of the ultraviolet (UV) spectrum of the DA-type white dwarf G191B2B,21 Zn iv lines were newly identified. Because of the lack of Zn iv data, transition probabilities of the isoelectronic Ge vi were adapted for a first, coarse determination of the photospheric Zn abundance.Aims. Reliable Zn iv and Zn v oscillator strengths are used to improve the Zn abundance determination and to identify more Zn lines in the spectra of G191B2B and the DO-type white dwarf RE 0503289. Methods. We performed new calculations of Zn iv and Zn v oscillator strengths to consider their radiative and collisional bound-bound transitions in detail in our NLTE stellar-atmosphere models for the analysis of the Zn iv v spectrum exhibited in high-resolution and high-SN UV observations of G191B2B and RE 0503289. Results. In the UV spectrum of G191B2B, we identify 31 Zn iv and 16 Zn v lines. Most of these are identified for the first time in any star. We can reproduce well almost all of them at log Zn 5.52 0.2 (mass fraction, about 1.7 times solar). In particular, the Zn iv Zn v ionization equilibrium, which is a very sensitive Teff indicator, is well reproduced with the previously determined Teff 60 000 2000 K and log g 7.60 0.05. In the spectrum of RE 0503289, we identified 128 Zn v lines for the first time and determined log Zn 3.57 0.2 (155 times solar). Conclusions. Reliable measurements and calculations of atomic data are a pre-requisite for stellar-atmosphere modeling. Observed Zn iv and Zn v line profiles in two white dwarf (G191B2B and RE 0503289) ultraviolet spectra were well reproduced with our newly calculated oscillator strengths. This allowed us to determine the photospheric Zn abundance of these two stars precisely.

strengths↗

Small Signal Point-to-Point Tracking of a Propellant Mixer

This paper addresses some theoretical modelling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties before it is fed into a test article. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions.

Barbieri, Enrique↗