Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Sensitivity analysis of a layered piezoelectric system using ZFEM

The complex variable finite element method (ZFEM) is a numerical technique which aims to find the partial derivatives of the independent variables with respect to variation in dependent parameters declared in the physics. This is done by combining the complex Taylor series expansion within the weak formulation of the governing equation in a coupled system of linear equations forming a complex valued block matrix given by the Cauchy–Riemann matrix representation. In this work, two-dimensional linear first-order elements have been implemented in ZFEM to predict the design derivatives of the mechanical displacement field and the voltage potential field for a layered piezoelectric system in a steady-state study with Dirichlet boundary condition applied at the top and bottom edges of the geometry. This approach allows the standard FEM solution to quantify the sensitivity of the mechanical displacement and voltage potential fields with respect to small variations in the material properties through the information obtained from the computation of the derivatives. The domain is formed by a layered body with PZT-4 and PZT-5 stacked together. For result verification, the numerical solution obtained with ZFEM was compared to results from a commercial FEM package and the solution from the imaginary part was compared to the exact solution of a well-known benchmark problem. In conclusion, comparison of the results showed good agreement for both the real and imaginary parts of the solution and the largest sensitivities were found in PZT-5 specifically in C 13 , C 33 , and ε 33 .

42 ENGINEERING↗
DOE OSTI2022

Numerical fluid dynamics for FRG flow equations: Zero-dimensional QFTs as numerical test cases. I. The O ( N ) model

The functional renormalization group (FRG) approach is a powerful tool for studies of a large variety of systems, ranging from statistical physics over the theory of the strong interaction to gravity. The practical application of this approach relies on the derivation of so-called flow equations, which describe the change of the quantum effective action under the variation of a coarse-graining parameter. In the present work, we discuss in detail a novel approach to solve such flow equations. This approach relies on the fact that RG equations can be rewritten such that they exhibit similarities with the conservation laws of fluid dynamics. This observation can be exploited in different ways. First of all, we show that this allows to employ powerful numerical techniques developed in the context of fluid dynamics to solve RG equations. In particular, it allows us to reliably treat the emergence of nonanalytic behavior in the RG flow of the effective action as it is expected to occur in studies of, e.g., spontaneous symmetry breaking. Second, the analogy between RG equations and fluid dynamics offers the opportunity to gain novel insights into RG flows and their interpretation in general, including the irreversibility of RG flows. Further, we work out this connection in practice by applying it to zero-dimensional quantum-field theoretical models. The generalization to higher-dimensional models is also discussed. Our findings are expected to help improving future FRG studies of quantum field theories in higher dimensions both on a qualitative and quantitative level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2019

Predicting the X-Ray Spectra of Stellar-mass Black Holes from Simulations

We describe results from a new technique for the prediction of complete, self-consistent X-ray spectra from three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations of black hole accretion flows. Density and cooling rate data from a harm3d GRMHD simulation are post-processed by an improved version of the Monte Carlo radiation transport code pandurata (in the corona) and the Feautrier solver PTRANSX (in the disk), with xstar subroutines. The codes are run in a sequential, iterative fashion to achieve globally energy-conserving and self-consistent radiation fields, temperature maps, and photoionization equilibria. The output is the X-ray spectrum as seen by a distant observer, including features, such as the Fe Kα emission line and corresponding K-edge absorption trough, due to disk reprocessing of coronal power. For the example cases we consider—a non-rotating 10 M(sub ⊙) black hole with solar abundances, accreting at 0.01, 0.03, 0.1, or 0.3 Eddington—we find spectra resembling actual observations of stellar-mass black holes in the soft or steep power-law state: broad thermal peaks (at 1–3 keV), steep power laws extending to high energy (Γ = 2.7–4.5), and prominent, asymmetric Fe Kα emission lines with equivalent widths in the range 40–400 eV (larger EW at lower accretion rates). By starting with simulation data, we obviate the need for parameterized descriptions of the accretion flow geometry—no a priori specification of the corona's shape or flux, or the disk temperature or density, etc., is needed. Instead, we apply the relevant physical principles to simulation output using appropriate numerical techniques; this procedure allows us to calculate inclination-dependent spectra after choosing only a small number of physically meaningful parameters: black hole mass and spin, accretion rate, and elemental abundances.

Kinch, Brooks E.↗
NASA NTRS1996

Theoretical and Numerical Investigation of Radiative Extinction of Diffusion Flames

