Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inverse Optimization”

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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Sparse-grid, reduced-basis Bayesian inversion: Nonaffine-parametric nonlinear equations

We extend the reduced basis (RB) accelerated Bayesian inversion methods for affine-parametric, linear operator equations which are considered in [16,17] to non-affine, nonlinear parametric operator equations. We generalize the analysis of sparsity of parametric forward solution maps in [20] and of Bayesian inversion in [48,49] to the fully discrete setting, including Petrov–Galerkin high-fidelity (“HiFi”) discretization of the forward maps. We develop adaptive, stochastic collocation based reduction methods for the efficient computation of reduced bases on the parametric solution manifold. The nonaffinity and nonlinearity with respect to (w.r.t.) the distributed, uncertain parameters and the unknown solution is collocated; specifically, by the so-called Empirical Interpolation Method (EIM). For the corresponding Bayesian inversion problems, computational efficiency is enhanced in two ways: first, expectations w.r.t. the posterior are computed by adaptive quadratures with dimension-independent convergence rates proposed in [49]; the present work generalizes [49] to account for the impact of the PG discretization in the forward maps on the convergence rates of the Quantities of Interest (QoI for short). Second, we propose to perform the Bayesian estimation only w.r.t. a parsimonious, RB approximation of the posterior density. Based on the approximation results in [49], the infinite-dimensional parametric, deterministic forward map and operator admit N-term RB and EIM approximations which converge at rates which depend only on the sparsity of the parametric forward map. In several numerical experiments, the proposed algorithms exhibit dimension-independent convergence rates which equal, at least, the currently known rate estimates for N-term approximation. We propose to accelerate Bayesian estimation by first offline construction of reduced basis surrogates of the Bayesian posterior density. The parsimonious surrogates can then be employed for online data assimilation and for Bayesian estimation. They also open a perspective for optimal experimental design.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Development and Implementation of an Automatic Continuous Online Monitoring and Control Platform for Polymerization Reactions to Sharply Boost Energy and Resource Efficiency in Polymer Manufacturing

