Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “normalization method”

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 955 records · Page 53

Progress in multi-dimensional upwind differencing

Multi-dimensional upwind-differencing schemes for the Euler equations are reviewed. On the basis of the first-order upwind scheme for a one-dimensional convection equation, the two approaches to upwind differencing are discussed: the fluctuation approach and the finite-volume approach. The usual extension of the finite-volume method to the multi-dimensional Euler equations is not entirely satisfactory, because the direction of wave propagation is always assumed to be normal to the cell faces. This leads to smearing of shock and shear waves when these are not grid-aligned. Multi-directional methods, in which upwind-biased fluxes are computed in a frame aligned with a dominant wave, overcome this problem, but at the expense of robustness. The same is true for the schemes incorporating a multi-dimensional wave model not based on multi-dimensional data but on an 'educated guess' of what they could be. The fluctuation approach offers the best possibilities for the development of genuinely multi-dimensional upwind schemes. Three building blocks are needed for such schemes: a wave model, a way to achieve conservation, and a compact convection scheme. Recent advances in each of these components are discussed; putting them all together is the present focus of a worldwide research effort. Some numerical results are presented, illustrating the potential of the new multi-dimensional schemes.

Vanleer, Bram↗

A simple method of calculating power-law velocity profile exponents from experimental data

Analytical expressions for the effects of compressibility and heat transfer on laminar and turbulent shape factors H have been developed. Solving the turbulent equation for the power law velocity profile exponent N has resulted in a simple technique by which the N values of experimental turbulent profiles can be calculated directly from the integral parameters. Thus the data plotting, curve fitting, and slope measuring, which is the normal technique of obtaining experimental N values, is eliminated. The N values obtained by this method should be within the accuracy with which they could be measured.

Allen, J. M.↗

Progression of damage and fracture in composites under dynamic loading

A new computational simulation method is presented to evaluate the dynamic aspects of composite structural response and durability that could not be simulated previously. Composite structural behavior under any loading, geometry, composite system, laminate configuration, and boundary conditions can now be simulated. Structural degradation, delamination, fracture, and damage propagation are included in the simulation. An angle-plied composite plate structure under normal impact loading is used as an example to demonstrate the versatility of the simulation method.

Minnetyan, Levon↗

Progression of damage and fracture in composites under dynamic loading

A new computational simulation method is presented to evaluate the dynamic aspects of composite structural response and durability that have not been simulated previously. Composite structural behavior under any loading condition, geometry, composite system, laminate configuration, and boundary conditions can now be simulated. Structural degradation, delamination, fracture, and damage propagation are included in the simulation. An angle-plied composite plate structure under normal impact loading is used as an example to demonstrate the versatility of the simulation method.

Minnetyan, Levon↗

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗

Discrete element model for powder grain interactions under high compressive stress

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. Finally, as validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

36 MATERIALS SCIENCE↗

Measurement and modeling of methane diffusion in hydrocarbon mixtures

Methane (CH 4 ) dissolution and diffusive mass transfer in liquid hydrocarbon mixtures is of key interest in the context of enhanced oil recovery from tight shales. In this paper, we have studied CH 4 dissolution and diffusion in normal alkane mixtures at a temperature of 50 ⁰C and at pressures of ~8 MPa. For the measurement of CH 4 dissolution/diffusion in bulk liquid hydrocarbon mixtures, we have utilized a high pressure and temperature, constant-volume diffusion (CV-D) setup. Here, we have studied CH 4 diffusion in mixtures with three long-chain normal alkanes: decane (C 10 ), dodecane (C 12 ) and hexadecane (C 16 ). We have measured CH 4 solubility and diffusion in its binary mixtures with each of these three normal alkanes, as well as in ternary and quaternary mixtures. During the experiments, the swelling of the liquid mixture due to CH 4 dissolution was measured in situ via a cathetometer and was subsequently integrated into the data analysis. A key conclusion from this study is that the solubility and transport properties of the multicomponent mixtures can be predicted accurately from binary mixture measurements using an appropriate Equation of State (EOS) and Wilke’s simplification of the classical Maxwell-Stefan (MS) diffusion theory. This observation can facilitate accurate prediction of diffusive mass transfer in more complex liquid hydrocarbon mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A data analytical approach for assessing the efficacy of Operational Technology active defenses against insider threats

