Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite difference 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,207 records · Page 67

On convergent finite difference schemes for variational–PDE-based image processing

We study an adaptive anisotropic Huber functional-based image restoration scheme. Using a combination of L2–L1 regularization functions, an adaptive Huber functional-based energy minimization model provides denoising with edge preservation in noisy digital images. We study a convergent finite difference scheme based on continuous piecewise linear functions and use a variable splitting scheme, namely the Split Bregman (In: Goldstein and Osher, SIAM J Imaging Sci 2(2):323–343, 2009) algorithm, to obtain the discrete minimizer. Experimental results are given in image denoising and comparison with additive operator splitting, dual fixed point, and projected gradient schemes illustrates that the best convergence rates are obtained for our algorithm.

97 MATHEMATICS AND COMPUTING↗

Symmetric spectrum current drive due to finite radial drift effects

We report the drift surfaces of minority heated ions differ from flux surfaces due to finite poloidal gyroradius effects. As the minority poloidal gyroradius approaches radial scale lengths in the plasma, the difference between drift and flux surfaces can modify the heating and lead to a symmetric spectrum minority counter-current being driven. In response, a corresponding overall net co-current of comparable size is driven. This beneficial symmetric spectrum current drive in a tokamak is due to the parallel velocity asymmetry in the drift departure from a flux surface. As this new source of driven current is a side effect of minority heating it comes without any additional economic cost to reactor power balance. The symmetric spectrum current driven for near Maxwellian minorities is evaluated by an adjoint method and found to be modest. However, minority heating typically results in strong non-Maxwellian features on minority distributions so it may be possible to drive a significantly larger co-current. A related evaluation is performed for alpha particles in a deuterium minority heated plasma with a tritium majority. The low density of the alphas tends to keep this driven symmetric spectrum current small, but at very high heating levels a significant co-current might be driven. Other mechanisms to drive co-current with a symmetric spectrum are discussed and estimated, including asymmetric electron drag and focusing of the applied minority heating radio frequency fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Seismic waves in medium with poroelastic/elastic interfaces: a two-dimensional P-SV finite-difference modelling

SUMMARY We present a new methodology of the finite-difference (FD) modelling of seismic wave propagation in a strongly heterogeneous medium composed of poroelastic (P) and (strictly) elastic (E) parts. The medium can include P/P, P/E and E/E material interfaces of arbitrary shapes. The poroelastic part can be with (i) zero resistive friction, (ii) non-zero constant resistive friction or (iii) JKD model of the frequency-dependent permeability and resistive friction. Our FD scheme is capable of subcell resolution: a material interface can have an arbitrary position in the spatial grid. The scheme keeps computational efficiency of the scheme for a smoothly and weakly heterogeneous medium (medium without material interfaces). Numerical tests against independent analytical, semi-analytical and spectral-element methods prove the efficiency and accuracy of our FD modelling. In numerical examples, we indicate effect of the P/E interfaces for the poroelastic medium with a constant resistive friction and medium with the JKD model of the frequency-dependent permeability and resistive friction. We address the 2-D P-SV problem. The approach can be readily extended to the 3-D problem.

58 GEOSCIENCES↗

Symmetric spectrum current drive due to finite radial drift effects

The drift surfaces of minority heated ions differ from flux surfaces due to finite poloidal gyroradius effects. As the minority poloidal gyroradius approaches radial scale lengths in the plasma the difference between drift and flux surfaces can modify the heating and lead to a symmetric spectrum minority counter-current being driven. In response, a corresponding overall net co-current of comparable size is driven. This beneficial symmetric spectrum current drive in a tokamak is due to the parallel velocity asymmetry in the drift departure from a flux surface. As this new source of driven current is a side effect of minority heating it comes without any additional economic cost to reactor power balance. The symmetric spectrum current driven for near Maxwellian minorities is evaluated by an adjoint method and found to be modest. However, minority heating typically results in strong non-Maxwellian features on minority distributions so it may be possible to drive significantly larger co-current. A related evaluation is performed for alpha particles in a deuterium minority heated plasma with a tritium majority. The low density of the alphas tends to keep this driven symmetric spectrum current small, but at very high heating levels a significant co-current might be driven. Other mechanisms to drive co-current with a symmetric spectrum are discussed and estimated, including asymmetric electron drag and focusing of the applied minority heating radio frequency fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterize traction–separation relation and interfacial imperfections by data-driven machine learning models

