Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mortar method”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Condensed Constrained Nonconforming Mortar-Based Approach for Preconditioning Finite Element Discretization Problems

This paper presents and studies an approach for constructing auxiliary space preconditioners for finite element problems using a constrained nonconforming reformulation that is based on a proposed modified version of the mortar method. The well-known mortar finite element discretization method is modified to admit a local structure, providing an element-by-element or subdomain-by-subdomain assembly property. This is achieved via the introduction of additional trace finite element spaces and degrees of freedom (unknowns) associated with the interfaces between adjacent elements or subdomains. The resulting nonconforming formulation and a reduced-via-static-condensation Schur complement form on the interfaces are used in the construction of auxiliary space preconditioners for a given conforming finite element discretization problem. Overall, the properties of these preconditioners are studied and their performance is illustrated on model second order scalar elliptic problems utilizing high order elements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Isogeometric Bézier dual mortaring: The Kirchhoff–Love shell problem

Here in this paper we develop an isogeometric Bézier dual mortar method for coupling multi-patch Kirchhoff–Love shell structures. The proposed approach weakly enforces the continuity of the solution at patch interfaces through a dual mortar method and can be applied to both conforming and non-conforming discretizations. As the employed dual basis functions have local supports and satisfy the biorthogonality property, the resulting stiffness matrix is sparse. In addition, the coupling accuracy is optimal because the dual basis possesses the polynomial reproduction property. We also formulate the continuity constraints through the Rodrigues’ rotation operator which gives a unified framework for coupling patches that are intersected with G 1 continuity as well as patches that meet at a kink. Several linear and nonlinear examples demonstrated the performance and robustness of the proposed coupling techniques.

42 ENGINEERING↗
DOE OSTI2019

Domain decomposition and multiscale mortar mixed finite element methods for linear elasticity with weak stress symmetry

Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier to impose interface continuity of normal stress or displacement, respectively. By eliminating the interior subdomain variables, the global problem is reduced to an interface problem, which is then solved by an iterative procedure. The condition number of the resulting algebraic interface problem is analyzed for both methods. A multiscale mortar mixed finite element method for the problem of interest on non-matching multiblock grids is also studied. It uses a coarse scale mortar finite element space on the non-matching interfaces to approximate the trace of the displacement and impose weakly the continuity of normal stress. A priori error analysis is performed. It is shown that, with appropriate choice of the mortar space, optimal convergence on the fine scale is obtained for the stress, displacement, and rotation, as well as some superconvergence for the displacement. Computational results are presented in confirmation of the theory of all proposed methods.

Mathematics↗
DOE OSTI2022

NEAMS Technical Area Support in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) framework is a foundational capability used by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to create over 15 different simulation tools for advanced nuclear reactors. Due to this ubiquity, improvements to the framework in support of modeling and simulation goals are critical to the program. These improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work transcribed in this report was conducted in direct support of the simulation tools and has already been deployed. The capabilities outlined in this report include a factor of 10^4 improvement in dependency resolution speed, sorting of user objects, ability to compute residuals and Jacobians together, transfer fixes, support for the mortar method in finite volume discretizations, addition of generalized advection schemes for fluid simulations, 10^2 speedup in some Griffin simulations due to a new matrix-only solve type, and much more.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Simulating the Fracture of Notched Mortar Beams through Extended Finite-Element Method and Peridynamics

This paper simulates fracture in notched mortar beams under three-point bending using extended finite element method (XFEM) and peridynamics. A three-phase microstructure (i.e., cement paste, aggregates, and paste-aggregate interface) is used for constitutive modeling of the mortar to obtain the elastic properties for simulation. In the XFEM approach, the simulated homogenized elastic modulus is used along with the total fracture energy of the cement mortar in a damage model to predict the fracture response of the mortar including crack propagation and its fracture parameters (Mode I stress intensity factor, K IC and critical crack tip opening displacement, CTOD C ). The damage model incorporates a maximum principal stress-based damage initiation criteria and a traction-separation law for damage evolution. In the peridynamics approach, a bond-based model involving a prototype microelastic brittle (PMB) material model is used. The elastic properties and fracture energy release rates are used as inputs in the PMB model, along with the choice of peridynamic horizon size. Comparison with experimental fracture properties ( K IC , CTOD C ) as well as crack propagation paths from digital image correlation show that both the approaches yield satisfactory results, particularly for K IC and crack extension. Furthermore, both these methods can be adopted for fracture simulation of cement-based materials.

