Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gradient information”

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 397 records · Page 22

Thermodynamic Consistent Neural Networks for Learning Material Interfacial Mechanics

For multilayer materials in thin substrate systems, interfacial failure is one of the most challenges. The traction-separation relations (TSR) quantitatively describe the mechanical behavior of a material interface undergoing openings, which is critical to understand and predict interfacial failures under complex loadings. However, existing theoretical models have limitations on enough complexity and flexibility to well learn the real-world TSR from experimental observations. A neural network can fit well along with the loading paths but often fails to obey the laws of physics, due to a lack of experimental data and understanding of the hidden physical mechanism. In this paper, we propose a thermodynamic consistent neural network (TCNN) approach to build a data-driven model of the TSR with sparse experimental data. The TCNN leverages recent advances in physics-informed neural networks (PINN) that encode prior physical information into the loss function and efficiently train the neural networks using automatic differentiation. We investigate three thermodynamic consistent principles, i.e., positive energy dissipation, steepest energy dissipation gradient, and energy conservative loading path. All of them are mathematically formulated and embedded into a neural network model with a novel defined loss function. A real-world experiment demonstrates the superior performance of TCNN, and we find that TCNN provides an accurate prediction of the whole TSR surface and significantly reduces the violated prediction against the laws of physics.

Zhang, Jiaxin↗

Leveraging gradient weighted class activation mapping to improve classification effectiveness: Case study in transportation infrastructure characterization

Roadway “corners†are common for pedestrian use, whether designated with markings or not. Different types of markings have been deployed, ranging from simple parallel lines to more complex designs. Understanding the impact of different types of crosswalks is important for public safety. In this work we explore methods to improve the logging of marked crosswalk types. We used the Roadway Information Database from the Second Strategic Highway Research Project and used active learning methods with transfer learning to identify the crosswalk types (marked or unmarked). Upon completion we found our classifiers were unable to perform above roughly 94% correct classifications. To improve their efficacy, we separated the crosswalks into their “fine grained†types and used Gradient-Weighted Class Activation Mapping to isolate and study the features that classified the crosswalks. We compared this with sampled manually marked crosswalks and present findings. We believe this use case can represent a process to improve the active learning method for some visual machine learning applications.

Karnowski, Thomas↗

Finding simplicity: unsupervised discovery of features, patterns, and order parameters via shift-invariant variational autoencoders *

Abstract Recent advances in scanning tunneling and transmission electron microscopies (STM and STEM) have allowed routine generation of large volumes of imaging data containing information on the structure and functionality of materials. The experimental data sets contain signatures of long-range phenomena such as physical order parameter fields, polarization, and strain gradients in STEM, or standing electronic waves and carrier-mediated exchange interactions in STM, all superimposed onto scanning system distortions and gradual changes of contrast due to drift and/or mis-tilt effects. Correspondingly, while the human eye can readily identify certain patterns in the images such as lattice periodicities, repeating structural elements, or microstructures, their automatic extraction and classification are highly non-trivial and universal pathways to accomplish such analyses are absent. We pose that the most distinctive elements of the patterns observed in STM and (S)TEM images are similarity and (almost-) periodicity, behaviors stemming directly from the parsimony of elementary atomic structures, superimposed on the gradual changes reflective of order parameter distributions. However, the discovery of these elements via global Fourier methods is non-trivial due to variability and lack of ideal discrete translation symmetry. To address this problem, we explore the shift-invariant variational autoencoders (shift-VAEs) that allow disentangling characteristic repeating features in the images, their variations, and shifts that inevitably occur when randomly sampling the image space. Shift-VAEs balance the uncertainty in the position of the object of interest with the uncertainty in shape reconstruction. This approach is illustrated for model 1D data, and further extended to synthetic and experimental STM and STEM 2D data. We further introduce an approach for training shift-VAEs that allows finding the latent variables that comport to known physical behavior. In this specific case, the condition is that the latent variable maps should be smooth on the length scale of the atomic lattice (as expected for physical order parameters), but other conditions can be imposed. The opportunities and limitations of the shift VAE analysis for pattern discovery are elucidated.

