SEARCH · Engineering Papers

Results for “Singularity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Split representation of adaptively compressed polarizability operator

The polarizability operator plays a central role in density functional perturbation theory and other perturbative treatment of first principle electronic structure theories. The cost of computing the polarizability operator generally scales as O(Ne4) where Ne is the number of electrons in the system. The recently developed adaptively compressed polarizability operator (ACP) formulation [L. Lin, Z. Xu and L. Ying, Multiscale Model. Simul. 2017] reduces such complexity to O(Ne 3 ) in the context of phonon calculations with a large basis set for the first time, and demonstrates its effectiveness for model problems. In this paper, we improve the performance of the ACP formulation by splitting the polarizability into a near singular component that is statically compressed, and a smooth component that is adaptively compressed. The new split representation maintains the O(Ne 3 ) complexity, and accelerates nearly all components of the ACP formulation, including Chebyshev interpolation of energy levels, iterative solution of Sternheimer equations, and convergence of the Dyson equations. For simulation of real materials, we discuss how to incorporate nonlocal pseudopotentials and finite temperature effects. In this work, we demonstrate the effectiveness of our method using one-dimensional model problem in insulating and metallic regimes, as well as its accuracy for real molecules and solids.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Sensitivity Analysis for Solutions to Heterogeneous Nonlocal Systems. Theoretical and Numerical Studies

The paper presents a collection of results on continuous dependence for solutions to nonlocal problems under perturbations of data and system parameters. The integral operators appearing in the systems capture interactions via heterogeneous kernels that exhibit different types of weak singularities, space dependence, even regions of zero-interaction. Here, the stability results showcase explicit bounds involving the measure of the domain and of the interaction collar size, nonlocal Poincaré constant, and other parameters. In the nonlinear setting, the bounds quantify in different L p norms the sensitivity of solutions under different nonlinearity profiles. The results are validated by numerical simulations showcasing discontinuous solutions, varying horizons of interactions, and symmetric and heterogeneous kernels.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Synchrotron X-ray microtomography and multifractal analysis for the characterization of pore structure and distribution in softwood pellet biochar

Biochar pores in the micrometer range (1–100 µm) derive from cellular structures of the plant biomass subjected to pyrolysis or can be the result of mechanical processing, such as pelleting. In this study, synchrotron X-ray microtomography was used to investigate the internal pore structure of softwood pellet biochar produced by slow pyrolysis at 550 and 700 °C. The microtomographic data sets consisted of 2025 images of 2560 × 2560 voxels with a voxel side length of 0.87 µm. The three-dimensional reconstructions revealed that pelleting and pyrolysis significantly altered the pore structures of the wood feedstock, creating a network of connected pores between fragments that resembled the wood morphology. While higher pyrolysis temperature increased the specific surface area (as determined by BET nitrogen adsorption), it did not affect the total observed porosity. Multifractal analysis was applied to assess the characteristics of the frequency distribution of pores along each of the three dimensions of reconstructed images of five softwood pellet biochar samples. The resulting singularity and Rényi spectra (generalized dimensions) indicated that the distribution of porosity had monofractal scaling behavior, was homogeneous within the analyzed volumes and consistent between replicate samples. Moreover, the pore distributions were isotropic (direction-independent), which is in strong contrast with the anisotropic pore structure of wood. As pores at the scale analyzed in this study are relevant, for example, for the supply of plant accessible water and habitable space for microorganisms, our findings combined with the ability to reproduce biochar with such pore distribution offer substantial advantages in various biochar applications.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Chapter 32 - Power Grid Resilience

