SEARCH · Search NASA
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.
Far-field patterns reconstructed. Nonuniform-sampling technique uses interpolation algorithm to obtain far-field radiation pattern of antenna at any point u, v based on measurements at few sampling points. Applicable to any components of electric field "E" at measurement locations in spherical coordinate system centered at antenna.
An algorithm for recovering a collection of linear cracks in a homogeneous electrical conductor from boundary measurements of voltages induced by specified current fluxes is described. The technique is a variation of Newton's method and is based on taking weighted averages of the boundary data. An apparatus that was constructed specifically for generating laboratory data on which to test the algorithm is also described. The algorithm is applied to a number of different test cases and the results are discussed.
We present results of simulations of the distribution of 1809 keV radiation from the decay of Al-26 in the Galaxy. Recent observations of this emission line using the Gamma Ray Imaging Spectrometer (GRIS) have indicated that the bulk of the AL-26 must have a velocity of approx. 500 km/ s. We have previously shown that a velocity this large could be maintained over the 10(exp 6) year lifetime of the Al-26 if it is trapped in dust grains that are reaccelerated periodically in the ISM. Here we investigate whether a dust grain velocity of approx. 500 km/ s will produce a distribution of 1809 keV emission in latitude that is consistent with the narrow distribution seen by COMPTEL. We find that dust grain velocities in the range 275 - 1000 km/ s are able to reproduce the COMPTEL 1809 keV emission maps reconstructed using the Richardson-Lucy and Maximum Entropy image reconstruction methods while the emission map reconstructed using the Multiresolution Regularized Expectation Maximization algorithm is not well fit by any of our models. The Al-26 production rate that is needed to reproduce the observed 1809 keV intensity yields in a Galactic mass of Al-26 of approx. 1.5 - 2 solar mass which is in good agreement with both other observations and theoretical production rates.
In addition to the awe it inspires, lightning can illuminate the microphysical processes hidden away within deep convection. The current generation of space-based lightning mapping uses mostly 2D optical imaging to connect overall flash characteristics to their parent storm dynamics, but are missing a dimension’s worth of information. With lightning now classified as an essential climate variable, future spaceborne mappers will need improved capabilities to take advantage of the 3D structure of lightning flashes to support meteorological and climate modeling. We report here on a study of the feasibility of implementing a radio frequency (RF)–based network of satellites for 3D lightning mapping in high resolution from low Earth orbit (LEO). Lightning sources are simulated using existing tools of the Lightning Mapping Array (LMA), modified for orbital detection, and spatially reconstructed using a Levenberg-Marquardt geolocation algorithm to assess sources of uncertainty in these solutions. We analyze the benefits and limitations of this approach compared to existing orbital and ground-based methods. Results of this study show that lightning can be mapped in 3D with a vertical location accuracy better than 2 km using as few as five satellites in LEO capable of measuring the time-of-arrival of impulsive RF signals in the very high frequency (VHF) band. The consequence of this study is that high-resolution, spaceborne 3D mapping of lightning is achievable across most of the globe, having crucial implications for our understanding of not only lightning, but also severe weather development, climate science, and more.
The multigrid properties of two data reconstruction methods used for achieving second-order spatial accuracy when solving the two-dimensional Euler equations are examined. The data reconstruction methods are used with an implicit upwind algorithm which uses linearized backward-Euler time-differencing. The solution of the resulting linear system is performed by an iterative procedure. In the present study only regular quadrilateral grids are considered, so a red-black Gauss-Seidel iteration is used. Although the Jacobian is approximated by first-order upwind extrapolation, two alternative data reconstruction techniques for the flux integral that yield higher-order spatial accuracy at steady state are examined. The first method, probably most popular for structured quadrilateral grids, is based on estimating the cell gradients using one-dimensional reconstruction along curvilinear coordinates. The second method is based on Green's theorem. Analysis and numerical results for the two dimensional Euler equations show that data reconstruction based on Green's theorem has superior multigrid properties as compared to the one-dimensional data reconstruction method.
We have developed a new Bayesian image reconstruction method that has been shown to be superior to the best implementations of other competing methods, including Goodness-of-Fit methods such as Least-Squares fitting and Lucy-Richardson reconstruction, as well as Maximum Entropy (ME) methods such as those embodied in the MEMSYS algorithms. Our new method is based on the concept of the pixon, the fundamental, indivisible unit of picture information. Use of the pixon concept provides an improved image model, resulting in an image prior which is superior to that of standard ME. Our past work has shown how uniform information content pixons can be used to develop a 'Super-ME' method in which entropy is maximized exactly. Recently, however, we have developed a superior pixon basis for the image, the Fractal Pixon Basis (FPB). Unlike the Uniform Pixon Basis (UPB) of our 'Super-ME' method, the FPB basis is selected by employing fractal dimensional concepts to assess the inherent structure in the image. The Fractal Pixon Basis results in the best image reconstructions to date, superior to both UPB and the best ME reconstructions. In this paper, we review the theory of the UPB and FPB pixon and apply our methodology to the reconstruction of far-infrared imaging of the galaxy M51. The results of our reconstruction are compared to published reconstructions of the same data using the Lucy-Richardson algorithm, the Maximum Correlation Method developed at IPAC, and the MEMSYS ME algorithms. The results show that our reconstructed image has a spatial resolution a factor of two better than best previous methods (and a factor of 20 finer than the width of the point response function), and detects sources two orders of magnitude fainter than other methods.
A novel method to reconstruct the energy of hadronic showersin the CMS High Granularity Calorimeter (HGCAL) is presented. TheHGCAL is a sampling calorimeter with very fine transverse andlongitudinal granularity. The active media are silicon sensors andscintillator tiles readout by SiPMs and the absorbers are acombination of lead and Cu/CuW in the electromagnetic section, andsteel in the hadronic section. The shower reconstruction method isbased on graph neural networks and it makes use of a dynamicreduction network architecture. It is shown that the algorithm isable to capture and mitigate the main effects that normally hinderthe reconstruction of hadronic showers using classicalreconstruction methods, by compensating for fluctuations in themultiplicity, energy, and spatial distributions of the shower'sconstituents. The performance of the algorithm is evaluated usingtest beam data collected in 2018 prototype of the CMS HGCALaccompanied by a section of the CALICE AHCAL prototype. Thecapability of the method to mitigate the impact of energy leakagefrom the calorimeter is also demonstrated.
The Short-Baseline Neutrino (SBN) program features three Liquid Argon Time Projection Chamber (LAr-TPC) detectors positioned along the Booster Neutrino Beam (BNB) axis: the Short Baseline Neutrino Near Detector, MicroBooNE, and the ICARUS T600. As the largest operational LAr-TPC, ICARUS T600 serves as the far detector, located 600 m from the BNB target. While its primary goal is to record neutrino events, it also detects other ionizing events, including cosmic rays. This work focuses on analyzing and modeling detector-specific reconstruction uncertainties in LAr-TPC. These inefficiencies, identified during the Pattern Recognition phase handled by the PANDORA algorithm, impact subsequent Particle Fits and Offline Analysis. Specifically, inaccuracies in track reconstruction can lead to significant physical consequences, such as erroneous particle energy estimates and poor Particle Identification (PID), reducing the efficiency of neutrino event characterization. A key issue addressed is split tracks, caused by missing hits or incomplete track stitching by PANDORA. The aim of this internship is to characterize, model, and quantify the impact of split tracks on track reconstruction.
Position vectors are not uniformly distributed in the plane following the height reconstruction process and thus need to be resampled to a uniform output grid. This process is call regridding.
Reconstruction of flight aerothermal environments often requires the solution of an inverse heat transfer problem, which is an ill-posed problem of specifying boundary conditions from discrete measurements in the interior of the domain. This paper will present the algorithms implemented in the CHAR code for use in reconstruction of EFT-1 flight data and future testing activities. Implementation nuances will be discussed, and alternative hybrid-methods that are permitted by the implementation will be described. Results will be presented for a number of one-dimensional and multi-dimensional problems
Reconstruction of flight aerothermal environments often requires the solution of an inverse heat transfer problem, which is an ill-posed problem of determining boundary conditions from discrete measurements in the interior of the domain. This paper will present the algorithms implemented in the CHAR code for use in reconstruction of EFT-1 flight data and future testing activities. Implementation details will be discussed, and alternative hybrid-methods that are permitted by the implementation will be described. Results will be presented for a number of problems.
In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.
The ultimate goal of this research is to create an anatomically accurate three-dimensional (3D) simulation model of the effects of microgravity in gastrointestinal physiology and to explore the role that such changes may have in the pharmacokinetics of drugs given to the space crews for prevention or therapy. To accomplish this goal the specific aims of this research are: 1) To generate a complete 3-D reconstructions of the human GastroIntestinal (GI) tract of the male and female Visible Humans. 2) To develop and implement time-dependent computer algorithms to simulate the GI motility using the above 3-D reconstruction.
It is recognized that mean square error (MSE) is not a sufficient criterion for determining the acceptability of an image reconstructed from data that has been compressed and decompressed using an encoding algorithm. In the case of Synthetic Aperture Radar (SAR) data, it is also deemed to be insufficient to display the reconstructed image (and perhaps error image) alongside the original and make a (subjective) judgment as to the quality of the reconstructed data. In this paper we suggest a number of additional evaluation criteria which we feel should be included as evaluation metrics in SAR data encoding experiments. These criteria have been specifically chosen to provide a means of ensuring that the important information in the SAR data is preserved. The paper also presents the results of an investigation into the effects of coding on SAR data fidelity when the coding is applied in (1) the signal data domain, and (2) the image domain. An analysis of the results highlights the shortcomings of the MSE criterion, and shows which of the suggested additional criterion have been found to be most important.
A new procedure is presented for the simultaneous coarsening and refinement of three-dimensional unstructured tetrahedral meshes. This algorithm allows for localized grid adaption that is used to capture aerodynamic flow features such as vortices and shock waves in helicopter flowfield simulations. The mesh-adaption algorithm is implemented in the C programming language and uses a data structure consisting of a series of dynamically-allocated linked lists. These lists allow the mesh connectivity to be rapidly reconstructed when individual mesh points are added and/or deleted. The algorithm allows the mesh to change in an anisotropic manner in order to efficiently resolve directional flow features. The procedure has been successfully implemented on a single processor of a Cray Y-MP computer. Two sample cases are presented involving three-dimensional transonic flow. Computed results show good agreement with conventional structured-grid solutions for the Euler equations.
This paper presents a procedure for estimating an accurate model of the CT imaging process including spectral effects. As raw projection data are typically unavailable to the end-user, we adopt a post-processing approach that utilizes the reconstructed images themselves. This approach includes errors from x-ray scatter and the nonidealities of the built-in soft tissue correction into the beam characteristics, which is crucial to beam hardening correction algorithms that are designed to be applied directly to CT reconstructed images. We formulate this approach as a quadratic programming problem and propose two different methods, dimension reduction and regularization, to overcome ill conditioning in the model. For the regularization method we use a statistical procedure, Cross Validation, to select the regularization parameter. We have constructed step-wedge phantoms to estimate the effective beam spectrum of a GE CT-I scanner. Using the derived spectrum, we computed the attenuation ratios for the wedge phantoms and found that the worst case modeling error is less than 3% of the corresponding attenuation ratio. We have also built two test (hybrid) phantoms to evaluate the effective spectrum. Based on these test phantoms, we have shown that the effective beam spectrum provides an accurate model for the CT imaging process. Last, we used a simple beam hardening correction experiment to demonstrate the effectiveness of the estimated beam profile for removing beam hardening artifacts. We hope that this estimation procedure will encourage more independent research on beam hardening corrections and will lead to the development of application-specific beam hardening correction algorithms.
Neutron diffraction based strain mapping is a useful technique for measuring residual strains in additively manufactured (AM) metal parts. The measurement is traditionally done by scanning the sample in a point-wise raster pattern to extract the strain at each position. Since the overall scan can span several hours, adaptive sampling approaches using Bayesian optimization based on Gaussian process (BO-GP) regression have been introduced—demonstrating that even with a fraction of the typically made measurements the dominant strain patterns in the sample can be reconstructed. However, the parameters of the BO-GP algorithm have to be carefully chosen for best performance, and the movement time between arbitrary points can offset the time savings from a reduced number of measurement locations. In this paper, we propose algorithms to refine the BO-GP based methods by using simulations of strain patterns in AM parts based on the materials and the process used to print them. We demonstrate that the simulated strain patterns can be used to help choose better parameters for the BO-GP based framework—leading to low reconstruction error for the final strain pattern. Furthermore, we show that the strain mapping experiment can be initialized with a sampling pattern learnt from the simulation data and ordered to reduce movement time, dramatically enabling reduction in the overall time required to run the baseline BO-GP method.
Image reconstruction from astronomical speckle data using the Knox-Thompson algorithm has now been extensively demonstrated as feasible through analysis and computer simulations. Experimental verification of the technique and its implementation for astronomy places stringent constraints on the recording system and requires complex processing algorithms to handle the many experimental details. In this paper we describe a video recording system specifically designed for the requirements of speckle imaging. This system has been built and tested, yielding high quality image reconstructions from a laboratory-atmospheric simulator for a wide range of input light levels and test objects.