Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian process”

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,711 records · Page 95

Effects of correlated noise on the full-spectrum combining and complex-symbol combining arraying techniques

The process of combining telemetry signals received at multiple antennas, commonly referred to as arraying, can be used to improve communication link performance in the Deep Space Network (DSN). By coherently adding telemetry from multiple receiving sites, arraying produces an enhancement in signal-to-noise ratio (SNR) over that achievable with any single antenna in the array. A number of different techniques for arraying have been proposed and their performances analyzed in past literature. These analyses have compared different arraying schemes under the assumption that the signals contain additive white Gaussian noise (AWGN) and that the noise observed at distinct antennas is independent. In situations where an unwanted background body is visible to multiple antennas in the array, however, the assumption of independent noises is no longer applicable. A planet with significant radiation emissions in the frequency band of interest can be one such source of correlated noise. For example, during much of Galileo's tour of Jupiter, the planet will contribute significantly to the total system noise at various ground stations. This article analyzes the effects of correlated noise on two arraying schemes currently being considered for DSN applications: full-spectrum combining (FSC) and complex-symbol combining (CSC). A framework is presented for characterizing the correlated noise based on physical parameters, and the impact of the noise correlation on the array performance is assessed for each scheme.

Vazirani, P.↗

Surface area and particle size distribution of plutonium oxides derived from the direct oxidation of Pu metal (LA-UR-18-22581)

The intent of this paper is to interrogate the physical characteristics of PuO 2 powders and aggregates produced by direct metal oxidation to examine possible relationships between process parameters and the physical characteristics of the product. The majority of the PuO 2 samples considered were generated during the pre-production and production phases of the Advanced Recovery and Integrated Extraction System (ARIES) program at Los Alamos National Laboratory (LANL). Oxide formed passively on Pu metal surfaces at ambient atmospheric conditions is also characterized. Plutonium oxide powders from the ARIES Direct Metal Oxidation (DMO) furnace and muffle furnace process lines are very consistent in terms of their bulk and tapped density, specific surface area (SSA) and particle size distribution (PSD). The SSAs of calcined PuO 2 generated in the ARIES DMO-2 and muffle furnaces are between 0.1 and 0.5 m2/g. Oxides from both the DMO-2 and muffle furnaces also had similar tri-modal PSDs with maxima at 0.73–3.16 μm, 10.9–18.5 μm, and 41.3–63.0 μm. The PSDs of the oxide formed passively on Pu metal surfaces at ambient conditions were unimodal, with near-Gaussian distributions and D50 close to 1 μm. In conclusion, the oxide formed passively at ambient conditions also had significantly higher SSAs (6–8 m 2 /g).

36 MATERIALS SCIENCE↗

Flower-shaped 2D crystals grown in curved fluid vesicle membranes

Abstract The morphologies of two-dimensional (2D) crystals, nucleated, grown, and integrated within 2D elastic fluids, for instance in giant vesicle membranes, are dictated by an interplay of mechanics, permeability, and thermal contraction. Mitigation of solid strain drives the formation of crystals with vanishing Gaussian curvature (i.e., developable domain shapes) and, correspondingly, enhanced Gaussian curvature in the surrounding 2D fluid. However, upon cooling to grow the crystals, large vesicles sustain greater inflation and tension because their small area-to-volume ratio slows water permeation. As a result, more elaborate shapes, for instance, flowers with bendable but inextensible petals, form on large vesicles despite their more gradual curvature, while small vesicles harbor compact planar crystals. This size dependence runs counter to the known cumulative growth of strain energy of 2D colloidal crystals on rigid spherical templates. This interplay of intra-membrane mechanics and processing points to the scalable production of flexible molecular crystals of controllable complex shape.

36 MATERIALS SCIENCE↗

Results of an integrated structure/control law design sensitivity analysis

A design sensitivity analysis method for Linear Quadratic Cost, Gaussian (LQG) optimal control laws, which predicts change in the optimal control law due to changes in fixed problem parameters using analytical sensitivity equations is discussed. Numerical results of a design sensitivity analysis for a realistic aeroservoelastic aircraft example are presented. In this example, the sensitivity of the optimally controlled aircraft's response to various problem formulation and physical aircraft parameters is determined. These results are used to predict the aircraft's new optimally controlled response if the parameter was to have some other nominal value during the control law design process. The sensitivity results are validated by recomputing the optimal control law for discrete variations in parameters, computing the new actual aircraft response, and comparing with the predicted response. These results show an improvement in sensitivity accuracy for integrated design purposes over methods which do not include changes in the optimal control law. Use of the analytical LQG sensitivity expressions is also shown to be more efficient than finite difference methods for the computation of the equivalent sensitivity information.