97 MATHEMATICS AND COMPUTING↗

Date Release Report for Large Surface Explosion Coupling Experiment (LSECE) Nevada National Security Site

The DTRA sponsored 2020-2022 Large Surface Explosion Coupling Experiment (LSECE) consists of two large ground surface chemical explosions, data collection, analysis and modeling carried out in 2020-2022. The LSECE chemical explosions were carried out at the site of prior NNSA sponsored buried chemical explosions that were part of the Source Physics Experiment (SPE) Phase II in Dry Alluvium Geology (DAG). This report described the data collected under LSECE that is being made publicly available. The prior buried explosion DAG data are described in a separate report (https://www.osti.gov/biblio/1825534) and are also publicly available. The LSECE explosions, data analysis and numerical modeling work sponsored by DTRA were intended to address three objectives: 1) Generate seismo-acoustic data to test and improve numerical models of explosion energy coupling in dry alluvium geology; 2) To improve seismo-acoustic yield estimation techniques as a function of depth and medium properties; 3) To study and improve acoustic propagation modeling under two different atmospheric conditions. The LSECE data that were collected have expanded the prior DAG buried explosion dataset to include surface chemical explosions recorded on a common set of seismo-acoustic stations. They have allowed a more detailed study of above and below surface seismo-acoustic energy coupling with applications to explosion monitoring and assessment. The two approximately 1-ton TNT equivalent yield LSECE explosions were detonated at different times of the day to explore the effects of the different atmospheric conditions. The first chemical explosion “Artemis” was conducted before dawn when temperature inversions were present. The second chemical explosion “Apollo” was conducted on a sunny afternoon when the temperature gradient was more linear. The LSECE chemical explosion data were also collected across a variety of different sensor types to allow evaluation of their effectiveness in recovering useful information. The LSECE instrumentation included fiber optic or Distributed Acoustic Sensing (DAS), a dense array or Large-N array of seismometers, borehole accelerometers, DAS and a velocity meter, airborne acoustic instruments and a variety of visual and remote sensing data.

58 GEOSCIENCES↗

Harvesting Reactor Pressure Vessel Beltline Material from the Decommissioned Zion Nuclear Power Plant Unit 1

The decommissioning of the Zion Nuclear Power Plant (NPP) provided a unique opportunity to harvest and study service-aged reactor pressure vessel (RPV) beltline materials. This work, conducted through the U.S. Department of Energy’s Light Water Reactor Sustainability (LWRS) Program, aims to improve the understanding of radiation-induced embrittlement to support extended nuclear plant operations. Material segments containing the Linde 80 flux, wire heat 72105 (WF-70) beltline weld and the A533B Heat B7835-1 base metal, obtained from the intermediate shell region with a peak fluence of 0.7 × 10 19 n/cm 2 (E > 1.0 MeV), were extracted, cut into blocks, and machined into test specimens for mechanical and microstructural characterization. The segmentation process involved oxy-propane torch-cutting, followed by precision machining using wire saws and electrical discharge machining (EDM). A chemical composition analysis confirmed the expected variations in alloying elements, with copper levels being notably higher in the weld metal. The harvested specimens enable a detailed evaluation of through-wall embrittlement gradients, a comparison with the existing surveillance data, and the validation of predictive embrittlement models. This study provides critical data for assessing long-term reactor vessel integrity, informing aging-management strategies, and supporting regulatory decisions to extend the life of nuclear plants. This article is a revised and expanded version of a paper entitled, “Current Status of the Characterization of RPV Materials Harvested from the Decommissioned Zion Unit 1 Nuclear Power Plant”, PVP2017-65090, which was accepted and presented at the ASME 2017 Pressure Vessels and Piping Conference, Waikoloa, HI, USA, 16–20 July 2017.

harvesting beltline material↗

Computational aspects of maximum likelihood estimation and reduction in sensitivity function calculations

This paper discusses numerical aspects of computing maximum likelihood estimates for linear dynamical systems in state-vector form. Different gradient-based nonlinear programming methods are discussed in a unified framework and their applicability to maximum likelihood estimation is examined. The problems due to singular Hessian or singular information matrix that are common in practice are discussed in detail and methods for their solution are proposed. New results on the calculation of state sensitivity functions via reduced order models are given. Several methods for speeding convergence and reducing computation time are also discussed.

Gupta, N. K.↗

Properties of echo spectra observed by MST radars

Turbulent scatter and Fresnel reflection are the fundamental echoing mechanisms to interpret the signals observed by Mesosphere-stratosphere-troposphere (MST) radars. Turbulent scattered echoes provide information about the turbulence structure and mean flow of the atmosphere. Observational results with VHF MST radars, however, show the importance of Fresnel reflection due to the infinite gradient of reflectivity at the edges of a scattering layer. This condition is excluded for the weak fluctuation models but it is still possible to include the observed aspect sensitivity by assuming an anisotropic structure of fluctuations. Another explanation of the aspect sensitivity observed by MST radars is advanced. Spectral estimates by the widely used periodogram were related to a four-dimensional spectrum of atmospheric fluctuations with anisotropic structure. Effects of the radar system such as antenna beam width, beam direction and Fast Fourier Transformations (FFT) data length were discussed for the anisotropic turbulent atmosphere. Echo parameters were also estimated.

Wakasugi, K.↗

Activities of the Solid State Physics Research Institute

Three research programs are reviewed. These programs are muon spin rotation, studies of annealing in gallium arsenide and Hall effect studies in semiconductors. The muon spin rotation work centers around the development of a facility at the Alternating Gradient Synchrotron of BNL. Studies of annealing in GaAs concerns itself with the measurement of depolarization in GaAs. The Hall effect studies of proton damaged semiconductors provide new information on the nature of defects and dislocations in GaAs.

Source record↗

The Vertical Error Characteristics of GOES-derived Winds: Description and Impact on Numerical Weather Prediction

Errors in the height assignment of some satellite-derived winds exist because the satellites sense radiation emitted from a finite layer of the atmosphere rather than a specific level. Potential problems in data assimilation may arise because the motion of a measured layer is often represented by a single-level value. In this research, cloud and water vapor motion winds that are derived from the Geostationary Operational Environmental Satellites (GOES winds) are compared to collocated rawinsonde observations (RAOBs). An important aspect of this work is that in addition to comparisons at each assigned height, the GOES winds are compared to the entire profile of the collocated RAOB data to determine the vertical error characteristics of the GOES winds. The impact of these results on numerical weather prediction is then investigated. The comparisons at individual vector height assignments indicate that the error of the GOES winds range from approx. 3 to 10 m/s and generally increase with height. However, if taken as a percentage of the total wind speed, accuracy is better at upper levels. As expected, comparisons with the entire profile of the collocated RAOBs indicate that clear-air water vapor winds represent deeper layers than do either infrared or water vapor cloud-tracked winds. This is because in cloud-free regions the signal from water vapor features may result from emittance over a thicker layer. To further investigate characteristics of the clear-air water vapor winds, they are stratified into two categories that are dependent on the depth of the layer represented by the vector. It is found that if the vertical gradient of moisture is smooth and uniform from near the height assignment upwards, the clear-air water vapor wind tends to represent a relatively deep layer. The information from the comparisons is then used in numerical model simulations of two separate events to determine the forecast impacts. Four simulations are performed for each case: 1) A control simulation that assimilates no satellite wind data, 2) assimilation of all GOES winds according to their assigned single level height, 3) assimilation of all GOES winds spread over multiple levels, and 4) assimilation of all GOES winds spread over multiple levels, but with variations in the vertical influence of clear-air water vapor winds based on the moisture profile in the model. In the first case, a strong mid-latitude cyclone is present and the use of the satellite data results in improved storm tracks during the initial approx. 36 h forecast period. This is because the satellite data improves the analysis of the environment into which the storm progresses. Statistics for mean wind vector and height differences show that, with the exception of the height field at later times in the first case, the use of GOES winds improves the simulation with time. The simulation results suggest that it is beneficial to spread the GOES wind information over multiple levels, particularly when the moisture profile is used to define the vertical influence.

