Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonconforming”

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.

33 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗
DOE OSTI2021

Scale-fixed predictions for γ + η c production in electron-positron collisions at NNLO in perturbative QCD

In the paper, we present QCD predictions for γ + η c production at an electron-positron collider up to next-to-next-to-leading order (NNLO) accuracy without renormalization scale ambiguities. The NNLO total cross-section for e + + e – → γ + η c using the conventional scale-setting approach has large renormalization scale ambiguities, usually estimated by choosing the renormalization scale to be the e + e – center-of-mass collision energy $\sqrt{s}$. The Principle of Maximum Conformality (PMC) provides a systematic way to eliminate such renormalization scale ambiguities by summing the nonconformal β contributions into the QCD coupling α s ( Q 2 ). The renormalization group equation then sets the value of α s for the process. The PMC renormalization scale reflects the virtuality of the underlying process, and the resulting predictions satisfy all of the requirements of renormalization group invariance, including renormalization scheme invariance. After applying the PMC, we obtain a renormalization scale-and-scheme independent prediction, σ | NNLO,PMC ≃ 41.18 fb for $\sqrt{s}$=10.6 GeV. The resulting pQCD series matches the series for conformal theory and thus has no divergent renormalon contributions. The large K factor which contributes to this process reinforces the importance of uncalculated NNNLO and higher-order terms. Using the PMC scale-and-scheme independent conformal series and the Padé approximation approach, we predict σ | NNNLO,PMC+Pade ≃ 18.99 fb, which is consistent with the recent BELLE measurement ${\sigma}^{\mathrm{obs}}={16.58}_{-9.93}^{+10.51}$ fb at $\sqrt{s}$ ≃ 10.6 GeV. This procedure also provides a first estimate of the NNNLO contribution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Uniform Subspace Correction Preconditioners for Discontinuous Galerkin Methods with hp-Refinement

In this paper, we develop subspace correction preconditioners for discontinuous Galerkin (DG) discretizations of elliptic problems with hp-refinement. These preconditioners are based on the decomposition of the DG finite element space into a conforming subspace, and a set of small nonconforming edge spaces. The conforming subspace is preconditioned using a matrix-free low-order refined technique, which in this work, we extend to the hp-refinement context using a variational restriction approach. The condition number of the resulting linear system is independent of the granularity of the mesh h, and the degree of the polynomial approximation p. The method is amenable to use with meshes of any degree of irregularity and arbitrary distribution of polynomial degrees. Furthermore, numerical examples are shown on several test cases involving adaptively and randomly refined meshes, using both the symmetric interior penalty method and the second method of Bassi and Rebay (BR2).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Supported Oxygen Evolution Catalysts by Design: Toward Lower Precious Metal Loading and Improved Conductivity in Proton Exchange Membrane Water Electrolyzers

Reducing the precious metal content of water oxidation catalysts for proton-exchange-membrane water electrolyzers remains a critical barrier to their large-scale deployment. Herein, we present an engineered architecture for supported iridium catalysts, which enables decreased precious metal content and improved activity and conductivity. The improvement in performance at lower precious metal loading is realized by the deposition of a conformal layer of platinum nanoparticles on titanium dioxide (TiO 2 ) using a facile photoreduction method to prepare conductive layer coated supports (CCSs). Platinum nanoparticles are homogeneously dispersed on TiO 2 , and the conductivity of the subsequent catalysts with 39 wt % precious group metal loadings is significantly higher than the commercial 75 wt % loaded IrO 2 -TiO 2 catalysts. The conformal conductive layer also maintains an enhanced conductivity and electrochemical activity upon thermal annealing when compared to catalysts without the conductive layer and nonconformal heterogeneous conductive layer. The iridium mass activity from half-cell studies shows a 141% improvement for CCS supported catalysts at 42% lower loadings compared to the commercial catalysts. The conductive layer also improves the single cell electrolyzer performance at a similar catalyst loading in comparison to a commercial state-of-the-art catalyst. Here, we correlate the physical properties of the engineered catalysts with their electrochemical performance in electrolyzers to understand structure-activity relationships, and we anticipate further performance improvements upon synthesis and materials optimizations.

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

