Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regularity”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2022

A neural network‐enhanced reproducing kernel particle method for modeling strain localization

Abstract Modeling the localized intensive deformation in a damaged solid requires highly refined discretization for accurate prediction, which significantly increases the computational cost. Although adaptive model refinement can be employed for enhanced effectiveness, it is cumbersome for the traditional mesh‐based methods to perform while modeling the evolving localizations. In this work, neural network‐enhanced reproducing kernel particle method (NN‐RKPM) is proposed, where the location, orientation, and shape of the solution transition near a localization is automatically captured by the NN approximation via a block‐level neural network (NN) optimization. The weights and biases in the blocked parameterization network control the location and orientation of the localization. The designed basic four‐kernel NN block is capable of capturing a triple junction or a quadruple junction topological pattern, while more complicated localization topological patters are captured by the superposition of multiple four‐kernel NN blocks. The standard RK approximation is then utilized to approximate the smooth part of the solution, which permits a much coarser discretization than the high‐resolution discretization needed to capture sharp solution transitions with the conventional methods. A regularization of the NN approximation is additionally introduced for discretization‐independent material responses. The effectiveness of the proposed NN‐RKPM is verified by a series of numerical verifications.

Baek, Jonghyuk↗
DOE OSTI2023

Local element operations for curved simplex meshes

Summary Mesh optimization procedures are generally a combination of node smoothing and discrete operations which affect a small number of elements to improve the quality of the overall mesh. These procedures are useful as a post‐processing step in mesh generation procedures and in applications such as fluid simulations with severely deforming domains. In order to perform high‐order mesh optimization, these ingredients must also be extended to high‐order (curved) meshes. In this work, we present a method to perform local element operations on curved meshes. The mesh operations discussed in this work are edge/face swaps, edge collapses, and edge splitting (more generally refinement) for triangular and tetrahedral meshes. These local operations are performed by first identifying the patch of elements which contain the edge/face being acted on, performing the operation as a “straight‐sided one" by placing the high‐order nodes via an isoparametric mapping from the master element, and smoothing the high‐order nodes on the elements in the patch by minimizing a Jacobian‐based high‐order mesh distortion measure. Since the initial “straight‐sided guess” from the placement of the nodes via the isoparametric mapping frequently results in invalid elements, the distortion measure must be regularized which allows for mesh untangling for the optimization to succeed. We present several examples in 2D and 3D to demonstrate these local operations and how they can be combined with a high‐order node smoothing procedure to maintain mesh quality when faced with severe deformations.

Shi, A.↗
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2024

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗
DOE OSTI2020

A least-squares finite element method based on the Helmholtz decomposition for hyperbolic balance laws

In this paper, a least-squares finite element method for scalar nonlinear hyperbolic balance laws is proposed and studied. The approach is based on a formulation that utilizes an appropriate Helmholtz decomposition of the flux vector and is related to the standard notion of a weak solution. This relationship, together with a corresponding connection to negative-norm least-squares, is described in detail. As a consequence, an important numerical conservation theorem is obtained, similar to the famous Lax–Wendroff theorem. The numerical conservation properties of the method in this paper do not fall precisely in the framework introduced by Lax and Wendroff, but they are similar in spirit as they guarantee that when L 2 convergence holds, the resulting approximations approach a weak solution to the hyperbolic problem. The least-squares functional is continuous and coercive in an H -1 -type norm, but not L 2 -coercive. Nevertheless, the L 2 convergence properties of the method are discussed. Convergence can be obtained either by an explicit regularization of the functional, that provides control of the L 2 norm, or by properly choosing the finite element spaces, providing implicit control of the L 2 norm. Numerical results for the inviscid Burgers equation with discontinuous source terms are shown, demonstrating the L 2 convergence of the obtained approximations to the physically admissible solution. The numerical method utilizes a least-squares functional, minimized on finite element spaces, and a Gauss–Newton technique with nested iteration. Finally, we believe that the linear systems encountered with this formulation are amenable to multigrid techniques and combining the method with adaptive mesh refinement would make this approach an efficient tool for solving balance laws (this is the focus of a future study).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Efficient quadrature rules for finite element discretizations of nonlocal equations

Here, in this paper, we design efficient quadrature rules for finite element (FE) discretizations of nonlocal diffusion problems with compactly supported kernel functions. Two of the main challenges in nonlocal modeling and simulations are the prohibitive computational cost and the nontrivial implementation of discretization schemes, especially in three-dimensional settings. In this work, we circumvent both challenges by introducing a parametrized mollifying function that improves the regularity of the integrand, utilizing an adaptive integration technique, and exploiting parallelization. We first show that the “mollified” solution converges to the exact one as the mollifying parameter vanishes, then we illustrate the consistency and accuracy of the proposed method on several two- and three-dimensional test cases. Furthermore, we demonstrate the good scaling properties of the parallel implementation of the adaptive algorithm and we compare the proposed method with recently developed techniques for efficient FE assembly.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A fractional model for anomalous diffusion with increased variability: Analysis, algorithms and applications to interface problems

