Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “observational methods”

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,747 records · Page 97

Shock interference studies on a circular cylinder at Mach 16

This paper presents the results of a computational study on the shock interference problems on a cylindrical body typical of an engine inlet cowl leading edge at a nominal Mach number of 16. The two-dimensional Navier-Stokes equations are solved assuming the flow to be in chemical and thermal equilibrium and using a finite element method. The algorithm employs a cell-centered fully implicit upwind scheme. Adaptively generated unstructured meshes of triangles and quadrilaterals are employed. Under certain conditions the finite element code resulted in oscillatory solutions for shock interference at Mach 16. Some of the causes of the unsteady behavior are identified, and to the extent possible, such situations are avoided in the present application. Two shock interference conditions involving a Type IV (supersonic jet) and a Type III (attaching shear layer) are solved. The results are compared with available experimental data and reasonable agreement is observed. A semi-empirical method is also used to estimate the maximum surface heat flux and static pressure.

Prabhu, Ramadas K.↗

Foldamer architectures of triazine-based sequence-defined polymers investigated with molecular dynamics simulations and enhanced sampling methods

Triazine-based sequence-defined polymers have recently been developed that are biomimetic and robust (Grate \latin{et al.} Angewandte Chemie 2016). In molecular dynamics (MD) simulations, the triazine polymers were shown to form linear nanorod foldamers through hydrogen bonding and $\pi$-$\pi$ interactions. The nanorod foldamers have motifs resembling those of DNA, $\alpha$-helices, and $\beta$-sheets, and have potential to be useful building blocks for new macromolecules and materials. To understand the formation of nanorod foldamers, we investigate how linker structures in the middle of the triazine polymers lead to folding using MD simulations. We found that a variety of linkers can participate in folding, but that specific linker structures are more favorable than others, depending on the polymer length. Folding of hexamers into well-defined nanorod foldamers was most favorable with pentanediamine and \emph{ortho}-xylenediamine linkers in the center of the polymers. Foldamers with \emph{ortho}-xylenediamine linkers in the center were investigated for longer polymers, i.e., octamers and decamers, using two different enhanced sampling methods, since regular MD simulations had failed to show any folding for these longer polymers. In particular, the recently developed concurrent adaptive sampling (CAS) algorithm (Ahn \latin{et al.} J. Chem. Phys. 2017) and replica exchange molecular dynamics (REMD) were used. We found that the two enhanced sampling methods did lead to the observation of foldamers, and that REMD revealed new foldamer architectures where \emph{cis}-\emph{trans} isomerizations had occurred. Foldamer formation, diversity, and the strengths and limitations of simulation techniques are discussed. These findings provide new insights into the diversity of foldamer architectures for a new type of biomimetic synthetic polymer.

Foldamers, Nanorods, Free energy landscapes, Molec↗

Turbulence in molecular clouds - A new diagnostic tool to probe their origin

A method is presented to uncover the instability responsible for the type of turbulence observed in molecular clouds and the value of the physical parameters of the 'placental medium' from which turbulence originated. The method utilizes the observational relation between velocities and sizes of molecular clouds, together with a recent model for large-scale turbulence (constructed by Canuto and Goldman, 1985).

Canuto, V. M.↗

Structural properties of optically clear bacterial cellulose produced by Komagataeibacter hansenii using arabitol

Here, bacterial cellulose (BC) exhibits beneficial properties for use in biomedical applications but is limited by its lack of tunable transparency capabilities. To overcome this deficiency, a novel method to synthesize transparent BC materials using an alternative carbon source, namely arabitol, was developed. Characterization of the BC pellicles was performed for yield, transparency, surface morphology, and molecular assembly. Transparent BC was produced using mixtures of glucose and arabitol. Zero percent arabitol pellicles exhibited 25% light transmittance, which increased with increasing arabitol concentration through to 75% light transmittance. While transparency increased, overall BC yield was maintained indicating that the altered transparency may be induced on a micro-scale rather than a macro-scale. Significant differences in fiber diameter and the presence of aromatic signatures were observed. Overall, this research outlines methods for producing BC with tunable optical transparency, while also bringing new insight to insoluble components of exopolymers produced by Komagataeibacter hansenii.