Simulating Sea-Ice Deformation in Viscous-Plastic Sea-Ice Models With CD-Grids

Linear kinematic features (LKFs) are found everywhere in the Arctic sea-ice cover. They are strongly localized deformations often associated with the formation of leads and pressure ridges. In viscous-plastic (VP) sea-ice models, the simulation of LKFs depends on several factors such as the grid resolution, the numerical solver convergence, and the placement of the variables on the mesh. In this study, we compare two recently proposed discretization with a CD-grid placement with respect to their ability to reproduce LKFs. The first (CD1) is based on a nonconforming finite element discretization, whereas the second (CD2) uses a conforming subgrid discretization. To analyze their resolution properties, we evaluate runs from different models (e.g., FESOM, MPAS) on a benchmark problem using quadrilateral, hexagonal and triangular meshes. Our findings show that the CD1 setup simulates more deformation structure than the CD2 setup. This highlights the importance of the type of spatial discretization for the simulation of LKFs. Due to the higher number of degrees of freedom, both CD-grids resolve more LKFs than traditional A, B, and C-grids at fixed mesh level. This is an advantage of the CD-grid approach, as high spatial mesh resolution is needed in VP sea-ice models to simulate LKFs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Reanalysis of the top-quark pair hadroproduction and a precise determination of the top-quark pole mass at the LHC *

In this study, we calculate the $t\bar{t}$ pQCD production cross-section at the NNLO and determine the top-quark pole mass from recent measurements at the LHC at the center-of-mass energy $\sqrt{S}=13$ TeV to a high precision by applying the principle of maximum conformality (PMC). The PMC provides a systematic method that rigorously eliminates QCD renormalization scale ambiguities by summing the nonconformal $\beta$ contributions into the QCD coupling constant. The PMC predictions satisfy the requirements of renormalization group invariance, including renormalization scheme independence, and the PMC scales accurately reflect the virtuality of the underlying production subprocesses. By using the PMC, an improved prediction for the $t\bar{t}$ production cross-section is obtained without scale ambiguities, which in turn provides a precise value for the top-quark pole mass. Moreover, the prediction of PMC calculations that the magnitudes of higher-order PMC predictions are well within the error bars predicted from the known lower-order has been demonstrated for the top-quark pair production. The resulting determination of the top-quark pole mass, $m_t^{\rm{pole}}=172.5\pm1.4$ GeV, from the LHC measurement at $\sqrt{S}=13$ TeV agrees with the current world average cited by the Particle Data Group (PDG). The PMC prediction provides an important high-precision test of the consistency of pQCD and the SM at $\sqrt{S}=13$ TeV with previous LHC measurements at lower CM energies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Self-consistent analysis for the η c → γ γ process

The next-to-next-to-leading-order (NNLO) perturbative QCD (pQCD) predictions for both the decay width and the transition form factor in the η c → γ γ process, based on nonrelativistic QCD (NRQCD), deviate from precise experimental measurements. These significant discrepancies have cast doubt on the applicability of NRQCD to charmonium processes. In this paper, we analyze the η c → γ γ process by applying the principle of maximum conformality (PMC), a systematic method for eliminating renormalization scheme-and-scale ambiguities. The PMC renormalization scales are determined by absorbing the nonconformal β -terms that govern the behavior of the QCD running coupling via the renormalization group equation. We obtain the PMC scale Q ⋆ = 4.49 m c for the η c → γ γ decay width. Even after using the PMC method, the convergence of the pQCD series is still poor, which indicates the importance of uncalculated next-to-next-to-next-to-leading-order and higher-order terms. The resulting value for Γ η c → γ γ is in agreement with the Particle Data Group’s reported value of Γ η c → γ γ = 5.1 ± 0.4 keV within the bounds of uncertainties. Moreover, the transition form factor obtained using the PMC is also in good agreement with precise experimental measurements. The application of the PMC suggests a potential resolution to η c → γ γ puzzle and supports the applicability of NRQCD to charmonium processes. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Partitioned Coupling for Structural Acoustics

