Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic parameter estimation”

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,441 records · Page 80

The Bursting Behavior of 4U 1728-34: Parameters of a Neutron Star and Geometry of a NS-disk system

We analyze a set of Type I X-ray bursts from the low mass X-ray binary 4U 1728-34, observed with Rossi X-ray Timing Explorer (RXTE). We infer the dependence of the neutron star (NS) mass and radius with respect to the assumed distance to the system using an analytical model of X-ray burst spectral formation. The model behavior clearly indicates that the burster atmosphere is helium-dominated. Our results strongly favor the soft equation of state (EOS) of NS for 4U 1728-34. We find that distance to the source should be within 4.5-5.0 kpc range. We obtain rather narrow constrains for the NS radius in 8.7-9.7 km range and interval 1.2-1.6 Ma for NS mass for this particular distance range. We uncover a temporal behavior of red-shift corrected burst flux for the radial expansion episodes and we put forth a dynamical evolution scenario for the NS accretion disk geometry during which an expanded envelope affects the accretion disk and increases the area of the neutron star exposed to the Earth observer. In the framework of this scenario we provide a new method for the estimation of the inclination angle which leads to the value of approximately 50 degrees for 4U 1728-34.

Shaposhnikov, Nickolai↗

NASA Tech Briefs, May 2011

Topics covered include: 1) Method to Estimate the Dissolved Air Content in Hydraulic Fluid; 2) Method for Measuring Collimator-Pointing Sensitivity to Temperature Changes; 3) High-Temperature Thermometer Using Cr-Doped GdAlO3 Broadband Luminescence; 4)Metrology Arrangement for Measuring the Positions of Mirrors of a Submillimeter Telescope; 5) On-Wafer S-Parameter Measurements in the 325-508-GHz Band; 6) Reconfigurable Microwave Phase Delay Element for Frequency Reference and Phase-Shifter Applications; 7) High-Speed Isolation Board for Flight Hardware Testing; 8) High-Throughput, Adaptive FFT Architecture for FPGA-Based Spaceborne Data Processors; 9) 3D Orbit Visualization for Earth-Observing Missions; 10) MaROS: Web Visualization of Mars Orbiting and Landed Assets; 11) RAPID: Collaborative Commanding and Monitoring of Lunar Assets; 12) Image Segmentation, Registration, Compression, and Matching; 13) Image Calibration; 14) Rapid ISS Power Availability Simulator; 15) A Method of Strengthening Composite/Metal Joints; 16) Pre-Finishing of SiC for Optical Applications; 17) Optimization of Indium Bump Morphology for Improved Flip Chip Devices; 18) Measuring Moisture Levels in Graphite Epoxy Composite Sandwich Structures; 19) Marshall Convergent Spray Formulation Improvement for High Temperatures; 20) Real-Time Deposition Monitor for Ultrathin Conductive Films; 21) Optimized Li-Ion Electrolytes Containing Triphenyl Phosphate as a Flame-Retardant Additive; 22) Radiation-Resistant Hybrid Lotus Effect for Achieving Photoelectrocatalytic Self-Cleaning Anticontamination Coatings; 23) Improved, Low-Stress Economical Submerged Pipeline; 24) Optical Fiber Array Assemblies for Space Flight on the Lunar Reconnaissance Orbiter; 25) Local Leak Detection and Health Monitoring of Pressurized Tanks; 26) Dielectric Covered Planar Antennas at Submillimeter Wavelengths for Terahertz Imaging; 27) Automated Cryocooler Monitor and Control System; 28) Broadband Achromatic Phase Shifter for a Nulling Interferometer; 29) Super Dwarf Wheat for Growth in Confined Spaces; 30) Fine Guidance Sensing for Coronagraphic Observatories; 31) Single-Antenna Temperature- and Humidity-Sounding Microwave Receiver; 32) Multi-Wavelength, Multi-Beam, and Polarization-Sensitive Laser Transmitter for Surface Mapping; 33) Optical Communications Link to Airborne Transceiver; 34) Ascent Heating Thermal Analysis on Spacecraft Adaptor Fairings; 35) Entanglement in Self-Supervised Dynamics; 36) Prioritized LT Codes; 37) Fast Image Texture Classification Using Decision Trees; 38) Constraint Embedding Technique for Multibody System Dynamics; 39) Improved Systematic Pointing Error Model for the DSN Antennas; 40) Observability and Estimation of Distributed Space Systems via Local Information-Exchange Networks; 41) More-Accurate Model of Flows in Rocket Injectors; 42) In-Orbit Instrument-Pointing Calibration Using the Moon as a Target; 43) Reliability of Ceramic Column Grid Array Interconnect Packages Under Extreme Temperatures; 44) Six Degrees-of-Freedom Ascent Control for Small-Body Touch and Go; and 45) Optical-Path-Difference Linear Mechanism for the Panchromatic Fourier Transform Spectrometer.