Rao, P. Anil↗

Valley Network Morphology and Topographic Gradients on Mars

Data returned from the Mars Orbiter Laser Altimeter allows construction of a high precision digital elevation model. Quantitative investigations into the geomorphic properties of drainage features, similar to ones carried out on Earth, are now possible Additional information is contained in the original extended abstract.

Aharonson, Oded↗

Physics-Driven Convolutional Autoencoder Approach for CFD Data Compressions: Preprint

With the growing size and complexity of turbulent flow models, data compression approaches are of the utmost importance to analyze, visualize, or restart the simulations. Recently, in-situ autoencoder-based compression approaches have been proposed and shown to be effective at producing reduced representations of turbulent flow data. However, these approaches focus solely on training the model using point-wise sample reconstruction losses that do not take advantage of the physical properties of turbulent flows. In this paper, we show that training autoencoders with additional physics-informed regularizations, e.g., enforcing incompressibility and preserving enstrophy, improves the compression model in three ways: (i) the compressed data better conform to known physics for homogeneous isotropic turbulence without negatively impacting point-wise reconstruction quality, (ii) inspection of the gradients of the trained model uncovers changes to the learned compression mapping that can facilitate the use of explainability techniques, and (iii) as a performance byproduct, training losses are shown to converge up to 12x faster than the baseline model.