Fractional equations have become the model of choice in several applications where heterogeneities at the microstructure result in anomalous diffusive behavior at the macroscale. Here, we introduce a new fractional operator characterized by a doubly-variable fractional order and possibly truncated interactions. Under certain conditions on the model parameters and on the regularity of the fractional order we show that the corresponding Poisson problem is well-posed. Additionally, we introduce a finite element discretization and describe an efficient implementation of the finite-element matrix assembly in the case of piecewise constant fractional order. Through several numerical tests, we illustrate the improved descriptive power of this new operator across media interfaces. Furthermore, we present one-dimensional and two-dimensional h-convergence results that show that the variable-order model has the same convergence behavior as the constant-order model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Photovoltaic module antireflection coating degradation survey using color microscopy and spectral reflectance

Abstract Commercial antireflective coatings (ARCs) on photovoltaic (PV) module glass can improve module power by 2.5%–3.0%, but their long‐term field performance requires additional study. In this paper, we investigate ARC performance on fielded modules using two nondestructive techniques: reflectance spectroscopy and RGB microscopy, finding a large variation in coating durability and performance. For new coatings, the fleet average nominal power enhancement is 2.8%. This power enhancement is observed to degrade at approximately −0.05%/year absolute over the first 8 years of fielding in arrays that are not regularly cleaned. Interferometry and imaging results imply that coating loss for these modules is due to slow chemical thinning of the initial 125 nm thick coating by 1.4 to 5.0 nm/year. Due to the physics of interference coatings, the performance loss is projected to accelerate as the coating is further thinned, resulting in a coating lifetime (time to 80% of initial performance) of 7.5 to 25 years. At the higher 5.0 nm/year coating loss rates, we estimate that −0.14%/year power degradation of the module can be attributed solely to ARC degradation over the first 20 years. Extreme coating loss is observed on some modules where after 8 years fielding, the coating can be completely removed with a single wet wipe with a lens tissue. A case study is also presented comparing ARC and noncoated modules installed at a single site in 2013. We find that the ARC modules have a 2.6% higher soiling power loss than the noncoated modules; this exceeds the ARC power enhancement of 2.4% and leads to the surprising conclusion that the ARC lowers production for modules in this location due to increased soiling. These results demonstrate that RGB microscopy is a powerful, field‐capable, and quantitative characterization technique for assessing degradation of PV module ARCs.

Karin, Todd↗
DOE OSTI2021

PDB ‐101: Educational resources supporting molecular explorations through biology and medicine

Abstract The Protein Data Bank (PDB) archive is a rich source of information in the form of atomic‐level three‐dimensional (3D) structures of biomolecules experimentally determined using macromolecular crystallography, nuclear magnetic resonance (NMR) spectroscopy, and electron microscopy (3DEM). Originally established in 1971 as a resource for protein crystallographers to freely exchange data, today PDB data drive research and education across scientific disciplines. In 2011, the online portal PDB‐101 was launched to support teachers, students, and the general public in PDB archive exploration ( pdb101.rcsb.org ). Maintained by the Research Collaboratory for Structural Bioinformatics PDB, PDB‐101 aims to help train the next generation of PDB users and to promote the overall importance of structural biology and protein science to nonexperts. Regularly published features include the highly popular Molecule of the Month series, 3D model activities, molecular animation videos, and educational curricula. Materials are organized into various categories (Health and Disease, Molecules of Life, Biotech and Nanotech, and Structures and Structure Determination) and searchable by keyword. A biennial health focus frames new resource creation and provides topics for annual video challenges for high school students. Web analytics document that PDB‐101 materials relating to fundamental topics (e.g., hemoglobin, catalase) are highly accessed year‐on‐year. In addition, PDB‐101 materials created in response to topical health matters (e.g., Zika, measles, coronavirus) are well received. PDB‐101 shows how learning about the diverse shapes and functions of PDB structures promotes understanding of all aspects of biology, from the central dogma of biology to health and disease to biological energy.

Zardecki, Christine↗
DOE OSTI2020

The Table of Standard Atomic Weights—an exercise in consensus

