Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary element 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 1,045 records · Page 58

Displacement Models for THUNDER Actuators having General Loads and Boundary Conditions

This paper summarizes techniques for quantifying the displacements generated in THUNDER actuators in response to applied voltages for a variety of boundary conditions and exogenous loads. The PDE (partial differential equations) models for the actuators are constructed in two steps. In the first, previously developed theory quantifying thermal and electrostatic strains is employed to model the actuator shapes which result from the manufacturing process and subsequent repoling. Newtonian principles are then employed to develop PDE models which quantify displacements in the actuator due to voltage inputs to the piezoceramic patch. For this analysis, drive levels are assumed to be moderate so that linear piezoelectric relations can be employed. Finite element methods for discretizing the models are developed and the performance of the discretized models are illustrated through comparison with experimental data.

Wieman, Robert↗

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗

Explicit modeling of composite plates and beams in the dynamics of multibody systems

The state of the art dynamic response analysis of flexible multibody systems is currently restricted to elastic bodies with homogeneous materials. The requirements for high speed operation has made it necessary to use lightweight multi layered composite bodies in robotic systems and space structure applications. Dynamic modeling and analysis of such systems are particularly important since the effects of body flexibility to the performance are likely to be more pronounced. The eight-noded isoperimetric quadrilateral element with independent rotational and displacement degrees of freedom is extended to laminated composite elements. The element includes an arbitrary number of bonded layers, each of which may have a different thickness. The transverse shear deformation which is a predominant factor in the analysis of laminated composite structures is taken into account in developing the stiffness and mass matrices. The corresponding 3-D mode shapes are then incorporated to the multibody system dynamical equations. Floating body reference frames allow the selection of different boundary conditions, and the dynamical equations contain all the nonlinear interactions between the rigid and elastic motion. Example simulations are presented to illustrate the methods proposed.

Amirouche, F. M. L.↗

Design Exploration of a Transonic Cruise Slotted Airfoil

A knowledge-based aerodynamic design method, CDISC, has been leveraged for the computational design exploration of the cruise slotted airfoil as a drag-saving technology for single aisle, transonic transport aircraft. Aerodynamic predictions were generated using the NASA USM3D Reynolds-averaged Navier-Stokes flow solver, and laminar flow assessments were conducted using the NASA BLSTA3D boundary layer profile solver paired with the LASTRAC stability analysis and transition prediction software. The CDISC design method was used to parametrically vary several slotted-airfoil design variables to understand their impact on aerodynamic performance. In addition to establishing best practices for slotted airfoil design, it was observed that the historically noted skin-friction penalty for cruise slotted airfoils was attributed to the low-Reynolds-number boundary layer of the flap. At cruise, this skin-friction penalty negated any potential decrease in pressure drag enabled by the cruise slotted airfoil architecture. This observation led to the development of an Aft-Laminar Multi-element Airfoil concept that uses airfoil shaping to promote natural laminar flow on the flap. At cruise conditions, the concept is predicted to offer a 2.7% reduction in sectional drag relative to a fully turbulent supercritical airfoil. Off-design analyses showed that laminar flow could be maintained with near-cruise variations in angle of attack and Mach number, resulting in sustained aerodynamic efficiency improvements and a delay in the drag rise Mach number. A low-speed, high-lift analysis at takeoff conditions predicted a maximum lift coefficient of 2.4 could be achieved by use of a variable-camber leading edge on the main element and simple deflection of the flap. These computationally-predicted benefits have motivated future research toward an aft-laminar cruise slotted wing design as an incremental step toward the application of natural laminar flow technology on transonic commercial transports.

Brett R Hiller↗

Comparison of chemical lateral boundary conditions for air quality predictions over the contiguous United States during pollutant intrusion events