Gilbert, Michael G.↗

SHADOZ (Southern Hemisphere Additional Ozonesondes): A Project Overview and New Insights on Tropical Tropospheric Ozone

The SHADOZ (Southern Hemisphere Additional Ozonesondes) ozone sounding network was initiated in 1998 to improve the coverage of tropical in-situ ozone measurements for satellite validation, algorithm development and related process studies. Over 2000 soundings have been archived at the website, http://croc.gsfc.nasa.gov/shadoz, for 12 stations: Ascension Island; Nairobi and Malindi, Kenya; Irene, South Africa; Reunion Island; Watukosek, Java; Fiji; Tahiti; American Samoa; San Cristobal, Galapagos; Natal, Brazil; Paramaribo, Surinam. Key results from SHADOZ will be described from among the following: 1) By using ECC sondes with similar procedures, 5-10% accuracy and precision (1-sigma) of the sonde total ozone measurement is achieved; 2) Week-to-week variability in tropospheric ozone is so great that statistics are frequently not Gaussian; most stations vary up to a factor of 3 in tropospheric column over the course of a year; 3) Longitudinal variability in tropospheric ozone profiles is a consistent feature, with a 10-15 DU column-integrated difference between Atlantic and Pacific sites; this causes a "zonal wave-one" feature in total ozone; 4) The ozone record from Paramaribo, Surinam (6N, 55W) is a marked contrast to southern tropical ozone because Surinam is often north of the Intertropical Convergence Zone.

Thompson, Anne M.↗

Triphenylamine based organic dyes for dye sensitized solar cells: A theoretical approach

The geometry, electronic structure and absorption spectra for newly designed triphenylamine based organic dyes were investigated by density functional theory (DFT) and time dependent density functional theory (TD-DFT) with the Becke 3-Parameter-Lee-Yang-parr(B3LYP) functional, where the 6-31G(d,p) basis set was employed. All calculations were performed using the Gaussian 09 software package. The calculated HOMO and LUMO energies show that charge transfer occurs in the molecule. Ultraviolet–visible (UV–vis) spectrum was simulated by TD-DFT in gas phase. The calculation shows that all of the dyes can potentially be good sensitizers for DSSC. The LUMOs are just above the conduction band of TiO{sub 2} and their HOMOs are under the reduction potential energy of the electrolytes (I{sup −}/I{sub 3}{sup −}) which can facilitate electron transfer from the excited dye to TiO{sub 2} and charge regeneration process after photo oxidation respectively. The simulated absorption spectrum of dyes match with solar spectrum. Frontier molecular orbital results show that among all the three dyes, the “dye 3” can be used as potential sensitizer for DSSC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Iterative reconstruction for CT perfusion with a prior-image induced hybrid nonlocal means regularization: Phantom studies

Purpose: In computed tomography perfusion (CTP) imaging, an initial phase CT acquired with a high-dose protocol can be used to improve the image quality of later phase CT acquired with a low-dose protocol. For dynamic regions, signals in the later low-dose CT may not be completely recovered if the initial CT heavily regularizes the iterative reconstruction process. The authors propose a hybrid nonlocal means (hNLM) regularization model for iterative reconstruction of low-dose CTP to overcome the limitation of the conventional prior-image induced penalty. Methods: The hybrid penalty was constructed by combining the NLM of the initial phase high-dose CT in the stationary region and later phase low-dose CT in the dynamic region. The stationary and dynamic regions were determined by the similarity between the initial high-dose scan and later low-dose scan. The similarity was defined as a Gaussian kernel-based distance between the patch-window of the same pixel in the two scans, and its measurement was then used to weigh the influence of the initial high-dose CT. For regions with high similarity (e.g., stationary region), initial high-dose CT played a dominant role for regularizing the solution. For regions with low similarity (e.g., dynamic region), the regularization relied on a low-dose scan itself. This new hNLM penalty was incorporated into the penalized weighted least-squares (PWLS) for CTP reconstruction. Digital and physical phantom studies were performed to evaluate the PWLS-hNLM algorithm. Results: Both phantom studies showed that the PWLS-hNLM algorithm is superior to the conventional prior-image induced penalty term without considering the signal changes within the dynamic region. In the dynamic region of the Catphan phantom, the reconstruction error measured by root mean square error was reduced by 42.9% in PWLS-hNLM reconstructed image. Conclusions: The PWLS-hNLM algorithm can effectively use the initial high-dose CT to reconstruct low-dose CTP in the stationary region while reducing its influence in the dynamic region.

