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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A hybrid adaptive multiresolution approach for the efficient simulation of reactive flows

Computational studies that use block-structured adaptive mesh refinement (AMR) approaches suffer from unnecessarily high mesh resolution in regions adjacent to important solution features. This deficiency limits the performance of AMR codes. In this work a novel hybrid adaptive multiresolution (HAMR) approach to AMR-based calculations is introduced to address this issue. The multiresolution (MR) smoothness indicators are used to identify regions of smoothness on the mesh where the computational cost of individual physics solvers may be decreased by replacing direct calculations with interpolation. We suggest an approach to balance the errors due to the adaptive discretization and the interpolation of physics quantities such that the overall accuracy of the HAMR solution is consistent with that of the MR-driven AMR solution. The performance of the HAMR scheme is evaluated for a range of test problems, from pure hydrodynamics to turbulent combustion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Extremized nonlinear and linearized responses in soft metamaterials enabled by gradient-based design and grayscale digital light processing

In this study, we develop a gradient-based design approach that exploits grayscale digital light processing (DLP) 3D printing for extremizing the nonlinear and linearized response of soft metamaterials — materials that harness engineered geometric instabilities to undergo large and programmable changes in configuration. Grayscale DLP approaches modulate local mechanical properties at the pixel scale by tuning the light intensity within a single grayscale image, unlocking an exceptionally large design space. To effectively navigate this space, we develop smooth mappings between local light intensity values and global quantities of interest that characterize the behavior of soft metamaterials. Enabling these smooth mappings are robust and differentiable nonlinear finite element simulations powered by a trust region solver. A PDE-constrained optimization problem is then solved to invert these mappings and produce light intensity distributions that endow the printed part with varying stiffness and flexibility in distinctive regions. It is shown that optimizing the distribution of soft and stiff phases throughout a metamaterial structure results in markedly different buckling and self-contact configurations to drive extremized nonlinear compression and linearized vibration responses. Optimized light intensity distributions are translated to grayscale images and directly used to print soft metamaterial samples, showing remarkable agreement between the buckling and self-contact response in simulated and measured deformed configurations.

Additive manufacturing↗
DOE OSTI2019

Experimental and theoretical investigations on the adherent behaviors of high viscosity liquid: The effect of surface topography

Wetting test and classical molecular dynamics (MD) simulations were employed to investigate the adherent behaviors of hydroxyl-terminated polybutadiene (HTPB), a typical high viscosity liquid widely used as a binder in rocket motor propellants, onto the copper (Cu) surfaces with different topographies. Experimental data show that the microstructures protruding above the substrate surfaces can be ‘sensed’ by the HTPB droplet, which has been indicated by the smaller contact angles and higher interface tension for this specific liquid on the rough surfaces (featured with sinusoidal structures or square pillars) than that on the smooth surface. Simulation results reveal the atomic details of spreading process of HTPB chains, with diverse diffusion propensities presented on different surfaces. The HTPB droplet spreads widely on the smooth surface, exhibiting a roughly isotropic trend. On the contrary, the state of this high viscosity liquid on the sinusoidal-structured surface is featured by wide contact diffusion along the run of groove, but relatively restricted extension in the direction perpendicular to the groove, as a result of the barriers exerted by the nearest neighbouring peaks. The existence of nanoscale pillars increases the exterior surface area of Cu substrate, therefore the HTPB chains can penetrate deep to envelop the peripheries of the pillars, forming abundant direct contacts with Cu surface and producing a highest HTPB–Cu interaction energy, which well accords with the experimental results about the adherent strength of HTPB.

Wu, Chunya↗
DOE OSTI2017

Formation of an ultracarbonaceous Antarctic micrometeorite through minimal aqueous alteration in a small porous icy body