Abstract Interfacial mechanical properties are important in composite materials and their applications, including vehicle structures, soft robotics, and aerospace. Determination of traction–separation (T–S) relations at interfaces in composites can lead to evaluations of structural reliability, mechanical robustness, and failures criteria. Accurate measurements on T–S relations remain challenging, since the interface interaction generally happens at microscale. With the emergence of machine learning (ML), data-driven model becomes an efficient method to predict the interfacial behaviors of composite materials and establish their mechanical models. Here, we combine ML, finite element analysis (FEA), and empirical experiments to develop data-driven models that characterize interfacial mechanical properties precisely. Specifically, eXtreme Gradient Boosting (XGBoost) multi-output regressions and classifier models are harnessed to investigate T–S relations and identify the imperfection locations at interface, respectively. The ML models are trained by macroscale force–displacement curves, which can be obtained from FEA and standard mechanical tests. The results show accurate predictions of T–S relations ( R 2 = 0.988) and identification of imperfection locations with 81% accuracy. Our models are experimentally validated by 3D printed double cantilever beam specimens from different materials. Furthermore, we provide a code package containing trained ML models, allowing other researchers to establish T–S relations for different material interfaces.

97 MATHEMATICS AND COMPUTING↗

Laser-based ultrasound interrogation of surface and sub-surface features in advanced manufacturing materials

Abstract Structures formed by advanced manufacturing methods increasingly require nondestructive characterization to enable efficient fabrication and to ensure performance targets are met. This is especially important for aerospace, military, and high precision applications. Surface acoustic waves (SAW) generated by laser-based ultrasound can detect surface and sub-surface defects relevant for a broad range of advanced manufacturing processes, including laser powder bed fusion (LPBF). In particular, an all-optical SAW generation and detection configuration can effectively interrogate laser melt lines. Here we report on scattered acoustic energy from melt lines, voids, and surface features. Sub-surface voids are also characterized using X-ray Computed Tomography (CT). High resolution CT results are presented and compared with SAW measurements. Finite difference simulations inform experimental measurements and analysis.

36 MATERIALS SCIENCE↗

Non-intrusive reduced order modeling of natural convection in porous media using convolutional autoencoders: Comparison with linear subspace techniques

Natural convection in porous media is a highly nonlinear multiphysical problem relevant to many engineering applications (e.g., the process of CO 2 sequestration). Here, we extend and present a non-intrusive reduced order model of natural convection in porous media employing deep convolutional autoencoders for the compression and reconstruction and either radial basis function (RBF) interpolation or artificial neural networks (ANNs) for mapping parameters of partial differential equations (PDEs) on the corresponding nonlinear manifolds. To benchmark our approach, we also describe linear compression and reconstruction processes relying on proper orthogonal decomposition (POD) and ANNs. Further, we present comprehensive comparisons among different models through three benchmark problems. The reduced order models, linear and nonlinear approaches, are much faster than the finite element model, obtaining a maximum speed-up of 7 × 10 6 because our framework is not bound by the Courant–Friedrichs–Lewy condition; hence, it could deliver quantities of interest at any given time contrary to the finite element model. Our model’s accuracy still lies within a relative error of 7% in the worst-case scenario. We illustrate that, in specific settings, the nonlinear approach outperforms its linear counterpart and vice versa. We hypothesize that a visual comparison between principal component analysis (PCA) and t-Distributed Stochastic Neighbor Embedding (t-SNE) could indicate which method will perform better prior to employing any specific compression strategy.