In recent years, the need for Operational Technology (OT) defenses has been recognized, serving as an additional line of defense when Information Technology (IT) defenses are bypassed. This is no longer considered an uncommon possibility when dealing with advanced persistent threat (APT) actors expected to be state-sponsored and receiving insider assistance. In these extreme adversarial situations, OT defenses aim to provide another layer of defense for the system, introduced directly at the physical process level, as described by the sensors data, the system model, and control actions. Just like IT defenses, two schools of thought, i.e., passive and active defenses, have emerged to address this challenge. In active defenses, representing the focus of this paper, known signatures, synthesized based on the system's unique characteristics, are inserted into the system. In contradistinction, passive methods rely solely on observing system behavior in search of patterns of normal behavior with deviations thereof representing abnormal behavior. In their most sophisticated implementations, both passive and active defenses rely on the use of data analytics to identify the patterns and synthesize the observed and/or inserted signatures. Past research has shown that passive defenses may be bypassed by APT actors relying on data analytics and their intimate knowledge of the system to evade detection by respecting the patterns identified by the defenders. Here this manuscript explores the use of active defenses under the assumption that the attacker has privileged access to the system, including access to the system's model and sensors data. Specifically, this manuscript assesses the ability of active defenses to remain invisible to the attackers, and discusses the associated challenges that must be addressed to ensure their resiliency against APT actors.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Machine learning control for disruption and tearing mode avoidance

