Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “smoothing”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Microstrata

Abstract Microstrata are the non-extremal analogues of superstrata: they are smooth, non-extremal (non-BPS) solitonic solutions to IIB supergravity whose deep-throat limits approximate black holes. Using perturbation theory and numerical methods, we construct families of solutions using a consistent truncation to three-dimensional supergravity. The most general families presented here involve two continuous parameters, or amplitudes, and four quantized parameters that set the angular momenta and energy levels. Our solutions are asymptotic to the vacuum of the D1-D5 system: AdS 3 ×S 3 × 𝕋 4 . Using holography, we show that the they are dual to multi-particle states in the D1-D5 CFT involving a large number of mutually non-BPS supergravitons and we determine the anomalous dimensions of these states from the binding energies in supergravity. These binding energies are uniformly negative and depend non-linearly on the amplitudes of the states. In one family of solutions, smoothness restricts some of the fields to lie on a special locus of the parameter space. Using precision holography we show that this special locus can be identified with the multi-particle states constructed via the standard OPE of the single-particle constituents. Our numerical analysis shows that microstrata are robust at large amplitudes and the solutions can be obtained to very high precision.

Physics↗
DOE OSTI2020

Strings in irrelevant deformations of AdS 3 /CFT 2

We generalize our recent analysis [1] of probe string dynamics to the case of general single-trace TT¯,JT¯ and TJ¯ deformations. We show that in regions in coupling space where the bulk geometry is smooth, the classical trajectories of such strings are smooth and approach the linear dilaton boundary at either the far past or the far future. These trajectories give rise to quantum scattering states with arbitrarily high energies. When the bulk geometry has closed timelike curves (CTC’s), the trajectories are singular for energies above a critical value E c . This singularity occurs in the region with CTC’s, and the value of E c agrees with that read off from the dual boundary theory for all values of the couplings and charges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A proximal trust-region method for nonsmooth optimization with inexact function and gradient evaluations

Many applications require minimizing the sum of smooth and nonsmooth functions. For example, basis pursuit denoising problems in data science require minimizing a measure of data misfit plus an $\ell^1$-regularizer. Similar problems arise in the optimal control of partial differential equations (PDEs) when sparsity of the control is desired. Here, we develop a novel trust-region method to minimize the sum of a smooth nonconvex function and a nonsmooth convex function. Our method is unique in that it permits and systematically controls the use of inexact objective function and derivative evaluations. When using a quadratic Taylor model for the trust-region subproblem, our algorithm is an inexact, matrix-free proximal Newton-type method that permits indefinite Hessians. We prove global convergence of our method in Hilbert space and demonstrate its efficacy on three examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Efficient proximal subproblem solvers for a nonsmooth trust-region method

In [R. J. Baraldi and D. P. Kouri, Mathematical Programming, (2022), pp. 1-40], we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex and nonsmooth convex function. The principle expense of this method is in computing a trial iterate that satisfies the so-called fraction of Cauchy decrease condition—a bound that ensures the trial iterate produces sufficient decrease of the subproblem model. In this paper, we expound on various proximal trust-region subproblem solvers that generalize traditional trust-region methods for smooth unconstrained and convex-constrained problems. We introduce a simplified spectral proximal gradient solver, a truncated nonlinear conjugate gradient solver, and a dogleg method. Finally, we compare algorithm performance on examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Accelerated Sparse Recovery via Gradient Descent with Nonlinear Conjugate Gradient Momentum

This paper applies an idea of adaptive momentum for the nonlinear conjugate gradient to accelerate optimization problems in sparse recovery. Specifically, we consider two types of minimization problems: a (single) differentiable function and the sum of a non-smooth function and a differentiable function. In the first case, we adopt a fixed step size to avoid the traditional line search and establish the convergence analysis of the proposed algorithm for a quadratic problem. This acceleration is further incorporated with an operator splitting technique to deal with the non-smooth function in the second case. As a result, we use the convex ι 1 and the nonconvex ι 1 – ι 2 functionals as two case studies to demonstrate the efficiency of the proposed approaches over traditional methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Nonlinear dynamics, bifurcations, and multi-stability in a vibro-impact system with geometric and multi-segmented freeplay nonlinearities