The new energy paradigm is altering the current power grid trend from synchronous generator reliant system toward power-electronics-based distributed energy resources (DERs). The resilient operation of modern power grid is highly dependent upon cyber and physical reliable operation of DERs. Power grid resilience broadly refers to the ability of the grid to be robust to eventualities and singularities that may disrupt the continuity of reliable power flow. This could involve resilience related to the physical system but more recently, the focus on cyber resilience has been evolving rapidly especially given the burgeoning growth of distributed generation and power-electronics-based DERs under smart grid and given that such threats to reliable operation do not have to evolve localized to where the physical asset is. Commonly, in power grids dominated by DERS, control and energy management are structured at a multitime scale multilayer fashion: (i) primary control layer at microseconds time scale, (ii) secondary control layer at millisecond time scale, and (iii) tertiary control layer at seconds to minutes time scale. The cyber-related issues that may affect the power grid resilience can be introduced in through primary, secondary, and tertiary control layers. The failure of each of these control layers may impact the reliable and resilient operation of the overall network and cause widespread failures which leads to unintended blackouts. As such, this chapter focuses on cyber-related resilience issues in primary control layer, secondary control layer, tertiary control layer, wide area/utility control, and anomaly detection and resilient communication.

anomaly detection↗

DOE OSTI2026

Emerging hierarchical dislocation structures: Insights from scanning electron microscopy-electron backscatter diffraction in situ tensile testing and multifractal analysis

Understanding the evolution of dislocation structures during plastic deformation is critical for predicting the mechanical performance of metallic materials. In this work, we applied in situ scanning electron microscopy/electron backscatter diffraction tensile testing combined with multifractal (MF) analysis to assess deformation-induced dislocation structure evolution in solution-annealed 304 L stainless steel, both in its as-received and neutron-irradiated states (5.4 displacements per atom). The analysis of kernel average misorientation patterns revealed the formation of hierarchical dislocation arrangements that exhibit clear MF scaling behavior. Despite pronounced visual differences between nonirradiated and irradiated specimens—most notably, the appearance of dislocation channels after irradiation—the singularity spectra suggest that both conditions give rise to similar underlying hierarchical structures. MF analysis provides a quantitative measure of the spatial complexity and self-organization of dislocation patterns, highlighting the accelerated emergence and evolution of the dislocation structures in irradiated polycrystalline materials, as well as the limitation of their spatial extent. The findings indicate that irradiation not only modifies microstructure but also alters correlation-driven dislocation organization. More generally, they demonstrate that MF analysis is a powerful tool for probing mesoscale deformation mechanisms.

Dislocation structures↗

DOE OSTI2024

Impacts of feeding three strains of microalgae alone or in combination on growth performance, protein metabolism, and meat quality of broiler chickens

Variations in nutrient compositions, especially amino acid (AA) profiles, among microalgal species may enable a superior feeding outcome from a combined than singular supplementation in poultry diets. Therefore, a feeding trial was conducted to compare the effects of three strains of microalgal biomass supplemented alone or in combination to replace 5 % (starter) and 10 % (grower) soybean meal (on weight-to-weight basis) on growth performance, protein metabolism, and meat quality of broiler chickens. Day-old Cornish Cross male chicks (total = 180) were divided into 5 groups (6 cages/treatment, 6 birds/cage) and fed a corn-soybean meal basal diet (BD), BD + H117 (Chlorella sp., H117), BD + C985 (Tetraselmis sp., C985), BD + Nannochloropsis oceanica (NO), and BD + H117 + C985 + NO (Combination). Feeding any of the microalgae diets did not alter growth performance nor meat quality including texture, pH, color, and water holding capacity of breast and thigh meats. However, the breast weight percentages were decreased (P < 0.05) by feeding the C985, NO, and Combination diets. Compared with the BD, the 4 microalgal diets led to higher (P < 0.05) plasma uric acid and protein concentrations at weeks 3 and (or) 6. The mRNA levels of MAFbx, MURF1, FOXO1, and calpastatin in the breast and thigh muscles were altered by the microalgal diets but not those of genes associated with other quality traits. In conclusion, replacing 5 % or 10 % soybean meal with three sources of microalgae in broiler diets decreased breast weights percentage but not absolute weight. Furthermore, feeding chickens with the combination of three microalgae did not restore the breast loss and induced different expressions of genes related to muscle hypertrophy or atrophy.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗

DOE OSTI2022

Parametric model-order reduction for radiation transport using multi-resolution proper orthogonal decomposition