97 MATHEMATICS AND COMPUTING↗

A tight-binding molecular dynamics study of the noble metals Cu, Ag and Au

Here we have used the Naval Research Lab (NRL) tight-binding (TB) method to study the electronic and mechanical properties of the noble metals. In order to perform molecular dynamics simulations, we used new TB parameters that work well at smaller interatomic distances. The TB parameters were fitted to the fcc, bcc and sc periodic structures and were demonstrated to be transferable and robust for calculating additional dynamical properties which they had not been fitted to. We calculated the phonon frequencies and density of states at finite temperature and we also performed simulations to determine the temperature dependence of the coefficient of thermal expansion and the mean squared displacement. The energy for vacancy formation as well as energy for fcc-based, bcc-based clusters and icosahedral clusters of different sizes were also calculated. The results compared very well with experimental observations and independent (non-fitted) first-principles density functional calculations.

36 MATERIALS SCIENCE↗

Quantifying the Uncertainty in Deterministic Phonon Transport Calculations of Thermal Conductivity using Polynomial Chaos Expansions

The nature of computational simulations requires the inclusion of an uncertainty analysis, as we have limited knowledge of all physically determined input parameters for a computationally simulated problem. We rely on uncertainty quantification (UQ) to characterize our confidence in the outcomes. Quantifying uncertainty can provide a basis for certifications in high-consequence decisions, such as nuclear reactor design, and is a fundamental component of model validation. We employ a previously developed method of uncertainty quantification, polynomial chaos expansion with stochastic collocation (PCE-SC), applied to deterministic phonon transport simulations. In these simulations, we use the neutron transport code Rattlesnake, which solves the Self-Adjoint Angular Flux (SAAF) formulation of the transport equation with a continuous finite element (CFEM) spatial discretization and discrete ordinates, spherical harmonics angular discretizations. Rattlesnake was developed in the multi-physics object oriented simulation environment (MOOSE) framework. We have previously shown Rattlesnake to be effective in simulating phonon transport. We aim to provide a deterministic phonon transport framework for heterogeneous nuclear fuel with fission product defects to predict thermal conductivity (?). A first principles, physics-based calculation of thermal conductivity must involve factors such as the microstructure of nuclear fuel, which constantly changes during the fission process through the formation of isotopic decay products. Heat transport in oxide nuclear fuels is dominated by phonon transport. Impurities in the bulk material influence the transport of energy at the fundamental level, altering the scattering behavior of phonons and electrons. Conventionally, heat transport follows classical physics based on the heat equation derived from Fourier's law q = -k∇T, where q is heat flux, k is thermal conductivity and ∇T is a temperature gradient. However, Fourier's law is a macroscopic empirical law in which the thermal conductivity? does not have a mechanistic connection to the underpinning heat transport processes. Thermal conductivity of a material depends both on the material's intrinsic ability to transport heat and a variety of resistive effects caused by defects in the material. Thus ab initio prediction of the macroscopic conductivity - the property of interest for safe reactor operation - requires simulating detailed processes of heat transport, and then determining the effective thermal conductivity of the simulated material from the resulting heat flux under the imposed temperature difference. The thermal conductivity of a bulk, homogeneous, dielectric material can be well estimated by the mechanistically derived expression: k{sub bulk} = 1/3 C{sub ν}νgλ (2), where C{sub ν}, νg, and λ are the volumetric specific heat, phonon speed, and phonon mean free path, respectively. In these problems the relevant parameter is the material's acoustic thickness, its characteristic distance L relative to the phonon mean free path λ. In Eq. (2) the largest source of uncertainty is λ. While propagating uncertainty in λ to uncertainty in k through Eq. (2) is trivial, systems with extrinsic scattering require more sophisticated approaches. At the microscopic level, uncertainty in λ changes both a material's intrinsic thermal conductivity. Our goal is to propagate the uncertainty in λ through the deterministic transport computation of k. Because there are a small number of uncertain input quantities, we use the method of polynomial chaos expansion (PCE), which expresses solutions in the form of spectral expansions of an uncertain variable. This approach combines both intrusive and nonintrusive methods of uncertainty propagation techniques and results in a unique formulation which is very effective and efficient for problems with few uncertain parameters. Uncertainty Quantification and Sensitivity Analysis We use PCE-SC to measure propagation of uncertainty in a 3-D phonon transport problem of homogeneous silicon, to compute the mean and variance in temperature, heat flux and thermal conductivity. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Schwarz alternating method for transient solid dynamics