Freeplay is a common type of piecewise-smooth nonlinearity in dynamical systems, and it can cause discontinuity-induced bifurcations and other behaviors that may bring about undesirable and potentially damaging responses. Prior research has focused on piecewise-smooth systems with two or three distinct regions, but less attention is devoted to systems with more regions (i.e., multi-segmented systems). In this work, numerical analysis is performed on a dynamical system with multi-segmented freeplay, in which there are four stiffness transitions and five distinct regions in the phase space. Here, the effects of the multi-segmented parameters are studied through bifurcation diagram evolution along with induced multi-stable behavior and different bifurcations. These phenomena are interrogated through various tools, such as harmonic balance, basins of attraction, phase planes, and Poincaré section analysis. Results show that among the three multi-segmented parameters, the asymmetry has the strongest effect on the response of the system.

42 ENGINEERING↗
DOE OSTI2026

Memory-efficient nonsmooth dynamic optimization using adaptive randomized compression

Dynamic optimization problems arise in many applications including flow control, full waveform inversion, and medical imaging. These problems are plagued by significant computational challenges. One such challenge — and the focus of this work — is the memory limitation induced by the size of the underlying dynamical system. In particular, the entire dynamic trajectory is required for derivative computation and therefore must be stored or recomputed using, e.g., checkpointing. Although recent work demonstrated the use of adaptive randomized sketching to overcome the memory challenge, that work only applies to smooth unconstrained problems, prohibiting its use for nonsmooth regularized and constrained problems. The inclusion of nonsmooth regularizers and constraints is critical as they often arise in an attempt to preserve certain physical properties or to promote sparsity. To solve these problems, we introduce a trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function that leverages randomized sketching to compress the dynamical system trajectories and adaptively adjust the sketch rank to satisfy a gradient inexactness condition. We prove convergence of this algorithm and demonstrate that it achieves substantial memory reduction on three discretized PDE-constrained optimization applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Local convergence analysis of an inexact trust-region method for nonsmooth optimization

In Baraldi, we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex function and a nonsmooth convex function in Hilbert space—a class of problems that is ubiquitous in data science, learning, optimal control, and inverse problems. Furthermore, this algorithm has demonstrated excellent performance and scalability with problem size. In this paper, we enrich the convergence analysis for this algorithm, proving strong convergence of the iterates with guaranteed rates. In particular, we demonstrate that the trust-region algorithm recovers superlinear, even quadratic, convergence rates when using a second-order Taylor approximation of the smooth objective function term.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Enhancing Surface Finish of Additively Manufactured 316L Stainless Steel with Pulse/Pulse Reverse Electropolishing

To improve the reliable integration of additively manufactured (AM) metal components, there is an urgent need for surface finishing methods capable of processing rough, “as-printed” surfaces. Here, we present a new and scalable electropolishing (EP) technique for use on AM 316L stainless steel (SS) that is capable of achieving a surface roughness equivalent to a machined finish, with arithmetic mean heights (S a ) as low as 1.0 µm. This technique utilizes the sequential application of a new, NaCl/H 2 O pulse/pulse reverse (P/PR) EP step along with a conventional DC (direct current) EP step, where P/PR expands the capability of EP by removing roughness features > 1 µm, which has been previously unattainable with DC EP alone without substantial material removal. This is achieved by independently varying pulse parameters, such as pulse width and height, allowing more control over electrodynamic processes at the part-electrolyte interface. Additionally, P/PR EP resulted in higher smoothing efficiency, where a higher degree of smoothing was achieved with less material removal when compared with conventional DC EP alone. This minimizes unwanted change in geometry from excessive metal dissolution.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mesoscopic-scale complexity in macroscopically-uniform plastic flow of an Al 0.3 CoCrFeNi high-entropy alloy