For parametric high-fidelity simulations, it is often desirable to utilize a reduced-order model (ROM) to emulate, at a reduced computational cost, parametric solutions of the governing partial differential equations (PDEs) for unseen parameter values. One commonly employed option is to utilize a data-driven, projection-based ROM supplemented with subspace identification via proper orthogonal decomposition (POD). POD discovers the ROM subspace by computing the singular value decomposition (SVD) of a set of training data from the full-order model (FOM). In streaming-dominated radiation transport simulations with localized sources, solutions often greatly vary over the spatial domain by many orders of magnitude. In such cases, machine-precision arithmetic can be insufficient to obtain an accurate SVD, resulting in a poorly performing ROM. We present a method called multiresolution POD (mrPOD) that mitigates these inaccuracies. The mrPOD method works by decomposing the spatial domain into regions and performing proper orthogonal decomposition on the training dataset separately in each region. In conclusion, mrPOD is tested on single energy group and multigroup atmospheric shielding transport problems and is shown to outperform classic POD.

42 ENGINEERING↗

DOE OSTI2020

Resurgence of the renormalization group equation

We show how the renormalons emerge from the renormalization group equation with a priori no reference to any Feynman diagrams. The proof is rather given by recasting the renormalization group equation as a resurgent equation studied in the mathematical literature, which describes a function with an infinite number of singularities in the positive axis of the Borel plane. Consistency requires a one-to-one correspondence between the existence of such kind of equation and the actual (generalized) Borel resummation of the renormalons through a one-parameter transseries. Our finding suggests how non-perturbative contributions can affect the running couplings. We also discuss these concepts within the context of gauge theories, making use of the large number of flavor expansion.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Geometry and entanglement in the scattering matrix

A formulation of nucleon–nucleon scattering is developed in which the S-matrix, rather than an effective-field theory (EFT) action, is the fundamental object. Spacetime plays no role in this description: the S-matrix is a trajectory that moves between RG fixed points in a compact theory space defined by unitarity. This theory space has a natural operator definition, and a geometric embedding of the unitarity constraints in four-dimensional Euclidean space yields a flat torus, which serves as the stage on which the S-matrix propagates. Trajectories with vanishing entanglement are special geodesics between RG fixed points on the flat torus, while entanglement is driven by an external potential. The system of equations describing S-matrix trajectories is in general complicated, however the very-low-energy S-matrix –that appears at leading-order in the EFT description– possesses a UV/IR conformal invariance which renders the system of equations integrable, and completely determines the potential. In this geometric viewpoint, inelasticity is in correspondence with the radius of a three-dimensional hyperbolic space whose two-dimensional boundary is the flat torus. This space has a singularity at vanishing radius, corresponding to maximal violation of unitarity. The trajectory on the flat torus boundary can be explicitly constructed from a bulk trajectory with a quantifiable error, providing a simple example of a holographic quantum error correcting code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

New static solutions of symmetric $\phi^4$ equation

In this paper, we provide new exact solutions of nonlinear Klein–Gordon ( $\phi^4$) equation in 1 + 1-dimension. For simplicity, we focus on the static equation and ignore the time-dependence. The symmetric $\phi^4$ equation has played an important role in several areas of physics. We obtain several novel non-singular solutions of the symmetric $\phi^4$ model in terms of the Jacobi elliptic functions and compare them with the well-known solutions. Finally, we categorize these solutions in terms of the potential parameters.

36 MATERIALS SCIENCE↗

DOE OSTI2024

HHL algorithm with mapping function and enhanced sampling for model predictive control in microgrids