The project goal was to create an energy saving paradigm shift in how polymers are manufactured in the 21st century. It used Automatic Continuous Online Monitoring of Polymerization reactions (ACOMP) integrated for the first time with automatic active control to create the innovative ‘ACOMP/Control Interface’, or ‘ACOMP/CI’. ACOMP/CI will begin the transformation from old, inefficient processes into highly evolved, energy and resource efficient ones. The ACOMP platform is broadly applicable to many types of reactions and processes throughout the vast polymer industry. The industry provides materials for sectors such as automotive, aerospace, oil recovery, agriculture, paints, resins, adhesives, pharmaceuticals and therapeutic proteins, optics, electronics, lightweight building materials, and many more. The U.S. chemical industry is one of the last major sectors in which the U.S. has top global stature. It consumes 24% of all U.S. manufacturing energy, produces over $800B of product annually, supports 25% of the U.S. GDP and employs over 6 million people. It is also a major source of GHG emissions. Polymers make up approximately 30% of this sector. It is estimated that annually 60 TBtu of energy could be saved and 3 million tons less of GHG emissions produced by optimizing production in the polyolefin manufacturing sector alone. The project scope included first time design and prototyping of an ACOMP/CI, creation of active reaction controllers, and demonstration of control capabilities on ideal, low concentration polymerization reactions. All these elements of the scope were met, including advances and findings not originally anticipated. Extensions to more complex reactions, beyond the reactor capabilities of the current project ACOMP/CI, such as polyolefins and other high pressure/high temperature reactions, are being proposed in Fall 2017 to CESMII, a DoE based NNMI. The initial proposal was for a three year funded project, but this was reduced to a two year project and budget due to funding constraints. Hence, some of the original plans, such as adaptation of the ACOMP/CI to more relevant industrial processes, such as emulsion and dispersion technologies, could not be carried out. A third year of funding was requested at the end of the project, but DoE did not have resources to grant this. The sub-contractor Fluence Analytics (previously Advanced Polymer Monitoring Technologies, Inc) designed, prototyped, and commissioned a working ACOMP/CI by June 2015. The reaction characteristics to be automatically controlled were i) conversion kinetics, ii) molecular weight, iii) copolymer composition, and iv) simultaneous molecular weight and composition. A two pronged control strategy was used. The Tulane/Fluence group took a basic principles approach that did not rely on kinetic models. The LSU group took a more complex, non-linear model-oriented approach involving complete kinetic descriptions of the reaction system Each of these approaches proved successful in their own way. By April 2016 fully automatic control of conversion and weight average molecular weight, Mw, trajectories was achieved using the Tulane/Fluence (TF) basic principles controller. Similar results were obtained by the LSU non-linear model controller by August 2016. The demonstration system was aqueous free radical polymerization of acrylamide, Am. The control variables were temperature and semi-batch feed to the reactor of Am monomer and initiator. A demonstration of active manual conversion control in an industrial process using high solids in inverse emulsion polymerization of Am was achieved. During Summer 2016 the TF controller was used in conjunction with a chain transfer agent, another control variable, to automatically produce multi-modal molecular weight distributions, MWD, in a single reactor. Industrially, multi-modal MWD are produced by mixing products made in separate reactors, requiring significant extra time, energy, and reactor resources. Recognizing the industrial potential a patent on automatic production of multi-modal polymers was filed, and DoE acknowledged. In Fall 2016 the TF team developed a basic principles controller for copolymer composition and demonstrated it on aqueous free radical copolymerization of comonomers Am and styrene sulfonate, SS. TF then fused the Mw and composition controllers to achieve simultaneous control of both Mw and copolymer composition trajectories. Numerous simultaneous trajectories were demonstrated, including a trimodal composition distribution with constant Mw. Meanwhile, the LSU group developed a Kalman filter to improve the results of their automatic Mw and conversion controller and successful tests were made. During the project the TF team developed a means of computing full MWD during polymer synthesis without need for any chromatographic separation, based on model distributions. This means the polymer product is ‘born characterized’ and this can eliminate post-manufacture analytical laboratory quality control. TF filed a joint patent application on this new approach to chromatography-free determination of MWD with acknowledgment to DoE. The Tulane group obtained a 60MHz NMR during the project and recently completed the first work on separating three comonomers, Am, SS, and Na-acrylate, with a first demonstration of terpolymer composition control with the TF basic principles controller. Widespread dissemination of ACOMP/CI in the polymer manufacturing sector will bolster DoE goals of energy efficiency and reduced GHG emissions: The ability to monitor and actively control polymerization reactions will lead to more efficient use of energy and non-renewable resources, plant and labor time, increase the safety of manufacturing personnel, and will enhance product quality and lead to feasibility of manufacturing of polymers currently too complex for industrial scale production, while leading to less GHG emissions per kilo of product, and allowing for increased U.S. competitiveness in this enormous manufacturing sector. When ACOMP/CI is expanded to the polyolefin industry it is estimated that 60 TeraBTU/year of energy can be saved. Much of this saving is anticipated to come from optimized control of grade changeovers in steady state reactors and maintenance of steady states. Conclusions: ACOMP’s ability to provide continuous realtime data streams of measured polymer and reaction characteristics made it possible, for the first time, to directly and automatically control free radical polymerization reactions. An industrial client of Fluence Analytics has requisitioned the first ACOMP/CI which uses the TF basic principles controller. This sets the stage for FA to add control features to the ACOMP systems it has begun to install on the industrial scale beginning in 2014. Recommendations: This successful project was mainly limited to ideal polymerizations not of an industrial sort. The most energy intensive portion of polymer manufacturing is polyolefins. Adoption of ACOMP/CI to this enormous industrial sector faces the enormous challenges of high temperature, high pressure continuous sampling and high temperature sensor operation to obtain the continuous data needed for direct reaction control. The project team has a strategy for achieving this ambitious goal and will present it in Fall 2017 as a proposal to CESMII/DoE. It is recommended that this upcoming proposal be funded in order to make full use of the achievements of this just ended DoE project as the next step towards making polyolefin ACOMP/CI an energy saving reality. It is projected that ACOMP/CI can have its first polyolefin testbed demonstrations within two years of beginning the proposed project.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Investigation of the Dzyaloshinskii-Moriya interaction and room temperature skyrmions in W/CoFeB/MgO thin films and microwires