60 APPLIED LIFE SCIENCES↗

Laser transit anemometer software development program

Algorithms were developed for the extraction of two components of mean velocity, standard deviation, and the associated correlation coefficient from laser transit anemometry (LTA) data ensembles. The solution method is based on an assumed two-dimensional Gaussian probability density function (PDF) model of the flow field under investigation. The procedure consists of transforming the data ensembles from the data acquisition domain (consisting of time and angle information) to the velocity space domain (consisting of velocity component information). The mean velocity results are obtained from the data ensemble centroid. Through a least squares fitting of the transformed data to an ellipse representing the intersection of a plane with the PDF, the standard deviations and correlation coefficient are obtained. A data set simulation method is presented to test the data reduction process. Results of using the simulation system with a limited test matrix of input values is also given.

Abbiss, John B.↗

Digital TV processing system

Two digital video data compression systems directly applicable to the Space Shuttle TV Communication System were described: (1) For the uplink, a low rate monochrome data compressor is used. The compression is achieved by using a motion detection technique in the Hadamard domain. To transform the variable source rate into a fixed rate, an adaptive rate buffer is provided. (2) For the downlink, a color data compressor is considered. The compression is achieved first by intra-color transformation of the original signal vector, into a vector which has lower information entropy. Then two-dimensional data compression techniques are applied to the Hadamard transformed components of this last vector. Mathematical models and data reliability analyses were also provided for the above video data compression techniques transmitted over a channel encoded Gaussian channel. It was shown that substantial gains can be achieved by the combination of video source and channel coding.

Source record↗

Sculpting 2D Crystals via Membrane Contractions before and during Solidification

When phospholipids crystallize within the otherwise fluid membranes of giant unilamellar vesicles, the resulting molecularly thin “2D” solids exhibit great variety in their morphology evolution. For example, within membranes containing moderate amounts of the crystallizing component, crystals grow with a fixed morphology depending on vesicle size. Conversely for membranes containing large amounts of the crystallizing species, we find small compact crystals on vesicles of all sizes. However, on large vesicles, growing crystals sprout flower petals that lengthen progressively. These behaviors result from two combined mechanisms: first, like other 2D solids, the shear rigidity of phospholipid crystals renders them intolerant to morphologies with nonzero Gaussian curvature. As a result and especially at elevated membrane tension, the cost of bending elasticity is reduced at the expense of line energy by the formation of flowers as opposed to compact crystals. Second, the composition-dependent tension rise during cooling relaxes via water permeation of the membrane with a time constant scaling as R2. The amount of crystal formed for a small decrease in temperature determines this composition-dependent increase in stress from thermal contractions versus solidification. Surface Evolver computations were motivated using the predicted tension evolution to develop a processing space that maps to experimental observations for initial and growing crystal morphology. Important variable groups are identified, including a scaled ratio of bending to line energy, a vesicle-size-independent group for membrane contractions, and a time constant for stress relaxation. Though processing stresses ultimately relax, the crystal morphology persists well beyond the processing window.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Design and analysis of a wake steering controller with wind direction variability

Wind farm control strategies are being developed to mitigate wake losses in wind farms, increasing energy production. Wake steering is a type of wind farm control in which a wind turbine's yaw position is misaligned from the wind direction, causing its wake to deflect away from downstream turbines. Current modeling tools used to optimize and estimate energy gains from wake steering are designed to represent wakes for fixed wind directions. However, wake steering controllers must operate in dynamic wind conditions and a turbine's yaw position cannot perfectly track changing wind directions. Research has been conducted on robust wake steering control optimized for variable wind directions. In this paper, the design and analysis of a wake steering controller with wind direction variability is presented for a two-turbine array using the FLOw Redirection and Induction in Steady State (FLORIS) control-oriented wake model. First, the authors propose a method for modeling the turbulent and low-frequency components of the wind direction, where the slowly varying wind direction serves as the relevant input to the wake model. Next, we explain a procedure for finding optimal yaw offsets for dynamic wind conditions considering both wind direction and yaw position uncertainty. We then performed simulations with the optimal yaw offsets applied using a realistic yaw offset controller in conjunction with a baseline yaw controller, showing good agreement with the predicted energy gain using the probabilistic model. Using the Gaussian wake model in FLORIS as an example, we compared the performance of yaw offset controllers optimized for static and dynamic wind conditions for different turbine spacings and turbulence intensity values, assuming uniformly distributed wind directions. For a spacing of five rotor diameters and a turbulence intensity of 10 %, robust yaw offsets optimized for variable wind directions yielded an energy gain equivalent to 3.24 % of wake losses recovered, compared to 1.42 % of wake losses recovered with yaw offsets optimized for static wind directions. In general, accounting for wind direction variability in the yaw offset optimization process was found to improve energy production more as the separation distance increased, whereas the relative improvement remained roughly the same for the range of turbulence intensity values considered.