Here, this paper presents a refined quantum Harrow Hassidim Lloyd (HHL) algorithm for microgrid control. The first novelty of the developed method is that a mapping shift function enables the original HHL algorithm to handle general linear equations with non-singular and indefinite matrix. Second, a method of Matrix Extension for Amplifying Sampling Probabilities of Intended Solution (ME-ASPI) is proposed to design the reformulated linear algebraic equations, allowing for improved sampling efficiency of the quantum tomography in the refined HHL algorithm. Then, we applied the method to solve the model predictive control (MPC) problem in nonlinear dynamical microgrids. Specifically, with the ME-ASPI method, the refined HHL algorithm can effectively obtain the intended partial optimal control inputs for MPC. The optimization of quadratic programming problem in each time step of MPC is transformed into a linear system problem, which is addressed by the proposed quantum solver through using only partial information, with the time complexity improved from $\mathscr{O}(\mathscr{N}^{2.37286})$ classically to $\mathscr{O}(\mathscr{N}^{2} log \mathscr{N}$ x $p$ log $p)$ in quantum. Numerical examples have validated the effectiveness of the refined HHL algorithm with the proposed mapping function and the ME-ASPI method. By leveraging quantum properties, the proposed method provides a hybrid quantum–classical framework for microgrid control. This generic method can also potentially tackle many other challenges in analyzing and controlling general complex engineered systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Towards a reliable assessment of charging effects during surface analysis: Accurate spectral shapes of ZrO 2 and Pd/ZrO 2 via X-ray Photoelectron Spectroscopy

X-ray Photoelectron Spectroscopy of large bandgap or insulating material surfaces relies on an effective mechanism that compensates for the emission (loss) of electrons by maintaining the material surface at a steady-state uniform potential. While a steady-state may be attained by utilizing an active compensation, such as low power electron emitting filament, there is the possibility that the surface potential is not uniform over the area analyzed, leading to peak shifts and incorrect spectral interpretation. Here, in this work, a spectral data processing method based on mapping the ZrO 2 and Pd/ZrO 2 surfaces utilizing photoemission peak binding energy is proposed, which provides information about the response of specific material surfaces to charge compensation. Spectromicroscopy of ZrO 2 and Pd/ZrO 2 surfaces without spatial information is used to monitor the efficacy of charge compensation. Exploiting counts distributed over many bins require the use of procedures and algorithms essential to practical mapping peak positions. Iterative singular value decomposition is therefore introduced and utilized as a means of efficiently delivering spatially resolved spectra from which binding energy for peaks is computed. The concepts developed in this work result in robust and accurate peak models of ZrO 2 and Pd/ZrO 2 that can be applied in XPS analysis of not only ZrO 2 but other large bandgap or insulating material surfaces. Supporting arguments for a peak model representing signal from Zr 3p and Pd 3d are developed within this work are presented.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Information theory optimization of signals from small-angle scattering measurements

Small-angle X-ray scattering (SAXS) of particles in solution informs on the conformational states and assemblies of biological macromolecules (bioSAXS) outside of cryo- and solid-state conditions. In bioSAXS, the SAXS measurement under dilute conditions is resolution limited, and through an inverse Fourier transform, the measured SAXS intensities directly relate to the physical space occupied by the particles via the P (r)-distribution. Yet, this inverse transform of SAXS data has been historically cast as an ill-posed, ill-conditioned problem requiring an indirect approach. Here, we show that through the applications of matrix and information theories, the inverse transform of SAXS intensity data is a well-conditioned problem. The so-called ill-conditioning of the inverse problem is directly related to the Shannon number. By exploiting the oversampling enabled by modern detectors, a direct inverse Fourier transform of the SAXS data is possible, provided the recovered information does not exceed the Shannon number. The Shannon limit corresponds to the maximum number of significant singular values that can be recovered in a SAXS experiment, suggesting this relationship is a fundamental property of band-limited inverse integral transform problems. This correspondence reduces the complexity of the inverse problem to the Shannon limit and maximum dimension. We propose a hybrid scoring function using an information theory framework that assesses both the quality of the model-data fit as well as the quality of the recovered P (r)-distribution. The hybrid score utilizes the Akaike information criteria and Durbin-Watson statistic that considers parameter-model complexity, i.e., degrees of freedom, and the randomness of the model-data residuals. The described tests and findings extend the boundaries for bioSAXS by completing the information theory formalism initiated by Peter B. Moore to enable a quantitative measure of resolution in SAXS, robustly determine maximum dimension, and more precisely define the best parameter model appropriately representing the observed scattering data.