42 ENGINEERING↗
DOE OSTI2022

Implement and Test 3D Mortar Contact in BISON

We leverage the extension of the generation of mortar segment meshes to three dimensions in MOOSE’s framework to extend thermomechanical modeling capabilities to problems with three dimensions. A modular approach to gap heat transfer physics using the mortar finite element method was created and documented, mechanical contact was extended to three dimensions—including frictional behavior, performance and ease of use were improved, and steps towards scalability of solid mechanics problems involving contact were taken. Many of these new developments are demonstrated in the simulation of 3D light-water reactor (LWR) problems, where the thermomechanical interface problem is solved using the mortar finite element method. Usage of the mortar framework has improved convergence in 2D problems and has enabled employing friction in 3D problems, of which we show results of a short, local stack of 3D pellets. Consequently, the benefits of mortar in terms of solution convergence and quality are extended to three dimensions. Section 2 discusses fundamental developments that enabled the simulation of practical mortar problems in three dimensions and other general improvements, including the reduction of the derivative container size, the modification of dual basis computations when edge dropping (lack of secondary element projection) takes place, the improvement of conditioning when employing the VCP in-edge dropping conditions, and code usability and quality improvements. These latter code enhancements include the migration of tests using “old” mortar contact constraints to using dual mortar with a semi-smooth Newton solution strategy and the reuse of lower dimensional domains for straightforwardly setting up a mortar thermomechanical LWR problem, i.e. the MOOSE action is employed for mechanical contact and the thermal LWR action is employed to capture the gas conductance, contact, and radiation components of gap heat transfer physics. Independently of the mortar LWR thermal action, we developed a modular approach to gap heat transfer that resides in MOOSE and can be leveraged, e.g., in metallic fuel problems. This approach, whose code design based on MOOSE’s user objects to model specific physics was proposed by the maintenance activity, is detailed in Section 3. Based on the dual mortar finite element method, the frictional contact constraints were extended to three dimensions. A block sheared in two directions in and out of contact with a rigid plane is employed in Section 4 to show the way the approach handles changes in frictional states (e.g. stick to slip) within a competitive number of Newton iterations. Equations and numerical results on the use of the VCP with Cartesian Lagrange multipliers, whose combination enables their direct condensation, are described in Section 5.3. Two-dimensional and three-dimensional BISON LWR simulations are discussed in Section 6. Particularly, a stack of five eccentric pellets with a surface defect is simulated and the effect of pellet-cladding friction is assessed. Finally, conclusions are outlined in Section 7.

42 ENGINEERING↗
DOE OSTI2021

Efficacy of SCMs to mitigate ASR in systems with higher alkali contents assessed by pore solution method

This study demonstrates the efficacy of supplementary cementitious materials (SCMs) to mitigate alkali-silica reaction (ASR) even when used in conjunction with cement of higher alkali contents (up to 1% Na{sub 2}O{sub eq}). The expansion of concrete prisms was studied immersed in simulated pore solution in order to address the limitations of conventional ASR testing methods, accelerated mortar bar test (AMBT) and concrete prism test (CPT). Expansion results demonstrate that 25% fly ash and 50% slag are both sufficient to mitigate ASR even with cements with alkali content up to 1% Na{sub 2}O{sub eq} and that the pore solution method is a viable alternative ASR testing method. Massive amounts of ASR products (~20 μm thickness) were observed in concretes without SCMs consistent with high degree of expansion and extensive cracking. Small amounts of ASR products (≤5 μm thickness) were also observed in concrete with SCMs despite absence of significant expansion.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Stochastic multiscale flux basis for Stokes-Darcy flows

Three algorithms are developed for uncertainty quantification in modeling coupled Stokes and Darcy flows. The porous media may consist of multiple regions with different properties. The permeability is modeled as a non-stationary stochastic variable, with its log represented as a sum of local Karhunen-Loève (KL) expansions. The problem is approximated by stochastic collocation on either tensor-product or sparse grids, coupled with a multiscale mortar mixed finite element method for the spatial discretization. A non-overlapping domain decomposition algorithm reduces the global problem to a coarse scale mortar interface problem, which is solved by an iterative solver, for each stochastic realization. In the traditional domain decomposition implementation, each subdomain solves a local Dirichlet or Neumann problem in every interface iteration. To reduce this cost, two additional algorithms based on deterministic or stochastic multiscale flux basis are introduced. The basis consists of the local flux (or velocity trace) responses from each mortar degree of freedom. It is computed for each subdomain independently before the interface iteration begins. The use of the multiscale flux basis avoids the need for subdomain solves on each iteration. The deterministic basis is computed at each stochastic collocation and used only at this realization. The stochastic basis is formed by further looping over all local realizations of a subdomain's KL region before the stochastic collocation begins. It is used over multiple realizations. Numerical tests are presented to illustrate the performance of the three algorithms, with the stochastic multiscale flux basis showing significant savings in computational cost compared to the other two algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Meso-scale modelling of compressive fracture in concrete with irregularly shaped aggregates