A comprehensive study of the organic chemistry and mineralogy of an ultracarbonaceous micrometeorite (UCAMM D05IB80) collected from near the Dome Fuji Station, Antarctica, was carried out to understand the genetic relationship among organic materials, silicates, and water. The micrometeorite is composed of a dense aggregate of ~5 µm-sized hollow ellipsoidal organic material containing submicrometer-sized phases such as glass with embedded metal and sulfides (GEMS) and mineral grains. There is a wide area of organic material (~15 × 15 μm) in its interior. Low-Ca pyroxene is much more abundant than olivine and shows various Mg/(Mg + Fe) ratios ranging from ~1.0 to 0.78, which is common to previous works on UCAMMs. By contrast, GEMS grains in this UCAMM have unusual chemical compositions. They are depleted in both Mg and S, which suggests that these elements were leached out from the GEMS grains during very weak aqueous alteration, without the formation of phyllosilicates. The organic materials have two textures—smooth and globular with an irregular outline—and these are composed of imine, nitrile and/or aromatic nitrogen heterocycles, and amide. The ratio of nitrogen to carbon (N/C) in the smooth region of the organics is ~0.15, which is five times higher than that of insoluble organic macromolecules in types 1 and 2 carbonaceous chondritic meteorites. In addition, the UCAMM organic materials are soluble in epoxy and are thus hydrophilic; this polar nature indicates that they are very primitive. The surface of the material is coated with an inorganic layer, a few nanometers thick, that consists of C, O, Si, S, and Fe. Sulfur is also contained in the interior, implying the presence of organosulfur moieties. There are no isotopic anomalies of D, 13 C, or 15 N in the organic material. Interstellar photochemistry alone would not be sufficient to explain the N/C ratio of the UCAMM organics; therefore, we imply that a very small amount of fluid on a comet must have been necessary for the formation of the UCAMM. The GEMS grains depleted in Mg and S in the UCAMM prove a very weak degree of aqueous alteration; weaker than that of carbonaceous chondrites. Short-duration weak alteration probably caused by planetesimal shock locally melted cometary ice grains and released water that dissolved the organics; the fluid would likely have not mobilized because of the very low thermal conductivity of the porous icy body. This event enabled the formation of the large organic puddle of the UCAMM, as well as organic matter sulfurization, formation of thin membrane-like layers of minerals, and deformation of organic nanoglobules.

58 GEOSCIENCES↗
DOE OSTI2019

Effect of helical micro-fins on two-phase flow behavior of R410A evaporating in horizontal round tubes obtained through visualization

This paper presents a new approach in visualization of the evaporation mechanism and the two phase flow. The visualization test section is a transparent micro-fin tube made by 3D printing. The 3D printing technology helps reproduce the real geometry in commercial metal micro-fin tubes. Even though the wall material is different from the real metal, this new method provides a more complete understanding of the mechanisms of flow boiling and the heat transfer enhancement due to micro-fin geometry. To investigate the effect of micro-fin geometries on two-phase flow behavior, R410A flow boiling experiments are conducted at 1°C saturation temperature and visualized through a clear smooth tube and micro-fin tubes of 0°, 10° and 18° helix angles. The results showed that the annular flow pattern in the helical micro-fin tube occurred at a lower vapor quality than the smooth and axial micro-fin tubes. Helix angle affects the transition boundary of stratified wavy flow to wavy annular flow, which shifts downward (to lower mass flux and vapor quality conditions) as the helix angle increases. Other transitions are not significantly affected by the helix angle. Helical micro-fin tubes with different groove directions were made to understand the flow behavior on the upward side and downward side of the tube.

3D printing↗
DOE OSTI2023

Enhanced internal condensation of R1233zd(E) on micro- and nanostructured copper and aluminum surfaces