Recent studies have shown that material structures, which lack structural inversion symmetry and have high spin-orbit coupling can exhibit chiral magnetic textures and skyrmions which could be a key component for next generation storage devices. The Dzyaloshinskii-Moriya Interaction (DMI) that stabilizes skyrmions is an anti-symmetric exchange interaction favoring non-collinear orientation of neighboring spins. It has been shown that materials systems with high DMI can lead to very efficient domain wall and skyrmion motion by spin-orbit torques. To engineer such devices, it is important to quantify the DMI for a given material system. Here, we extract the DMI at the Heavy Metal/Ferromagnet interface using two complementary measurement schemes, namely, asymmetric domain wall motion and the magnetic stripe annihilation. By using the two different measurement schemes, we find for W(5 nm)/Co 20 Fe 60 B 20 (0.6 nm)/MgO(2 nm) the DMI to be 0.68 ± 0.05 mJ/m 2 and 0.73 ± 0.5 mJ/m 2 , respectively. Furthermore, we show that this DMI stabilizes skyrmions at room temperature and that there is a strong dependence of the DMI on the relative composition of the CoFeB alloy. Finally, we optimize the layers and the interfaces using different growth conditions and demonstrate that a higher deposition rate leads to a more uniform film with reduced pinning and skyrmions that can be manipulated by spin orbit torques.

36 MATERIALS SCIENCE↗
NASA NTRS2012

NASA's Carbon Monitoring System Flux-Pilot Project: A Multi-Component Analysis System for Carbon-Cycle Research and Monitoring

The importance of greenhouse gas increases for climate motivates NASA s observing strategy for CO2 from space, including the forthcoming Orbiting Carbon Observatory (OCO-2) mission. Carbon cycle monitoring, including attribution of atmospheric concentrations to regional emissions and uptake, requires a robust modeling and analysis infrastructure to optimally extract information from the observations. NASA's Carbon-Monitoring System Flux-Pilot Project (FPP) is a prototype for such analysis, combining a set of unique tools to facilitate analysis of atmospheric CO2 along with fluxes between the atmosphere and the terrestrial biosphere or ocean. NASA's analysis system is unique, in that it combines information and expertise from the land, oceanic, and atmospheric branches of the carbon cycle and includes some estimates of uncertainty. Numerous existing space-based missions provide information of relevance to the carbon cycle. This study describes the components of the FPP framework, assessing the realism of computed fluxes, thus providing the basis for research and monitoring applications. Fluxes are computed using data-constrained terrestrial biosphere models and physical ocean models, driven by atmospheric observations and assimilating ocean-color information. Use of two estimates provides a measure of uncertainty in the fluxes. Along with inventories of other emissions, these data-derived fluxes are used in transport models to assess their consistency with atmospheric CO2 observations. Closure is achieved by using a four-dimensional data assimilation (inverse) approach that adjusts the terrestrial biosphere fluxes to make them consistent with the atmospheric CO2 observations. Results will be shown, illustrating the year-to-year variations in land biospheric and oceanic fluxes computed in the FPP. The signals of these surface-flux variations on atmospheric CO2 will be isolated using forward modeling tools, which also incorporate estimates of transport error. The results will be discussed in the context of interannual variability of observed atmospheric CO2 distributions.

Pawson, S.↗
NASA NTRSDate not supplied

Toward an IMU-based Space Suit Motion Capture System

Spacesuits are complex engineering systems that sustain human health and enable performance outside of Earth-like environments. These systems must support human mobility and physical workload demands while minimizing injury risk during extravehicular activity (EVA). Future EVA on the lunar surface during the Artemis program is expected to be more frequent and require higher physical workloads than previous EVAs during the ISS, Shuttle, or Apollo programs. Hence it is important to optimize future as well as current spacesuits to be efficient and comfortable for the success of space and planetary missions. To enable this, an efficient method is needed to test these spacesuits on the ground.When testing spacesuits in ground environments, it is often necessary to understand the kinematics of the suit to validate the design against relevant requirements or characterize the physical workload necessary to operate the suit. This is a challenging task for traditional optical motion capture (OMC) approaches: suit-mounted OMC markers are easily occluded by the subject or environment and may become detached during testing. Controlling lighting and reflectivity of objects in the motion capture volume is also difficult. Fixed-position OMC cameras also constrain testing to a small and contrived laboratory environment, disallowing kinematics capture in field environments.One promising alternative is the use of suit-mounted inertial measurement units (IMUs). These sensors are small, unobtrusive, and portable, but come at the cost of increased sensor noise and complexity of the software and mathematics to analyze the collected data. To this end, engineers at NASA are developing the Augmented Suit Inverse Kinematics (ASIK) system, a complete motion capture methodand inverse kinematics solver which relies solely on a network of wireless IMUs attached to the major kinematic segments of the spacesuit. The ASIK modeling language allows for the simple inclusion of probabilistic priors such as suit size and shape or IMU positions and rotations. Furthermore, to increase accuracy and reduce operational overhead to use this motion capture approach, the developed inverse kinematics solver exploits so-called self-calibratingalgorithmic techniques, which reduce the need for precise alignment of the sensors on the segments or scripted functional calibration procedures. The ASIK system was tested in a 7-subject pilot study. Each subject donned NASA’s new prototype exploration spacesuit in the Active Response Gravity Offload System (ARGOS) facility at the NASA Johnson Space Center. The subjects were outfitted with a set of 14 APDM (Portland, OR, USA) Opal IMUs, 12 of which were used in the ASIK model to estimate lower body and trunk kinematics. The subjects were also outfitted with a set of reflective OMC markers and traditional OMC data was collected and processed. Presented results will include characterization of ASIK-derived suit joint angles accuracy against an optical motion capture datum. Discussion of these results, as well as discussion of system calibration and nuances of mathematical observability, will be included.If successful, IMU-based motion capture will enable testing and validation of spacesuits more frequently, with less overhead, in more extreme environments. Future work will apply these techniques to common spacesuit testing tasks, such as gait, mobility, and balance assessment, physical workload characterization, and ergonomics evaluations.