The influence of soot radiation on diffusion flames was investigated using both analytical and numerical techniques. Soot generated in diffusion flames dominate the flame radiation over gaseous combustion products and can significantly lower the temperature of the flame. In low gravity situations there can be significant accumulation of soot and combustion products in the vicinity of the primary reaction zone owing to the absence of any convective buoyant flow. Such situations may result in substantial suppression of chemical activities in a flame, and the possibility of a radiative extinction may also be anticipated. The purpose of this work was to not only investigate the possibility of radiative extinction of a diffusion flame but also to qualitatively and quantitatively analyze the influence of soot radiation on a diffusion flame. In this study, first a hypothetical radiative loss profile of the form of a sech(sup 2) was assumed to influence a pure diffusion flame. It was observed that the reaction zone can, under certain circumstances, move through the radiative loss zone and locate itself on the fuel side of the loss zone contrary to our initial postulate. On increasing the intensity and/or width of the loss zone it was possible to extinguish the flame, and extinction plots were generated. In the presence of a convective flow, however, the movement of the temperature and reaction rate peaks indicated that the flame behavior is more complicated compared to a pure diffusional flame. A comprehensive model of soot formation, oxidation and radiation was used in a more involved analysis. The soot model of Syed, Stewart and Moss was used for soot nucleation and growth and the model of Nagle and Strickland-Constable was used for soot oxidation. The soot radiation was considered in the optically thin limit. An analysis of the flame structure revealed that the radiative loss term is countered both by the reaction term and the diffusion term. The essential balance for the soot volume fraction was found to be between the processes of soot convection and soot growth. Such a balance yielded to analytical treatment and the soot volume fraction could be expressed in the form of an integral. The integral was evaluated using two approximate methods and the results agreed very well with the numerical solutions for all cases examined.

Ray, Anjan↗
DOE OSTI2025

Simulation Models for Exploring Magnetic Reconnection

Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This article reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

Atomistic Method Applied to Computational Modeling of Surface Alloys

The formation of surface alloys is a growing research field that, in terms of the surface structure of multicomponent systems, defines the frontier both for experimental and theoretical techniques. Because of the impact that the formation of surface alloys has on surface properties, researchers need reliable methods to predict new surface alloys and to help interpret unknown structures. The structure of surface alloys and when, and even if, they form are largely unpredictable from the known properties of the participating elements. No unified theory or model to date can infer surface alloy structures from the constituents properties or their bulk alloy characteristics. In spite of these severe limitations, a growing catalogue of such systems has been developed during the last decade, and only recently are global theories being advanced to fully understand the phenomenon. None of the methods used in other areas of surface science can properly model even the already known cases. Aware of these limitations, the Computational Materials Group at the NASA Glenn Research Center at Lewis Field has developed a useful, computationally economical, and physically sound methodology to enable the systematic study of surface alloy formation in metals. This tool has been tested successfully on several known systems for which hard experimental evidence exists and has been used to predict ternary surface alloy formation (results to be published: Garces, J.E.; Bozzolo, G.; and Mosca, H.: Atomistic Modeling of Pd/Cu(100) Surface Alloy Formation. Surf. Sci., 2000 (in press); Mosca, H.; Garces J.E.; and Bozzolo, G.: Surface Ternary Alloys of (Cu,Au)/Ni(110). (Accepted for publication in Surf. Sci., 2000.); and Garces, J.E.; Bozzolo, G.; Mosca, H.; and Abel, P.: A New Approach for Atomistic Modeling of Pd/Cu(110) Surface Alloy Formation. (Submitted to Appl. Surf. Sci.)). Ternary alloy formation is a field yet to be fully explored experimentally. The computational tool, which is based on the BFS (Bozzolo, Ferrante, and Smith) method for the calculation of the energetics, consists of a small number of simple PCbased computer codes that deal with the different aspects of surface alloy formation. Two analysis modes are available within this package. The first mode provides an atom-by-atom description of real and virtual stages 1. during the process of surface alloying, based on the construction of catalogues of configurations where each configuration describes one possible atomic distribution. BFS analysis of this catalogue provides information on accessible states, possible ordering patterns, and details of island formation or film growth. More importantly, it provides insight into the evolution of the system. Software developed by the Computational Materials Group allows for the study of an arbitrary number of elements forming surface alloys, including an arbitrary number of surface atomic layers. The second mode involves large-scale temperature-dependent computer 2. simulations that use the BFS method for the energetics and provide information on the dynamic processes during surface alloying. These simulations require the implementation of Monte-Carlo-based codes with high efficiency within current workstation environments. This methodology capitalizes on the advantages of the BFS method: there are no restrictions on the number or type of elements or on the type of crystallographic structure considered. This removes any restrictions in the definition of the configuration catalogues used in the analytical calculations, thus allowing for the study of arbitrary ordering patterns, ultimately leading to the actual surface alloy structure. Moreover, the Monte Carlo numerical technique used for the large-scale simulations allows for a detailed visualization of the simulated process, the main advantage of this type of analysis being the ability to understand the underlying features that drive these processes. Because of the simplicity of the BFS method for e energetics used in these calculations, a detailed atom-by-atom analysis can be performed at any point in the simulation, providing necessary insight on the details of the process. The main objective of this research program is to develop a tool to guide experimenters in understanding and interpreting often unexpected results in alloy formation experiments. By reducing the computational effort without losing physical accuracy, we expect that powerful simulation tools will be developed in the immediate future, which will allow material scientists to easily visualize and analyze processes at a level not achievable experimentally.