Here, we expand the second-order fluid–structure coupling scheme of Farhat et al. (1998, “Load and Motion Transfer Algorithms for 19 Fluid/Structure Interaction Problems With Non-Matching Discrete Interfaces: Momentum and Energy Conservation, Optimal Discretization and Application to Aeroelasticity,” Comput. Methods Appl. Mech. Eng., 157(1–2), pp. 95–114; 2006, “Provably Second-Order Time-Accurate Loosely-Coupled Solution Algorithms for Transient Nonlinear Computational Aeroelasticity,” Comput. Methods Appl. Mech. Eng., 195(17), pp. 1973–2001) to structural acoustics. The staggered structural acoustics solution method is demonstrated to be second-order accurate in time, and numerical results are compared to a monolithically coupled system. The partitioned coupling method is implemented in the Sierra Mechanics software suite, allowing for the loose coupling of time domain acoustics in SIERRA/SD to structural dynamics ( SIERRA/SD ) or solid mechanics ( SIERRA/SD ). The coupling is demonstrated to work for nonconforming meshes. Results are verified for a one-dimensional piston, and the staggered and monolithic results are compared to an exact solution. Huang, H. (1969, “Transient Interaction of Plane Acoustic Waves With a Spherical Elastic Shell,” J. Acoust. Soc. Am., 45(3), pp. 661–670) sphere scattering problem with a spherically spreading acoustic load demonstrates parallel capability on a complex problem. Our numerical results compare well for a bronze plate submerged in water and sinusoidally excited (Fahnline and Shepherd, 2017, “Transient Finite Element/Equivalent Sources Using Direct Coupling and Treating the Acoustic Coupling Matrix as Sparse,” J. Acoust. Soc. Am., 142(2), pp. 1011–1024).

42 ENGINEERING↗
DOE OSTI2024

Energy Exascale Computational Fluid Dynamics Simulations With the Spectral Element Method

Development and application of the open-source GPU-based fluid-thermal simulation code, NekRS, are described. Time advancement is based on an efficient kth-order accurate timesplit formulation coupled with scalable iterative solvers. Spatial discretization is based on the high-order spectral element method (SEM), which affords the use of fast, low-memory, matrix-free operator evaluation. Further, recent developments include support for nonconforming meshes using overset grids and for GPU-based Lagrangian particle tracking. Results of large-eddy simulations of atmospheric boundary layers for wind-energy applications as well as extensive nuclear energy applications are presented.

42 ENGINEERING↗
DOE OSTI2023

End-to-end GPU acceleration of low-order-refined preconditioning for high-order finite element discretizations

In this article, we present algorithms and implementations for the end-to-end GPU acceleration of matrix-free low-order-refined preconditioning of high-order finite element problems. The methods described here allow for the construction of effective preconditioners for high-order problems with optimal memory usage and computational complexity. The preconditioners are based on the construction of a spectrally equivalent low-order discretization on a refined mesh, which is then amenable to, for example, algebraic multigrid preconditioning. The constants of equivalence are independent of mesh size and polynomial degree. For vector finite element problems in H(curl) and H(div) (e.g., for electromagnetic or radiation diffusion problems), a specially constructed interpolation–histopolation basis is used to ensure fast convergence. Detailed performance studies are carried out to analyze the efficiency of the GPU algorithms. The kernel throughput of each of the main algorithmic components is measured, and the strong and weak parallel scalability of the methods is demonstrated. The different relative weighting and significance of the algorithmic components on GPUs and CPUs is discussed. Results on problems involving adaptively refined nonconforming meshes are shown, and the use of the preconditioners on a large-scale magnetic diffusion problem using all spaces of the finite element de Rham complex is illustrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A deep learning-enhanced framework for multiphysics joint inversion