Timothy Mcgrath↗
DOE OSTI2016

Emitter/absorber interface of CdTe solar cells

The performance of CdTe solar cells can be very sensitive to the emitter/absorber interface, especially for high-efficiency cells with high bulk lifetime. Performance losses from acceptor-type interface defects can be significant when interface defect states are located near mid-gap energies. Numerical simulations show that the emitter/absorber band alignment, the emitter doping and thickness, and the defect properties of the interface (i.e., defect density, defect type, and defect energy) can all play significant roles in the interface recombination. In particular, a type I heterojunction with small conduction-band offset (0.1 eV ≤ ΔE{sub C} ≤ 0.3 eV) can help maintain good cell efficiency in spite of high interface defect density, much like with Cu(In,Ga)Se{sub 2} (CIGS) cells. The basic principle is that positive ΔE{sub C}, often referred to as a “spike,” creates an absorber inversion and hence a large hole barrier adjacent to the interface. As a result, the electron-hole recombination is suppressed due to an insufficient hole supply at the interface. A large spike (ΔE{sub C} ≥ 0.4 eV), however, can impede electron transport and lead to a reduction of photocurrent and fill-factor. In contrast to the spike, a “cliff” (ΔE{sub C} < 0 eV) allows high hole concentration in the vicinity of the interface, which will assist interface recombination and result in a reduced open-circuit voltage. Another way to mitigate performance losses due to interface defects is to use a thin and highly doped emitter, which can invert the absorber and form a large hole barrier at the interface. CdS is the most common emitter material used in CdTe solar cells, but the CdS/CdTe interface is in the cliff category and is not favorable from the band-offset perspective. The ΔE{sub C} of other n-type emitter choices, such as (Mg,Zn)O, Cd(S,O), or (Cd,Mg)Te, can be tuned by varying the elemental ratio for an optimal positive value of ΔE{sub C}. These materials are predicted to yield higher voltages and would therefore be better candidates for the CdTe-cell emitter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

BLOC Site - ASSIST Thermodynamic Retrievals TROPoe v0.18 / Derived Data

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). This is a post-processed dataset and recommended for use. The profiles are retrieved every 10 minutes from instantaneous radiances observed with an Atmospheric Sounder Spectrometer by Infrared Spectral Technology (ASSIST, Michaud-Belleau et al. 2025) operated by NOAA Physical Sciences Laboratory (PSL) on Block Island for WFIP3. The spectral bands used in the retrieval are in the wavenumber range from 612 - 905.4 cm-1 and are specified in Turner and Löhnert (2021). Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a collocated surface tower operated by NOAA PSL. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY. The TROPoe docker container (version 0.18) is available from Docker Hub at https://hub.docker.com/r/davidturner53/tropoe/tags, and the source code code is available in the GitHub repository https://github.com/OAR-atmospheric-observations/TROPoe.

17 WIND ENERGY↗
DOE OSTI2026