The concept of a smooth and homogeneous plastic flow of solids is nowadays constantly challenged by various observations of the self-organization of crystal defects on mesoscopic scales pertaining, e.g., to acoustic emission or the evolution of the local strain field. Such investigations would be of particular interest for High-Entropy Alloys (HEAs) characterized by extremely complex microstructures. However, the complexity of their deformation has been only studied in relation to the jerky flow. The present work explores a mesoscale complexity of the macroscopically-smooth plastic flow of an Al 0.3 CoCrFeNi alloy. Here, an inclusive approach involving various experimental methods and complementary nonlinear analyses allowed to reveal transitions between distinct dynamical regimes, which may be assimilated with noises of different colors, in particular the blue noise that is rarely observed in complex systems, and various kinds of reddened behavior. Such findings provide significant new insights into the micro/macro transition in the deformation behavior of HEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Rules for the crystallite size and dislocation density evolution in phases during α-ω transformation in Zr under high-pressure and severe plastic flow

The first in-situ X-ray diffraction (XRD) study of the evolution of the crystallite size and dislocation density in phases during plastic strain-induced phase transformation (PT) is performed utilizing α-ω PT in strongly pre-deformed commercially pure Zr as an important example. Rough diamond anvils (rough-DA) are introduced to intensify all occurring processes during heterogeneous compression of Zr in a diamond anvil cell (DAC). The main rule is found that during α-ω PT the crystallite size and dislocation density in ω-Zr depend solely on the volume fraction of ω-Zr and are independent of pressure, plastic strain tensor, its path, and initial nanostructure. Crystallite size in ω-Zr increases from 10 to 60 nm during the PT, while dislocation density reduces from 1.83×10 15 /m 2 to 0.6×10 15 /m 2 . Rough-DA produce a steady nanostructure in α-Zr before PT with smaller crystallite size and larger dislocation density than smooth-DA, leading to a reduction of the minimum pressure for α-ω PT to a record value 0.67 GPa, 9 times smaller than under hydrostatic loading and 5.1 times lower than the phase equilibrium pressure. In addition to strain, the kinetics of strain-induced PT unexpectedly depends on time. Also, strain-controlled part of kinetics is zero order, in contrast to the first-order kinetics with smooth-DA. The obtained results open a new window for understanding the mutual effects of nanostructure evolution and PT during severe plastic flow in various technological and natural processes. In conclusion, they may bring up economic strategies of producing nanocomposites and single-phase nanostructured materials with optimal properties.

Crystalline size↗
DOE OSTI2025

Microstructure characteristics of LPBF&HIP fabricated graded composite transition joint between ferritic steel and austenitic stainless steel

Graded composite transition joints (GCTJs) offer a promising alternative to conventional dissimilar metal welds (DMWs) by enabling smooth compositional and microstructural transitions. However, GCTJs fabricated solely through additive manufacturing (AM) face challenges such as heat accumulation, complex parameter control, and elemental segregation. In this study, we propose a novel approach that relies on AM to design a spatially graded structure in one alloy and then employs hot isostatic pressing (HIP) as a diffusion bonding method to join it with a second alloy. Here, this method combines the flexibility of AM with the powder net-shaping advantage of HIP. Specifically, a series of closely packed austenitic stainless steel 304 conical structures were printed using laser powder bed fusion (LPBF) and then combined with ferritic steel P91 powder via HIP. By using electron backscatter diffraction (EBSD), electron probe microanalysis (EPMA), and transmission electron microscopy (TEM) techniques, the microstructure characteristics of the GCTJ of 304&P91, especially the interdiffusion zone (IDZ), have been systematically investigated. The microstructure at the interface transitions from austenite-ferrite (A+F) to austenite-martensite-ferrite (A+M+F), and finally to martensite-ferrite (M+F) due to diffusion. Additionally, the diffusion width between 304 and P91 increases with the volume fraction of P91. This unique design also ensures a gradual transition in both hardness and thermal expansion coefficient from 304 to P91, thereby enabling a smooth gradient in functional properties. Overall, this study proposes a novel approach for fabricating GCTJs and contributes to advancing design concepts in the field of dissimilar metal joining.