Source record↗

Process design tools and techno-economic analysis for capacitive deionization

Capacitive deionization (CDI) devices use cyclical electrosorption on porous electrode surfaces to achieve water desalination. Process modeling and design of CDI systems requires accurate treatment of the coupling among input electrical forcing, input flow rates, and system responses including salt removal dynamics, water recovery, energy storage, and dissipation. Techno-economic analyses of CDI further require a method to calculate and compare between a produced commodity (e.g. desalted water) versus capital and operational costs of the system. In this work, we demonstrate a new modeling and analysis tool for CDI developed as an installable Matlab program that allows direct numerical simulation of CDI dynamics and calculation of key performance and cost parameters. The program is provided for free and is used to run open-source Simulink models. The Simulink environment sends information to the program and allows for a drag and drop design space where users can connect CDI cells to relevant periphery blocks such as grid energy, battery, solar panel, waste disposal, and maintenance/labor cost streams. The program allows for simulation of arbitrary current forcing and arbitrary flow rate forcing of one or more CDI cells. We employ validated well-mixed reactor formulations together with a non-linear circuit model formulation that can accommodate a variety of electric double layer sub-models (e.g. for charge efficiency). The program includes a graphical user interface (GUI) to specify CDI plant parameters, specify operating conditions, run individual tests or parameter batch-mode simulations, and plot relevant results. The techno-economic models convert among dimensional streams of species (e.g. feed, desalted water, and brine), energy, and cost and enable a variety of economic estimates including levelized water costs.

42 ENGINEERING↗

Dielectron production in central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02TeV

The first measurement of the 𝑒 + ⁢𝑒 − pair production at midrapidity and low invariant mass in central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV at the Large Hadron Collider is presented. The yield of 𝑒 + ⁢𝑒 − pairs is compared with a cocktail of expected hadronic decay contributions in the invariant mass (𝑚𝑒⁢𝑒) and pair transverse momentum (𝑝 T,𝑒⁢𝑒 ) ranges 𝑚 𝑒⁢𝑒 < 3.5 GeV/𝑐 2 and 𝑝 T,𝑒⁢𝑒 < 8 GeV/𝑐. For 0.18 < 𝑚 𝑒⁢𝑒 < 0.5 GeV/𝑐 2 the ratio of data to the cocktail of hadronic contributions amounts to 1.40 ± 0.11⁢(stat.) ± 0.23 ⁢(syst.) ± 0.16⁢ (cocktail) and 1.42 ± 0.11⁢ (stat.) ± 0.23 (syst.)$^{+0.24}_{−0.29}$ ⁢(cocktail), including or not including medium effects in the estimation of the heavy-flavor background, respectively. It is consistent with predictions from two different models for an additional contribution of thermal 𝑒+⁢𝑒− pairs from the hadronic and partonic phases. In the intermediate-mass range (1.2 < 𝑚 𝑒⁢𝑒 < 2.6 GeV/𝑐 2 ), the pair transverse impact parameter of the 𝑒 + ⁢𝑒 − pairs (DCA𝑒⁢𝑒, where “DCA” denotes “distance of closest approach”) is used for the first time in Pb-Pb collisions to separate displaced dielectrons from heavy-flavor hadron decays from a possible (thermal) contribution produced at the interaction point. The data are consistent with a suppression of 𝑒 + ⁢𝑒 − pairs from $𝑐\bar{𝑐}$ and an additional prompt component. Finally, the first direct-photon measurement in the 10% most central Pb-Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV is reported via the study of virtual direct photons in the transverse momentum range 1 < 𝑝 T < 5 GeV/𝑐. A model including prompt photons, as well as photons from the preequilibrium and fluid-dynamic phases, can reproduce the result, while being at the upper edge of the data uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗

A Sustainable Approach to Decision-Making in CCUS Systems

This study proposes a methodology for operational decision-making leading to the sustainability CO 2 -EOR systems. The methodology is a simple, yet comprehensive, integration of environmental and economic carbon capture, utilization and storage (CCUS) performance. Understanding the interplay between the environmental and economic conditions allows the EOR operator to manage the CO 2 -EOR flood so it remains within the limits of carbon balance neutrality and favorable economics under given prices and costs. The environmental performance and the main operational inputs (fluid production, CO 2 injection, and purchase requirements) were taken from a dynamic carbon Life Cycle Assessment (LCA) study previously performed for a CO 2 -EOR site in Cranfield, Mississippi. The LCA study accounts for greenhouse gases (GHG) emissions within a gate-to-grave CCUS system boundary, considering four different CO 2 injection strategies and four gas separation processes. The economic performance assessment was made through a marginalist analysis (Marginalist Production Theory), which determines important technical and economic relations through basic differential calculus derived from production, cost, and income functions. This method allows to define the operational economic optimum (optimum profit). The assessment considers four CO 2 price scenarios and uses World Bank estimated oil prices. In a first step of the proposed methodology, environmental limits of the operation are estimated. These limits are then compared with operative maximum productivity as well as with the economic optimum. According to the methodology, the EOR operation is sustainable as long as the environmental optimum is greater, or lasts longer, than the economic optimum along the project’s life. Results show that for this operation and assumed parameters, CO 2 -EOR could be classified as a sustainable activity depending on the operator’s strategic decision-making. In the case study, all injection strategies reach the environmental optimum after maximum productivity is reached. As a result, only two strategies (continuous gas injection and water alternating gas) find their economic optimum within the environmental limits, fulfilling the defined necessary and sufficient condition for CCUS sustainability.

02 PETROLEUM↗

Photometry of Outer Solar System Objects from the Dark Energy Survey. II. A Joint Analysis of Trans-Neptunian Absolute Magnitudes, Colors, Light Curves and Dynamics

For the 696 trans-Neptunian objects (TNOs) with absolute magnitudes 5.5 < H r < 8.2 detected in the Dark Energy Survey, we characterize the relationships between their dynamical state and physical properties—namely H r , indicating size; colors, indicating surface composition; and flux variation semiamplitude A, indicating asphericity and surface inhomogeneity. We seek “birth” physical distributions that can recreate these parameters in every dynamical class. We show that the observed colors of these TNOs are consistent with two Gaussian distributions in griz space, “near-infrared bright” (NIRB) and “near-infrared faint” (NIRF), presumably an inner and outer birth population, respectively. We find a model in which both the NIRB and NIRF H r and A distributions are independent of current dynamical states, supporting their assignment as birth populations. All objects are consistent with a common rolling p(H r ), but NIRF objects are significantly more variable. Cold classicals (CCs) are purely NIRF, while hot classical (HC), scattered, and detached TNOs are consistent with ≈ 70% NIRB and the resonance NIRB fractions show significant variation. The NIRB components of the HCs and of some resonances have broader inclination distributions than the NIRFs, i.e. their current dynamics retains information about birth location. We find evidence for radial stratification within the birth NIRB population, in that HC NIRBs are on average redder than detached or scattered NIRBs; a similar effect distinguishes CCs from other NIRFs. We estimate total object counts and masses of each class within our H r range. These results will strongly constrain models of the outer solar system.

79 ASTRONOMY AND ASTROPHYSICS↗

A Multi-Satellite Assimilation and Modeling Platform to Construct a Global and Complete view of the Hydrologic Cycle