59 BASIC BIOLOGICAL SCIENCES↗

Cognitive simulation models for inertial confinement fusion: Combining simulation and experimental data

The design space for inertial confinement fusion (ICF) experiments is vast, and experiments are extremely expensive. Researchers rely heavily on computer simulations to explore the design space in search of high-performing implosions. However, ICF multiphysics codes must make simplifying assumptions, and thus deviate from experimental measurements for complex implosions. For more effective design and investigation, simulations require input from past experimental data to better predict future performance. In this work, we describe a cognitive simulation method for combining simulation and experimental data into a common, predictive model. This method leverages a machine learning technique called “transfer learning,” the process of taking a model trained to solve one task, and partially retraining it on a sparse dataset to solve a different, but related task. In the context of ICF design, neural network models are trained on large simulation databases and partially retrained on experimental data, producing models that are far more accurate than simulations alone. Here, we demonstrate improved model performance for a range of ICF experiments at the National Ignition Facility and predict the outcome of recent experiments with less than 10% error for several key observables. We discuss how the methods might be used to carry out a data-driven experimental campaign to optimize performance, illustrating the key product—models that become increasingly accurate as data are acquired.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Bilinear modeling and nonlinear estimation

New methods are illustrated for online nonlinear estimation applied to the lateral deflection of an elastic beam on board measurements of angular rates and angular accelerations. The development of the filter equations, together with practical issues of their numerical solution as developed from global linearization by nonlinear output injection are contrasted with the usual method of the extended Kalman filter (EKF). It is shown how nonlinear estimation due to gyroscopic coupling can be implemented as an adaptive covariance filter using off-the-shelf Kalman filter algorithms. The effect of the global linearization by nonlinear output injection is to introduce a change of coordinates in which only the process noise covariance is to be updated in online implementation. This is in contrast to the computational approach which arises in EKF methods arising by local linearization with respect to the current conditional mean. Processing refinements for nonlinear estimation based on optimal, nonlinear interpolation between observations are also highlighted. In these methods the extrapolation of the process dynamics between measurement updates is obtained by replacing a transition matrix with an operator spline that is optimized off-line from responses to selected test inputs.

Dwyer, Thomas A. W., III↗

Variation of stratospheric NO2 during the solar eclipse

The numerical modeling of the changes of the concentration of trace gases in the atmosphere during the eclipse shows that the NO2 total content in the vertical column increases approximately by 80 percent. The first observations of the NO2 total content during the eclipse of 1981 have given 60 plus or minus 20 percent. In the observations of 1990 the more precise methods and instruments for stratospheric NO2 measurements were used. The surface ozone, NO, and NO2 concentrations were under control. The results of the observations give the increasing of the stratospheric NO2 during the eclipse by 55 plus or minus 6 percent. The maximum increasing of the NO2 content is observed at the moment of the maximum phase.

Elansky, Nikolay F.↗

Conditional Karhunen-Loève expansion for uncertainty quantification and active learning in partial differential equation models