The National Air Quality Forecast Capability (NAQFC) operated in U.S.’s National Oceanic and Atmospheric Administration (NOAA) provides the operational forecast guidance for ozone and fine particulate matters with aerodynamic diameters less than 2.5μm (PM2.5) over the contiguous 48 U.S. states (CONUS) using the Community Multi-scale Air Quality (CMAQ) model. The existing NAQFC uses climatological chemical lateral boundary conditions ( CLBCs), which cannot capture pollutant intrusion events originating outside of the model domain. In this study, we developed a model framework to use dynamic CLBCs from the Goddard Earth Observing System Model, version 5 (GEOS) to drive NAQFC. A mapping of the GEOS chemical species to the CMAQ’s Carbon Bond 5 (CB05)-Aero6 species was developed. The utilization of the GEOS dynamic CLBCs in NAQFC showed the best overall performance in simulating the surface observations during the Saharan dust intrusion and Canadian wildfire events in summer 2015.The simulated PM2.5 was improved from 0.18 to 0.37 and the mean bias was reduced from -6.74 μg/m3 to -2.96 μg/m3 over CONUS. Although the effect of CLBCs on the PM2.5 correlation was mainly near the inflow boundary, its impact on the background concentrations reached further inside the domain. The CLBCs could affect background ozone concentrations through the inflows of ozone itself and its precursors, such as CO. It was further found that the aerosol optical thickness (AOT) from satellite retrievals correlated well with the column CO and elemental carbon from GEOS. The satellite-derived AOT CLBCs generally improved the model performance for the wildfire intrusion events during a summer 2018 case study, and demonstrated how satellite observations of atmospheric composition could be used as an alternative method to capture the air quality effects of intrusions when the global model CLBCs, such as GEOS CLBCs, are not available.

Community Multi-scale Air Quality (CMAQ) model↗

On The Incorporation of Both Function-Driven Requirements and Topology Optimization in The Development of Lunar Launch & Landing Pads

In the design of additively manufactured structures, there can be a constant battle between focusing on function-driven design and topology optimization in an attempt to improve a products ability to meet all requirements. Function-driven optimization is defined as the process of identifying the best design parameters that satisfy project functionality requirements. Function-driven optimization is typically implemented by a designer(s) who, instead of computational limitations, has human bias for the input or inputs determined. On the other hand, topology optimization is defined as a mathematical method that optimizes material layout within a given design space, for a given set of loads, boundary conditions, and constraints with the goal of maximizing the performance of the system. In topology optimization, the work is defined by the limited mathematical inputs and computationally available goals/restraints. If not evaluated simultaneously, these two approaches can result in fundamentally different design concepts meant to address the same requirements. NASA’s Artemis program includes the development of a base camp near the lunar South Pole. After site selection is made, one of the first infrastructure elements needed will be a launch and landing pad (LLP). The construction of such a pad will prevent damage to nearby structures from plume-surface interaction (PSI) while also providing a stable platform for the lander. These benefits reduce the required distance between the established landing zones and other settlement infrastructure. In this paper, the authors evaluate the requirements associated with lunar LLP development, application of both functional and topological optimization techniques, methodologies to combine the two techniques, and potential impacts on final design of such a pad.

In-situ↗

Grid generation by elliptic partial differential equations for a tri-element Augmentor-Wing airfoil

Two efforts to numerically simulate the flow about the Augmentor-Wing airfoil in the cruise configuration using the GRAPE elliptic partial differential equation grid generator algorithm are discussed. The Augmentor-Wing consists of a main airfoil with a slotted trailing edge for blowing and two smaller airfoils shrouding the blowing jet. The airfoil and the algorithm are described, and the application of GRAPE to an unsteady viscous flow simulation and a transonic full-potential approach is considered. The procedure involves dividing a complicated flow region into an arbitrary number of zones and ensuring continuity of grid lines, their slopes, and their point distributions across the zonal boundaries. The method for distributing the body-surface grid points is discussed.

Sorenson, R. L.↗

Transfinite element methodology for nonlinear/linear transient thermal modelling/analysis - Progress and recent advances

The 'transfinite element method' (TFEM) proposed by Tamma and Railkar (1987 and 1988) for the analysis of linear and nonlinear heat-transfer problems is described and demonstrated. The TFEM combines classical Galerkin and transform approaches with state-of-the-art FEMs to obtain a flexible hybrid modeling scheme. The fundamental principles of the TFEM and the derivation of the governing equations are reviewed, and numerical results for sample problems are presented in extensive graphs and briefly characterized. Problems analyzed include a square plate with a hole, a rectangular plate with natural and essential boundary conditions and varying thermal conductivity, the Space Shuttle thermal protection system, a bimaterial plate subjected to step temperature variations, and solidification in a semiinfinite liquid slab.

Tamma, Kumar K.↗

Nonlinear mesomechanics of composites with periodic microstructure