Bozzolo, Guillermo H.↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2022

Probabilistic constrained Bayesian inversion for transpiration cooling

Abstract To enable safe operations in applications such as rocket combustion chambers, the materials require cooling to avoid material damage. Here, transpiration cooling is a promising cooling technique. Numerous studies investigate possibilities to simulate and evaluate the complex cooling mechanism. One naturally arising question is the amount of coolant required to ensure a safe operation. To study this, we introduce an approach that determines the posterior probability distribution of the Reynolds number using an inverse problem and constraining the maximum temperature of the system under parameter uncertainties. Mathematically, this chance inequality constraint is dealt with by a generalized polynomial chaos expansion of the system. The posterior distribution will be evaluated by different Markov chain Monte Carlo based methods. A novel method for the constrained case is proposed and tested among others on two‐dimensional transpiration cooling models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

An Optothermal Raman Model Accounting for Convection, Radiation, and Substrate Effects for Evaluation of the Thermal Conductivity of Two‐Dimensional Coatings

Thermal conductivity of two‐dimensional (2D) materials is important to characterize as this will govern thermal transport physics in the various thermal and energy applications benefitting from the unique properties of these materials. Here, analytical and numerical techniques are presented to determine the thermal conductivity and interfacial conductance of a 2D material coating using confocal micro‐Raman spectroscopy as a noncontact diagnostic. While several methods have been proposed for calculating the values of thermal conductivity, interfacial conductance, and spatial temperature profile, these models often do not consider effects of radiation, convection, and substrate resistance on the temperature and heat flux profiles. We present a model to calculate the thermal conductivity and interfacial conductance which accounts for convection, radiation, and substrate effects to characterize a variety of 2D material coatings, which is demonstrated using large area graphene transferred onto copper and nickel substrates. Convection and radiation effects are found to have a negligible effect on the temperature profile of supported coatings, while the substrate effects have a considerable effect and therefore impact the quantification of thermal conductivity and interfacial conductance in addition to limiting the spatial resolution of the technique. These findings will allow for more accurate and reproducible extraction of the thermal conductivity of graphene and other 2D material coatings when employing noncontact optothermal Raman spectroscopy methods.

2D coatings↗
DOE OSTI2025

Defect Engineering in Large‐Scale CVD‐Grown Hexagonal Boron Nitride: Formation, Spectroscopy, and Spin Relaxation Dynamics

Recently, numerous techniques have been reported for generating optically active defects in exfoliated hexagonal boron nitride (hBN), which hold transformative potential for quantum photonic devices. However, achieving on-demand generation of desirable defect types in scalable hBN films remains a significant challenge. Here, it is demonstrated that formation of negative boron vacancy defects, V B − , in suspended, large-area CVD-grown hBN is strongly dependent on the type of bombarding particles (ions, neutrons, and electrons) and irradiation conditions. In contrast to suspended hBN, defect formation in substrate-supported hBN is more complex due to the uncontrollable generation of secondary particles from the substrate, and the outcome strongly depends on the thickness of the hBN. Different defect types are identified by correlating spectroscopic and optically detected magnetic resonance features, distinguishing boron vacancies (formed by light ions and neutrons and emitting at 800 nm) from other optically active defects emitting at 650 nm assigned to anti-site nitrogen vacancy (N B V N ) and reveal the presence of additional “dark” paramagnetic defects that influence spin-lattice relaxation time (T 1 ) and zero-field splitting parameters, all of which strongly depend on the defect density. These results underscore the potential for precisely engineered defect formation in large-scale CVD-grown hBN, paving the way for the scalable fabrication of quantum photonic devices.

CVD↗
DOE OSTI2023

Calculation of Dynamical Response Functions Using a Bound-State Method