auto-encoders↗

Gaining Perspective on Unconventional Well Design Choices through Play-level Application of Machine Learning Modeling

The recent development of unconventional oil and gas (O&G) reservoirs has led to an abundant hydrocarbon supply, both domestically and globally. However, there is a continued push to develop new and innovative approaches to improve exploration and extraction efficiencies and overall well productivity moving forward. Substantial improvements in unconventional O&G development are expected through optimized well completion and stimulation strategies aimed at maximizing well productivity. Optimizing well designs will require tailoring to the distinctive geologic conditions present for any newly placed well. To better evaluate the impact of well design attributes and their associated interactions on productivity in a major unconventional play, multivariate machine learning-based models that use empirical datasets were developed. A gradient boosted regression tree (GBRT) algorithm was applied. GBRT has been narrowly investigated for O&G applications but enables straightforward parametric importance and influence evaluation, as well as assessment of parameter interaction effects. Models were trained on well design and locational parameters that serve as a proxy for variable geologic conditions to estimate two types of productivity indicator response variables strongly correlated to estimated ultimate recovery (EUR). The dataset utilized consists of over 7,000 well observations that cover the majority of the productive region of the Marcellus Shale. Model performance was evaluated and algorithm parameters tuned by analyzing the goodness-of-fit for simulated results against observed data in a cross-validation approach. Models were found capable of 73–79 percent prediction accuracy on held out testing data of gas equivalent production and can be used to inform future well design and placement decisions for increasing EUR per well and improving overall field-level recovery. Study results indicate that Marcellus well performance improves most with upscaling perforated interval lengths and water and proppant volumes per foot; but relative productivity improvements are spatially dependent across the play. Finally, optimal combinations of water and proppant on well performance were found to vary depending on well location, emphasizing the utility of data-driven models capable of broad application across a play of interest for informing tailored well design approaches prior to their field deployment.

04 OIL SHALES AND TAR SANDS↗

In situ synchrotron investigation of degenerate graphite nodule evolution in ductile cast iron

