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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Data Assimilation for Robust UQ Within Agent-Based Simulation on HPC Systems

Agent-based simulation provides a powerful tool for in silico system modeling. However, these simulations do not provide built-in methods for uncertainty quantification (UQ). Within these types of models a typical approach to UQ is to run multiple realizations of the model then compute aggregate statistics. This approach is limited due to the compute time required for a solution. When faced with an emerging biothreat, public health decisions need to be made quickly and solutions for integrating near real-time data with analytic tools are needed. We propose an integrated Bayesian UQ framework for agent-based models based on sequential Monte Carlo sampling. Given streaming or static data about the evolution of an emerging pathogen this Bayesian framework provides a distribution over the parameters governing the spread of a disease through a population. These estimates of the spread of a disease may be provided to public health agencies seeking to abate the spread. By coupling agent-based simulations with Bayesian modeling in a data assimilation, our proposed framework provides a powerful tool for modeling dynamical systems in silico. We propose a method which reduces model error and provides a range of realistic possible outcomes. Moreover, our method addresses two primary limitations of ABMs: the lack of UQ and an inability to assimilate data. Our proposed framework combines the flexibility of an agent-based model with UQ provided by the Bayesian paradigm in a workflow which scales well to HPC systems. We provide algorithmic details and results on a simulated outbreak with both static and streaming data.

Spannaus, Adam [ORNL] (ORCID:0000000225213657)↗

DOE OSTI2015

Consequences of bounds on longitudinal emittance growth for the design of recirculating linear accelerators

Recirculating linear accelerators (RLAs) are a cost-effective method for the acceleration of muons for a muon collider in energy ranges from a couple GeV to a few 10s of GeV. Muon beams generally have longitudinal emittances that are large for the RF frequency that is used, and it is important to limit the growth of that longitudinal emittance. This has particular consequences for the arc design of the RLAs. I estimate the longitudinal emittance growth in an RLA arising from the RF nonlinearity. Given an emittance growth limitation and other design parameters, one can then compute the maximum momentum compaction in the arcs. I describe how to obtain an approximate arc design satisfying these requirements based on the deisgn in [1]. Longitudinal dynamics also determine the energy spread in the beam, and this has consequences on the transverse phase advance in the linac. This in turn has consequences for the arc design due to the need to match beta functions. I combine these considerations to discuss design parameters for the acceleration of muons for a collider in an RLA from 5 to 63 GeV.

43 PARTICLE ACCELERATORS↗

DOE OSTI2016

Determination of crystal growth rates during rapid solidification of polycrystalline aluminum by nano-scale spatio-temporal resolution in situ transmission electron microscopy

In situ investigations of rapid solidification in polycrystalline Al thin films were conducted using nano-scale spatio-temporal resolution dynamic transmission electron microscopy. Differences in crystal growth rates and asymmetries in melt pool development were observed as the heat extraction geometry was varied by controlling the proximity of the laser-pulse irradiation and the associated induced melt pools to the edge of the transmission electron microscopy support grid, which acts as a large heat sink. Experimental parameters have been established to maximize the reproducibility of the material response to the laser-pulse-related heating and to ensure that observations of the dynamical behavior of the metal are free from artifacts, leading to accurate interpretations and quantifiable measurements with improved precision. Interface migration rate measurements revealed solidification velocities that increased consistently from ∼1.3 m s{sup −1} to ∼2.5 m s{sup −1} during the rapid solidification process of the Al thin films. Under the influence of an additional large heat sink, increased crystal growth rates as high as 3.3 m s{sup −1} have been measured. The in situ experiments also provided evidence for development of a partially melted, two-phase region prior to the onset of rapid solidification facilitated crystal growth. Using the experimental observations and associated measurements as benchmarks, finite-element modeling based calculations of the melt pool evolution after pulsed laser irradiation have been performed to obtain estimates of the temperature evolution in the thin films.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2012

Dynamic Behavior of Sand: Annual Report FY 11