In-tube condensation of refrigerants is an important process which affects thermal efficiency in many applications, ranging from refrigeration and air conditioning to electronics thermal management. In-tube heat transfer and pressure drop are important to heat exchanger sizing and design. Here, in this work, micro- and nanostructured surfaces are applied to the internal wetted areas of copper and aluminum mini-channels to enhance the condensation heat transfer coefficient of hydrofluorocarbon R1233zd(E) refrigerant. To achieve scalable nanomanufacturing, surfaces were uniformly structured by relying on hydrochloric acid etching of aluminum and chemical oxidation of copper. The etched aluminum surfaces exhibited a 150% increase in heat transfer coefficient compared to smooth aluminum channels at specific qualities, with a 66% heat transfer coefficient improvement for complete phase change from saturated vapor to saturated liquid. Copper oxide structures showed no discernable difference in thermal-hydraulic performance when compared to smooth copper channels. Critical dimensionless parameters governing the heat transfer enhancement were identified by varying the tube internal diameter (2.3 mm to 4.7 mm), refrigerant mass flux (50 to 300 kg/(m 2 ·s)), and refrigerant quality (0 to 1). The dimensionless parameters include the Bond number normalized to the condensate film thickness, and the Weber number modified by the vapor friction factor. The relatively small increase in pressure drop (< 10%) associated with these surface enhancements further supports the promise of this method. The scalable and cost-effective techniques used to create these aluminum microstructures may reduce manufacturing cost when compared with current enhancement approaches such as extrusion, drawing, and welding.

42 ENGINEERING↗
DOE OSTI2018

Glass-particle adhesion-force-distribution on clean (laboratory) and contaminated (outdoor) surfaces

This article explores whether laboratory-clean samples are representative of real outdoor surfaces – often covered in dust and other contaminants – with respect to adhesion properties, which ultimately affect particle resuspension. In this study, Atomic Force Microscopy (AFM) was used to measure the adhesion force between a clean, smooth glass particle (5 µm in diameter) and two glass slides, one prepared according to common laboratory practices (hereafter, clean glass) and a second one left outdoors for six weeks (hereafter, contaminated glass). It was found that the adhesion force distribution was bimodal for the clean glass while it showed three peaks for the contaminated glass. For both surfaces, the adhesion force showed a peak at very low adhesion values (< 100 nN) representing 1.2% of the measurements for the clean glass and 58% for the contaminated glass. This low-adhesion peak is important because it represents the particles that would be resuspended first even in light winds. The low adhesion-force measurements were associated with the presence of contaminants on the surface. Another notable difference between the adhesion force distributions was a long tail at high adhesion forces for clean glass (up to ~1900 nN) and the absence of the tail for the contaminated glass, which had a maximum adhesion force of 950 nN. The mean adhesion force on the contaminated glass was also lower than on the clean glass (206 ± 240 nN vs. 516 ± 354 nN). These observations indicate that the adhesion force distribution between a 5-µm smooth and clean glass sphere and the outdoor contaminated glass surface is quite different from that of the clean glass surface. Hence, one should be cautious when using adhesion force measurements from clean surfaces to estimate particle resuspension from outdoor surfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Novel nonwetting solid-infused surfaces for superior fouling mitigation

Fouling is a ubiquitous issue in several environmental and energy applications. Here we introduce novel nonwetting solid-infused surfaces (SIS) with superior anti-fouling characteristics that are durable than conventional nonwetting surfaces in a dynamic flow environment. A systematic study is presented to elucidate the fouling mitigation performance of SIS in comparison to lubricant-infused surface (LIS) and conventional smooth surface. Copper tubes with SIS, LIS or smooth inner walls are fabricated and subjected to accelerated calcium sulfate fouling in a flow fouling experimental setup. Fouling on the various surface types is quantified in terms of asymptotic fouling resistance, and the fundamental morphological differences in the interactions of the foulant and the various surface types are analyzed. Based on a systematic sweep of the parameter combinations using design of experiments and Taguchi analysis, an analytical dependence of asymptotic fouling resistance on the governing parameters namely, Reynolds number, foulant concentration and temperature is derived. The analytical model is shown to predict the asymptotic fouling resistance to within 20% accuracy with a 95% confidence. In addition, for the first time, the effects of shear durability on the fouling mitigation performance of LIS vis-à-vis SIS are studied. It is shown that the novel nonwetting SIS offers a robust option for superior fouling mitigation over LIS in the long run.

42 ENGINEERING↗
DOE OSTI2020

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

A frozen Gaussian approximation-based multi-level particle swarm optimization for seismic inversion