Ductile cast irons (DCIs) are of increasing importance in the renewable energy and transportation sectors. The distribution and morphology of the graphite nodules, in particular the formation of degenerate features during solidification, dictate the mechanical performance of DCIs. In situ high-speed synchrotron X-ray tomography was used to capture the evolution of graphite nodules during solidification of DCI, including degenerate features and the effect of the carbon concentration field. The degeneration of nodules is observed to increase with re-melting cycles, which is attributed to Mg-loss. The dendritic primary austenite and carbon concentration gradients in the surrounding liquid phase were found to control nodule morphology by locally restricting and promoting growth. A coupled diffusion-mechanical model was developed, confirming the experimentally informed hypothesis that protrusions form through liquation cracking of the austenite shell and subsequent localised growth. In conclusion, these results provide valuable insights into the solidification kinetics of cast irons, supporting the design of advanced alloys.

36 MATERIALS SCIENCE↗

An autoencoder compression approach for accelerating large-scale inverse problems

Abstract Partial differential equation (PDE)-constrained inverse problems are some of the most challenging and computationally demanding problems in computational science today. Fine meshes required to accurately compute the PDE solution introduce an enormous number of parameters and require large-scale computing resources such as more processors and more memory to solve such systems in a reasonable time. For inverse problems constrained by time-dependent PDEs, the adjoint method often employed to compute gradients and higher order derivatives efficiently requires solving a time-reversed, so-called adjoint PDE that depends on the forward PDE solution at each timestep. This necessitates the storage of a high-dimensional forward solution vector at every timestep. Such a procedure quickly exhausts the available memory resources. Several approaches that trade additional computation for reduced memory footprint have been proposed to mitigate the memory bottleneck, including checkpointing and compression strategies. In this work, we propose a close-to-ideal scalable compression approach using autoencoders to eliminate the need for checkpointing and substantial memory storage, thereby reducing the time-to-solution and memory requirements. We compare our approach with checkpointing and an off-the-shelf compression approach on an earth-scale ill-posed seismic inverse problem. The results verify the expected close-to-ideal speedup for the gradient and Hessian-vector product using the proposed autoencoder compression approach. To highlight the usefulness of the proposed approach, we combine the autoencoder compression with the data-informed active subspace (DIAS) prior showing how the DIAS method can be affordably extended to large-scale problems without the need for checkpointing and large memory.

Mathematics↗

The Circular Velocity Curve of the Milky Way from 5–25 kpc Using Luminous Red Giant Branch Stars

We present a sample of 254,882 luminous red giant branch (LRGB) stars selected from the APOGEE and LAMOST surveys. By combining photometric and astrometric information from the Two Micron All Sky Survey and Gaia survey, the precise distances of the sample stars are determined by a supervised machine-learning algorithm: the gradient-boosted decision trees. To test the accuracy of the derived distances, member stars of globular clusters (GCs) and open clusters are used. The tests by cluster member stars show a precision of about 10% with negligible zero-point offsets, for the derived distances of our sample stars. The final sample covers a large volume of the Galactic disk(s) and halo of 0 < R < 30 kpc and |Z| ≤ 15 kpc. The rotation curve (RC) of the Milky Way across the radius of 5 ≲ R ≲ 25 kpc has been accurately measured with ~54,000 stars of the thin disk population selected from the LRGB sample. The derived RC shows a weak decline along R with a gradient of -1.83 ± 0.02 (stat.) ± 0.07 (sys.) km s -1 kpc -1 , in excellent agreement with the results measured by previous studies. The circular velocity at the solar position, yielded by our RC is 234.04 ± 0.08 (stat.) ± 1.36 (sys.) km s -1 , again in great consistency with other independent determinations. From the newly constructed RC, as well as constraints from other data, we have constructed a mass model for our Galaxy, yielding a mass of the dark matter halo of M 200 = (8.05 ± 1.15) × 10 11 M ⊙ with a corresponding radius of R 200 = 192.37 ± 9.24 kpc and a local dark matter density of 0.39 ± 0.03 GeV cm -3 .