The present Table of Standard Atomic Weights (TSAW) of the elements is perhaps one of the most familiar datasets in science. Unlike most parameters in physical science whose values and uncertainties are evaluated using the “Guidelines for Uncertainty in Measurement” (GUM), the majority of standard atomic-weight values and their uncertainties are consensus values, not GUM-evaluated values. The Commission on Isotopic Abundances and Atomic Weights of the International Union of Pure and Applied Chemistry (IUPAC) regularly evaluates the literature for new isotopic-abundance measurements that can lead to revised standard atomic-weight values, A r (E) for element E. The Commission strives to provide utmost clarity in products it disseminates, namely the TSAW and the Table of Isotopic Compositions of the Elements (TICE). In 2016, it was recognized by the Commission that a guideline recommending expression of uncertainty listed in parentheses following the standard atomic-weight value, e.g. A r (Se) = 78.971(8), is not in accord with the GUM, which suggests that this parenthetic notation be reserved to express standard uncertainty, not the expanded uncertainty used in the TSAW and TICE. In 2017, to eliminate this noncompliance with the GUM, a new format was adopted in which the uncertainty value is specified by the plus-minus symbol, e.g. A r (Se) = 78.971 ± 0.008. To clarify the definition of uncertainty, a new footnote has been added to the TSAW. This footnote emphasizes that an atomic-weight uncertainty is a consensus (decisional) uncertainty. Not only has the Commission shielded users of the TSAW and TICE from unreliable measurements that appear in the literature as a result of unduly small uncertainties, but the aim of IUPAC has been fulfilled by which any scientist, taking any natural sample from commerce or research, can expect their sample atomic weight to lie within A r (E) ± its uncertainty almost all of the time.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Absolute quantitation of propranolol from 200-μm regions of mouse brain and liver thin tissues using laser ablation-dropletProbe-mass spectrometry

Rationale: The ability to quantify drugs and metabolites in tissue with sub-mm resolution is a challenging but much needed capability in pharmaceutical research. To fill this void, a novel surface sampling approach combining laser ablation with the commercial dropletProbe automated liquid surface sampling system (LA-dropletProbe) was developed and presented here.Methods: Parylene C-coated 200 × 200 μm tissue regions of mouse brain and kidney thin tissue sections were analyzed for propranolol by laser ablation of tissue directly into a preformed liquid junction. Propranolol was detected by high performance liquid chromatography with positive ion mode electrospray ionization tandem mass spectrometry. Quantitation was achieved via application of a stable isotope-labeled internal standard and an external calibration curve.Results: Absolute concentration of propranolol determined from 200 × 200 µm tissue regions were compared to propranolol concentrations obtained from 2.3-mm-diameter tissue punches of adjacent, non-coated sections using standard bulk tissue extraction protocols followed by regular HPLC-MS/MS analysis. Average concentration of propranolol in both organs determined by the two employed methods agreed within ±12%. Furthermore, the relative abundance of phase II hydroxypropranolol glucuronide metabolites were recorded and found to be consistent with previous results.Conclusions: This work illustrates that depositing a thin parylene C layer onto thin tissue prior to analysis, which seals the surface and prevents direct liquid extraction of the drug from the tissue, coupled to the novel LA-dropletProbe surface sampling system is a viable approach for sub-mm resolution quantitative drug distribution analysis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

An approach to characterizing spatial aspects of image system blur

Quantitative X-ray radiographic imaging systems that utilize a charged couple device (CCD) camera connected to a thick, monolithic scintillator can exhibit blur that varies spatially across the field of view, especially for thick scintillators used in pulse-power radiography of dynamically compressed objects. Here, a three-point approach to estimating and accounting for this effect is demonstrated by (a) using a local estimation technique to measure the effect of blurring a calibration object at key locations across the field of view, (b) combining each of the local estimates into a spatially varying blurring function via partitions of unity interpolation, and (c) resolving the effects of that blur on the image by solving an ill-posed inverse problem using a spatially varying regularization term. The technique is demonstrated on synthetic examples and actual radiographs collected at the Naval Research Laboratory’s (NRL)Mercury pulsed power facility.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Amorphous K–Co–Mo–S x Chalcogel: A Synergy of Surface Sorption and Ion-Exchange

Chalcogel represents a unique class of meso- to macroporous nanomaterials that offer applications in energy and environmental pursuits. We report the synthesis of an ion-exchangeable amorphous chalcogel using a nominal composition of K 2 CoMo 2 S 10 (KCMS) at room temperature is reported. Synchrotron X-ray pair distribution function (PDF), X-ray absorption near-edge structure (XANES), and extended X-ray absorption fine structure (EXAFS) reveal a plausible local structure of KCMS gel consisting of Mo 5+ 2 and Mo 4+ 3 clusters in the vicinity of di/polysulfides which are covalently linked by Co 2+ ions. The ionically bound K + ions remain in the percolating pores of the Co–Mo–S covalent network. XANES of Co K-edge shows multiple electronic transitions, including quadrupole (1s→3d), shakedown (1s→4p + MLCT), and dipole allowed 1s→4p transitions. Remarkably, despite a lack of regular channels as in some crystalline solids, the amorphous KCMS gel shows ion-exchange properties with UO 2 2+ ions. Additionally, it also presents surface sorption via [S∙∙∙∙UO 2 2+ ] covalent interactions. Overall, this study underscores the synthesis of quaternary chalcogels incorporating alkali metals and their potential to advance separation science for cations and oxo-cationic species by integrating a synergy of surface sorption and ion-exchange.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Probing Distinctive Redox Mechanism in Ni-Rich Cathode Via Real-Time Quick X-Ray Absorption Spectroscopy