This work is concerned with modeling the mechanical deformation or constitutive behavior of composites comprised of a periodic microstructure under small displacement conditions at elevated temperature. A mesomechanics approach is adopted which relates the microimechanical behavior of the heterogeneous composite with its in-service macroscopic behavior. Two different methods, one based on a Fourier series approach and the other on a Green's function approach, are used in modeling the micromechanical behavior of the composite material. Although the constitutive formulations are based on a micromechanical approach, it should be stressed that the resulting equations are volume averaged to produce overall effective constitutive relations which relate the bulk, volume averaged, stress increment to the bulk, volume averaged, strain increment. As such, they are macromodels which can be used directly in nonlinear finite element programs such as MARC, ANSYS and ABAQUS or in boundary element programs such as BEST3D. In developing the volume averaged or efective macromodels from the micromechanical models, both approaches will require the evaluation of volume integrals containing the spatially varying strain distributions throughout the composite material. By assuming that the strain distributions are spatially constant within each constituent phase-or within a given subvolume within each constituent phase-of the composite material, the volume integrals can be obtained in closed form. This simplified micromodel can then be volume averaged to obtain an effective macromodel suitable for use in the MARC, ANSYS and ABAQUS nonlinear finite element programs via user constitutive subroutines such as HYPELA and CMUSER. This effective macromodel can be used in a nonlinear finite element structural analysis to obtain the strain-temperature history at those points in the structure where thermomechanical cracking and damage are expected to occur, the so called damage critical points of the structure.

Walker, Kevin P.↗

Crack-mouth displacements for semielliptical surface cracks subjected to remote tension and bending loads

The exact analytical solution for an embedded elliptical crack in an infinite body subjected to arbitrary loading was used in conjunction with the finite element alternating method to obtain crack-mouth-opening displacements (CMOD) for surface cracks in finite plates subjected to remote tension. Identical surface-crack configurations were also analyzed with the finite element method using 20-noded element for plates subjected to both remote tension and bending. The CMODs from these two methods generally agreed within a few percent of each other. Comparisons made with experimental results obtained from surface cracks in welded aluminum alloy specimens subjected to tension also showed good agreement. Empirical equations were developed for CMOD for a wide range of surface-crack shapes and sizes subjected to tension and bending loads. These equations were obtained by modifying the Green-Sneddon exact solution for an elliptical crack in an infinite body to account for finite boundary effects. These equations should be useful in monitoring surface-crack growth in tests and in developing complete crack-face-displacement equations for use in three-dimensional weight-function methods.

Raju, Ivatury S.↗

The buckling response of symmetrically laminated composite plates having a trapezoidal planform area

The focus of this work is the buckling response of symmetrically laminated composite plates having a planform area in the shape of an isosceles trapezoid. The loading is assumed to be inplane and applied perpendicular to the parallel ends of the plate. The tapered edges of the plate are assumed to have simply supported boundary conditions, while the parallel ends are assumed to have either simply supported or clamped boundary conditions. A semi-analytic closed-form solution based on energy principles and the Trefftz stability criterion is derived and solutions are obtained using the Rayleigh-Ritz method. Intrinsic in this solution is a simplified prebuckling analysis which approximates the inplane force resultant distributions by the forms Nx=P/W(x) and Ny=Nxy=0, where P is the applied load and W(x) is the plate width which, for the trapezoidal planform, varies linearly with the lengthwise coordinate x. The out-of-plane displacement is approximated by a double trigonometric series. This analysis is posed in terms of four nondimensional parameters representing orthotropic and anisotropic material properties, and two nondimensional parameters representing geometric properties. For comparison purposes, a number of specific plate geometry, ply orientation, and stacking sequence combinations are investigated using the general purpose finite element code ABAQUS. Comparison of buckling coefficients calculated using the semi-analytical model and the finite element model show agreement within 5 percent, in general, and within 15 percent for the worst cases. In order to verify both the finite element and semi-analytical analyses, buckling loads are measured for graphite/epoxy plates having a wide range of plate geometries and stacking sequences. Test fixtures, instrumentation system, and experimental technique are described. Experimental results for the buckling load, the buckled mode shape, and the prebuckling plate stiffness are presented and show good agreement with the analytical results regarding the buckling load and the prebuckling plate stiffness. However, the experimental results show that for some cases the analysis underpredicts the number of halfwaves in the buckled mode shape. In the context of the definitions of taper ratio and aspect ratio used in this study, it is concluded that the buckling load always increases as taper ratio increases for a given aspect ratio for plates having simply supported boundary conditions on the parallel ends. There are combinations of plate geometry and ply stackling sequences, however, that reverse this trend for plates having clamped boundary conditions on the parallel ends such that an increase in the taper ratio causes a decrease in the buckling load. The clamped boundary conditions on the parallel ends of the plate are shown to increase the buckling load compared to simply supported boundary conditions. Also, anisotropy (the D16 and D26 terms) is shown to decrease the buckling load and skew the buckled mode shape for both the simply supported and clamped boundary conditions.