In this paper, we propose a frozen Gaussian approximation (FGA)-based multi-level particle swarm optimization (MLPSO) method for seismic inversion of high-frequency wave data. The method addresses two challenges in it: First, the optimization problem is highly non-convex, which makes hard for gradient-based methods to reach global minima. This is tackled by MLPSO which can escape from undesired local minima. Second, the character of high-frequency of seismic waves requires a large number of grid points in direct computational methods, and thus renders an extremely high computational demand on the simulation of each sample in MLPSO. We overcome this difficulty by three steps: First, we use FGA to compute high-frequency wave propagation based on asymptotic analysis on phase plane; Then we design a constrained full waveform inversion problem to prevent the optimization search getting into regions of velocity where FGA is not accurate; Last, we solve the constrained optimization problem by MLPSO that employs FGA solvers with different fidelity. The performance of the proposed method is demonstrated by a two-dimensional full-waveform inversion example of the smoothed Marmousi model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Multiscale Universal Interface: A concurrent framework for coupling heterogeneous solvers

Graphical abstract: - Abstract: Concurrently coupled numerical simulations using heterogeneous solvers are powerful tools for modeling multiscale phenomena. However, major modifications to existing codes are often required to enable such simulations, posing significant difficulties in practice. In this paper we present a C++ library, i.e. the Multiscale Universal Interface (MUI), which is capable of facilitating the coupling effort for a wide range of multiscale simulations. The library adopts a header-only form with minimal external dependency and hence can be easily dropped into existing codes. A data sampler concept is introduced, combined with a hybrid dynamic/static typing mechanism, to create an easily customizable framework for solver-independent data interpretation. The library integrates MPI MPMD support and an asynchronous communication protocol to handle inter-solver information exchange irrespective of the solvers' own MPI awareness. Template metaprogramming is heavily employed to simultaneously improve runtime performance and code flexibility. We validated the library by solving three different multiscale problems, which also serve to demonstrate the flexibility of the framework in handling heterogeneous models and solvers. In the first example, a Couette flow was simulated using two concurrently coupled Smoothed Particle Hydrodynamics (SPH) simulations of different spatial resolutions. In the second example, we coupled the deterministic SPH method with the stochastic Dissipative Particle Dynamics (DPD) method to study the effect of surface grafting on the hydrodynamics properties on the surface. In the third example, we consider conjugate heat transfer between a solid domain and a fluid domain by coupling the particle-based energy-conserving DPD (eDPD) method with the Finite Element Method (FEM)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A novel coupling of noise reduction algorithms for particle flow simulations

Proper orthogonal decomposition (POD) and its extension based on time-windows have been shown to greatly improve the effectiveness of recovering smooth ensemble solutions from noisy particle data. However, to successfully de-noise any molecular system, a large number of measurements still need to be provided. In order to achieve a better efficiency in processing time-dependent fields, we have combined POD with a well-established signal processing technique, wavelet-based thresholding. In this novel hybrid procedure, the wavelet filtering is applied within the POD domain and referred to as WAVinPOD. The algorithm exhibits promising results when applied to both synthetically generated signals and particle data. In this work, the simulations compare the performance of our new approach with standard POD or wavelet analysis in extracting smooth profiles from noisy velocity and density fields. Numerical examples include molecular dynamics and dissipative particle dynamics simulations of unsteady force- and shear-driven liquid flows, as well as phase separation phenomenon. Simulation results confirm that WAVinPOD preserves the dimensionality reduction obtained using POD, while improving its filtering properties through the sparse representation of data in wavelet basis. This paper shows that WAVinPOD outperforms the other estimators for both synthetically generated signals and particle-based measurements, achieving a higher signal-to-noise ratio from a smaller number of samples. The new filtering methodology offers significant computational savings, particularly for multi-scale applications seeking to couple continuum informations with atomistic models. It is the first time that a rigorous analysis has compared de-noising techniques for particle-based fluid simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Seismic modeling with radial basis function-generated finite differences (RBF-FD) – a simplified treatment of interfaces