Currently, design of earth-penetrating munitions relies heavily on empirical relationships to estimate behavior, making it difficult to design novel munitions or address novel target situations without expensive and time-consuming full-scale testing with relevant system and target characteristics. Enhancing design through numerical studies and modeling could help reduce the extent and duration of full-scale testing if the models have enough fidelity to capture all of the relevant parameters. This can be separated into three distinct problems: that of the penetrator structural and component response, that of the target response, and that of the coupling between the two. This project focuses on enhancing understanding of the target response, specifically granular geomaterials, where the temporal and spatial multi-scale nature of the material controls its response. As part of the overarching goal of developing computational capabilities to predict the performance of conventional earth-penetrating weapons, this project focuses specifically on developing new models and numerical capabilities for modeling sand response in ALE3D. There is general recognition that granular materials behave in a manner that defies conventional continuum approaches which rely on response locality and which degrade in the presence of strong response nonlinearities, localization, and phase gradients. There are many numerical tools available to address parts of the problem. However, to enhance modeling capability, this project is pursuing a bottom-up approach of building constitutive models from higher fidelity, smaller spatial scale simulations (rather than from macro-scale observations of physical behavior as is traditionally employed) that are being augmented to address the unique challenges of mesoscale modeling of dynamically loaded granular materials. Through understanding response and sensitivity at the grain-scale, it is expected that better reduced order representations of response can be formulated at the continuum scale as illustrated in Figure 1 and Figure 2. The final result of this project is to implement such reduced order models in the ALE3D material library for general use.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Data Synthesis and Data Assimilation at Global Change Experiments and Fluxnet Toward Improving Land Process Models

The project was conducted during the period from 7/1/2012 to 6/30/2017 with three major tasks: (1) data synthesis and development of data assimilation (DA) techniques to constrain modeled ecosystem feedback to climate change; (2) applications of DA techniques to improve process models at different scales from ecosystem to regions and the globe; and 3) improvements of modeling soil carbon (C) dynamics by land surface models. During this period, we have synthesized published data from soil incubation experiments (e.g., Chen et al., 2016; Xu et al., 2016; Feng et al., 2016), global change experiments (e.g., Li et al., 2013; Shi et al., 2015, 2016; Liang et al., 2016) and fluxnet (e.g., Niu et al., 2012., Xia et al., 2015; Li et al., 2016). These data have been organized into multiple data products and have been used to identify general mechanisms and estimate parameters for model improvement. We used the data sets that we collected and the DA techniques to improve model performance of both ecosystem models and global land models. The objectives are: 1) to improve model simulations of litter and soil carbon storage (e.g., Schädel et al., 2013; Hararuk and Luo, 2014; Hararuk et al., 2014; Liang et al., 2015); 2) to explore the effects of CO 2 , warming and precipitation on ecosystem processes (e.g., van Groenigen et al., 2014; Shi et al., 2015, 2016; Feng et al., 2017); and 3) to estimate parameters variability in different ecosystems (e.g., Li et al., 2016). We developed a traceability framework, which was based on matrix approaches and decomposed the modeled steady-state terrestrial ecosystem carbon storage capacity into four can trace the difference in ecosystem carbon storage capacity among different biomes to four traceable components: net primary productivity (NPP), baseline C residence times, environmental scalars and climate forcing (Xia et al., 2013). With this framework, we can diagnose the differences in modeled carbon storage across ecosystems, biomes, and models. This framework has been successfully implemented by several global land models, such as CABLE (Xia et al., 2013), LPJ-GUESS (Ahlström et al., 2015), CLM (Hararuk et al., 2014; Huang et al., 2017, submitted; Shi et al., 2017, submitted), and ORCHIDEE (Huang et al., 2017, unpublished). Moreover, we have identified the theoretical foundation of the determinants of transient C storage dynamics by adding another term, C storage potential, to the steady-state traceability framework (Luo et al., 2017). The theoretical foundation of transient C storage dynamics has been applied to develop a transient traceability framework to explore the traceable components of transient C storage dynamics responded to the rising CO 2 and climate change in the two contrasting ecosystem types Duke needleleaved forest and Harvard deciduous broadleaved forest (Jiang et al., 2017, in revision). Overall, with the data synthesis, data assimilation techniques, and the steady-state and transient traceability frameworks, we have greatly improved land process models for predicting responses and feedback of terrestrial C dynamics to global change. The matrix approaches has the potential to be applied in theoretical research on nitrogen and phosphorus cycle, and therefore, the coupling of carbon-nitrogen-phosphorus.

58 GEOSCIENCES↗

DOE OSTI2024

Aboveground Biomass Estimation Using NISAR Simulated ALOS-2 Time Series Data