Joint inversion has drawn considerable attention due to the availability of multiple geophysical data sets, ever-increasing computational resources, the development of advanced algorithms, and its ability to reduce inversion uncertainty. A key issue of joint inversion is to develop effective strategies to link different geophysical data in a unified mathematical framework, in which the information obtained from different models can complement each other. We have developed a deep learning-enhanced joint inversion framework to simultaneously reconstruct different physical models by fusing different types of geophysical data. Traditionally, structure similarity constraints are pursued by joint inversion algorithms using manually crafted formulations (e.g., cross gradient). The constraint is constructed by a deep neural network (DNN) during the learning process. The framework is designed to combine the DNN and a traditional independent inversion workflow and improve the joint inversion result iteratively. The network can be easily extended to incorporate multiphysics without structural changes. Numerical experiments on the joint inversion of 2D DC resistivity data and seismic traveltime are used to validate our method. In addition, this learning-based framework demonstrates excellent generalization abilities when tested on data sets using different geologic structures. It also can handle different sensing configurations and nonconforming discretization.

Geochemistry & Geophysics↗
DOE OSTI2019

Site Environmental Report for 2018

Lawrence Berkeley National Laboratory (LBNL, Berkeley Lab) is a multi-program scientific facility operated by the University of California (UC) for the U.S. Department of Energy (DOE). Berkeley Lab’s research is focused on the physical, biological, environmental, and computational sciences, with the objective of delivering scientific knowledge and discoveries pertinent to DOE’s mission. This annual report describes environmental protection activities and potential impacts resulting from LBNL operations conducted in 2018. The format and content of this report satisfy the requirements of both DOE Order 231.1B, Administrative Change 1 (Environment, Safety, and Health Reporting) and the operating contract between UC and DOE (DOE Contract No. DE-AC02-05CH11231, also known as Contract 31). LBNL activities are planned and conducted with full regard to protecting employees, the public, and the environment, as well as complying with all applicable environmental, safety, and health laws and regulations. As presented in Chapter 2, Berkeley Lab implements an Environmental Management System (EMS) to oversee environmental compliance activities and continually improve overall environmental performance while maintaining operational capability and sustaining its overall mission. The effectiveness of the EMS and environmental programs is reviewed annually as part of the performance evaluation process of Contract 31. For fiscal year (FY) 2018, which began October 1, 2017, and ended September 30, 2018, the EMS was given a performance rating of B plus for its management of environmental activities (on a scale from A plus as the highest grade, to F as the lowest). The measurement and rating system was developed jointly by Berkeley Lab, UC, and DOE. The FY 2018 rating was based on how Berkeley Lab met the objective in DOE’s FY 2018 Performance Evaluation and Measurement Plan (Appendix B in Contract 31, Section J) of providing an efficient and effective EMS. Six significant accomplishments, which ranged from effective teaming on projects and with regulatory agencies to improved program assessment approaches, were factors in the rating. The EMS was also graded through the federal Office of Management and Budget’s annual EMS performance metrics, in which a reporting scorecard rates elements of the International Organization for Standardization (ISO) 14001 standard and the degree of integration between the EMS and Berkeley Lab’s sustainable practices. Overall scores fall into one of three categories: green (highest), yellow (middle), or red (lowest). For FY 2018, Berkeley Lab received a score of green. A week-long formal conformity audit of the EMS was conducted at the end of July by a qualified external party, as is required every three years. The audit team found that Berkeley Lab’s EMS continues to conform to the ISO standard with no major nonconformance. An overview of environmental protection and restoration programs is provided in Chapter 3, including information about compliance activities, operating permits, and regulatory agency inspections and audits that occurred during 2018. Twenty-one minor violations and two Class II violations were issued during agency inspections of programs governing aboveground and underground storage tanks, hazardous waste treatment units, hazardous waste storage areas, hazardous materials management, and wastewater, stormwater, and air emission sources. This report also includes information on environmental monitoring performed in 2018 (Chapter 4). The results of these monitoring activities confirmed that groundwater cleanup actions have been effective in reducing concentrations of volatile organic compounds (VOCs) in the groundwater, although concentrations appear to be reaching asymptotic levels in some areas. Site groundwater plumes are stable or are attenuating, and VOCs are not migrating off site. All emissions and discharges from LBNL operations were within environmental compliance release limits except oxides of nitrogen and carbon monoxide from boilers (under 5 million BTUs) from Buildings 2, 66, and 88. The radiological dose assessments (Chapter 5) performed in 2018 concluded that the maximum potential dose to a hypothetical resident from Berkeley Lab’s airborne radionuclide releases was approximately 0.04% of the DOE and U.S. Environmental Protection Agency annual limit of 10 millirem per year (mrem/yr); the potential dose from all radiation sources at Berkeley Lab was approximately 0.15% of the average natural background radiation dose of 310 mrem/yr in the United States, and about 0.5% of the DOE annual limit of 100 mrem/yr from all sources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Fiberboard Properties and Degradation Rates for Storage of 9975 Shipping Packages in KAC (2021 Status Report)