What captivates me about the research field of satellite remote sensing of the earth system science is the opportunity to help human beings by addressing challenging questions related to the current and future availability of natural resources for domestic, agricultural, and industrial needs. For example, a continuously varying climate poses a threat to the water resource, not only in developing countries, but also in economically well-developed regions, creating the urgent need to characterize and predict its spatial and temporal availability. Specifically, my curiosity is driven by the challenge of merging cutting-edge space technology and earth observations (i.e., remote sensing) with state-of-the-art models for the purpose of improving our scientific knowledge about the variability and the change of the hydrologic cycle.The only practical way to observe the land surface processes (e.g., hydrologic cycle) on continental to global scales is via satellite remote sensing. Though remote sensing can make spatially comprehensive measurements of various components of the land surface system, it cannot provide direct information on the entire system, and the measurements represent only a snapshot in time. Physical based models may be used to continuously predict the temporal and spatial earth processes, but these predictions are often poor, due to model initialization, parameter and meteorological forcing errors, and inadequate model physics and/or resolution. Thus, an attractive prospect is to combine the strengths of land surface models and observations (and minimize the weaknesses of both) to provide a superior land surface state estimate. This is the goal of land surface data assimilation. Data assimilation provides a better estimate of the environmental states than either models or observations could individually do. The broad concept of data assimilation can be applied to various disciplines such as hydrology, ecology, environmental hazards, agriculture and economy.My research is targeted to integrating multiple satellite information having multi-sensor and multi-resolution assimilation within today's state-of-the-art hydrologic models. Multi-sensor and multi-resolution assimilation techniques represent a necessary milestone in future earth science applications because they offer the capability of comprehensively integrating disparate types of earth observations via deep-learning techniques. A multi-satellite assimilation and modeling platform will ultimately provide a robust and complete dataset to more fully understand earth's dynamics. In fact, the multi-satellite assimilation platforms build a comprehensive description of the earth processes that is geared toward accurately representing the complexity of natural and anthropogenic interactions in land surface processes. Most importantly, these platforms are suitable as decision support tools for applications across different aspects of earth system science.My recent work has focused on recent multi-sensor data assimilation techniques targeted at improving snow, soil moisture, groundwater, and terrestrial water storage hydrological states.

Girotto, Manuela↗

Estimating ionization states and continuum lowering from ab initio path integral Monte Carlo simulations for warm dense hydrogen

Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under these conditions, but can be difficult to diagnose since experimental campaigns are limited and often require model-dependent analysis. This is especially true for hydrogen, which has a comparably low scattering cross section, making high-quality data particularly difficult to obtain. Consequently, building equation of state tables often relies on simulations in combination with untested approximations to extract properties from experiments. Here, we investigate an approach for extracting the ionization potential depression and ionization degree—quantities which are otherwise not directly accessible from the physical model—from first-principles path integral Monte Carlo (PIMC) simulations utilizing a chemical model. In contrast to experimental measurements, where noise and nonequilibrium effects add to the uncertainty of the inferred parameters, PIMC simulations provide a clean signal with well-defined thermodynamic conditions. Comparisons against commonly used models show a qualitative agreement, but we find deviations primarily for the high-density and high-temperature cases. We also demonstrate the decreasing sensitivity of the dynamic structure factor with respect to both ionization and continuum lowering for increasing scattering angles in x-ray Thomson scattering experiments. Our work has important implications for the design of future experiments, but also offers qualitative understanding of structure factors and the imaginary-time correlation function obtained from first-principles quantum Monte Carlo simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bayesian Parameter Estimation of the k-ω Shear Stress Transport Model for Accurate Simulations of Impinging-Jet Heat Transfer