Radloff, H. D., II↗

Identifying Challenges in Safeguards for Metallic Fuel Fabrication Facilities

As new advanced reactors gain popularity, there is an increasing interest in metallic fuel fabrication for fast reactors. While metallic fuels themselves are not a new idea, as many of the first reactors employed metallic fuels, new designs, compositions, and fabrication methods are appearing throughout the nuclear community. As the interest grows and facilities are constructed, both domestic and international safeguards will need to be heavily involved to support safeguards-by-design (SBD) measures from the start. This work compiles a review of historical and modern fuel types and fabrication methods, fabrication processes, safeguards gaps, and potential safeguards solutions. Metallic nuclear fuel types have been around for many decades and were included in some of the first reactors including the Experimental Breeder Reactor (EBR)-I and -II, the Fermi 1 reactor, the Integral Fast Reactor (IFR), and the Dounreay Fast Reactor (DFR). These reactors used various compositions including pure uranium (U) metal, U-zirconium (Zr) alloys, plutonium (Pu)-aluminum (Al) alloys, U-fissium (Fs) alloys, U-Pu-Zr alloys, and U-molybdenum (Mo) alloys [1, 2, 3, 4, 5]. These small alloying additions are included to improve the material properties of the pure U metal. The alpha-phase U (stable below 661C) suffers elongation in one direction causing grain boundary cracking and increasing creep rate due to irradiation growth, thermal cycling, and preferential crystal orientation. It is ideal to utilize the gamma-phase U (typically stable above 769C) by adding small amounts of alloying elements such as Zr or Mo to stabilize this phase down to room temperature [3]. Additionally, some research has been focused on U with transuranic (TRU) elements present, typically coming from the used fuel recycling process. Including these elements in fast reactor fuel can aid in the reduction of nuclear waste by burning minor long-lived actinides. However, the additions of TRU elements can cause concerns to arise when trying to fabrication or safeguard metallic fuels. A typical metallic fuel element is shown in Figure 1. Sodium is added into the cladding to create a thermal bond between the fuel slug and cladding wall. The fuel slug is then inserted and the end plug is welded on to the top of the fuel element. A gas plenum is left to create a headspace for gaseous fission products to escape rather than continue to build in the fuel itself [1, 5]. Other fuel element geometries exist as well, such as the Lightbridge twisted cruciform geometry shown in Figure 2 [6]. This design allows for better cooling performance and provides room for fuel rod swelling without impacting the fuel rod diameter. There are many different fabrication methods for metallic fuels, which is one of the many benefits of these fuel types. Many of these fabrication methods are relatively easy and cost-efficient. The most popular fabrication method is injection casting, sometimes called vacuum induction melting (VIM), shown in Figure 3 [4, 8, 9, 7, 10]. This method was largely used for EBR-II fuel fabrication. The injection casting system is contained inside of a vessel consisting of a Y2O3-coated graphite crucible surrounded by an induction coil with ZrO2-coated quartz molds suspended above the crucible. The fuel feedstock is placed inside of the graphite crucible and melted using the induction furnace. The induction furnace utilizes a dual frequency with the high frequency melting the feedstock and the low frequency causing stirring of the melted feedstock to form a homogeneous mixture. The mixture is heated to approximately 1600C in an argon environment. The vessel is evacuated and then the quartz molds are lowered into the graphite crucible containing the molten metal and the vessel is repressurized to inject the metal fuel upwards into the molds. The molds are removed and then shattered to release the fuel slugs. This fabrication method was used to fabricate 39,000 metallic fuel pins for EBR-II. While injection casting has been the most common metallic fuel fabrication method throughout the decades, many other methods have been explored including low-pressure gravity casting, microwave casting, continuous casting, centrifugal casting, coextrusion, and many others [11, 12, 8, 13, 14, 15]. Some of these methods aim to mitigate challenges that arise with americium (Am) volatilization during the casting process for TRU-containing fuel feedstocks, an issue with injection casting. Coextrusion is one of the methods explored at the Idaho National Laboratory (INL) and has been utilized for the initial fabrication tests of Lightbridge's unique fuels, as well as other metallic fuels with cladding coextruded. In this process, large billets are formed and machined and then inserted into a molten salt bath for approximately 30 minutes. The billets are then loaded into the extrusion press and extruded. This process can be seen in Figure 4 [15].

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unstructured and adaptive mesh generation for high Reynolds number viscous flows