NANT Site - ASSIST Thermodynamic Retrievals TROPoe v0.18 / Derived Data

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). This is a post-processed dataset and recommended for use. The profiles are retrieved every 10 minutes from instantaneous radiances observed with an Atmospheric Sounder Spectrometer by Infrared Spectral Technology (ASSIST, Michaud-Belleau et al. 2025) operated by NOAA Physical Sciences Laboratory (PSL) on Nantucket Island for WFIP3. The spectral bands used in the retrieval are in the wavenumber range from 612 - 905.4 cm-1 and are specified in Turner and Löhnert (2021). Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a collocated surface tower operated by NOAA PSL. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY. The TROPoe docker container (version 0.18) is available from Docker Hub at https://hub.docker.com/r/davidturner53/tropoe/tags, and the source code code is available in the GitHub repository https://github.com/OAR-atmospheric-observations/TROPoe.

17 WIND ENERGY↗
DOE OSTI2017

Identifying attributes of CO 2 leakage zones in shallow aquifers using a parametric level set method: Original Research Article: Identifying attributes of CO 2 leakage zones in shallow aquifers using a parametric level set method

Leakage through abandoned wells and geologic faults poses the greatest risk to CO 2 storage permanence. For shallow aquifers, secondary CO 2 plumes emanating from the leak zones may go undetected for a sustained period of time and has the greatest potential to cause large-scale and long-term environmental impacts. Identification of the attributes of leak zones, including their shape, location, and strength, is required for environmental risk assessment. This study applies a parametric level set (PaLS) method to characterize leak zones. Level set methods are appealing for tracking topological changes and recovering unknown shapes of objects. However, level set evolution using the conventional level set methods is computationally challenging. In PaLS, the level set function is approximated by using a weighted sum of basis functions and the level set evolution problem is replaced by an optimization problem. The efficacy of PaLS is demonstrated by reconstructing attributes of a CO 2 leak zone in a carbonate aquifer. Our results show that PaLS is a robust source identification method that can recover the approximate leak zone locations in the presence of measurement errors, model parameter uncertainty, and inaccurate guesses of source flux strengths. Accurate delineation of source geometry is achievable when a relatively dense observation network is used. The PaLS-based, source recovery framework introduced in this work is generic and can be adapted for any reactive transport model by simply switching the pre- and post-processing routines.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Analysis of Alfven eigenmodes destabilization by energetic particles in TJ-II using a Landau-closure model

Alfvén Eigenmodes (AE) can be destabilized by energetic particles in neutral beam injection (NBI) heated plasmas through inverse Landau damping and couplings with gap modes in the shear Alfvén continua. In this paper, we describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system using the reduced MHD equations, density and parallel velocity moments for the energetic particles as well as the geodesic acoustic wave dynamics. A closure relation adds the Landau damping and resonant destabilization effects in the model. We apply the model to study the Alfvén modes stability in TJ-II, performing a parametric analysis in a range of realistic values of energetic particle β ($\beta_{f}$ ), ratios of thermal/Alfvén velocities ($V_{\rm th}/V_{\rm A0}$ ), energetic particle density profiles and toroidal modes (n) including toroidal and helical couplings. The study predicts a large helical coupling between different toroidal modes and the destabilization of helical Alfvén eigenmodes (HAE) with frequencies similar to the AE activity measured in TJ-II, between 50–400 kHz. The analysis has also revealed the destabilization of GAE (global Alfvén eigenmodes), TAE (toroidal Alfvén eigenmodes) and EPM (energetic particle modes). For the modes considered here, optimized TJ-II operations require a $\rlap{-} \iota$ profile in the range of $[0.845, 0.979]$ to stabilize AEs in the inner and middle plasma. AEs in the plasma periphery cannot be fully stabilized, although for a configuration with $\rlap{-} \iota = [0.945, 1.079]$ , only $n=7, 11, 15$ AE are unstable with a growth rate 4 times smaller compared to the standard $\rlap{-} \iota = [1.54, 1.68]$ case and a frequency of 100 kHz. Finally, we reproduce the frequency sweeping evolution of the AE frequency observed in TJ-II as the $\rlap{-} \iota$ profile is varied. The AE frequency sweeping is caused by consecutive changes of the instability dominant modes between different helical families.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2003

Mathematical Inversion of Lightning Data: Techniques and Applications