Given the lack of fusion-relevant component test facilities, current estimates of the thermo-fluid performance of plasma-facing components are based for the most part on numerical simulations. A major source of uncertainty in these simulations is the semiempirical turbulence (closure) models for the Reynolds stresses appearing in the governing Reynolds-averaged Navier-Stokes equations, which involve a set of constants that depend upon the flow. The objective of this study is to evaluate Bayesian parameter estimation of turbulence closure constants in ANSYS Fluent to model heat transfer in impinging jets. The Bayesian statistical calibration produces a probability distribution for these constants from experimental data; the maximum a posteriori estimates are then taken to be the calibrated constants, or parameters. The turbulence model constants are calibrated using an experimental study of a submerged jet of air impinging on a flat heated surface at Reynolds numbers Re = O(10 4 ) and impingement distance in jet diameters H/d = 2. Numerical predictions using the calibrated model parameters are then compared with those generated using the default constants. Predictions obtained with model parameters calibrated on datasets of two different sizes are compared to evaluate the effect of the number of calibration samples. Lastly, the extrapolative ability of the calibrated model is examined by predictions at a Re beyond the calibration values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantitative agreement between dynamical rocking curves in ultrafast electron diffraction for x-ray lasers

Electron diffraction through a thin patterned silicon membrane can be used to create complex spatial modulations in electron distributions. By precisely varying parameters such as crystallographic orientation and wafer thickness, the intensity of reflections in the diffraction plane can be controlled and by placing an aperture to block all but one spot, we can form an image with different parts of the patterned membrane, as is done for bright-field imaging in microscopy. The patterned electron beams can then be used to control phase and amplitude of subsequent x-ray emission, enabling novel coherent x-ray methods. Furthermore, the electrons themselves can also be used for femtosecond time resolved diffraction and microscopy. As a first step toward patterned beams, we demonstrate experimentally and through simulation the ability to accurately predict and control diffraction spot intensities. We simulate MeV transmission electron diffraction patterns using the multislice method for various crystallographic orientations of a single crystal Si(001) membrane near beam normal. The resulting intensity maps of the Bragg reflections are compared to experimental results obtained at the Accelerator Structure Test Area Ultrafast Electron Diffraction (ASTA UED) facility at SLAC. Furthermore, the fraction of inelastic and elastic scattering of the initial charge is estimated along with the absorption of the membrane to determine the contrast that would be seen in a patterned version of the Si(001) membrane.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Repulsive Tomonaga-Luttinger liquid in the quasi-one-dimensional alternating spin-$\frac{1}{2}$ antiferromagnet NaVOPO 4

Here, we probe the magnetic field-induced Tomonaga-Luttinger liquid (TLL) state in the bond-alternating spin1/2 antiferromagnetic (AFM) chain compound NaVOPO 4 using thermodynamic as well as local μSR and 31 P NMR probes down to mK temperatures in magnetic fields up to 14 T. The μSR and NMR relaxation rates in the gapless TLL regime decay slowly following characteristic power-law behavior, enabling us to directly determine the interaction parameter K as a function of the magnetic field. These estimates are crosschecked using magnetization and specific heat data. The field-dependent K lies in the range of 0.4 < K < 1 and indicates the repulsive nature of interactions between the spinless fermions, in line with the theoretical predictions. This renders NaVOPO 4 the first experimental realization of TLL with repulsive fermionic interactions in hitherto studied S = 1/2 bond-alternating AFM-AFM chain compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Physics and technology considerations for the deuterium–tritium fuel cycle and conditions for tritium fuel self sufficiency