79 ASTRONOMY AND ASTROPHYSICS↗

Oceanographic applications of color-enhanced satellite imageries.

Black and white infrared imageries obtained from satellites over the oceans were transformed into color presentations. Investigations in different regions (Persian Gulf, Arabian Coast, Somali Coast and the Northwest Coast of Australia) revealed that temperature gradients and temperature differences of two degrees Celsius can be displayed by the color process from the imageries. This data display can be used for a rapid analysis of information obtained with an APT station.

Szekielda, K.-H.↗

Cold Flow Determination of the Internal Flow Environment Around the Submerged TVC Nozzle for the Space Shuttle SRM

A series of subscale cold flow tests was performed to quantify the gas flow characteristics at the aft end of the Space Shuttle Solid Rocket Motor. This information was used to support the analyses of the redesigned nozzle/case joint. A portion of the thermal loads at the joint are due to the circumferential velocities and pressure gradients caused primarily by the gimbaling of the submerged nose TVC nozzle. When the nozzle centerline is vectored with respect to the motor centerline, asymmetries are set up in the flow field under the submerged nozzle and immediately adjacent to the nozzle/case joint. Specific program objectives included: determination of the effects of nozzle gimbal angle and propellant geometry on the circumferential flow field; measurement of the static pressure and gas velocities in the vicinity of the nozzle/case joint; use of scaling laws to apply the subscale cold flow data to the full scale SRM; and generation of data for use in validation of 3-D computational fluid dynamic, CFD, models of the SRM flow field. These tests were conducted in the NASA Marshall Space Flight Center Airflow Facility with a 7.5 percent scale model of the aft segment of the SRM. Static and dynamic pressures were measured in the model to quantify the flow field. Oil flow data was also acquired to obtain qualitative visual descriptions of the flow field. Nozzle gimbal angles of 0, 3.5, and 7 deg were used with propellant grain configurations corresponding to motor burn times of 0, 9, 19, and 114 seconds. This experimental program was successful in generating velocity and pressure gradient data for the flow field around the submerged nose nozzle of the Space Shuttle SRM at various burn times and gimbal angles. The nature of the flow field adjacent to the nozzle/case joint was determined with oil droplet streaks, and the velocity and pressure gradients were quantified with pitot probes and wall static pressure measurements. The data was applied to the full scale SRM thru a scaling analysis and the results compared well with the 3-D computational fluid dynamics computer model.

Whitesides, R. H.↗

A carbon isotope mass balance for an anoxic marine sediment: Isotopic signatures of diagenesis

A carbon isotope mass balance was determined for the sediments of Cape Lookout Bight, NC to constrain the carbon budgets published previously. The diffusive, ebullitive and burial fluxes of sigma CO2 and CH4, as well as the carbon isotope signatures of these fluxes, were measured. The flux-weighted isotopic signature of the remineralized carbon (-18.9 plus or minus 2.7 per mil) agreed with the isotopic composition of the remineralized organic carbon determined from the particulate organic carbon (POC) delta(C-13) profiles (-19.2 plus or minus 0.2), verifying the flux and isotopic signature estimates. The measured delta(C-13) values of the sigma CO2 and CH4 diffusive fluxes were significantly different from those calculated from porewater gradients. The differences appear to be influenced by methane oxidation at the sediment-water interface, although other potential processes cannot be excluded. The isotope mass balance provides important information concerning the locations of potential diagenetic isotope effects. Specifically, the absence of downcore change in the delta(C-13) value of the POC fraction and the identical isotopic composition of the POC and the products of remineralization indicate that no isotopic fractionation is expressed during the initial breakdown of the POC, despite its isotopically heterogeneous composition.

Boehme, Susan E.↗