Additively manufacturing (AM)↗
DOE OSTI2019

Applications of the worldline Monte Carlo formalism in quantum mechanics

Highlights: Universal approach to the numerical computation of the Euclidean path integral. Inspired by recent work in relativistic quantum field theory. Here adapted to non-relativistic quantum mechanics. Worked out for the computation of propagators and ground-state energies. Special smoothing procedure for singular potentials. -- Abstract: In recent years efficient algorithms have been developed for the numerical computation of relativistic single-particle path integrals in quantum field theory. Here, we adapt this “worldline Monte Carlo” approach to the standard problem of the numerical approximation of the non-relativistic path integral, resulting in a formalism whose characteristic feature is the fast, non-recursive generation of an ensemble of trajectories that is independent of the potential, and is thus universally applicable. The numerical implementation discretises the trajectories with respect to their time parametrisation but maintains a continuous spatial domain. In the case of singular potentials, the discretised action gets adapted to the singularity through a “smoothing” procedure. We show for a variety of examples (the harmonic oscillator in various dimensions, the modified Pöschl–Teller potential, delta-function potentials, the Coulomb and Yukawa potentials) that the method allows one to obtain fast and reliable estimates for the Euclidean propagator and use them in a certain time window suitable for extracting the ground state energy. Usually Monte Carlo calculations in quantum mechanics are performed using algorithms that involve a more specific adaptation to the potential, and we briefly compare our method with the most widely used such algorithms, which are of Metropolis type. As an aside, we apply it for studying the classical limit where nearly classical trajectories are expected to dominate in the path integral. We expect the advances made here to be useful also in the relativistic case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Wettability and CHF limits of Accident-Tolerant nuclear fuel cladding materials in light water reactor conditions

Here, we present the results of experimental investigations aimed at evaluating the thermal–hydraulic performance of chromium-coated zircaloy, i.e., one of the most promising accident tolerant fuel (ATF) cladding material for light water nuclear reactors. Precisely, we investigate the wettability and critical heat flux (CHF) limits of chromium-coated and conventional zircaloy surfaces in prototypical reactor conditions. For both surface types, we measure the contact angle in a vapor-saturated atmosphere from ambient pressure to the operating pressure of pressurized water reactors (PWRs), i.e., ~15 MPa. We measure the ambient-pressure steady-state flow boiling CHF with a spatially uniform heat flux. We measure the high-pressure steady-state flow boiling CHF with a cosine shape heat flux (up to 20 MPa) and with a uniform heat flux (up to 15 MPa), also exploring the effect of CRUD deposits on the chromium-coated surface. Our results reveal that the chromium surface and the bare zircaloy surface have similar wettability and both become super-hydrophilic in PWR conditions, and that there is practically no difference in the steady-state CHF limits, both at low-pressure and high-pressure conditions, also when the chromium-coated surface is covered by a CRUD deposit. However, while the chromium-coating does not improve the CHF compared to the bare zircaloy surface, it improves the post-CHF behavior. The chromium coating prevents the reaction between zircaloy and steam, which results in the formation of a brittle zirconium oxide through the surface of the cladding. We also measure the transient CHF under exponentially escalating heat flux inputs of a nano-smooth and a rough surface mimicking a chromium-coated zircaloy cladding. Interestingly, the results of the transient heat flux tests suggest that the CHF limit for very short periods (i.e., fast transients) is independent of the surface finish, being the same for a rough chromium surface or a nano-smooth surface.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Enhanced pool boiling heat transfer with metal foam tubes in inline tube bundle configuration