This paper presents a meso-scale modelling framework to investigate the fracture process in concrete subjected to uniaxial and biaxial compression accounting for its mesostructural characteristics. 3D mesostructure of concrete consisting of coarse aggregates, mortar and interfacial transition zone between them was developed using an in-house code based on the Voronoi tessellation and splining method, which enables to generate the realistic-look aggregates with controllable structural features such as content, location, size and shape. Based on the generated 3D mesostructure, the concrete damage plasticity approach was employed to simulate the compressive fracture behaviour of concrete in terms of crack morphology and stress-strain response against the shape parameters of aggregate. Results indicate that the shape of aggregate has a negligible effect on compressive strength of concrete, which is highly associated with the random location and size distribution of aggregate. The aggregate irregularity has a significant influence on crack initiation and growth of concrete.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Optimizing dicyandiamide pretreatment conditions for enhanced structure and electronic properties of polymeric graphitic carbon nitride

Graphitic carbon nitride, a polymeric semiconductor, possesses a distinctive electronic band structure and exceptional chemical stability, making it a highly promising material for various catalytic applications such as electrocatalysis, photocatalysis, and photo-electrocatalysis. However, its practical applications remain limited due to its low active site density and poor electrical conductivity. In this study, to overcome such limitations, we have conducted a thorough investigation to explore the impact of dicyandiamide (DCDA) precursor pretreatment prior to thermal polymerization to graphitic carbon nitride. The DCDA precursor was subjected to various pretreatment methods including grinding using mortar and pestle, recrystallization through stirring or probe sonication after dissolving in deionized water, and freeze drying, prior to thermal polymerization at 550 °C for 2 hours. The structural and morphological properties of the catalysts prepared were compared and characterized by X-ray diffraction (XRD), attenuated total reflectance–Fourier transform infrared (ATR-IR) spectroscopy, scanning electron microscopy (SEM), transmission electron microscopy (TEM), and X-ray photoelectron spectroscopy (XPS) analysis while the electronic band gap properties were determined based on density functional theory (DFT) calculations for a set of crystalline systems having C/N ratios similar to those identified experimentally. In comparison to direct thermal polymerization, pretreated samples rendered the same product yield, diverse morphologies with flat or wrinkle structures, and reduced electrochemical resistance, making them suitable for use in various catalytic processes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Rapid Assessment of Sulfate Resistance in Mortar and Concrete

Extensive research has been conducted on the sulfate attack of concrete structures; however, the need to adopt the use of more sustainable materials is driving a need for a quicker test method to assess sulfate resistance. This work presents accelerated methods that can reduce the time required for assessing the sulfate resistance of mixtures by 70%. Class F fly ash has historically been used in concrete mixtures to improve sulfate resistance. However, environmental considerations and the evolving energy industry have decreased its availability, requiring the identification of economically viable and environmentally friendly alternatives to fly ash. Another challenge in addressing sulfate attack durability issues in concrete is that the standard sulfate attack test (ASTM C1012) is time-consuming and designed for only standard mortars (not concrete mixtures). To expedite the testing process, accelerated testing methods for both mortar and concrete mixtures were adopted from previous work to further the development of the accelerated tests and to assess the feasibility of testing the sulfate resistance of mortar and concrete mixtures rapidly. This study also established criteria for interpreting sulfate resistance for each of the test methods used in this work. A total of 14 mortar mixtures and four concrete mixtures using two types of Portland cement (Type I and Type I/II) and various supplementary cementitious materials (SCMs) were evaluated in this study. The accelerated testing methods significantly reduced the evaluation time from 12 months to 21 days for mortar mixtures and from 6 months to 56 days for concrete mixtures. The proposed interpretation method for mortar accelerated test results showed acceptable consistency with the ACI 318-19 interpretations for ASTM C1012 results. The interpretation methods proposed for the two concrete sulfate attack tests demonstrated excellent consistency with the ASTM C1012 results from mortar mixtures with the same cementitious materials combinations. Metakaolin was shown to improve sulfate resistance for both mortar and concrete mixtures, while silica fume and natural pozzolan had a limited impact. Using 15% metakaolin in mortar or concrete mixtures with Type I/II cement provided the best sulfate resistance.