Abstract We investigate a method to extract response functions (dynamical polarisabilities) directly from a bound-state approach applied to calculations of perturbation-induced reactions. The use of a square-integrable basis leads to a response in the form of a sum of $$\delta $$ δ functions. We integrate this over energy and fit a smooth function to the resulting stepwise-continuous one. Its derivative gives the final approximation to the physical response function. We show that the method reproduces analytical results where known, and analyse the details for a variety of models. We apply it to some simple models, using the stochastic variational method as the numerical method. Albeit we find that this approach, and other numerical techniques, have some difficulties with the threshold behavior in coupled-channel problems with multiple thresholds, its stochastic nature allows us to extract robust results even for such cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Origin of the background absorption in carbon nanotubes: Phonon-assisted excitonic continuum

Excitonic effects in 1D semiconductors can be qualitatively different from those in higher dimensions. In particular, the Sommerfeld factor, the ratio of the above-band-edge excitonic continuum absorption to free electron-hole pair generation, has been shown to be less than 1 (suppressed) in 1D systems while it is larger than 1 (enhanced) in 2D and 3D systems. Strong continuum suppression indeed exists in semiconducting single-wall carbon nanotubes, a prototypical 1D semiconductor. However, absorption spectra for carbon nanotubes are typically fit with a combination of Lorentzians and a polynomial background baseline with little physical meaning. We performed absorption measurements in aligned single-chirality (6,5) carbon nanotube films. The obtained spectra were fit with our theoretical model obtained by solving the Boltzmann scattering equation involving fifty-nine different types of transitions among three different types of quasiparticles. Specifically, we took into account microscopic interactions between photons, phonons, and excitons, including their dispersions, which unambiguously demonstrated that the background absorption is due to phonon-assisted transitions from the semiconductor vacuum to finite-momentum continuum states of excitons. The excellent agreement we obtained between experiment and theory suggests that our numerical technique can be seamlessly extended to compute strongly out-of-equilibrium many-body dynamics and time-resolved spectra in low-dimensional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Recent progress in the phase-field dislocation dynamics method

The phase-field dislocation dynamics (PFDD) method, originated in 2002, is a continuum dislocation model that uses order parameters to describe dislocation slips in crystalline materials. In the past two decades, and especially since it was last reviewed in 2016, PFDD was advanced significantly in terms of the mathematical formulation, numerical implementation, and applicability. The main purpose of this short review is to summarize recent progress made to improve the energy functional formulation and numerical techniques of PFDD as well as its recent applications. Additionally, some recommendations for future work to further extend the PFDD method are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A predictive strain-gradient model with no undetermined constants or length scales

A general meso–scale (GM) crystal plasticity (CP) model was developed that accounts for lower-order (strain hardening) and higher-order (internal stress) effects of geometrically necessary dislocations (GNDs). It is predictive: no arbitrary parameters or length scales were invoked and no ad hoc numerical techniques were employed. It uses general stress field equations for GND content and a novel harmonization technique to enforce consistency of elastic long-range singular defect fields with applied elastic-plastic fields. The model facilitates implementation in commercial finite element programs without requiring special elements, special boundary conditions, or access to element shape functions. GM simulations confirmed, with improved accuracy, previously published predictions of the Hall-Petch effect, Bauschinger effect, and anelasticity. Previously unpredicted phenomena were also predicted: anelasticity and hysteresis for single Ta crystals and strain-hardening stagnation. The internal stresses (higher-order effect) dominate at large length scales, while at small length scales, the GND density hardening (lower-order effect) dominates. As a result, GM predicts that strain heterogeneity and consequent GND internal stresses are important factors in anelasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Emergence of 4 H $J^π$ = 1 - resonance in contact theories

We obtain the s- and p-wave low-energy scattering parameters for n 3 H elastic scattering and the position of the 4 H J π = 1 - resonance using the pionless effective field theory at leading order. Results are extracted with three numerical techniques: confining the system in a harmonic oscillator trap, solving the Faddeev-Yakubovsky equations in configuration space, and using an effective two-body cluster approach. The renormalization of the theory for the relevant amplitudes is assessed in a cutoff-regulator range between 1 fm -1 and 10 fm -1 . Most remarkably, we find a cutoff-stable/RG-invariant resonance in the 4 H J π = 1 - system. This p-wave resonance is a universal consequence of a shallow two-body state and the introduction of a three-body s-wave scale set by the triton binding energy. The stabilization of a resonant state in a few-fermion system through pure contact interactions has a significant consequence for the powercounting of the pionless theory. Specifically, it suggests the appearance of similar resonant states also in larger nuclei, like 16-oxygen, in which the theory’s leading order does not predict stable states. Those resonances would provide a starting state to be moved to the correct physical position by the perturbative insertion of sub-leading orders, possibly resolving the discrepancy between data and contact EFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cascading from $\mathscr{N}$ = 2 supersymmetric Yang–Mills theory to confinement and chiral symmetry breaking in adjoint QCD