Abstract In our earlier work, we formulated the Schwarz alternating method as a means for concurrent multiscale coupling in finite deformation solid mechanics for quasi‐static problems. Herein, we advance this method for the study of transient dynamic multiscale solid mechanics problems where information is exchanged back and forth between small and large scales. The extension to dynamics relies on the notion of a global time stepper. Within each global time step, the subdomains are coupled by the standard Schwarz iterative process. Remarkably, each subdomain can use its own time step or even its own time integrator to advance its solution in time, provided that they synchronize at each global time step. We study the performance of the Schwarz method on several examples designed for this purpose. Our numerical experiments demonstrate that the method is capable of coupling regions with different mesh resolutions, different element types, and different time integration schemes (e.g., implicit and explicit), all without introducing any artifacts that afflict other coupling methods for transient dynamics. Finally, we apply the dynamic Schwarz alternating method to the simulation of a bolted joint subjected to dynamic loading, as a demonstration of the performance of the method in a realistic scenario.

Mota, Alejandro↗

Variational Self-Consistent Theory for Beam-Loaded Cavities

A new variational theory is presented for beam loading in microwave cavities. The beam{field interaction is formulated as a dynamical interaction whose stationarity according to Hamilton’s principle will naturally lead to steady-state solutions that indicate how a cavity’s resonant frequency, Q and optimal coupling coefficient will detune as a result of the beam loading. A driven cavity Lagrangian is derived from first principles, including the effects of cavity wall losses, input power and beam interaction. The general formulation is applied to a typical klystron input cavity to predict the appropriate detuning parameters required to maximize the gain (or modulation depth) in the average Lorentz factor boost, (Δ γ). Numerical examples are presented, showing agreement with the general detuning trends previously observed in the literature. The developed formulation carries several advantages for beam-loaded cavity structures. It provides a self-consistent model for the dynamical (nonlinear) beam{field interaction, a procedure for maximizing gain under beam-loading conditions, and a useful set of parameters to guide cavity-shape optimization during the design of beam-loaded systems. Enhanced clarity of the physical picture underlying the problem seems to be gained using this approach, allowing straightforward inclusion or exclusion of different field configurations in the calculation and expressing the final results in terms of measurable quantities. Two field configurations are discussed for the klystron input cavity, using finite magnetic confinement or no confinement at all. Formulating the problem in a language that is directly accessible to the powerful techniques found in Hamiltonian dynamics and canonical transformations may potentially carry an additional advantage in terms of analytical computational gains, under suitable conditions

43 PARTICLE ACCELERATORS↗

An asymptotic preserving unified gas kinetic scheme for frequency-dependent radiative transfer equations

This paper presents an extension of previous work (Sun et al., 2015 [22]) of the unified gas kinetic scheme (UGKS) for the gray radiative transfer equations to the frequency-dependent (multi-group) radiative transfer system. Different from the gray radiative transfer equations, where the optical opacity is only a function of local material temperature, the simulation of frequency-dependent radiative transfer is associated with additional difficulties from the frequency-dependent opacity. For the multiple frequency radiation, the opacity depends on both the spatial location and the frequency. For example, the opacity is typically a decreasing function of frequency. At the same spatial region the transport physics can be optically thick for the low frequency photons, and optically thin for high frequency ones. Therefore, the optical thickness is not a simple function of space location. In this paper, the UGKS for frequency-dependent radiative system is developed. The UGKS is a finite volume method and the transport physics is modeled according to the ratio of the cell size to the photon's frequency-dependent mean free path. When the cell size is much larger than the photon's mean free path, a diffusion solution for such a frequency radiation will be obtained. On the other hand, when the cell size is much smaller than the photon's mean free path, a free transport mechanism will be recovered. In the regime between the above two limits, with the variation of the ratio between the local cell size and photon's mean free path, the UGKS provides a smooth transition in the physical and frequency space to capture the corresponding transport physics accurately. The seemingly straightforward extension of the UGKS from the gray to multiple frequency radiation system is due to its intrinsic consistent multiple scale transport modeling, but it still involves lots of work to properly discretize the multiple groups in order to design an asymptotic preserving (AP) scheme in all regimes. The current scheme is tested in a few frequency-dependent radiation problems, and the results are compared with the solutions from the well-defined implicit Monte Carlo (IMC) method. The UGKS is much more efficient than IMC, and the computational times of both schemes for all test cases are listed. The UGKS seems to be the first discrete ordinate method (DOM) for the accurate capturing of multiple frequency radiative transport physics from ballistic particle motion to the diffusive wave propagation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Novel Honeycomb-DD Multi Coil Design for Wireless Power Transfer Systems

Wireless power transfer (WPT) technology has received significant attention recently as an alternative charging method for batteries in a wide range of power levels. In the literature, many types of coil structures have been well-studied. Honeycomb coil arrays, which have been used mostly for low-power applications, have not been well-studied for high-power-level applications. In this paper, a novel honeycomb-DD coil design is proposed for high-power wireless battery charging systems. The proposed coil design with step by step design process is given and Finite Element Analysis (FEA) was performed to observe the full performance characteristics of the system for a 100 kW system. In addition, the misalignment tolerance of the proposed system was observed by shifting the secondary side charging pad at different misalignment positions, and the electromagnetic compatibility to the standards was investigated. The core and strand losses were obtained. The results show that the honeycomb coil array provides high coupling coefficient and better misalignment tolerance, making it a potential coil topology for WPT applications.

Aydin, Emrullah↗

Scalable and compact magnetocaloric heat pump technology

Magnetocaloric heat pumping (MCHP) promises to be more efficient than traditional vapor compression while also eliminating the deleterious effects of gaseous refrigerants. While MCHP devices have shown the temperature spans and efficiencies needed for different heating and cooling applications, they struggle to become commercially viable due to their large size and mass, and resultant high cost. This paper evaluates a baseline MCHP device and explores methods to boost its system power density (SPD). The key components of the baseline system are the gadolinium packed-particle bed active magnetic regenerator (AMR) and a magnetic source composed of permanent magnets and high permeability magnetic steel. To enhance the SPD, the paper evaluates maximizing the AMR volume, opting for first-order magnetocaloric materials, optimizing the magnet and AMR geometry, and reducing the size of magnets and magnetic steel parts. At larger thermal powers, increasing the AMR diameter and the number of magnetic poles were evaluated. Using finite element models, solid models, and estimates of magnetocaloric material performance, thermal powers ranging from 37 W to 44 kW at a nominal 10 K temperature span were projected, and SPD was estimated to improve from 6 W/kg to 81 W/kg. Neglecting end effects, an upper limit of 114 W/g is estimated. Compared to SPD of off-the-shelf compressors with similar environment temperatures, MCHP power density using gadolinium is competitive up to roughly 200 W of cooling power. This is extended to 1 kW when using LaFeSi alloys and up to 3 kW in the limiting case. In conclusion, these results indicate that the performance and mass of MCHP can match that of compressors, which is a critical step toward cost-competitive magnetocaloric technology.

42 ENGINEERING↗

A well-posed and stable stochastic Galerkin formulation of the incompressible Navier–Stokes equations with random data

We present a well-posed stochastic Galerkin formulation of the incompressible Navier–Stokes equations with uncertainty in model parameters or the initial and boundary conditions. The stochastic Galerkin method involves representation of the solution through generalized polynomial chaos expansion and projection of the governing equations onto stochastic basis functions, resulting in an extended system of equations. A relatively low-order generalized polynomial chaos expansion is sufficient to capture the stochastic solution for the problem considered. We derive boundary conditions for the continuous form of the stochastic Galerkin formulation of the velocity and pressure equations. The resulting problem formulation leads to an energy estimate for the divergence. With suitable boundary data on the pressure and velocity, the energy estimate implies zero divergence of the velocity field. Based on the analysis of the continuous equations, we present a semi-discretized system where the spatial derivatives are approximated using finite difference operators with a summation-by-parts property. With a suitable choice of dissipative boundary conditions imposed weakly through penalty terms, the semi-discrete scheme is shown to be stable. Numerical experiments in the laminar flow regime corroborate the theoretical results and we obtain high-order accurate results for the solution variables and the velocity divergence converges to zero as the mesh is refined.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reinvestigation into Closure Predictions of Room D at the Waste Isolation Pilot Plant

Room D was an in-situ, isothermal, underground experiment conducted at the Waste Isolation Pilot Plant between 1984 and 1991. The room was carefully instrumented to measure the horizontal and vertical closure immediately upon excavation and for several years thereafter. Early finite element simulations of salt creep around Room D under-predicted the vertical closure by 4.5×, causing investigators to explore a series of changes to the way Room D was modeled. Discrepancies between simulations and measurements were resolved through a series of adjustments to model parameters, which were openly acknowledged in published reports. Interest in Room D has been rekindled recently by the U.S./German Joint Project III and Project WEIMOS, which seek to improve the predictions of rock salt constitutive models. Joint Project participants calibrate their models solely against laboratory tests, and benchmark the models against underground experiments, such as room D. This report describes updating legacy Room D simulations to today’s computational standards by rectifying several numerical issues. Subsequently, the constitutive model used in previous modeling is recalibrated two different ways against a suite of new laboratory creep experiments on salt extracted from the repository horizon of the Waste Isolation Pilot Plant. Simulations with the new, laboratory-based, calibrations under-predict Room D vertical closure by 3.1×. A list of potential improvements is discussed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗

Uncertainty Quantification for Dissimilar Material Joints Under Corrosion Environment

Abstract Self-Piercing Riveting (SPR) is one of the most commonly used methods for joining dissimilar materials in the automotive industry. These joints are popular due to their adaptability, high performance and short cycle time. However, since these joints involve two dissimilar materials, they are susceptible to galvanic corrosion in the presence of an electrolyte which is driven by the difference in the equilibrium potential of the metals. This can affect the safety and resilience of these joints. In this paper, we focus on galvanic corrosion in Al-Fe SPR joints. A Machine learning (ML) based surrogate model, which is based off of FE simulations, for statistical corrosion analysis is developed. This model enables the resilience and reliability analysis of SPR joints under corrosion environment. In this study, first a physics-based finite element (FE) corrosion model has been developed to simulate the galvanic corrosion between a Fe cathode and an Al anode of a SPR joint. This model takes into account the effect of the crystal microstructure of the Al anode and the precipitation of the corrosion product. Several geometric and environmental factors including crevice gap, roughness of anode, conductivity, pH and the temperature of the electrolyte that effect corrosion are investigated. A thorough Uncertainty Quantification (UQ) analysis is conducted for the overall corrosion behavior of the Fe-Al SPR joints using a novelistic Probabilistic Confidence-Based Adaptive Sampling (PCAS) technique. PCAS is used to train the surrogate model by identifying the critical sampling points and thus reducing the overall computational costs. It is found that the electrolyte temperature has the largest effects on the material loss and needs to be managed closely for better corrosion control. By understanding the corrosion performance and resultant uncertainty impact on SPR joints, the reliability and resilience of these joints can be improved.

36 MATERIALS SCIENCE↗