Thermal, mechanical, and physical properties have been measured on fiberboard samples following accelerated aging for up to 16 years. The aging environments include elevated temperature up to 250 ºF (the maximum allowed service temperature for fiberboard in 9975 packages) and elevated humidity. Accelerated aging results have been analyzed and used to build aging models. Correlations relating several properties (thermal conductivity, energy absorption, weight, dimensions, and density) to their rate of change in potential storage environments have been developed. Combined with an estimate of the actual conditions the fiberboard experiences in K-Area Complex (KAC), these models allow development of service life predictions. KAC completed calculations that supported extending the service life of 9975 packages with 3013 and non3013 containers in storage from 15 years to 20 years. These calculations addressed the potential for degraded fiberboard properties following 20 years in storage using degradation models developed previously. The updated degradation models developed in this report have been compared to the assumptions in these calculations, and the current results for fiberboard dimensions, density, radial thermal conductivity, and energy absorption remain consistent with those calculations. The axial thermal conductivity degradation rate slightly exceeds calculations estimate of 0.5% per year between 130 and 160 °F. Additional data will continue to be collected to permit future refinements to the models and assumptions. Some of the predicted degradation rates presented in this report are extreme. However, these relate to environments that do not exist within KAC or would be postulated only as upset conditions that would not likely persist for an extended period. For a typical package stored in KAC with approximately 10-12 watts internal heat load or less, and ambient temperatures below 90 ºF, the fiberboard assemblies experience storage conditions less severe than any of the aging environments. Fiberboard in conforming packages with lower internal heat loads should experience little or no degradation and is expected to provide a service life beyond the currently approved 20 year storage period. Packages with higher internal heat loads may not continue to perform their required safety functions beyond 20 years. Ultimately, service life will be determined by the cumulative effect of degradation from all the conditions these packages might encounter. The results and model predictions presented in this report are applicable to 9975 packages with cane or softwood fiberboard overpack assemblies. These degradation models do not address the effects of nonconforming conditions such as the presence of excess moisture and mold or beetle infestations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

ISO 14001:2015 Environmental Management System Management Assessment Internal Review Assessing Organization: EPC-ES LANL-ASMT-2022-0289 (Spring 2023)

The Environmental Protection and Compliance Division-Environmental Stewardship Group of Los Alamos National Laboratory conducted a management assessment that supports the implementation of the Los Alamos National Laboratory Environmental Management System according to International Organization for Standardization (ISO) 14001:2015, Environmental Management System. The assessment resulted in zero findings of nonconformity, five opportunities for improvement, four noteworthy practices, and three observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces

In this work, we analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously unexplored relation between the two types of summation-by-parts operators is investigated. By combining them we develop a new fourth order accurate finite difference discretization with hanging nodes on the mesh refinement interface. We take the model problem as the two-dimensional acoustic wave equation in second order form in terms of acoustic pressure, and we prove energy stability for the proposed method. Compared to previous approaches using ghost points, the proposed method leads to a smaller system of linear equations that needs to be solved for the ghost point values. Another attractive feature of the proposed method is that the explicit time step does not need to be reduced relative to the corresponding periodic problem. Numerical experiments, both for smoothly varying and discontinuous material properties, demonstrate that the proposed method converges to fourth order accuracy. A detailed comparison of the accuracy and the time-step restriction with the simultaneous-approximation-term penalty method is also presented. (An erratum is attached.)

97 MATHEMATICS AND COMPUTING↗