We use a conditional Karhunen-Lo` eve (KL) model to quantify and reduce uncertainty in a stochastic partial differential equation (SPDE) problem with partially-known space-dependent coefficient, Y (x). We assume that a small number of Y (x) measurements are available and model Y (x) with a KL expansion. We achieve reduction in uncertainty by conditioning the KL expansion coefficients on measurements. We consider two approaches for conditioning the KL expansion: In Approach 1, we condition the KL model first and then truncate it. In Approach 2, we first truncate the KL expansion and then condition it. We employ the conditional KL expansion together with Monte Carlo and sparse grid collocation methods to compute the moments of the solution of the SPDE problem. Uncertainty of the problem is further reduced by adaptively selecting additional observation locations using two active learning methods. Method 1 minimizes the variance of the PDE coefficient, while Method 2 minimizes the variance of the solution of the PDE. We demonstrate that conditioning leads to dimension reduction of the KL representation of Y (x). For a linear diffusion SPDE with uncertain log-normal coefficient, we show that Approach 1 provides a more accurate approximation of the conditional log-normal coefficient and solution of the SPDE than Approach 2 for the same number of random dimensions in a conditional KL expansion. Furthermore, Approach 2 provides a good estimate for the number of terms of the truncated KL expansion of the conditional field of Approach 1. Finally, we demonstrate that active learning based on Method 2 is more efficient for uncertainty reduction in the SPDE’s states (i.e., it leads to a larger reduction of the variance) than active learning using Method 2.

Conditioned Karhunen-Lo` eve expanion, machine lea↗

Nimbus 7 CZCS - Reduction of its radiometric sensitivity with time

Preliminary results are described for an effort to quantify the sensitivity decay of a radiometry sensor (the Coastal Zone Color Scanner or CZCS aboard Nimbus 7). The method used in the study is to (1) compute the water-leaving radiance for imagery acquired in regions where this radiance is known or can be independently estimated, and (2) adjust the sensor calibration to force agreement between the two radiances. Decay factors for orbit numbers from 0 to 20,000 are plotted, and surface and space measurements are compared for the Gulf Stream and the Northern Sargasso Sea at different seasons. The fact that a seasonal variability in the chlorophyll a concentration in the Sargasso Sea was found in the sensor analysis (apparently the first such satellite observation) increases confidence in the method.

Gordon, H. R.↗

Modeling and data assimilation for Mars Observer

The Mars Observer mission will return a substantial set of meteorological data for the atmosphere of Mars, principally in the form of thermal infrared soundings from the Pressure Modulator Infrared Radiometer (PMIRR). The dataset will enable global studies of a range of dynamical and physical phenomena in the martian atmosphere. The solar orbital configuration of Mars Observer, however, means that the coverage of the martian atmosphere and surface by remote sounding instruments will be asynchronous. This leads to difficulties of interpretation when studying synoptic-scale phenomena with timescales of a few days. In an attempt to overcome such problems we propose to use data assimilation techniques, as currently employed for operational weather forecasting on the Earth, in conjunction with a Martian General Circulation Model (MGCM) that is under development at Oxford and Reading Universities. An overview of the model and assimilation scheme will be given, with reference to problems encountered so far in applying such methods to remote-sensing observations of the martian atmosphere.

Lewis, Stephen R.↗

Data for Clumping Index Estimation With 30°-tilted Cameras in Row Crops: Evaluation of Methods and Segment Size Effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

Modeling↗

Locating Surface Explosions by Combining Seismic and Infrasound Data

An outline of a Bayesian source location framework for using seismic and acoustic observations is developed and tested on synthetic and real data. Seismic and acoustic phenomena are both commonly used in detection and location of a variety of natural or man–made events, such as volcanic eruptions, quarry blasts, and military exercises. Typically, seismic and acoustic observations have been utilized independently of each other. Here, we outline a Bayesian formulation for combining the two observations in a single estimate of the location and origin time. Using realistic estimates of uncertainty, we subsequently explore how combining the different observation types can benefit event location at local to near–regional distances. We apply the method to synthetic data and to real observations from a mining blast in Bingham Mine in Utah. Our findings suggest that, for relatively sparse or azimuthally limited observations, the relative strengths of the two different phenomenologies enable more precise joint–event localization than either seismic or infrasonic measurements alone.

58 GEOSCIENCES↗

Maximum a posteriori classification of multifrequency, multilook, synthetic aperture radar intensity data