A method for generating and adaptively refining a highly stretched unstructured mesh suitable for the computation of high-Reynolds-number viscous flows about arbitrary two-dimensional geometries was developed. The method is based on the Delaunay triangulation of a predetermined set of points and employs a local mapping in order to achieve the high stretching rates required in the boundary-layer and wake regions. The initial mesh-point distribution is determined in a geometry-adaptive manner which clusters points in regions of high curvature and sharp corners. Adaptive mesh refinement is achieved by adding new points in regions of large flow gradients, and locally retriangulating; thus, obviating the need for global mesh regeneration. Initial and adapted meshes about complex multi-element airfoil geometries are shown and compressible flow solutions are computed on these meshes.

Mavriplis, Dimitri J.↗

Unstructured and adaptive mesh generation for high Reynolds number viscous flows

A method for generating and adaptively refining a highly stretched unstructured mesh suitable for the computation of high-Reynolds-number viscous flows about arbitrary two-dimensional geometries was developed. The method is based on the Delaunay triangulation of a predetermined set of points and employs a local mapping in order to achieve the high stretching rates required in the boundary-layer and wake regions. The initial mesh-point distribution is determined in a geometry-adaptive manner which clusters points in regions of high curvature and sharp corners. Adaptive mesh refinement is achieved by adding new points in regions of large flow gradients, and locally retriangulating; thus, obviating the need for global mesh regeneration. Initial and adapted meshes about complex multi-element airfoil geometries are shown and compressible flow solutions are computed on these meshes.

Mavriplis, D. J.↗

Half-closed discontinuous Galerkin discretisations

Here we introduce the concept of half-closed nodes for nodal discontinuous Galerkin (DG) discretisations. Unlike more commonly used closed nodes in DG, where on every element nodes are placed on all of its boundaries, half-closed nodes only require nodes to be placed on a subset of the element's boundaries. The effect of using different nodes on DG operator sparsity is studied and we find in particular for there to be no difference in the sparsity pattern of the Laplace operator whether closed or half-closed nodes are used. On quadrilateral/hexahedral elements we use the Gauss-Radau points as the half-closed nodes of choice, which we demonstrate is able to speed up DG operator assembly in addition to leverage previously known superconvergence results. We also discuss in this work some linear solver techniques commonly used for Finite Element or discontinuous Galerkin methods such as static condensation and block-based methods, and how they can be applied to half-closed DG discretisations.

97 MATHEMATICS AND COMPUTING↗

The estimation of Poisson’s ratio by time-averaging and Cornu’s method for isotropic beams

Prior research suggests that direct static (e.g., uniaxial testing) and dynamic (e.g., ultrasonic wave propagation analysis) measurements differ in the estimation of Poisson’s ratio because anisotropies and heterogeneities in the sample material affect the two types of tests differently. Even assuming isotropic and homogeneous material properties, prior research further suggests that discrepancies between static/dynamic test results will exist because the error of the diagnostic techniques for the measurand are inherently different. Finally, thermodynamic effects are not present in static tests but can significantly affect dynamic test results. Given the potential for all these variables to produce discrepancies, it would be helpful to have the measurement of Poisson’s ratio obtainable from the same theory and experimental measurements by either static or dynamic testing methods. Our finite element calculations show that by combining time-averaged scanning digital holography with Cornu’s method, it is theoretically possible to estimate the effective Poisson’s ratio from the anticlastic contours at the antinode of the first out-of-plane bending mode shape. This is true regardless of frequency and therefore applicable for both static and dynamic measurements. Our results show that the estimate of Poisson’s ratio by Cornu’s method using data from simulations of mode shapes approaches the true value of Poisson’s ratio. Additionally, our research suggests that beam geometry and boundary conditions are fundamental factors limiting the convergence of the estimate of Poisson’s ratio to the true value of Poisson’s ratio regardless of performing a static or dynamic test.

42 ENGINEERING↗

Meshfree Methods for Modeling Chemo-Mechanical Cathode Cracking in Li-Ion Batteries

Through repeated charging and discharging, the electrodes of a Li-ion battery experience cyclic loading causing swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries. Chemo-mechanical models are simulated mainly by two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a chemo-mechanical crack propagation model that achieves both intermediate accuracy and computational time between a CDM and CZM. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

cathode cracking↗