Combination of Discrete and Finite Element Methods for Coupled Electrochemical-Mechanical Simulations of Lithium-Ion Battery Electrodes.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
Explore the source record for details and available documents.
Abstract not provided.
A new transport solver option based on the hybrid FEM (HFEM) was implemented in GRIFFIN, the MOOSE-based reactor analysis code, as an effort to support routine core design calculations for advanced reactor applications. The HFEM formulation with P{sub N} (spherical harmonics expansion), akin to the variational nodal method, is effective for solving a spatially homogenized problem with strong transport effect. The residual and Jacobian evaluations of the HFEM weak form were derived and successfully implemented in GRIFFIN, having the diffusion and the PN options available in the new HFEM based transport solver. The performance was tested with the simplified ABTR benchmark problems. The results indicate that the HFEM-based transport solver is a feasible option for solving problems with spatially homogenized and strong streaming by providing superior accuracy with a proper p-refinement. (authors)
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The finite element method (FEM) is widely used to simulate a variety of physics phenomena. Approaches that integrate FEM with neural networks (NNs) are typically leveraged as an alternative to conducting expensive FEM simulations in order to reduce the computational cost without significantly sacrificing accuracy. However, these methods can produce biased predictions that deviate from those obtained with FEM, since these hybrid FEM-NN approaches rely on approximations trained using physically relevant quantities. In this work, an uncertainty estimation framework is introduced that leverages ensembles of Bayesian neural networks to produce diverse sets of predictions using a hybrid FEM-NN approach that approximates internal forces on a deforming solid body. The uncertainty estimator developed herein reliably infers upper bounds of bias/variance in the predictions for a wide range of interpolation and extrapolation cases using a three-element FEM-NN model of a bar undergoing plastic deformation. This proposed framework offers a powerful tool for assessing the reliability of physics-based surrogate models by establishing uncertainty estimates for predictions spanning a wide range of possible load cases.
Abstract not provided
This research focuses on electromagnetic interference (EMI) / electromagnetic compatibility (EMC) design and analysis in power electronics systems. To limit the EMI under the standards, different methods and strategies are investigated. Parasitic parameters of high frequency (HF) transformer are analyzed using a novel analytical method, finite element method (FEM), and experimental measurements for different structures and windings arrangements. Also, the magnetic field, electric field, electric displacement, and electric potential distribution are simulated and analyzed. Moreover, a high voltage system is considered and analyzed to improve the EMC. The EMI propagation paths are analyzed. The EMI noise level of the system is obtained and compared to the IEC61800-3 standard. To improve the EMC, the parasitic parameters of the transformer, as the main path of EMI circulation, are analyzed and optimized to block the propagation. Furthermore, the geometry structure of the HF transformer is optimized to lower the parasitics in the system. Three pareto-optimal techniques are investigated for the optimization. The models and results are verified by 3D-FEM and experimental results for several given scenarios. Furthermore, the EMC modeling and conducted EMI analysis are developed for a system including an AC-DC-DC power supply (rectifier and dual active bridge (DAB) converter). Moreover, the common mode (CM) EMI noise propagation through the system is discussed and the noise sources and effect of components on the noise are analyzed. Additionally, the CM impedance of different parts of the system and the noise levels are discussed. Finally, EMI attenuation techniques were applied to the system.
In this study, two modifications are introduced for improving the accuracy, versatility, and robustness of a class of hybrid methods for radiation transport. In general, such methods are constructed by splitting the radiative flux into collided and uncollided components to which low- and high-resolution angular approximations are applied, respectively. In this work we focus on discrete ordinates discretizations of high and low order. The first modification we introduce changes the way in which the collided component is mapped into the uncollided component at the end of each time step in a simulation. The new mapping is a Nyström-type reconstruction that is applicable to arbitrary discrete ordinates quadratures, is guaranteed to preserve positivity of the solution provided that all ordinate weights are positive, is significantly more accurate than previous methods, and can be readily extended to other discretizations such as moment methods, finite element methods, and diffusion approximations. The second modification leverages integral deferred correction (IDC) to iteratively correct for the splitting error introduced by the inconsistency in angular discretization between the collided and uncollided components, in addition to improving the accuracy of the low-order temporal error that is treated by traditional IDC methods. Numerical tests in one- and two-dimensional geometries are used to demonstrate the increased accuracy and efficiency of the proposed modifications. It is found that the two techniques combined yield methods with solution accuracy and memory requirements comparable to that of monolithic discrete ordinates methods while reducing runtime by as much as a factor of between two and ten, depending on the problem.
High-frequency (HF) transformers are an essential part of many power electronic devices. The performance and behavior of HF transformers can greatly affect the efficiency and performance of all systems, particularly, from a parasitic parameter point of view. In this article, HF transformers' parasitic parameters, such as leakage inductances and parasitic capacitances, are analyzed using a novel analytical method, finite element method (FEM), and experimental measurements of different structures and winding arrangements. Also, the magnetic field, electric field, electric displacement field, and electric potential distribution within the transformers are simulated and analyzed. Four different HF transformers with E and U cores with different windings are designed and analyzed. Investigation outcomes help to classify structures according to the trade-off between leakage inductances and series parasitic capacitances. This information can later be used for the optimal selection and design of transformers as a function of their operating frequency for any power rating and voltage level. Moreover, 3-D FEM and experimental results validate the proposed methodology to be used for designing HF transformers in high-voltage/power applications.
From its inception, the combined finite discrete element method has used a distributed potential contact force algorithm to resolve interaction between finite elements. The contact interaction algorithm relies on evaluation of the contact force potential field. The problem with existing algorithms is that the potential field introduces artificial numerical non-smoothness in the contact force. This work introduces a smooth potential field based on the finite element topology, and a generalized contact interaction law is constructed on top of the smooth potential field. Further, a number of validation cases for the proposed algorithm, considering different shapes of discrete elements, are presented, and detailed aspects of the proposed contact interaction law are tested with numerical examples.
Numerical simulation of nonlinear elastic wave propagation in solids with cracks is indispensable for decoding the complicated mechanisms associated with the nonlinear ultrasonic techniques in NonDestructive Testing (NDT). Here, we introduce a two-dimensional implementation of the combined finitediscrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), to explicitly simulate the crack induced nonlinear elasticity in solids with both horizontal and inclined cracks. In the FDEM model, the solid is discretized into finite elements to capture the wave propagation in the bulk material, and the finite elements along the two sides of the crack also behave as discrete elements to track the normal and tangential interactions between crack surfaces. The simulation results show that for cracked models, nonlinear elasticity is generated only when the excitation amplitude is large enough to trigger the contact between crack surfaces, and the nonlinear behavior is very sensitive to the crack surface contact. The simulations reveal the influence of normal and tangential contact on the nonlinear elasticity generation. Moreover, the results demonstrate the capabilities of FDEM for decoding the causality of nonlinear elasticity in cracked solid and its potential to assist in Non-Destructive Testing (NDT).
The purpose of this paper is to construct a new class of discrete generalized Korn’s inequalities for piecewise $H^1$ vector fields and piecewise $H^2$ vector fields in three-dimensional space. The resulting Korn’s inequalities are different from the standard Korn’s inequalities, as they involve the trace-free symmetric gradient operator, in place of the usual symmetric gradient operator. Furthermore, it is anticipated that the new generalized Korn’s inequalities will be useful for the analysis of a broad range of finite element methods, including mixed finite element methods and discontinuous Galerkin methods.
Two codes were tested for the first order neutron transport equation using finite element methods. The un-collided-flux solution is used as a preconditioner for each of these methods. These codes include a least squares finite element method and a discontinuous finite element method. The performance of each code is shown on problems in one and two dimensions. The un-collided-flux preconditioner shows good speedup on each of the given methods. The un-collided-flux preconditioner has been used on the second-order equation, and here we extend those results to the first order equation. (authors)
Here in this paper, two-level methods for solving transport problems in one-dimensional slab geometry based on the quasi-diffusion (QD) method are developed. A linear discontinuous finite element method (LDFEM) is derived for the spatial discretization of the low-order QD (LOQD) equations. It involves special interface conditions at the cell edges based on the idea of QD boundary conditions (BCs). We consider different kinds of QD BCs to formulate the necessary cell-interface conditions. We develop two-level methods with independent discretization of the high-order transport equation and LOQD equations, where the transport equation is discretized using the method of characteristics and the LDFEM is applied to the LOQD equations. We also formulate closures that lead to the discretization consistent with a LDFEM discretization of the transport equation. The proposed methods are studied by means of test problems formulated with the method of manufactured solutions. Numerical experiments are presented demonstrating the performance of the proposed methods. Lastly, we also show that the method with independent discretization has the asymptotic diffusion limit.
A simple and effective decoupled finite element analysis method was developed for simulating both the piezoelectric and flexoelectric effects of zinc oxide (ZnO) and barium titanate (BTO) nanowires (NWs). The piezoelectric potential distribution on a ZnO NW was calculated under three deformation conditions (cantilever, three-point, and four-point bending) and compared to the conventional fully coupled method. The discrepancies of the electric potential maximums from these two methods were found very small, validating the accuracy and effectiveness of the decoupled method. Both ZnO and BTO NWs yielded very similar potential distributions. Comparing the potential distributions induced by the piezoelectric and flexoelectric effects, we identified that the middle segment of a four-point bending NW beam is the ideal place for measuring the flexoelectric coefficient, because the uniform parallel plate capacitor-like potential distribution in this region is exclusively induced by the flexoelectric effect. This decoupled method could provide a valuable guideline for experimental measurements of the piezoelectric effects and flexoelectric effects in the nanometer scale.
This paper is concerned with a theoretical and numerical study for the initial-boundary value problem for a linear hyperbolic equation with variable coefficient and acoustic boundary conditions. On the theoretical results, we prove the existence and uniqueness of global solutions, and the uniform stability of the total energy. Numerical simulations using the finite element method associated with the finite difference method are employed, for one-dimensional and two-dimensional cases, to validate the theoretical results. In addition, numerically the uniform decay rate for energy and the order of convergence of the approximate solution are also shown.