Chemistry↗
DOE OSTI2021

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗
DOE OSTI2022

Sensitivity of the strength and toughness of concrete to the properties of the interfacial transition zone

We report Civil infrastructure is made primarily of concrete structures or components and therefore understanding durability and fracture behavior of concrete is of utmost importance. Concrete contains an interfacial transition zone (ITZ), a porous region surrounding the aggregates, that is often considered to be the weakest region in the concrete. The ITZ is poorly characterized and property estimates for the ITZ differ considerably. In this simulation study, representative concrete mesostructures are produced by packing coarse aggregates with realistic geometries into a mortar matrix. A meshless numerical method, peridynamics, is utilized to simulate the mechanical response including fracture under uniaxial compression and tension. The sensitivity of the stiffness and fracture toughness of the samples to the ITZ properties is computed, showing strong relationships between the ITZ properties and the effective modulus and effective yield strength of the concrete. These results provides insight into the influence of the poorly characterized ITZ on the stiffness and strength of concrete. This work showcases the applicability of peridynamics to concrete systems, matching experimental strength and modulus values. Additionally, relationships between the ITZ’s mechanical properties and the overall concrete strength and stiffness are presented to enable future design decisions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Interface Problem Formulation Improvements with Application to Nuclear Fuel Performance Analysis

The U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation Program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. This report documents new developments and robustness improvements in mechanical and thermal (gap heat transfer) contact formulations. The improvements range from the migration of industrial level (“assessment”) nuclear fuel model setups to the usage of mortar formulations, the addition of frictional contact to one-dimensional layered representations of fuel and cladding components, and the addition of the Petrov-Galerkin approach to dual mortar, which improves performance on curved, relatively coarse meshes. In addition, the Lagrange-multiplier enforcement of mechanical mortar contact constraints has been extended to two additional types of enforcement: penalty and augmented Lagrange-Uzawa. We show that the latter approach yields the same interface results as dual mortar in the Multiphysics Object-Oriented Simulation Environment, with the advantage of not worsening the condition number of the system matrix—thereby enabling the use of some general implementations of iterative preconditioners, at the expense of additional system evaluations (i.e., augmentations).

42 ENGINEERING↗
DOE OSTI2023

Residual estimation for grid modification in wall-modeled large eddy simulation using unstructured high-order methods

Here, the accuracy and computational cost of a large eddy simulation are highly dependent on the computational grid. Building optimal grids manually from a priori knowledge is not feasible in most practical use cases; instead, solution-adaptive strategies can provide a robust and cost-efficient method to generate a grid with the desired accuracy. We adapt the residual estimation algorithm developed by Toosi and Larsson for Discontinuous Galerkin Spectral Elements Methods (DGSEM) to guide the grid-adaptation process. The core of the method is the computation of the estimated modeling residual using the polynomial basis functions used in DGSEM, and the averaging of the estimated residual over each element. The final method is assessed in multiple channel flow test cases and for the transonic flow over an airfoil, in both cases making use of mortar interfaces between elements with hanging nodes. The method is found to be robust and reliable, and to provide solutions on grids with significantly fewer elements at comparable accuracy compared to when using human-generated grids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A mortar thermomechanical contact computational framework for nuclear fuel performance simulation

Nuclear fuel performance simulations involve the modeling of complex physical phenomena, ranging from fission gas release to fuel swelling and other temperature-induced effects. For light-water reactors (LWRs), swelling of the fuel and the pressure it imposes on the clad when they come into contact causes permanent clad deformation. Accurately characterizing the fuel-cladding interaction, which involves multiple physics, is essential to accurately simulate the fuel/cladding system. Thermomechanical modeling of this problem using a variationally consistent enforcement (e.g., a mortar approach) has been shown to improve the quality of results and facilitate convergence. Here, we present a general multiphysics computational framework for solving nuclear fuel problems using a mortar approach in BISON, a nuclear fuel performance code. In this study analyses show that using the mortar approach, which enables variationally consistent constraint enforcement, improves the quality of results as compared to the more commonly used node-on-face enforcement for representative LWR nuclear fuel simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