Aboveground biomass (AGB) is a critical parameter to better understand the global carbon cycle and to develop sustainable forest management. However, a large uncertainty prevails. L-band SAR data have demonstrated strong potential to accurately retrieve AGB over low-biomass regions (<100 Mg ha-1). The upcoming NASA-ISRO Synthetic Aperture Radar mission will collect data at L- and S-band over earth’s landmass with a repeat period of 12 days, allowing us to have ample data for monitoring biomass and its dynamics. One of the key science requirements of the mission is to produce annual AGB maps at 1-ha resolution with RMS accuracy of 20 Mg/ha for 80 percentage of area over low-biomass regions in Calibration/Validation sites. The NISAR biomass algorithm will generate AGB maps based on the parameterization of semi-empirical model along with NISAR time-series dual pol data (HH and HV). To calibrate and validate the model for mission requirements, the mission will use reference estimates of AGB produced from ground inventory plots and airborne LiDAR data collected over selected sites distributed across different global ecoregions. This paper presents the initial results of the calibration/validation of the NISAR AGB retrieval algorithm over the Lenoir Landing (LENO), Alabama, USA site using NISAR simulated ALOS-2 time series data. Five multi-temporal dual-pol HH and HV NISAR Simulated ALOS 2 data collections were used as input to assess the performance of the model. The model AGB retrieval results shows that the NISAR model was able to achieve RMS accuracy within 20 Mg/ha.