In a previous study of seismic modeling with radial basis function-generated finite differences (RBF-FD), we outlined a numerical method for solving 2-D wave equations in domains with material interfaces between different regions. The method was applicable on a mesh-free set of data nodes. It included all information about interfaces within the weights of the stencils (allowing the use of traditional time integrators), and was shown to solve problems of the 2-D elastic wave equation to 3rd-order accuracy. In the present paper, we discuss a refinement of that method that makes it simpler to implement. It can also improve accuracy for the case of smoothly-variable model parameter values near interfaces. We give several test cases that demonstrate the method solving 2-D elastic wave equation problems to 4th-order accuracy, even in the presence of smoothly-curved interfaces with jump discontinuities in the model parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Construction and simulation of a novel continuous traffic flow model

In this paper, we aim to propose a novel mathematical model for traffic flow and apply a newly developed characteristic particle method to solve the associate governing equations. As compared with the existing non-equilibrium higher-order traffic flow models, the present one is put forward to satisfy the following three conditions: 1.Preserve the equilibrium state in the smooth region. 2.Yield an anisotropic propagation of traffic flow information. 3.Expressed with a conservation law form for traffic momentum. These conditions will ensure a more practical simulation in traffic flow physics: The current traffic will not be influenced by the condition in the behind and result in unambiguous condition across a traffic shock. Through analyses of characteristics, stability condition and steady-state solution adherent to the equation system, it is shown that the proposed model actually conform to these conditions. Furthermore, this model can be cast into its characteristic form which, incorporated with the Rankine-Hugoniot relation, is appropriate to be simulated by the characteristic particle method to obtain accurate computational results. - Highlights: • The traffic model expressed with the momentum conservation law. • Traffic flow information propagate anisotropically and preserve the equilibrium state in the smooth region. • Computational particles of two families are invented to mimic forward-running and backward-running characteristics. • Formation of shocks will be naturally detected by the intersection of computational particles of same family. • A newly developed characteristic particle method is used to simulate traffic flow model equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Adaptive Activation Functions Accelerate Convergence in Deep and Physics-informed Neural Networks

We employ adaptive activation functions for regression in deep and physics-informed neural networks (PINNs) to approximate smooth and discontinuous functions as well as solutions of linear and nonlinear partial differential equations. In particular, we solve the nonlinear Klein-Gordon equation, which has smooth solutions, the nonlinear Burgers equation, which can admit high gradient solutions, and the Helmholtz equation. We introduce a scalable hyper-parameter in the activation function, which can be optimized to achieve best performance of the network as it changes dynamically the topology of the loss function involved in the optimization process. The adaptive activation function has better learning capabilities than the traditional one (fixed activation) as it improves greatly the convergence rate, especially at early training, as well as the solution accuracy. To better understand the learning process, we plot the neural network solution in the frequency domain to examine how the network captures successively different frequency bands present in the solution. We consider both forward problems, where the approximate solutions are obtained, as well as inverse problems, where parameters involved in the governing equation are identified. Our simulation results show that the proposed method is a very simple and effective approach to increase the efficiency, robustness and accuracy of the neural network approximation of nonlinear functions as well as solutions of partial differential equations, especially for forward problems. We theoretically prove that in the proposed method, gradient descent algorithms are not attracted to suboptimal critical points or local minima.

machine leaning, Bad minima, Inverse problems, Phy↗
DOE OSTI2021

Higher-order particle representation for particle-in-cell simulations

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A short note on the accuracy of the discontinuous Galerkin method with reentrant faces

In this work, we study the convergence of the discontinuous Galerkin (DG) method applied to the advection–reaction equation on meshes with reentrant faces. On such meshes, the upwind numerical flux is not smooth, and so the numerical integration of the resulting face terms can only be expected to be first-order accurate. Despite this inexact integration, we prove that the DG method converges with order $\mathscr{O}$(h p+1/2 ), which is the same rate as in the case of exact integration. Consequently, specialized quadrature rules that accurately integrate the non-smooth numerical fluxes are not required for high-order accuracy. These results are numerically corroborated on examples of linear advection and discrete ordinates transport equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