A survey of some interesting mathematical inversion studies dealing with radio, optical, and electrostatic measurements of lightning are presented. A discussion of why NASA is interested in lightning, what specific physical properties of lightning are retrieved, and what mathematical techniques are used to perform the retrievals are discussed. In particular, a relatively new multi-station VHF time-of-arrival (TOA) antenna network is now on-line in Northern Alabama and will be discussed. The network, called the Lightning Mapping Array (LMA), employs GPS timing and detects VHF radiation from discrete segments (effectively point emitters) that comprise the channel of lightning strokes within cloud and ground flashes. The LMA supports on-going ground-validation activities of the low Earth orbiting Lightning Imaging Sensor (LIS) satellite developed at NASA Marshall Space Flight Center (MSFC) in Huntsville, Alabama. The LMA also provides detailed studies of the distribution and evolution of thunderstorms and lightning in the Tennessee Valley, and offers interesting comparisons with other meteorological/geophysical datasets. In order to take full advantage of these benefits, it is essential that the LMA channel mapping accuracy (in both space and time) be fully characterized and optimized. A new channel mapping retrieval algorithm is introduced for this purpose. To characterize the spatial distribution of retrieval errors, the algorithm has been applied to analyze literally tens of millions of computer-simulated lightning VHF point sources that have been placed at various ranges, azimuths, and altitudes relative to the LMA network. Statistical results are conveniently summarized in high-resolution, color-coded, error maps.

Koshak, William↗
DOE OSTI2015

SU-C-BRD-05: Non-Invasive in Vivo Biodosimetry in Radiotherapy Patients Using Electron Paramagnetic Resonance (EPR) Spectroscopy

Purpose: Medical intervention following a major, unplanned radiation event can elevate the human whole body exposure LD50 from 3 to 7 Gy. On a large scale, intervention cannot be achieved effectively without accurate and efficient triage. Current methods of retrospective biodosimetry are restricted in capability and applicability; published human data is limited. We aim to further develop, validate, and optimize an automated field-deployable in vivo electron paramagnetic resonance (EPR) instrument that can fill this need. Methods: Ionizing radiation creates highly-stable, carbonate-based free radicals within tooth enamel. Using a process similar to nuclear magnetic resonance, EPR directly measures the presence of radiation-induced free radicals. We performed baseline EPR measurements on one of the upper central incisors of total body irradiation (TBI) and head and neck (H&N) radiotherapy patients before their first treatment. Additional measurements were performed between subsequent fractions to examine the EPR response with increasing radiation dose. Independent dosimetry measurements were performed with optically-stimulated luminescent dosimeters (OSLDs) and diodes to more accurately establish the relationship between EPR signal and delivered radiation dose. Results: 36 EPR measurements were performed over the course of four months on two TBI and four H & N radiotherapy patients. We observe a linear increase in EPR signal with increasing dose across the entirety of the tested range. A linear least squares-weighted fit of delivered dose versus measured signal amplitude yields an adjusted R-square of 0.966. The standard error of inverse prediction (SEIP) is 1.77 Gy. For doses up to 7 Gy, the range most relevant to triage, we calculate an SEIP of 1.29 Gy. Conclusion: EPR spectroscopy provides a promising method of retrospective, non-invasive, in vivo biodosimetry. Our preliminary data show an excellent correlation between predicted signal amplitude and delivered dose. With further development, a robust means of predicting delivered radiation dose from EPR measurements is expected. This project was funded by the Biomedical Advanced Research and Development Authority (BARDA) within the U.S. Department of Health and Human Services subcontracted through the Geisel School of Medicine at Dartmouth and by the Dartmouth Physically-Based Biodosimetry Center for Medical Countermeasures Against Radiation (Dart-Dose CMCR) Pilot Program.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1991

A unified viscous theory of lift and drag of 2-D thin airfoils and 3-D thin wings