Ramachandran, Naveen [Jet Propulsion Laboratory, C↗

DOE OSTI2015

Image-derived and arterial blood sampled input functions for quantitative PET imaging of the angiotensin II subtype 1 receptor in the kidney

Purpose: The radioligand {sup 11}C-KR31173 has been introduced for positron emission tomography (PET) imaging of the angiotensin II subtype 1 receptor in the kidney in vivo. To study the biokinetics of {sup 11}C-KR31173 with a compartmental model, the input function is needed. Collection and analysis of arterial blood samples are the established approach to obtain the input function but they are not feasible in patients with renal diseases. The goal of this study was to develop a quantitative technique that can provide an accurate image-derived input function (ID-IF) to replace the conventional invasive arterial sampling and test the method in pigs with the goal of translation into human studies. Methods: The experimental animals were injected with [{sup 11}C]KR31173 and scanned up to 90 min with dynamic PET. Arterial blood samples were collected for the artery derived input function (AD-IF) and used as a gold standard for ID-IF. Before PET, magnetic resonance angiography of the kidneys was obtained to provide the anatomical information required for derivation of the recovery coefficients in the abdominal aorta, a requirement for partial volume correction of the ID-IF. Different image reconstruction methods, filtered back projection (FBP) and ordered subset expectation maximization (OS-EM), were investigated for the best trade-off between bias and variance of the ID-IF. The effects of kidney uptakes on the quantitative accuracy of ID-IF were also studied. Biological variables such as red blood cell binding and radioligand metabolism were also taken into consideration. A single blood sample was used for calibration in the later phase of the input function. Results: In the first 2 min after injection, the OS-EM based ID-IF was found to be biased, and the bias was found to be induced by the kidney uptake. No such bias was found with the FBP based image reconstruction method. However, the OS-EM based image reconstruction was found to reduce variance in the subsequent phase of the ID-IF. The combined use of FBP and OS-EM resulted in reduced bias and noise. After performing all the necessary corrections, the areas under the curves (AUCs) of the AD-IF were close to that of the AD-IF (average AUC ratio =1 ± 0.08) during the early phase. When applied in a two-tissue-compartmental kinetic model, the average difference between the estimated model parameters from ID-IF and AD-IF was 10% which was within the error of the estimation method. Conclusions: The bias of radioligand concentration in the aorta from the OS-EM image reconstruction is significantly affected by radioligand uptake in the adjacent kidney and cannot be neglected for quantitative evaluation. With careful calibrations and corrections, the ID-IF derived from quantitative dynamic PET images can be used as the input function of the compartmental model to quantify the renal kinetics of {sup 11}C-KR31173 in experimental animals and the authors intend to evaluate this method in future human studies.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

Hot Deformation and Processing Maps of B{sub 4}C/6061Al Nanocomposites Fabricated by Spark Plasma Sintering

We developed 3 vol.% B{sub 4}C/6061Al nanocomposites by spark plasma sintering. We then characterized their flow and hot deformation behaviors using isothermal compression tests at deformation temperatures of 673-823 K and strain rates of 0.001-1 s{sup −1}. Results indicated that their steady flow characteristic can be described by the Zener–Hollomon parameter (Z) as a hyperbolic-sine function of flow stress. Based on pinning effect from B{sub 4}C nanoparticles, the estimated activation energy was calculated to be 274.80 kJ/mol. At the unstable deformation condition (723 K/1 s{sup −1}) with high Z, flow localization appeared; the main softening mechanisms involved are dynamic recovery and dynamic recrystallization (DRX). At the stable deformation condition (823 K/0.001 s{sup −1}) with low Z, adequate DRX occurred. Based on the processing map and microstructure evolution, the optimum hot working conditions for B{sub 4}C/6061Al nanocomposites were determined to be a deformation temperature of 773-823 K and a strain rate of 0.001-0.003 s{sup −1}. These working conditions were verified by a hot extrusion test.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Multiscale Modeling of a Direct Nonoxidative Methane Dehydroaromatization Reactor with a Validated Model for Catalyst Deactivation

Due to the recent boom in shale gas production, aromatics production using direct nonoxidative methane dehydroaromatization (DHA) is being investigated extensively. However, due to rapid coke formation, catalysts in the nonoxidative methane DHA reactors get deactivated, which is one of the critical issues for the commercial success of the methane DHA process. In this paper, a model for catalyst deactivation is developed. Rate models for other DHA reactions are developed by considering the decrease in the catalyst activity with time. Due to the very fast coke formation rate on the fresh catalyst, there is coke formation immediately upon the introduction of the feed. Therefore, an algorithm is developed for estimation of the initial state of the reactor and the kinetic parameters by coupling an iterative direct substitution approach with an optimization approach. Transient experimental data from an in-house reactor are first reconciled and then used for developing the kinetic model including the coke formation model. Using the rate model, a dynamic, heterogeneous, multiscale reactor model with embedded heating is developed. Here, the model couples the catalyst pellet level model with a reactor level model. Impacts of temperature, L/D ratio, and scheduling of reactors on variability in conversion and yield with time are studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Quantification of the Mass and Viscoelasticity of Interfacial Films on Tin Anodes Using EQCM-D

Electrochemical quartz crystal microbalance coupled with dissipation (EQCM-D) is employed to explore the solid electrolyte interphase (SEI) formation and Li insertion/deinsertion into thin film electrodes of tin. Due to the frequency change we find that the initial SEI formation process is rapid before Li insertion but varies significantly with increasing concentration of the additive fluoroethylene carbonate (FEC) in the electrolyte. The extent of dissipation, which represents the film rigidity, increases with cycle number, reflecting film thickening and softening. Dissipation values are almost twice as large in the baseline electrolyte (1.2 M LiPF 6 in 3:7 wt % ethylene carbonate:ethyl methyl carbonate), indicating the film in baseline electrolyte is roughly twice as soft as in the FEC-containing cells. More importantly, we detail how quantitative data about mass, thickness, shear elastic modulus, and shear viscosity in a time-resolved manner can be obtained from the EQCM-D response. These parameters were extracted from the frequency and dissipation results at multiple harmonics using the Sauerbrey and Voigt viscoelastic models. From these modeled results we show the dynamic mass changes for each half cycle. We also demonstrate that different amounts of FEC additive influence the SEI formation behavior and result in differences in the estimated mass, shear modulus and viscosity. After three cycles, the film in baseline electrolyte exhibits a 1.2 times larger mass change compared with the film in the FEC-containing electrolyte. The shear elastic modulus of films formed in the presence of FEC is larger than in the baseline electrolyte at early stages of lithiation. Also with lithiation is a marked increase in film viscosity, which together point to a much stiffer and more homogeneous SEI formed in the presence of FEC.

36 MATERIALS SCIENCE↗

DOE OSTI2019

PMU-Based Monitoring and Mitigation of Delayed Voltage Recovery Using Admittances

This project analyzes the delayed voltage recovery phenomenon by simplifying the Western Electricity Coordinating Council (WECC) composite load model using the load admittance and the thermal relay dynamics. From this analysis, a closed form expression approximating the recovery time is derived and the key load parameters impacting the behavior of voltage recovery are identified. A monitoring scheme based on the measured load admittance is then proposed to detect the onset of stalling, even in the presence of voltage oscillations, and estimate the duration of the delayed voltage recovery. A mitigation scheme utilizing smart thermostats and offline learning is also derived to ensure that the voltage recovers to the pre-contingency voltage within a specified time. Both the monitoring and mitigation schemes only need local measurements at a substation making them promising for online applications. Results for the monitoring and mitigation schemes are described in detail for the IEEE 162 bus system validating the various assumptions used for the analysis and establishing the connection between the delayed voltage recovery phenomenon and load admittance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2019

Methodology for Estimating Safe Charge Limits of Flammable Refrigerants in HVAC&R Applications - Part 2

This report provides a detailed description of the development of a reduced order model (ROM) that could be used for estimation of safe charge limits based on a computational fluid dynamics (CFD) parametric study. The ROM was developed to simulate leak events from a single side-wall duct penetration into a room and provide estimates of refrigerant concentration together with flammable volume fraction (FVF), accumulated fuel volume (AFV) and other outputs vs. time. The ROM is based on a range of seven parameters listed in Table ES1 and assumes that the AC unit air circulation blower is off. A total of 586 CFD “training” cases were run and used as the basis for development and training or calibration of the ROM. The regression analysis used to develop the ROM was performed in Tasmanian ( T oolkit for A daptive S tochastic M odeling A nd N on- I ntrusive A pproximatio N ). This is a tool developed at ORNL for high-dimensional integration, interpolation and parameter calibration. Comparison of the ROM vs. several CFD test cases (at points not used in the training set) showed maximum absolute and root mean square errors (RMSE) for FVF of 1.32% and 0.5%, respectively. Initial CFD runs were conducted to gauge the sensitivity of the results to each parameter. Based on the findings, the numbers of design of experiment (DoE) levels for each parameter was determined. Then, for each parameter, Tasmanian calculated the value of each level and defined a DoE of 586 unique CFD simulation cases.

42 ENGINEERING↗

DOE OSTI2020

Implementing an aerosol dynamic shape factor for cerium oxide powder

The shape factor (X B ) of bulk powder CeO 2 (cerium oxide) particles is estimated to be X B = 1.34 (dimensionless). This is defined in relation to the measured aerated bulk density (0.91 g cm -3 , Hosokawa Micron, Summit NJ, 2018) of the CeO 2 powder. This information is used to determine the RRF (respirable release factor), and other parameters during drop test experiments of nuclear material storage containers at Los Alamos National Laboratory. The calculation is based on an experimental shape factor (X P = 2.66 dimensionless) for CeO 2 particles of the same (chain aggregate) morphology studied by other researchers. In that work, the shape factor was defined in terms of the maximum possible (compacted) density of CeO 2 material (7.13 g cm -3 ).

61 RADIATION PROTECTION AND DOSIMETRY↗

DOE OSTI2012

Dynamic recrystallization in friction surfaced austenitic stainless steel coatings

Friction surfacing involves complex thermo-mechanical phenomena. In this study, the nature of dynamic recrystallization in friction surfaced austenitic stainless steel AISI 316L coatings was investigated using electron backscattered diffraction and transmission electron microscopy. The results show that the alloy 316L undergoes discontinuous dynamic recrystallization under conditions of moderate Zener-Hollomon parameter during friction surfacing. - Highlights: Black-Right-Pointing-Pointer Dynamic recrystallization in alloy 316L friction surfaced coatings is examined. Black-Right-Pointing-Pointer Friction surfacing leads to discontinuous dynamic recrystallization in alloy 316L. Black-Right-Pointing-Pointer Strain rates in friction surfacing exceed 400 s{sup -1}. Black-Right-Pointing-Pointer Estimated grain size matches well with experimental observations in 316L coatings.

36 MATERIALS SCIENCE↗

DOE OSTI2015

DECIPHERING THERMAL PHASE CURVES OF DRY, TIDALLY LOCKED TERRESTRIAL PLANETS

Next-generation space telescopes will allow us to characterize terrestrial exoplanets. To do so effectively it will be crucial to make use of all available data. We investigate which atmospheric properties can, and cannot, be inferred from the broadband thermal phase curve of a dry and tidally locked terrestrial planet. First, we use dimensional analysis to show that phase curves are controlled by six nondimensional parameters. Second, we use an idealized general circulation model to explore the relative sensitivity of phase curves to these parameters. We find that the feature of phase curves most sensitive to atmospheric parameters is the peak-to-trough amplitude. Moreover, except for hot and rapidly rotating planets, the phase amplitude is primarily sensitive to only two nondimensional parameters: (1) the ratio of dynamical to radiative timescales and (2) the longwave optical depth at the surface. As an application of this technique, we show how phase curve measurements can be combined with transit or emission spectroscopy to yield a new constraint for the surface pressure and atmospheric mass of terrestrial planets. We estimate that a single broadband phase curve, measured over half an orbit with the James Webb Space Telescope, could meaningfully constrain the atmospheric mass of a nearby super-Earth. Such constraints will be important for studying the atmospheric evolution of terrestrial exoplanets as well as characterizing the surface conditions on potentially habitable planets.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Adaptive pruning-based optimization of parameterized quantum circuits

Abstract Variational hybrid quantum–classical algorithms are powerful tools to maximize the use of noisy intermediate-scale quantum devices. While past studies have developed powerful and expressive ansatze, their near-term applications have been limited by the difficulty of optimizing in the vast parameter space. In this work, we propose a heuristic optimization strategy for such ansatze used in variational quantum algorithms, which we call ‘parameter-efficient circuit training (PECT)’. Instead of optimizing all of the ansatz parameters at once, PECT launches a sequence of variational algorithms, in which each iteration of the algorithm activates and optimizes a subset of the total parameter set. To update the parameter subset between iterations, we adapt the Dynamic Sparse Reparameterization scheme which was originally proposed for training deep convolutional neural networks. We demonstrate PECT for the Variational Quantum Eigensolver, in which we benchmark unitary coupled-cluster ansatze including UCCSD and k -UpCCGSD, as well as the Low-Depth Circuit Ansatz (LDCA), to estimate ground state energies of molecular systems. We additionally use a layerwise variant of PECT to optimize a hardware-efficient circuit for the Sycamore processor to estimate the ground state energy densities of the one-dimensional Fermi-Hubbard model. From our numerical data, we find that PECT can enable optimizations of certain ansatze that were previously difficult to converge and more generally can improve the performance of variational algorithms by reducing the optimization runtime and/or the depth of circuits that encode the solution candidate(s).

Physics↗

DOE OSTI2016

SU-F-J-158: Respiratory Motion Resolved, Self-Gated 4D-MRI Using Rotating Cartesian K-Space Sampling

Purpose: Dynamic MRI has been used to quantify respiratory motion of abdominal organs in radiation treatment planning. Many existing 4D-MRI methods based on 2D acquisitions suffer from limited slice resolution and additional stitching artifacts when evaluated in 3D{sup 1}. To address these issues, we developed a 4D-MRI (3D dynamic) technique with true 3D k-space encoding and respiratory motion self-gating. Methods: The 3D k-space was acquired using a Rotating Cartesian K-space (ROCK) pattern, where the Cartesian grid was reordered in a quasi-spiral fashion with each spiral arm rotated using golden angle{sup 2}. Each quasi-spiral arm started with the k-space center-line, which were used as self-gating{sup 3} signal for respiratory motion estimation. The acquired k-space data was then binned into 8 respiratory phases and the golden angle ensures a near-uniform k-space sampling in each phase. Finally, dynamic 3D images were reconstructed using the ESPIRiT technique{sup 4}. 4D-MRI was performed on 6 healthy volunteers, using the following parameters (bSSFP, Fat-Sat, TE/TR=2ms/4ms, matrix size=500×350×120, resolution=1×1×1.2mm, TA=5min, 8 respiratory phases). Supplemental 2D real-time images were acquired in 9 different planes. Dynamic locations of the diaphragm dome and left kidney were measured from both 4D and 2D images. The same protocol was also performed on a MRI-compatible motion phantom where the motion was programmed with different amplitude (10–30mm) and frequency (3–10/min). Results: High resolution 4D-MRI were obtained successfully in 5 minutes. Quantitative motion measurements from 4D-MRI agree with the ones from 2D CINE (<5% error). The 4D images are free of the stitching artifacts and their near-isotropic resolution facilitates 3D visualization and segmentation of abdominal organs such as the liver, kidney and pancreas. Conclusion: Our preliminary studies demonstrated a novel ROCK 4D-MRI technique with true 3D k-space encoding and respiratory motion self-gating. The technique leads to high-resolution and artifacts-free 4D images for improved abdominal organ motion studies. K.S acknowledges funding support from NIH R01CA188300.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2017

Spectrum of spin waves in cold polarized gases

The spin dynamics of cold polarized gases are investigated using the Boltzmann equation. The dispersion relation for spin waves (transverse component of the magnetic moment) and the spin diffusion coefficient of the longitudinal component of the magnetic moment are calculated without using fitting parameters. The spin wave frequency and the diffusion coefficient for rubidium atoms are estimated numerically.

74 ATOMIC AND MOLECULAR PHYSICS↗