17 WIND ENERGY↗

Evaluating Subgrid-Scale Properties of Low Clouds over the Beaufort Sea in Arctic System Reanalysis using ARISE Airborne In Situ Observations

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. To better understand the magnitude and possible causes of this problem, we analyze the disagreements between the Arctic System Reanalysis version 2 (ASR) and data taken during the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign conducted over the Beaufort Sea. Previous analysis of grid scale-averaged meteorological variables revealed a lack of cloud water produced in ASR, related to warm/dry biases, a failure to produce sufficient cloud water in conditions that are observed to be favorable for large cloud water values, and sampling issues with the ARISE flight paths. However, ARISE collected meteorological observations at one second intervals, offering spatial data resolution comparable to the subgrid scale of ASR. Subgrid scale properties and processes must be parameterized in conventional atmospheric models and reanalyses such as ASR, and often rely on assumed distributions of meteorological variables about the grid box mean state. To enable an examination of the realism of ASR’s subgrid scale properties, we use high resolution data from ARISE determine the “subgrid scale” distributions of cloud water and thermodynamic variables in the real atmosphere. Preliminary results show that the subgrid scale distribution of thermodynamic variables about the grid scale mean resembles the triangular distribution assumed in many bulk microphysical parameterization schemes, but also with tails resembling Gaussian distributions. The main goal is to understand the partitioning of the total available atmospheric water into water vapor and cloud water. Models typically covert some percentage of total water in a grid box to cloud water based on the estimated amount of supersaturation from the assumed vapor distribution. The ARISE measurements can reveal the relationship between subgrid scale supersaturation and the partition between cloud water and water vapor.

Jason Brant Dodson↗

An Integrated Signaling-Encryption Mechanism to Reduce Error Propagation in Wireless Communications: Performance Analyses

Traditional encryption techniques require packet overhead, produce processing time delay, and suffer from severe quality of service deterioration due to fades and interference in wireless channels. These issues reduce the effective transmission data rate (throughput) considerably in wireless communications, where data rate with limited bandwidth is the main constraint. In this paper, performance evaluation analyses are conducted for an integrated signaling-encryption mechanism that is secure and enables improved throughput and probability of bit-error in wireless channels. This mechanism eliminates the drawbacks stated herein by encrypting only a small portion of an entire transmitted frame, while the rest is not subject to traditional encryption but goes through a signaling process (designed transformation) with the plaintext of the portion selected for encryption. We also propose to incorporate error correction coding solely on the small encrypted portion of the data to drastically improve the overall bit-error rate performance while not noticeably increasing the required bit-rate. We focus on validating the signaling-encryption mechanism utilizing Hamming and convolutional error correction coding by conducting an end-to-end system-level simulation-based study. The average probability of bit-error and throughput of the encryption mechanism are evaluated over standard Gaussian and Rayleigh fading-type channels and compared to the ones of the conventional advanced encryption standard (AES).

advanced encryption standard (AES)↗

Combined optimization of image-gathering and image-processing systems for scene feature detection

The relationship between the image gathering and image processing systems for minimum mean squared error estimation of scene characteristics is investigated. A stochastic optimization problem is formulated where the objective is to determine a spatial characteristic of the scene rather than a feature of the already blurred, sampled and noisy image data. An analytical solution for the optimal characteristic image processor is developed. The Wiener filter for the sampled image case is obtained as a special case, where the desired characteristic is scene restoration. Optimal edge detection is investigated using the Laplacian operator x G as the desired characteristic, where G is a two dimensional Gaussian distribution function. It is shown that the optimal edge detector compensates for the blurring introduced by the image gathering optics, and notably, that it is not circularly symmetric. The lack of circular symmetry is largely due to the geometric effects of the sampling lattice used in image acquisition. The optimal image gathering optical transfer function is also investigated and the results of a sensitivity analysis are shown.

Halyo, Nesim↗

Stochastic generation of explicit pore structures by thresholding Gaussian random fields

We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone. -- Graphical abstract: -- Highlights: •An efficient method to stochastically generate realistic pore structures is provided. •Samples are generated by applying a level threshold to a Gaussian field realization. •Two user prescribed quantities determine the topology and geometry of the pore space. •Multiple pore structures and preferential flow directions can be produced. •A pore space based on Berea sandstone is generated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