A unified viscous theory of 2-D thin airfoils and 3-D thin wings is developed with numerical examples. The viscous theory of the load distribution is unique and tends to the classical inviscid result with Kutta condition in the high Reynolds number limit. A new theory of 2-D section induced drag is introduced with specific applications to three cases of interest: (1) constant angle of attack; (2) parabolic camber; and (3) a flapped airfoil. The first case is also extended to a profiled leading edge foil. The well-known drag due to absence of leading edge suction is derived from the viscous theory. It is independent of Reynolds number for zero thickness and varies inversely with the square root of the Reynolds number based on the leading edge radius for profiled sections. The role of turbulence in the section induced drag problem is discussed. A theory of minimum section induced drag is derived and applied. For low Reynolds number the minimum drag load tends to the constant angle of attack solution and for high Reynolds number to an approximation of the parabolic camber solution. The parabolic camber section induced drag is about 4 percent greater than the ideal minimum at high Reynolds number. Two new concepts, the viscous induced drag angle and the viscous induced separation potential are introduced. The separation potential is calculated for three 2-D cases and for a 3-D rectangular wing. The potential is calculated with input from a standard doublet lattice wing code without recourse to any boundary layer calculations. Separation is indicated in regions where it is observed experimentally. The classical induced drag is recovered in the 3-D high Reynolds number limit with an additional contribution that is Reynold number dependent. The 3-D viscous theory of minimum induced drag yields an equation for the optimal spanwise and chordwise load distribution. The design of optimal wing tip planforms and camber distributions is possible with the viscous 3-D wing theory.

Yates, John E.↗
NASA NTRS2013

Support Routines for In Situ Image Processing

This software consists of a set of application programs that support ground-based image processing for in situ missions. These programs represent a collection of utility routines that perform miscellaneous functions in the context of the ground data system. Each one fulfills some specific need as determined via operational experience. The most unique aspect to these programs is that they are integrated into the large, in situ image processing system via the PIG (Planetary Image Geometry) library. They work directly with space in situ data, understanding the appropriate image meta-data fields and updating them properly. The programs themselves are completely multimission; all mission dependencies are handled by PIG. This suite of programs consists of: (1)marscahv: Generates a linearized, epi-polar aligned image given a stereo pair of images. These images are optimized for 1-D stereo correlations, (2) marscheckcm: Compares the camera model in an image label with one derived via kinematics modeling on the ground, (3) marschkovl: Checks the overlaps between a list of images in order to determine which might be stereo pairs. This is useful for non-traditional stereo images like long-baseline or those from an articulating arm camera, (4) marscoordtrans: Translates mosaic coordinates from one form into another, (5) marsdispcompare: Checks a Left Right stereo disparity image against a Right Left disparity image to ensure they are consistent with each other, (6) marsdispwarp: Takes one image of a stereo pair and warps it through a disparity map to create a synthetic opposite- eye image. For example, a right eye image could be transformed to look like it was taken from the left eye via this program, (7) marsfidfinder: Finds fiducial markers in an image by projecting their approximate location and then using correlation to locate the markers to subpixel accuracy. These fiducial markets are small targets attached to the spacecraft surface. This helps verify, or improve, the pointing of in situ cameras, (8) marsinvrange: Inverse of marsrange . given a range file, re-computes an XYZ file that closely matches the original. . marsproj: Projects an XYZ coordinate through the camera model, and reports the line/sample coordinates of the point in the image, (9) marsprojfid: Given the output of marsfidfinder, projects the XYZ locations and compares them to the found locations, creating a report showing the fiducial errors in each image. marsrad: Radiometrically corrects an image, (10) marsrelabel: Updates coordinate system or camera model labels in an image, (11) marstiexyz: Given a stereo pair, allows the user to interactively pick a point in each image and reports the XYZ value corresponding to that pair of locations. marsunmosaic: Extracts a single frame from a mosaic, which will be created such that it could have been an input to the original mosaic. Useful for creating simulated input frames using different camera models than the original mosaic used, and (12) merinverter: Uses an inverse lookup table to convert 8-bit telemetered data to its 12-bit original form. Can be used in other missions despite the name.

Deen, Robert G.↗
DOE OSTI2025

BLOC Site - ASSIST Thermodynamic retrievals TROPoe / Derived Data

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous radiances observed with an Atmospheric Sounder Spectrometer by Infrared Spectral Technology (ASSIST, Rochette et al. 2009) operated by NOAA Physical Sciences Laboratory on Block Island for WFIP3. The spectral bands used in the retrieval are in the wavenumber range from 612 - 905.4 cm-1 and are specified in Turner and Löhnert (2021). Radiances are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a sensor attached to the ASSIST housing. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗
NASA NTRS1999

Cloud Thickness from Diffusion of Lidar Pulses in Clouds