We argue that adjoint QCD in 3 + 1 dimensions, with any SU(N) gauge group and two Weyl fermion flavors (i.e. one adjoint Dirac fermion), confines and spontaneously breaks its chiral symmetries via the condensation of a fermion bilinear. We flow to this theory from pure $\mathscr{N}$ = 2 SUSY Yang–Mills theory with the same gauge group, by giving a SUSY-breaking mass M to the scalars in the $\mathscr{N}$ = 2 vector multiplet. This flow can be analyzed rigorously at small M, where it leads to a deconfined vacuum at the origin of the $\mathscr{N}$ = 2 Coulomb branch. The analysis can be extended to all M using an Abelian dual description that arises from the N multi-monopole points of the $\mathscr{N}$ = 2 theory. At each such point, there are N −1 hypermultiplet Higgs fields h$^{i=1,2}_m$, which are SU(2) R doublets. We provide a detailed study of the phase diagram as a function of M, by analyzing the semi-classical phases of the dual using a combination of analytic and numerical techniques. The result is a cascade of first-order phase transitions, along which the Higgs fields h i m successively turn on, and which interpolates between the Coulomb branch at small M, where all h$^{i}_m$ = 0, and a maximal Higgs branch, where all h$^{i}_m$ ≠ 0, at sufficiently large M. We show that this maximal Higgs branch precisely matches the confining and chiral symmetry breaking phase of two-flavor adjoint QCD, including its broken and unbroken symmetries, its massless spectrum, and the expected large-N scaling of various observables. The spontaneous breaking pattern SU(2) R → U(1) R , consistent with the Vafa–Witten theorem, is ensured by an intricate alignment mechanism for the h$^{i}_m$ in the dual, and leads to a CP 1 sigma model of increasing radius along the cascade.

D’Hoker, Eric [Univ. of California, Los Angeles, C↗
DOE OSTI2022

Simulating self-powered neutron detector responses to infer burnup-induced power distribution perturbations in next-generation light water reactors

Understanding how 3D power distribution will be monitored throughout reactor core volumetric space in next-generation nuclear power reactors is crucial to the design, deployment, and licensing of these reactors. Although numerous techniques exist for 3D power distribution monitoring based on the response of both in situ and ex situ sensors currently implemented or proposed for use in the US reactor fleet, crucial details about these techniques are often unclear. The publicly available documentation does not include information such as how well these techniques are characterized and optimized in their implementations and the levels of uncertainty in the inferred 3D power distribution. The work described herein investigated a recently developed 3D power distribution inferencing method as applied to two next-generation reactor simulations: (1) the NuScale small modular reactor design and (2) the Westinghouse AP1000 design, both of which contain in-core strings of vanadium self-powered neutron detectors (SPNDs). This investigation considered a range of SPND string sensor densities, as well as a range of 3D power distribution axial segment sizes. In this work, SPND response simulation is informed by neutron flux calculations in representative homogenized cores. For the different sensor densities and power distribution axial segment sizes in these simulations, the average solution error, solver iterations, and run time were tracked to parameterize the sensor-core configuration.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Optimizing the location and configuration of disaster resilience hubs under transportation and electric power network failures

Natural disasters often result in failures of transportation network components and blackouts that imperil the wellbeing of vulnerable populations. In response to these events, resilience hubs have been proposed as a pre-disaster planning strategy to improve access to critical services. This paper introduces an optimization-based approach to locate and configure electric power-generating resilience hubs considering the possibility of failures in transportation and electric power systems. The model's objective is to identify hub locations and configurations that maximize transportation accessibility to the hubs and maximize the satisfaction of basic energy needs through hub-generated electric power. Besides a budget constraint, the model accounts for limits on the levels of hub energy generation vis-à-vis community energy demands, and on the transportation network distance of communities to hubs. Three heuristics are presented for the proposed planning problem. The first heuristic is a genetic algorithm (GA) with problem-specific solution generation procedures. The other two heuristics implement greedy search techniques. Numerical experiments were conducted, using data from rural Puerto Rico, to illustrate the application of the proposed model and heuristics, and examine their performance. In the numerical experiments, the GA heuristic found better solutions than the greedy heuristics. Additionally, design solutions consisting of spatially dispersed hubs with low energy generation capacity were better than solutions with spatially concentrated high-capacity hubs. Lastly, across a wide range of hub demand scenarios, only a small number of candidate hub locations consistently ranked among the best locations for establishing a hub.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