X-ray radiation damage on the measuring system has been a critical issue regularly for a long-time exposure to X-ray beam during the in operando characterizations, which is particularly severe when the applied X-ray energy is near the absorption edges (M, L, K, etc.) of the interest element. Here, to minimize the negative effects raised by beam radiation, we employ quick X-ray absorption spectroscopy (QXAS) to study the electrochemical reaction mechanism of a Ni-rich layered structure cathode for lithium-ion batteries. With the advanced QXAS technique, the electronic structure and local coordination environment of the transition metals (TMs) are monitored in-operando with limited radiation damage. Compared to the conventional step-mode X-ray absorption spectroscopy, the QXAS can provide more reliable oxidation state change and more detailed local structure evolutions surrounding TMs (Ni and Co) in Ni-rich layered oxides. By leveraging these advantages of QXAS, we demonstrated that the Ni dominates the electrochemical process with the Co being almost electrochemically inactive. Reversible Ni ions movement between TMs sites and Li sites is also revealed by the time-resolved QXAS technique.

25 ENERGY STORAGE↗
DOE OSTI2024

Implementation and Validation of a Generalized Actuator Disk Parameterization for Wind Turbine Simulations Within the FastEddy Model

ABSTRACT Fast and accurate large‐eddy simulation (LES) of the atmospheric boundary layer plays a crucial role in advancing wind energy research. Long‐duration wind farm studies at turbine‐resolving scales have become increasingly important to understand the intricate interactions between large wind farms and the atmospheric boundary layer. However, the prohibitive computational cost of these turbulence‐ and turbine‐resolving simulations has precluded such modeling to be exercised on a regular basis. To that end, we implement and validate the generalized actuator disk (GAD) model in the computationally efficient, graphics processing unit (GPU)–resident, LES model FastEddy. We perform single‐turbine simulations under three atmospheric stabilities (neutral, unstable, and stable) and compare them against observations from the Scaled Wind Farm Technology (SWiFT) facility and other LES codes from the recent Wakebench turbine wake model benchmark. Our idealized LES results agree well with observed wake velocity deficit and downstream recovery across stability regimes. Turbine response in terms of rotational speed, generated power, torque, and thrust coefficient are well predicted across stability regimes and are consistent with the LES results from the benchmark. The FastEddy simulations are found to be at least two orders of magnitude more efficient than the traditional CPU‐based LES models, opening the door for realistic LES simulations of full wind plants as a viable standard practice.

17 WIND ENERGY↗
DOE OSTI2022

MAPredict: Static Analysis Driven Memory Access Prediction Framework for Modern CPUs

Application memory access patterns are crucial in deciding how much traffic is served by the cache and forwarded to the dynamic random-access memory (DRAM). However, predicting such memory traffic is difficult because of the interplay of prefetchers, compilers, parallel execution, and innovations in manufacturer-specific micro-architectures. This research introduced MAPredict, a static analysis-driven framework that addresses these challenges to predict last-level cache (LLC)-DRAM traffic. By exploring and analyzing the behavior of modern Intel processors, MAPredict formulates cache-aware analytical models. MAPredict invokes these models to predict LLC-DRAM traffic by combining the application model, machine model, and user-provided hints to capture dynamic information. MAPredict successfully predicts LLC-DRAM traffic for different regular access patterns and provides the means to combine static and empirical observations for irregular access patterns. Evaluating 130 workloads from six applications on recent Intel micro-architectures, MAPredict yielded an average accuracy of 99% for streaming, 91% for strided, and 92% for stencil patterns. By coupling static and empirical methods, up to 97% average accuracy was obtained for random access patterns on different micro-architectures.

Monil, M. A. H.↗
DOE OSTI2024

Apparently superluminal superfluids

We consider the superfluid phase of a specific renormalizable relativistic quantum field theory. We prove that, within the regime of validity of perturbation theory and of the superfluid effective theory, there are consistent and regular vortex solutions where the superfluid’s velocity field as traditionally defined smoothly interpolates between zero and arbitrarily large superluminal values. We show that this solution is free of instabilities and of superluminal excitations. We show that, in contrast, a generic vortex solution for an ordinary fluid does develop an instability if the velocity field becomes superluminal. All this questions the characterization of a superfluid velocity field as the actual velocity of “something”.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