In this work, real-time feedback control based on machine learning algorithms (MLA) was successfully developed and tested on DIII-D plasmas to avoid tearing modes and disruptions while maximizing the plasma performance, which is measured by normalized plasma beta. Control uses MLAs that were trained with ensemble learning methods using only the data available to the real-time Plasma Control System (PCS) from several thousand DIII-D discharges. A `tearability' metric that quantifies the likelihood of the onset of 2/1 tearing modes (TM) in a given time window, and a 'disruptivity' metric that quantifies the likelihood of the onset of plasma disruptions were first tested o -line then implemented on the PCS. A real-time control system based on these MLAs was successfully tested on DIII-D discharges, using feedback algorithms to maximize β Ν while avoiding tearing modes and to dynamically adjust ramp down to avoid high-current disruptions in ramp down.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine learning methods for probabilistic locked-mode predictors in tokamak plasmas

A rotating tokamak plasma can interact resonantly with the external helical magnetic perturbations, also known as error fields. This can lead to locking and then to disruptions. We leverage machine learning (ML) methods to predict the locking events. We use a coupled third-order nonlinear ordinary differential equation model to represent the interaction of the magnetic perturbation and the plasma rotation with the error field. This model is sufficient to describe qualitatively the locking and unlocking bifurcations. Here, we explore using ML algorithms with the simulation data and experimental data, focusing on the methods that can be used with sparse datasets. These methods lead to the possibility of the avoidance of locking in real-time operations. We describe the operational space in terms of two control parameters: the magnitude of the error field and the rotation frequency associated with the momentum source that maintains the plasma rotation. The outcomes are quan- tified by order parameters that completely characterize the state, whether locked or unlocked. We use unsupervised ML methods to classify locked/unlocked states and note the usefulness of a certain normalization of the order parameters. Three supervised ML classifiers are used in suite to estimate the probability of locking in the region of control parameter space with hysteresis, i.e., the set of control parameters for which both locked and unlocked states can exist. The results show that a neural network gives the best estimate of the locking probability. An analogy of the present locking model with the van der Waals equation of state is also provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temporal and spatial laser intensification within nodular defects overcoated with multilayer dielectric mirrors over a wide range of defect geometries

The nodular defect shape and the laser incidence angle have a dramatic impact on the spatial distribution of light intensification within the nodule as well as how the laser light is drained from the defect. Nodular defect geometries unique to ion beam sputtering, ion-assisted deposition, and electron-beam (e-beam) deposition, respectively, are modeled in this parametric study over a wide range of nodular inclusion diameters and layer count for optical interference mirror coatings constructed with quarter-wave thicknesses and capped with a half wave of the low index material. It was found for hafnia ($n$=1.9) and silica ($n$=1.45) multilayer mirrors that the light intensification in nodular defects with a $C$ factor of 8, typical of e-beam deposited coatings deposited with a wide range of deposition angles, was maximized for a 24-layer design. For intermediate size inclusion diameters, increasing the layer count for normal incidence multilayer mirrors reduced the light intensification within the nodular defect. A second parametric study explored the impact of the nodule shape on the light intensification for a fixed number of layers. In this case, there is a strong temporal trend for the different nodule shapes. Narrow nodules tend to drain more laser energy through the bottom of the nodule into the substrate while wide nodules tend to drain more laser energy through the top of the nodule when irradiated at normal incidence. Further, at a 45° incidence angle, waveguiding is an additional method to drain laser energy from the nodular defect. Finally, laser light resonates within nodular defects longer than within the adjacent nondefective multilayer structure.

47 OTHER INSTRUMENTATION↗

XANES spectroscopy of sulfides stable under reducing conditions

X-ray absorption near-edge structure (XANES) spectroscopy is a powerful technique to quantitatively investigate sulfur speciation in geologically complex materials such as minerals, glasses, soils, organic compounds, industrial slags, and extraterrestrial materials. This technique allows non-destructive investigation of the coordination chemistry and oxidation state of sulfur species ranging from sulfide (2–oxidation state) to sulfate (6+ oxidation state). Each sulfur species has a unique spectral shape with a characteristic K-edge representing the s → p and d hybridization photoelectron transitions. As such, sulfur speciation is used to measure the oxidation state of samples by comparing the overall XANES spectra to that of reference compounds. Although many S XANES spectral standards exist for terrestrial applications under oxidized conditions, new sulfide standards are needed to investigate reduced (oxygen fugacity, f O 2 , below IW) silicate systems relevant for studies of extraterrestrial materials and systems. Sulfides found in certain meteorites (e.g., enstatite chondrites and aubrites) and predicted to exist on Mercury, such as CaS (oldhamite), MgS (niningerite), and FeCr 2 S 4 (daubréelite), are stable at f O 2 below IW-3 but rapidly oxidize to sulfate and/or produce sulfurous gases under terrestrial surface conditions. XANES spectra of these compounds collected to date have been of variable quality, possibly due to the unstable nature of certain sulfides under typical (e.g., oxidizing) laboratory conditions. A new set of compounds was prepared for this study and their XANES spectra are analyzed for comparison with potential extraterrestrial analogs. S K -edge XANES spectra were collected at Argonne National Lab for FeS (troilite), MnS (alabandite), CaS (oldhamite), MgS (niningerite), Ni 1–x S, NiS 2 , CaSO 4 (anhydrite), MgSO 4 , FeSO 4 , Fe 2 (SO 4 ) 3 , FeCr 2 S 4 (daubréelite), Na 2 S, Al 2 S 3 , Ni 7 S 6 , and Ni 3 S 2 ; the latter five were analyzed for the first time using XANES. These standards expand upon the existing S XANES end-member libraries at a higher spectral resolution (0.25 eV steps) near the S K -edge. Processed spectra, those that have been normalized and “flattened,” are compared to quantify uncertainties due to data processing methods. Future investigations that require well-characterized sulfide standards, such as the ones presented here, may have important implications for understanding sulfur speciation in reduced silicate glasses and minerals with applications for the early Earth, Moon, Mercury, and enstatite chondrites.

Anzures, Brendan A.↗

Transfer Learning for HVAC System Fault Detection

Faults in HVAC systems degrade thermal comfort and energy efficiency in buildings and have received significant attention from the research community, with data driven methods gaining in popularity. Yet the lack of labeled data, such as normal versus faulty operational status, has slowed the application of machine learning to HVAC systems. In addition, for any particular building, there may be an insufficient number of observed faults over a reasonable amount of time for training. To overcome these challenges, we present a transfer methodology for a novel Bayesian classifier designed to distinguish between normal operations and faulty operations. The key is to train this classifier on a building with a large amount of sensor and fault data (for example, via simulation or standard test data) then transfer the classifier to a new building using a small amount of normal operations data from the new building. We demonstrate a proof-of-concept for transferring a classifier between architecturally similar buildings in different climates and show few samples are required to maintain classification precision and recall.

transfer learning, Building HVAC, Bayesian framewo↗

An improved viscid/inviscid interaction procedure for transonic flow over airfoils

A new interacting boundary layer approach for computing the viscous transonic flow over airfoils is described. The theory includes a complete treatment of viscous interaction effects induced by the wake and accounts for normal pressure gradient effects across the boundary layer near trailing edges. The method is based on systematic expansions of the full Reynolds equation of turbulent flow in the limit of Reynolds numbers, Reynolds infinity. Procedures are developed for incorporating the local trailing edge solution into the numerical solution of the coupled full potential and integral boundary layer equations. Although the theory is strictly applicable to airfoils with cusped or nearly cusped trailing edges and to turbulent boundary layers that remain fully attached to the airfoil surface, the method was successfully applied to more general airfoils and to flows with small separation zones. Comparisons of theoretical solutions with wind tunnel data indicate the present method can accurately predict the section characteristics of airfoils including the absolute levels of drag.

Melnik, R. E.↗

Evidence of an infrared luminosity indicator for galaxies

To elucidate the nature of infrared-luminous galaxies discovered with the IRAS satellite, the optical and infrared luminosities of 1161 Markarian galaxies and 2146 'normal' galaxies from the CfA redshift survey are compared. Survival analysis statistical methods that take upper limits fully into account are used. It is found that L(IR)/L(B) is statistically correlated with L(60) in both samples, though they differ in the distribution at low luminosities. The derived correlation shows that L(IR)/L(B) provides an indicator for L(60). Since galaxies selected in unbiased IRAS surveys will have higher L(IR)/L(B) than optically selected galaxies, they are therefore also selected for high L(60).

Feigelson, Eric D.↗

Mass transport phenomena between bubbles and dissolved gases in liquids under reduced gravity conditions

This paper will describe the experimental and analytical work that has been done to establish justification and feasibility for a Shuttle mid-deck experiment involving mass transfer between a gas bubble and a liquid. The experiment involves the observation and measurement of the dissolution of an isolated, immobile gas bubble of specified size and composition in a thermostatted solvent liquid of known concentration in the reduced gravity environment of earth orbit. Methods to generate and deploy the bubble have been successful both in normal gravity using mutually buoyant fluids and under reduced gravity conditions in the NASA Lear Jet. Initialization of the experiment with a bubble of a prescribed size and composition in a liquid of known concentration has been accomplished using the concept of unstable equilibrium. Subsequent bubble dissolution or growth is obtained by a step increase or decrease in the liquid pressure. A numerical model has been developed which simulates the bubble dynamics and can be used to determine molecular parameters by comparison with the experimental data. The primary objective of the experiment is the elimination of convective effects that occur in normal gravity. The results will yield information on transport under conditions of pure diffusion.

Dewitt, Kenneth J.↗

Probabilities and statistics for backscatter estimates obtained by a scatterometer

Methods for the recovery of winds near the surface of the ocean from measurements of the normalized radar backscattering cross section must recognize and make use of the statistics (i.e., the sampling variability) of the backscatter measurements. Radar backscatter values from a scatterometer are random variables with expected values given by a model. A model relates backscatter to properties of the waves on the ocean, which are in turn generated by the winds in the atmospheric marine boundary layer. The effective wind speed and direction at a known height for a neutrally stratified atmosphere are the values to be recovered from the model. The probability density function for the backscatter values is a normal probability distribution with the notable feature that the variance is a known function of the expected value. The sources of signal variability, the effects of this variability on the wind speed estimation, and criteria for the acceptance or rejection of models are discussed. A modified maximum likelihood method for estimating wind vectors is described. Ways to make corrections for the kinds of errors found for the Seasat SASS model function are described, and applications to a new scatterometer are given.

Pierson, Willard J., Jr.↗