Rambo, Robert P. [Science and Technology Facilitie↗

DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗

DOE OSTI2023

A unified hp-HDG framework for Friedrichs' PDE systems

Here, this work proposes a unified hp-adaptivity framework for hybridized discontinuous Galerkin (HDG) method for a large class of partial differential equations (PDEs) of Friedrichs' type. In particular, we present unified hp-HDG formulations for abstract one-field and two-field structures and prove their well-posedness. In order to handle non-conforming interfaces we simply take advantage of HDG built-in mortar structures. With split-type mortars and the approximation space of trace, a numerical flux can be derived via Godunov approach and be naturally employed without any additional treatment. As a consequence, the proposed formulations are parameter-free. We perform several numerical experiments for time-independent and linear PDEs including elliptic, hyperbolic, and mixed-type to verify the proposed unified hp-formulations and demonstrate the effectiveness of hp-adaptation. Two adaptivity criteria are considered: one is based on a simple and fast error indicator, while the other is rigorous but more expensive using an adjoint-based error estimate. The numerical results show that these two approaches are comparable in terms of convergence rate even for problems with strong gradients, discontinuities, and singularities.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Mechanistic pathways and role of oxygen in oxidative coupling of methane derived from transient kinetic studies

Oxidative coupling of methane (OCM) is a promising industrial process to upgrade natural gas to high value chemicals. In this study, Temporal Analysis of Products (TAP) and steady-state experiments were conducted to distinguish how the composition of surface and gas phase oxygen influence mechanistic details of the selective conversion of CH 4 to C 2 H 4 over the Mn-Na 2 WO 4 /SiO 2 catalyst. The results from TAP studies indicate that methane activation on this catalyst proceeds predominantly via a short-lived, transient surface oxygen species and there is a competition for this species to form either CO or methyl radicals on the surface. This active species has a total lifetime of 3 s and is identified to have a dioxygen (e.g. O 2 2- or O 2 -) form. We show that the concentration of the transient surface oxygen species significantly impacts the OCM performance. Oxygen attributed to the catalyst lattice (in a singular form e.g., O - ), is found to activate methane to a lesser degree, but exclusively forms CO 2 . Evidence for these surface pathways for methyl radical, CO and CO 2 formation identified by TAP are also validated through steady-state experiments. Finally, by distinguishing different catalyst oxygen species and their role in selective/nonselective pathways, important screening criteria have been identified for the advancement of superior catalyst formulations.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Hydro-chemo-mechanical phase field formulation for corrosion induced cracking in reinforced concrete

Steel corrosion induced cracking in reinforced concrete structures is often caused by volumetric straining at the steel-concrete interface. This volumetric strain is mainly controlled by the rate of generation of corrosion product, which depends on the diffusion of Cl{sup −} and dissolved O{sub 2} as well as the chemical reaction at the steel-concrete interface. While many numerical modeling techniques have been proposed to tackle this complex phenomenon, most have been limited to using staggered solution schemes, where the diffusion, chemical reaction, mechanical strain, and crack fields are solved for separately. Unfortunately, such approaches introduce numerical errors that can lead to unrealistic predictions. Therefore, a coupled solution, is needed to remove the formulation induced errors and understand this complex phenomenon. To provide a unifying solution framework for corrosion induced cracking in reinforced concrete, a phase field formulation is presented and implemented with FEM algorithms. The governing equations are derived from the system's free energy and Fick's law of diffusion. The proposed formulation was first validated in 1D by comparing the modeling results with the semi-analytical solution approximated via Laurent series to avoid singularities in the complex domain. 2D plane strain solutions to the embedded rebar undergoing corrosion induced cracking with pre-cracks are also presented and compared with the approximated solutions from linear elastic fracture mechanics theory (LEFM). The modeling of crack propagation enabled by the proposed formulation is also presented along with parametric studies to reveal the roles of the fracture toughness and permeability of the concrete on the corrosion induced cracking. A case study was also conducted on H-pile steel with concrete jackets to demonstrate the feasibility of modeling corrosion induced fracture in reinforced concrete structures with complex steel-concrete interfaces.

36 MATERIALS SCIENCE↗