The tritium aspects of the DT fuel cycle embody some of the most challenging feasibility and attractiveness issues in the development of fusion systems. The review and analyses in this paper provide important information to understand and quantify these challenges and to define the phase space of plasma physics and fusion technology parameters and features that must guide a serious R&D in the world fusion program. We focus in particular on components, issues and R&D necessary to satisfy three 'principal requirements': (1) achieving tritium self-sufficiency within the fusion system, (2) providing a tritium inventory for the initial start-up of a fusion facility, and (3) managing the safety and biological hazards of tritium. A primary conclusion is that the physics and technology state-of-the-art will not enable DEMO and future power plants to satisfy these principal requirements. We quantify goals and define specific areas and ideas for physics and technology R&D to meet these requirements. A powerful fuel cycle dynamics model was developed to calculate time-dependent tritium inventories and flow rates in all parts and components of the fuel cycle for different ranges of parameters and physics and technology conditions. Dynamics modeling analyses show that the key parameters affecting tritium inventories, tritium start-up inventory, and tritium self-sufficiency are the tritium burn fraction in the plasma (f b ), fueling efficiency (η f ), processing time of plasma exhaust in the inner fuel cycle (t p ), reactor availability factor (AF), reserve time (tr) which determines the reserve tritium inventory needed in the storage system in order to keep the plant operational for time t r in case of any malfunction of any part of the tritium processing system, and the doubling time (t d ). Results show that η f f b > 2% and processing time of 1–4 h are required to achieve tritium self-sufficiency with reasonable confidence. For η f f b = 2% and processing time of 4 h, the tritium start-up inventory required for a 3 GW fusion reactor is ~11 kg, while it is <5 kg if η f f b = 5% and the processing time is 1 h. To achieve these stringent requirements, a serious R&D program in physics and technology is necessary. The EU-DEMO direct internal recycling concept that carries fuel directly from the plasma exhaust gas to the fueling systems without going through the isotope separation system reduces the overall processing time and tritium inventories and has positive effects on the required tritium breeding ratio (TBR R ). A significant finding is the strong dependence of tritium self-sufficiency on the reactor availability factor. Simulations show that tritium self-sufficiency is: impossible if AF < 10% for any η f f b , possible if AF > 30% and 1% ≤ η f f b ≤ 2%, and achievable with reasonable confidence if AF > 50% and η f f b > 2%. These results are of particular concern in light of the low availability factor predicted for the near-term plasma-based experimental facilities (e.g. FNSF, VNS, CTF), and can have repercussions on tritium economy in DEMO reactors as well, unless significant advancements in RAMI are made. There is a linear dependency between the tritium start-up inventory and the fusion power. The required tritium start-up inventory for a fusion facility of 100 MW fusion power is as small as 1 kg. Since fusion power plants will have large powers for better economics, it is important to maintain a 'reserve' tritium inventory in the tritium storage system to continue to fuel the plasma and avoid plant shutdown in case of malfunctions of some parts of the tritium processing lines. But our results show that a reserve time as short as 24 h leads to unacceptable reserve and start-up inventory requirements. Therefore, high reliability and fast maintainability of all components in the fuel cycle are necessary in order to avoid the need for storing reserve tritium inventory sufficient for continued fusion facility operation for more than a few hours. The physics aspects of plasma fueling, tritium burn fraction, and particle and power exhaust are highly interrelated and complex, and predictions for DEMO and power reactors are highly uncertain because of lack of experiments with burning plasma. Fueling by pellet injection on the high field side of tokamak has evolved to be the preferred method to fuel a burning plasma. Extrapolation from the DIII-D penetration scaling shows fueling efficiency expected in DEMO to be <25%, but such extrapolations are highly uncertain. The fueling efficiency of gas in a reactor relevant regime is expected to be extremely poor and not very useful for getting tritium into the core plasma efficiently. Gas fueling will nonetheless be useful for feedback control of the divertor operating parameters. Extensive modeling has been carried out to predict burn fraction, fueling requirements, and fueling efficiency for ITER, DEMO, and beyond. The fueling rate required to operate Q = 10 ITER plasmas in order to provide the required core fueling, helium exhaust and radiative divertor plasma conditions for acceptable divertor power loads was calculated. If this fueling is performed with a 50–50 DT mix, the tritium burn fraction in ITER would be ~0.36%, which is too low to satisfy the self-sufficiency conditions derived from the dynamics modeling for fusion reactors. Extrapolation to DEMO using this approach would also yield similarly low burn fraction. Extensive analysis presented shows that specific features of edge neutral dynamics in ITER and fusion reactors, which are different from present experiments, open possibilities for optimization of tritium fueling and thus to improve the burn fraction. Using only tritium in pellet fueling of the plasma core, and only deuterium for edge density, divertor power load and ELM control results in significant increase of the burn fraction to 1.8–3.6%. These estimates are performed with physics models whose results cannot be fully validated for ITER and DEMO plasma conditions since these cannot be achieved in present tokamak experiments. Thus, several uncertainties remain regarding particle transport and scenario requirements in ITER and DEMO. The safety standard requirements for protection of the public and release guidelines for tritium have been reviewed. General safety approaches including minimizing tritium inventories, reducing tritium permeation through materials, and decontaminating material for waste disposal have been suggested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Bayesian Approach to Estimating Coupling Between Neural Components: Evaluation of the Multiple Component, Event-Related Potential (mcERP) Algorithm