We present a maximum a posteriori (MAP) classifier for classifying multifrequency, multilook, single polarization SAR intensity data into regions or ensembles of pixels of homogeneous and similar radar backscatter characteristics. A model for the prior joint distribution of the multifrequency SAR intensity data is combined with a Markov random field for representing the interactions between region labels to obtain an expression for the posterior distribution of the region labels given the multifrequency SAR observations. The maximization of the posterior distribution yields Bayes's optimum region labeling or classification of the SAR data or its MAP estimate. The performance of the MAP classifier is evaluated by using computer-simulated multilook SAR intensity data as a function of the parameters in the classification process. Multilook SAR intensity data are shown to yield higher classification accuracies than one-look SAR complex amplitude data. The MAP classifier is extended to the case in which the radar backscatter from the remotely sensed surface varies within the SAR image because of incidence angle effects. The results obtained illustrate the practicality of the method for combining SAR intensity observations acquired at two different frequencies and for improving classification accuracy of SAR data.

Rignot, E.↗

On the determination of the number of O stars in H II regions and starburst galaxies

The hot star population in H II regions, H II galaxies, and starburst galaxies is often described in terms of the number of 'equivalent' O stars of a single representative subtype and luminosity class needed to produce the ionizing luminosity deduced from the nebular recombination lines in the optical spectra. In this paper we define conversion factors eta(sub 0), eta(sub 1), and zeta(sub 5000) with which the total number of O V stars and their flux contribution at 5000 A can be derived from the number of these 'equivalent' stars. These quantities depend primarily on three parameters: the slope and upper mass limit of the stellar mass function and the metallicity of the region. Using the latest stellar atmosphere and evolution models, we calculate eta(sub 0), eta(sub 1), and zeta(sub 5000) for a large number of values of these parameters. The results are presented in tabular as well as graphical form. We apply our results to two H II regions for which the hot star population are known and find that the predicted numbers of O stars agree well the observed counts. In addition, we describe a method by which the values of eta(sub 0) and eta(sub 1) and the observed emission-line fluxes can be used to place constraints on the allowed values of the slope and upper mass limit of the stellar mass function in a region.

Vacca, William D.↗

Microstructural Development of Directionally Solidified Hg(1-x)Zn(x)Se Alloys

Hg(1-x)Zn(x)Se alloys have been studied as an alternative to HgCdTe for the detection of electromagnetic radiation, because the shorter ZnSe and HgSe bonds were predicted to improve lattice stability. Previous studies showed that optical properties were stable with time and that the microhardness was higher than for HgTe based alloys. However, for this material to be commercially viable, compositional variations must be reduced, and the distribution of defects induced by the growth process must be understood. The basic structural properties of bulk materials were characterized, and dislocation etch pit densities were measured to determine whether the addition of Zn to the HgSe lattice reduced the potential for dislocation formation. Hg(0.9)Zn(0.1)Se alloys were directionally solidified using a modified Bridgman-Stockbarger method and in an applied magnetic field. Radial compositional variations were greatly reduced when solidification occurred in an applied magnetic field. The combination of a convex liquid-solid interface shape and a method for reduced wetting produced boules that were a single crystal after reaching full diameter. Observed surface features indicated ampoule wetting could be eliminated using a graphite getter. Microstructural characteristics were greatly improved over HgCdTe alloys, In six boules, a total of only one twin was observed, and no inclusions were found. A method for polishing and producing dislocation etch pits was developed for these alloys, revealing dislocation etch pit densities one to two orders of magnitude less than HgTe based alloys. A kink in the thermal profile during processing of one boule generated more dislocations than did lattice mismatch due to compositional variations. The results of this investigation indicate this alloy has improved microstructural properties and resistance to dislocation formation compared with similar II-VI alloys, and should be further investigated for the detection of electromagnetic radiation.

Cobb, S. D.↗

Reducing Errors in Velocity–Azimuth Display (VAD) Wind and Deformation Retrievals from Airborne Doppler Radars in Convective Environments