Measurements of the distribution of reflected light from a laser beam incident on an aqueous suspension of particles or "cloud" with known thickness and particle size distribution are reported. The distribution is referred to as the "cloud radiative Green's function", G. In the diffusion domain, G is sensitive to cloud thickness, allowing that important quantity to be retrieved. The goal of the laboratory simulation is to provide preliminary estimates of sensitivity of G to cloud thickness,for use in the optimal design of an offbeam Lidar instrument for remote sensing of cloud thickness (THOR, Thickness from Offbeam Returns). These clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The microsphere size distribution is roughly lognormal, from 0.5 microns to 25 microns, similar to real clouds. Density of suspended spheres is adjusted so mean-free-path of visible photons is about 10 cm, approximately 1000 times smaller than in real clouds. The light source is a ND:YAG laser at 530 nm. Detectors are flux and photon-counting Photomultiplier Tube (PMTS), with a glass probe for precise positioning. A Labview 5 VI controls positioning, and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider, and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns, and step size is selectable from the VI software. Far from the incident beam, the rate of exponential increase as the direction of the incident beam is approached scales as expected from diffusion theory, linearly with the cloud thickness, and inversely as the square root of the reduced optical thickness, and is independent of particle size. Near the beam the signal begins to increase faster than exponential, due to single and low-order scattering near the backward direction, and here the distribution depends on particle size. Results are being used to verify 3D Monte Carlo radiative transfer simulations, used to estimate signal-to-noise ratios for remotely sensed off beam returns, for both homogeneous and inhomogeneous clouds. Signal-to-noise estimates show that unfiltered observations are straight forward at night, while narrow band pass filters are being studied for day.

Cahalan, Robert F.↗
NASA NTRS1998

Combined Numerical/Analytical Perturbation Solutions of the Navier-Stokes Equations for Aerodynamic Ejector/Mixer Nozzle Flows

In spite of rapid advances in both scalar and parallel computational tools, the large number of variables involved in both design and inverse problems make the use of sophisticated fluid flow models impractical, With this restriction, it is concluded that an important family of methods for mathematical/computational development are reduced or approximate fluid flow models. In this study a combined perturbation/numerical modeling methodology is developed which provides a rigorously derived family of solutions. The mathematical model is computationally more efficient than classical boundary layer but provides important two-dimensional information not available using quasi-1-d approaches. An additional strength of the current methodology is its ability to locally predict static pressure fields in a manner analogous to more sophisticated parabolized Navier Stokes (PNS) formulations. To resolve singular behavior, the model utilizes classical analytical solution techniques. Hence, analytical methods have been combined with efficient numerical methods to yield an efficient hybrid fluid flow model. In particular, the main objective of this research has been to develop a system of analytical and numerical ejector/mixer nozzle models, which require minimal empirical input. A computer code, DREA Differential Reduced Ejector/mixer Analysis has been developed with the ability to run sufficiently fast so that it may be used either as a subroutine or called by an design optimization routine. Models are of direct use to the High Speed Civil Transport Program (a joint government/industry project seeking to develop an economically.viable U.S. commercial supersonic transport vehicle) and are currently being adopted by both NASA and industry. Experimental validation of these models is provided by comparison to results obtained from open literature and Limited Exclusive Right Distribution (LERD) sources, as well as dedicated experiments performed at Texas A&M. These experiments have been performed using a hydraulic/gas flow analog. Results of comparisons of DREA computations with experimental data, which include entrainment, thrust, and local profile information, are overall good. Computational time studies indicate that DREA provides considerably more information at a lower computational cost than contemporary ejector nozzle design models. Finally. physical limitations of the method, deviations from experimental data, potential improvements and alternative formulations are described. This report represents closure to the NASA Graduate Researchers Program. Versions of the DREA code and a user's guide may be obtained from the NASA Lewis Research Center.

DeChant, Lawrence Justin↗
DOE OSTI2025

Chapter 4 - Recent Advances in Identification of Differential Equations from Noisy Data: IDENT Review

Differential equations and numerical methods are extensively used to model various real-world phenomena in science and engineering. With modern developments, we aim to find the underlying differential equation from a single observation of time-dependent data. If we assume that the differential equation is a linear combination of various linear and nonlinear differential terms, then the identification problem can be formulated as solving a linear system. The goal then reduces to finding the optimal coefficient vector that best represents the time derivative of the given data. We review some recent works on the identification of differential equations. We find some common themes for the improved accuracy: (i) The formulation of linear system with proper denoising is important, (ii) how to utilize sparsity and model selection to find the correct coefficient support needs careful attention, and (iii) there are ways to improve the coefficient recovery. We present an overview and analysis of recent developments on the topic.

97 MATHEMATICS AND COMPUTING↗