Accurate measurement of single-trial responses is key to a definitive use of complex electromagnetic and hemodynamic measurements in the investigation of brain dynamics. We developed the multiple component, Event-Related Potential (mcERP) approach to single-trial response estimation. To improve our resolution of dynamic interactions between neuronal ensembles located in different layers within a cortical region and/or in different cortical regions. The mcERP model assets that multiple components defined as stereotypic waveforms comprise the stimulus-evoked response and that these components may vary in amplitude and latency from trial to trial. Maximum a posteriori (MAP) solutions for the model are obtained by iterating a set of equations derived from the posterior probability. Our first goal was to use the ANTWERP algorithm to analyze interactions (specifically latency and amplitude correlation) between responses in different layers within a cortical region. Thus, we evaluated the model by applying the algorithm to synthetic data containing two correlated local components and one independent far-field component. Three cases were considered: the local components were correlated by an interaction in their single-trial amplitudes, by an interaction in their single-trial latencies, or by an interaction in both amplitude and latency. We then analyzed the accuracy with which the algorithm estimated the component waveshapes and the single-trial parameters as a function of the linearity of each of these relationships. Extensions of these analyses to real data are discussed as well as ongoing work to incorporate more detailed prior information.

Shah, Ankoor S.↗

SMOS/SMAP Synergy for SMAP Level 2 Soil Moisture Algorithm Evaluation

Soil Moisture Active Passive (SMAP) satellite has been proposed to provide global measurements of soil moisture and land freeze/thaw state at 10 km and 3 km resolutions, respectively. SMAP would also provide a radiometer-only soil moisture product at 40-km spatial resolution. This product and the supporting brightness temperature observations are common to both SMAP and European Space Agency's Soil Moisture and Ocean Salinity (SMOS) mission. As a result, there are opportunities for synergies between the two missions. These include exploiting the data for calibration and validation and establishing longer term L-band brightness temperature and derived soil moisture products. In this investigation we will be using SMOS brightness temperature, ancillary data, and soil moisture products to develop and evaluate a candidate SMAP L2 passive soil moisture retrieval algorithm. This work will begin with evaluations based on the SMOS product grids and ancillary data sets and transition to those that will be used by SMAP. An important step in this analysis is reprocessing the multiple incidence angle observations provided by SMOS to a global brightness temperature product that simulates the constant 40 degree incidence angle observations that SMAP will provide. The reprocessed brightness temperature data provide a basis for evaluating different SMAP algorithm alternatives. Several algorithms are being considered for the SMAP radiometer-only soil moisture retrieval. In this first phase, we utilized only the Single Channel Algorithm (SCA), which is based on the radiative transfer equation and uses the channel that is most sensitive to soil moisture (H-pol). Brightness temperature is corrected sequentially for the effects of temperature, vegetation, roughness (dynamic ancillary data sets) and soil texture (static ancillary data set). European Centre for Medium-Range Weather Forecasts (ECMWF) estimates of soil temperature for the top layer (as provided as part of the SMOS ancillary data) were used to correct for surface temperature effects and to derive microwave emissivity. ECMWF data were also used for precipitation forecasts, presence of snow, and frozen ground. Vegetation options are described below. One year of soil moisture observations from a set of four watersheds in the U.S. were used to evaluate four different retrieval methodologies: (1) SMOS soil moisture estimates (version 400), (2) SeA soil moisture estimates using the SMOS/SMAP data with SMOS estimated vegetation optical depth, which is part of the SMOS level 2 product, (3) SeA soil moisture estimates using the SMOS/SMAP data and the MODIS-based vegetation climatology data, and (4) SeA soil moisture estimates using the SMOS/SMAP data and actual MODIS observations. The use of SMOS real-world global microwave observations and the analyses described here will help in the development and selection of different land surface parameters and ancillary observations needed for the SMAP soil moisture algorithms. These investigations will greatly improve the quality and reliability of this SMAP product at launch.