Flooded evaporators, commonly known as shell-and-tube heat exchangers, are widely used in large-scale industrial refrigeration and air-conditioning systems due to their high pool boiling heat transfer efficiency. However, their bulky configuration requires a large refrigerant inventory, which poses safety challenges—particularly when employing next-generation A2L refrigerants that are mildly flammable. Here, to address this issue, the development of compact heat exchangers with enhanced heat transfer performance is essential. In this study, a novel metal foam tube design is proposed to augment pool boiling heat transfer. Experiments were conducted at a saturation temperature of 20 °C and across a heat flux range of 7–60 kW m −2 . The tests were conducted for both smooth and metal foam tubes with an inline tube bundle configuration and for a pitch–to–diameter (P/D) ratios of 1.3 and 1.5. Pool boiling performance was evaluated for next-generation hydrofluoroolefin (HFO) refrigerants, R1234yf and R1234ze(E), with R134a serving as the baseline fluid. The results indicate that the heat transfer coefficient (HTC) of R1234yf is comparable to R134a, while R1234ze(E) shows slightly lower performance. Importantly, the metal foam tubes achieved up to a 217% enhancement in HTC compared to smooth tubes.

Metal foam tubes↗
DOE OSTI2020

Exploring driving forces for length growth in graphene nanoribbons during chemical vapor deposition of hydrocarbons on Ge(0 0 1) via kinetic Monte Carlo simulations

Graphene grown slowly on Ge(0 0 1) using chemical vapor deposition of hydrocarbons leads to high-aspect ratio graphene nanoribbons with smooth edges and a technologically relevant band gap at room temperature; however, the driving forces leading to one-dimensional growth of such graphene crystals are not well understood. Here, we combine a lattice kinetic Monte Carlo approach based on steps in graphene growth and experimental measurements to study the growth of graphene nanoribbons via chemical vapor deposition on Ge(0 0 1). To identify potential reasons for growth of graphene as anisotropic ribbons, we study the impact of anisotropy in various growth parameters on the resulting graphene crystals. Comparing our model with experimental measurements indicates that anisotropy in the stabilization of a graphene precursor species bound to the graphene edge is the most likely reason why high aspect ratio graphene ribbons with smooth edges grow on Ge(0 0 1). Using the growth model developed here, we reproduce experimental trends in the synthesis of graphene nanoribbons on Ge(001) and arrive at an intuitive picture for their growth. These insights shed light on the driving forces governing this highly anisotropic regime of crystal growth.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tuning the fluid wetting dynamics on gold microstructures using photoactive compounds

This work demonstrates how photoactive compounds can be used to tune the surface chemistry, surface free energy, and the wetting velocity of fluids on Spiropyran functionalized surfaces with different surface microstructures. Evidence of this photowetting effect is based on data showing: (1) the cyclic changes in the static, advancing, and receding contact angles $\theta$ CA for multiple UV$\rightleftarrows$vis photoswitching cycles with both smooth and microstructured surfaces and (2) the changes in the fluid wicking velocity (via UV$\rightleftarrows$vis photoswitching) on different Spiropyran functionalized hemiwicking surfaces. X-ray photoelectron spectroscopy is used to determine the efficiency of photoswitching caused by the quality of the Spiropyran functionalization. Reversible photoisomerization of Spiropyran generates a gradient in surface free energy, resulting in a similar change in contact angle ($\theta$ CA ) in the photoswitchable surfaces. In conclusion, by incorporating these changes into the Owens-Wendt and Van Oss surface energy models, the energy of smooth Au surfaces can also be predicted for both UV and visible light irradiations, where reversible contact angle variations of $\Delta$$\theta$ CA ≈ 5–10° are achieved with both water and water-ethanol mixtures due to the conversion of SP to merocyanine, which corresponds to a total free energy change of ~13% on the photoswitchable surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