From asymmetric to symmetric fission in the fermium isotopes within the time-dependent generator-coordinate-method formalism

Background: Predicting the properties of neutron-rich nuclei far from the valley of stability is one of the major challenges of modern nuclear theory. In heavy and superheavy nuclei, a difference of only a few neutrons is sufficient to change the dominant fission mode. A theoretical approach capable of predicting such rapid transitions for neutron-rich systems would be a valuable tool to better understand $r$-process nucleosynthesis or the decay of superheavy elements. Purpose: In this work, we investigate for the first time the transition from asymmetric to symmetric fission through the calculation of primary fission yields with the time-dependent generator coordinate method (TDGCM). We choose here the transition in neutron-rich fermium isotopes, which was the first to be observed experimentally in the late 1970s and is often used as a benchmark for theoretical studies. Methods: We compute the primary fission fragment mass and charge yields for 254 Fm, 256 Fm, and 258 Fm from the TDGCM under the Gaussian overlap approximation. The static part of the calculation (generation of a potential energy surface) consists of a series of constrained Hartree-Fock-Bogoliubov calculations based on the D1S, D1M, or D1N parametrization of the Gogny effective interaction in a two-center harmonic oscillator basis. The two-dimensional dynamics in the collective space spanned by the quadrupole and octupole moments ( $\hat{Q}_{20}$, $\hat{Q}_{30}$) is then computed with the finite element solver Felix-2.0. Results: The available experimental data and the TDGCM post-dictions are consistent and agree especially on the position in the fermium isotopic chain at which the transition occurs. In addition, the TDGCM predicts two distinct asymmetric modes for the fission of 254 Fm. Conclusions: Thanks to its intrinsic accounting of shell effects and to its ability to describe the dynamics of the system up to configurations close to scission, the TDGCM is able to describe qualitatively the fission yield transition in the neutron-rich fermium isotopes. This makes it a promising tool to study the evolution of the fission yields far from the valley of stability. The main limitation of the method lies in the presence of discontinuities in the two-dimensional manifold of generator states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Range data description based on multiple characteristics

An algorithm for describing range images based on Mean curvature (H) and Gaussian curvature (K) is presented. Range images are unique in that they directly approximate the physical surfaces of a real world 3-D scene. The curvature parameters are derived from the fundamental theorems of differential geometry and provides visible invariant pixel labels that can be used to characterize the scene. The sign of H and K can be used to classify each pixel into one of eight possible surface types. Due to the sensitivity of these parameters to noise the resulting HK-sing map does not directly identify surfaces in the range images and must be further processed. A region growing algorithm based on modeling the scene points with a Markov Random Field (MRF) of variable neighborhood size and edge models is suggested. This approach allows the integration of information from multiple characteristics in an efficient way. The performance of the proposed algorithm on a number of synthetic and real range images is discussed.

Al-Hujazi, Ezzet↗

Forces and stress in second order Møller-Plesset perturbation theory for condensed phase systems within the resolution-of-identity Gaussian and plane waves approach

The forces acting on the atoms as well as the stress tensor are crucial ingredients for calculating the structural and dynamical properties of systems in the condensed phase. Here, these derivatives of the total energy are evaluated for the second-order Møller-Plesset perturbation energy (MP2) in the framework of the resolution of identity Gaussian and plane waves method, in a way that is fully consistent with how the total energy is computed. This consistency is non-trivial, given the different ways employed to compute Coulomb, exchange, and canonical four center integrals, and allows, for example, for energy conserving dynamics in various ensembles. Based on this formalism, a massively parallel algorithm has been developed for finite and extended system. The designed parallel algorithm displays, with respect to the system size, cubic, quartic, and quintic requirements, respectively, for the memory, communication, and computation. All these requirements are reduced with an increasing number of processes, and the measured performance shows excellent parallel scalability and efficiency up to thousands of nodes. Additionally, the computationally more demanding quintic scaling steps can be accelerated by employing graphics processing units (GPU’s) showing, for large systems, a gain of almost a factor two compared to the standard central processing unit-only case. In this way, the evaluation of the derivatives of the RI-MP2 energy can be performed within a few minutes for systems containing hundreds of atoms and thousands of basis functions. With good time to solution, the implementation thus opens the possibility to perform molecular dynamics (MD) simulations in various ensembles (microcanonical ensemble and isobaric-isothermal ensemble) at the MP2 level of theory. Geometry optimization, full cell relaxation, and energy conserving MD simulations have been performed for a variety of molecular crystals including NH{sub 3}, CO{sub 2}, formic acid, and benzene.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