Bindlish, Rajat↗

Kalman Filtering with Inequality Constraints for Turbofan Engine Health Estimation

Kalman filters are often used to estimate the state variables of a dynamic system. However, in the application of Kalman filters some known signal information is often either ignored or dealt with heuristically. For instance, state variable constraints (which may be based on physical considerations) are often neglected because they do not fit easily into the structure of the Kalman filter. This paper develops two analytic methods of incorporating state variable inequality constraints in the Kalman filter. The first method is a general technique of using hard constraints to enforce inequalities on the state variable estimates. The resultant filter is a combination of a standard Kalman filter and a quadratic programming problem. The second method uses soft constraints to estimate state variables that are known to vary slowly with time. (Soft constraints are constraints that are required to be approximately satisfied rather than exactly satisfied.) The incorporation of state variable constraints increases the computational effort of the filter but significantly improves its estimation accuracy. The improvement is proven theoretically and shown via simulation results. The use of the algorithm is demonstrated on a linearized simulation of a turbofan engine to estimate health parameters. The turbofan engine model contains 16 state variables, 12 measurements, and 8 component health parameters. It is shown that the new algorithms provide improved performance in this example over unconstrained Kalman filtering.

Simon, Dan↗

Response properties of pigeon otolith afferents to linear acceleration

In the present study, the sensitivity to sinusoidal linear accelerations in the plane of the utricular macula was tested in afferents. The head orientation relative to the translation axis was varied in order to determine the head position that elicited the maximal and minimal responses for each afferent. The response gain and phase values obtained to 0.5-Hz and 2-Hz linear acceleration stimuli were then plotted as a function of head orientation and a modified cosine function was fit to the data. From the best-fit cosine function, the predicted head orientations that would produce the maximal and minimal response gains were estimated. The estimated maximum response gains to linear acceleration in the utricular plane for the afferents varied between 75 and 1420 spikes s-1 g-1. The mean maximal gains for all afferents to 0.5-Hz and 2-Hz sinusoidal linear acceleration stimuli were 282 and 367 spikes s-1 g-1, respectively. The minimal response gains were essentially zero for most units. The response phases always led linear acceleration and remained constant for each afferent, regardless of head orientation. These response characteristics indicate that otolith afferents are cosine tuned and behave as one-dimensional linear accelerometers. The directions of maximal sensitivity to linear acceleration for the afferents varied throughout the plane of the utricle; however, most vectors were directed out of the opposite ear near the interaural axis. The response dynamics of the afferents were tested using stimulus frequencies ranging between 0.25 Hz and 10 Hz (0.1 g peak acceleration). Across stimulus frequencies, most afferents had increasing gains and constant phase values. These dynamic properties for individual afferents were fit with a simple transfer function that included three parameters: a mechanical time constant, a gain constant, and a fractional order distributed adaptation operator.

Non-NASA Center↗

Spectroscopy of GR1 centers in synthetic diamonds

We report on the spectroscopic characterization, absorption saturation, and excited-state dynamics of GR1 centers in synthetic diamonds. The non-linear optical measurements reveal an efficient bleaching of the GR1 center’s ground level under ns-pulsed 633 nm excitation. The maxima of absorption and emission cross sections were estimated to be 4.5 × 10 −17 cm 2 and 9 × 10 −17 cm 2 at 630 nm and 780 nm, respectively. The radiative lifetime of the excited level was estimated to be 8.5 ns. The 658 nm probe kinetics uncovered relaxation of 1 T 2 excited level going predominantly to a metastable state with a lifetime of 220 μs. An induced absorption detected with the use of a highly concentrated diamond sample could be due to up-conversion and photoionization processes in additional impurity-vacancy center with absorption at the short-wavelength tail of GR1. The results presented here indicate that synthetic diamonds with GR1 centers could serve as near infrared gain media or passive Q-switchers for laser cavities over the 633–750 nm spectral range. Optimization of the center concentration and the crystal parameters are required to minimize the induced absorption at the lasing wavelengths.

36 MATERIALS SCIENCE↗