The present study describes methods to reduce the uncertainty of velocity–azimuth display (VAD) wind and deformation retrievals from downward-pointing, conically scanning, airborne Doppler radars. These retrievals have important applications in data assimilation and real-time data processing. Several error sources for VAD retrievals are considered here, including violations to the underlying wind field assumptions, Doppler velocity noise, data gaps, temporal variability, and the spatial weighting function of the VAD retrieval. Specific to airborne VAD retrievals, we also consider errors produced due to the radar scans occurring while the instrument platform is in motion. While VAD retrievals are typically performed using data from a single antenna revolution, other strategies for selecting data can be used to reduce retrieval errors. Four such data selection strategies for airborne VAD retrievals are evaluated here with respect to their effects on the errors. These methods are evaluated using the second hurricane nature run numerical simulation, analytic wind fields, and observed Doppler radar radial velocities. The proposed methods are shown to reduce the median absolute error of the VAD wind retrievals, especially in the vicinity of deep convection embedded in stratiform precipitation. The median absolute error due to wind field assumption violations for the along-track and for the across-track wind is reduced from 0.36 to 0.08 m s−1 and from 0.35 to 0.24 m s−1, respectively. Although the study focuses on Doppler radars, the results are equally applicable to conically scanning Doppler lidars as well.

Charles N Helms↗

SECARB Annual Report on In-Situ Monitoring Program

The Southeast Carbon Sequestration Partnership Phase III Anthropogenic Test at Citronelle Field is the United States’ largest prototype of commercial CO 2 capture from a coal-fired power station and transportation, injection, storage and monitoring. This report details the results of the second and third year of injection monitoring results for the Anthropogenic Test. Monitoring methods include both direct and indirect measurements. The Class V Experimental Permit requires extensive monitoring of injected CO 2 at the surface and at deep and shallow depths. The safety, integrity and information objectives of the CO 2 monitoring plan are in place to: (1) create and sustain wellbore integrity; (2) assure safe CO 2 injection operations; 3) verify the location and migration of the injected CO 2 plume; and (4) monitor for any CO 2 leakage. A robust suite of monitoring approaches are utilized for the insitu monitoring and include, pressure monitoring, various seismic surveys, pulsed neutron logging, reservoir fluid sampling and injection surveys. Data and observations collected using these monitoring methods will be used in the numerical model to inform the post-injection site care monitoring plan. This report describes actions pertaining to in situ monitoring at the Anthropogenic Test Site that were undertaken during the period of October 2013 to January 2015. This report details the results of the second and third year of injection monitoring results for the Anthropogenic Test. Work Product 9.2.2.b details the shallow monitoring results for 2014. CO 2 injection operations began on August 20, 2012 and terminated on September 1st 2014.

01 COAL, LIGNITE, AND PEAT↗

Evaluation of glass model precipitation forecasts for North America during SOP-1

A study has been made of the accuracy of the GLAS fourth order forecast model's precipitation forecasts over North America during the first Special Observing Period (SOP-1) of FGGE and of the impact of the FGGE special observing systems on these forecasts. Fourteen 120 h predictions were generated using the coarse 4 deg latitude by 5 deg longitude version of the GLAS model from both the FGGE and NOSAT assimilation cycles. These forecasts were then verified against a detailed set of precipitation observations. Separate verifications were performed for precipitation accumulations or = .01 in., or = .1 in., and or = 1 in. occuring in 6, 12 and 24 time periods. In addition, three different methods of verification were applied. In the first method, model precipitation forecasts are verified against the greatest precipitation amount observed within a 4 deg latitude by 5 deg longitude gridbox centered on each gridpoint. For the prediction at a gridpoint to be verified only one observation within the gridbox is required. In the second method, a weighted mean of the closest observations to the gridpoint is used, provided that at least two observations on opposite sides of a gridpoint are available. For the third method, the model precipitation forecasts were interpolated to the observation locations.

Atlas, R.↗