Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iteration 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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗
DOE OSTI2021

LoTToR: An Algorithm for Missing-Wedge Correction of the Low-Tilt Tomographic 3D Reconstruction of a Single-Molecule Structure (LoTToR) v1.0

The low-tilt tomographic reconstruction (LoTToR) method is a post-processing method to correct the missing-wedge artifact, which contains a model-free iteration process under a set of constraints in real and reciprocal spaces. A proof of concept is conducted by using the LoTToR on a phantom, i.e., a simulated 3D reconstruction from a low-tilt series of images, including that within a tilt range of ±15°. The method is validated by using both negative-staining (NS) and cryo-electron tomography (cryo-ET) experimental data. A significantly reduced missing-wedge artifact verifies the capability of LoTToR, suggesting a new tool to support the future study of macromolecular dynamics, fluctuation and chemical activity from the viewpoint of single-molecule 3D structure determination. Moreover, the software can be used for hospital computer tomography (CT) to reduce the radiation damage to patients without reducing the CT image quality.

Ren, Gang↗
DOE OSTI2022

Randomized probe imaging through deep k-learning

Randomized probe imaging (RPI) is a single-frame diffractive imaging method that uses highly randomized light to reconstruct the spatial features of a scattering object. The reconstruction process, known as phase retrieval, aims to recover a unique solution for the object without measuring the far-field phase information. Typically, reconstruction is done via time-consuming iterative algorithms. In this work, we propose a fast and efficient deep learning based method to reconstruct phase objects from RPI data. The method, which we call deep k-learning, applies the physical propagation operator to generate an approximation of the object as an input to the neural network. This way, the network no longer needs to parametrize the far-field diffraction physics, dramatically improving the results. Deep k-learning is shown to be computationally efficient and robust to Poisson noise. The advantages provided by our method may enable the analysis of far larger datasets in photon starved conditions, with important applications to the study of dynamic phenomena in physical science and biological engineering.

Guo, Zhen (ORCID:0000000213473451)↗
DOE OSTI2020

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method—the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]—we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. Overall, the resulting qDSRG(2)+(T) approach matches the accuracy of the “gold-standard” coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Parallel derivative-free optimization for simulation-based design of behind-the-meter energy systems

In this work, the integrated design and dispatch of behind-the-meter or distributed resources (e.g. stationary battery storage and solar PV generation) is considered. A simulation-based framework is employed, generating high-fidelity results with closed-loop predictive control at a fine resolution, at the expense of high computational cost (several minutes to a few hours per design point). To address this challenge, parallel derivative-free design methods are considered. Four methods are compared, including state-of-the-art surrogate-based methods (Radial-Basis Functions and Gaussian processes) and sampling strategies, an evolutionary-based method, and a simple sequential grid refinement method. As a case study, two types of design problem with increasing complexity are considered, namely, the design of behind-the-meter resources (three design variables) and the inclusion of grid capacity (four design variables). The second yields a constrained design problem for which violations can only be determined after solving the computationally expensive simulation. For the three-dimensional case, all methods present a good performance, achieving a solution within 1% of the optimum after the first iteration, with the sequential grid refinement exhibiting the fastest convergence and achieving the best final objective value. This indicates that the parallel evaluation of multiple sampling points may be more important than the choice of method for small decision spaces. For the four-dimensional constrained case, the Genetic Algorithm presents the best tradeoff between performance and computational effort, while the rough objective function terrain generated by constraint violation penalties reduces the performance of surrogate-based methods. Contour plots with flat regions indicate flexibility in the optimal design and highlight the importance of characterizing the solution space.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

A Faster-Than-Real-Time Framework for Reliability-Oriented Simulation of PV Inverters

Physics-of-Failure (PoF) based reliability assessment for photovoltaic (PV) inverters requires long-duration electrical and electrothermal stress histories, yet generating such stress histories with high-fidelity switching models over year long mission profiles is computationally prohibitive. Conventional methods either sacrifice modeling fidelity for speed or require runtimes that are impractical for design iteration and uncertainty studies. To address this bottleneck, this paper presents a High-Performance Computing (HPC) based simulation frame work for faster-than-real-time reliability-oriented simulation. The proposed framework integrates the Average-to-Switching (A2S) method with parallel computing techniques to accelerate switching-level waveform reconstruction. We further introduce optimization strategies, including cluster merging and sensitivity based mission profile screening, to reduce the computational burden. Evaluated using real-world mission profile inputs and a MATLAB/Simulink switching-model reference, the framework reduces the simulation time for a one-year mission from an intractable multi-year duration to approximately 7.3 minutes while maintaining low waveform error. This acceleration provides a practical reliability-oriented simulation engine that can be coupled with component-specific aging models for subsequent PV inverter PoF assessment.

High-performance Computing↗
DOE OSTI2022

Joint routing of conventional and range-extended electric vehicles in a large metropolitan network

Range-extended electric vehicles combine the higher efficiency and environmental benefits of battery-powered electric motors with the longer mileage and autonomy of conventional internal combustion engines. This combination is particularly advantageous for time-constrained delivery routing in dense urban areas, where battery recharging along routes can be too time-consuming to economically justify the use of all-electric vehicles. However, switching from electric to conventional fossil fuel modes also results in higher costs and emissions and lower efficiency. This paper analyzes this heterogeneous vehicle routing problem and describes two solution methods: an exact branch-price-and-cut algorithm and an iterated tabu search metaheuristic. From a methodological perspective, we find that the exact algorithm consistently obtains tight lower bounds that also serve to certify the metaheuristic solutions as near-optimal. From a policy standpoint, we examine a large-scale real-world case study concerning parcel deliveries in the Chicago metropolitan area and quantify various operational metrics including energy costs and vehicle miles traveled. We find that by deploying roughly 20% of range -extended vehicles with a modest all-electric range of 33 miles, parcel distributors can save energy costs by up to 17% while incurring less than 0.5% increase in vehicle miles traveled. Increasing the range to 60 miles further reduces costs by only 4%, which can alternatively be achieved by decreasing the average service time by 1 minute or increasing driver working time by 1 hour. In conclusion, our study reveals several key areas of improvement on which vehicle manufacturers, distributors, and policy makers can focus their attention.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Reducing Artifacts in BF and HAADF-STEM Images of Pt/C Fuel Cells using MBIR-ARAR

Electron tomography is a powerful tool for characterization of three-dimensional (3D) nanoscale materials and devices. Bright-field (BF) and high-angle annular dark-field (HAADF) scanning transmission electron microscopy (STEM) are two widely used imaging modes in electron tomography. These imaging modalities have proven to be useful for characterizing the structure of carbon-supported platinum (Pt/C) electrocatalysts used in fuel cells, which are an important class of clean energy conversion systems. BF and HAADF-STEM images are typically simultaneously acquired due to the complementary information they contain, as BF-STEM is more suitable for characterizing carbon due to its ability to detect lighter elements, while HAADF-STEM is more suitable for characterizing platinum due to its sensitivity to atomic number. However, the quality of BF- and HAADF-STEM images is often compromised by various artifacts, such as local blurring due to abrupt contrast changes as well as missing wedge artifacts which can severely impact the accuracy of the reconstructed 3D images. The traditional methods of reconstruction, such as Filtered Back-projection (FBP), Simultaneous Iterative Reconstruction Technique (SIRT), and Model-Based Iterative Reconstruction (MBIR), are often not suitable for removing these artifacts effectively. Here, to address this challenge, we present a novel approach called MBIR with Artifact Reduction and Adaptive Regularization (MBIR-ARAR).

25 ENERGY STORAGE↗
DOE OSTI2026

Impact of Resonant Magnetic Perturbations on the Toroidal Location of the Runaway Electron Final Loss Strike Point

It is demonstrated that the peak toroidal impact phase of the post-disruption runaway electrons (REs) can be varied shot-to-shot by means of applied static non-axisymmetric (3D) magnetic fields, or resonant magnetic perturbations (RMPs). In the experiments, high-current (500 kA), purged RE plateaus (i.e. low impurity post-disruption plasmas with current carried completely by REs) were terminated on the DIII-D tokamak center post (CP). The global toroidal impact peak of the REs appears to be dominantly set by the toroidal phase of a wall-locked tearing mode, typically (m/n) = (2/1), but sometimes (m/n) = (3/2). In the absence of external RMP, the locked mode reliably locks close to a preferred toroidal phase Φ ≃ 75◦ probably set by intrinsic error fields. With applied RMP, the locked mode and RE impact phase is pushed to a new location. When the tearing mode is (m/n) = (2/1), the phase is well controlled and tracks the applied RMP phase, but in the case of an (m/n) = (3/2) tearing mode, the new phase location is not well controlled. These results could have application in ITER and other fusion power plant (FPP) - sized tokamaks as a method for spreading disruption first wall damage from REs

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Solar Thermochemical Carbon Dioxide Splitting Using Ceria and Iron Aluminate Foam Devices and Simulation of a Plant System for Demonstration

An international research project has been undertaken to integrate a unique solar thermal processing reactor system with ceria and iron aluminate as active redox materials for CO2 splitting. Experimental investigations for CO2 splitting were conducted using a solar simulator and tube furnace at Niigata University, followed by demonstrations using a high-flux solar furnace (HFSF) at the National Renewable Energy Laboratory (NREL) in Golden, CO. Each experimental setup consisted of foam devices composed of reticulated porous ceramic (RPC). The RPC has a full ceria or iron aluminate body. It fabricated using the replica method and subjected to a two-step redox reaction, which iteratively separated a stream of CO2 into O2 and CO. Reactivity was evaluated using CO production per mass of the reactive material. The tubular furnace yielded a CO production of 6.41 mL/g at a reduction temperature of 1600degrees C, showcasing a higher CO production rate and total amount than those obtained from experiments conducted with solar simulators and solar furnace setups. For iron aluminate RPC, the productivity was measured as 3.57 mL/g using HFSF at a reduction temperature of 1450degrees C. These results are somewhat higher than those of the previous experiment at lower reduction temperatures of 1400degrees C-1500degrees C. Additionally, the production of CO in the case of ceria RPC was compared with the steady flow model simulation, which assumed chemical equilibrium at various levels of oxygen partial pressure during the reduction process. On the basis of these results, this study proposes a solar fuel system with an open receiver that uses a high-temperature heat transfer fluid.

carbon dioxide thermochemical splitting↗
DOE OSTI2022

NF-κB perturbation reveals unique immunomodulatory functions in Prx1 + fibroblasts that promote development of atopic dermatitis

Skin is composed of diverse cell populations that cooperatively maintain homeostasis. Up-regulation of the nuclear factor κB (NF-κB) pathway may lead to the development of chronic inflammatory disorders of the skin, but its role during the early events remains unclear. Here, through analysis of single-cell RNA sequencing data via iterative random forest leave one out prediction, an explainable artificial intelligence method, we identified an immunoregulatory role for a unique paired related homeobox-1 (Prx1) + fibroblast subpopulation. Disruption of Ikkb–NF-κB under homeostatic conditions in these fibroblasts paradoxically induced skin inflammation due to the overexpression of C-C motif chemokine ligand 11 (CCL11; or eotaxin-1) characterized by eosinophil infiltration and a subsequent T H 2 immune response. Because the inflammatory phenotype resembled that seen in human atopic dermatitis (AD), we examined human AD skin samples and found that human AD fibroblasts also overexpressed CCL11 and that perturbation of Ikkb–NF-κB in primary human dermal fibroblasts up-regulated CCL11. Monoclonal antibody treatment against CCL11 was effective in reducing the eosinophilia and T H 2 inflammation in a mouse model. Together, the murine model and human AD specimens point to dysregulated Prx1 + fibroblasts as a previously unrecognized etiologic factor that may contribute to the pathogenesis of AD and suggest that targeting CCL11 may be a way to treat AD-like skin lesions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Direct Feed High-Level Waste APPS Model Glass Testing (DFHLW APPS) Matrix, Phase 2

This report summarizes the data collected during the batching and melting of a second matrix of Direct Feed High-Level Waste (DFHLW) glasses generated using the preliminary enhanced waste glass models (EWG2.5) and the Britton and Anderson (2024) preliminary DFHLW feed vector. The purpose of these glasses is two-fold: 1. Validate EWG2.5 glass calculations being used in the Aspen Process Performance Simulation (APPS) model. 2. Evaluate and ultimately improve the glass property models and formulation methods used for design of DFHLW glasses as part of an iterative process of data collection and model refinement. Some of the 16 APPS2 glasses tested did not satisfy all target property constraints due to the limited data on DFHLW glass supporting the EWG2.5 models. • One glass, APPS2-10, formed nepheline on canister centerline cooling (CCC) heat-treatment and failed the product consistency test (PCT) response limits. This glass also had high B and Cr release rates for the toxicity characteristic leaching procedure (TCLP). All other glasses were found to satisfy the PCT and TCLP constraints for both quenched and CCC samples. • One glass, APPS2-08, had higher than acceptable viscosity due to magnetite crystallization. • One glass, APPS2-09, formed greater than 2 vol% crystals at 950 °C. As the glass design criterion was that the temperature at 2 vol% crystal (T 2% ) be less than 950 °C, only one glass failed the criteria. However, this criterion is being reevaluated. Four additional glasses formed crystal fractions between 1 and 2 vol% at 950 °C (APPS2-03, -08, -12, and -14). • Four glasses – APPS2-01, -02, -04, and -16 – failed the Monofrax K-3 refractory neck corrosion (k neck ) design limit of 0.04 in. at 1208 °C for 6 d. This is another criterion being reevaluated. Four additional glasses (APPS2-05, -06, -11, and -13) exhibited 0.025 = k neck = 0.04 in. • All 16 glasses passed the sulfur solubility and TCLP constraints. The measured property values were compared to predicted values using EWG2.5 and a selection of other existing models. A few models (e.g., electrical conductivity, TCLP) were found to be adequate for designing DFHLW glasses in the near future, while others require refits or offsets. It is recommended that new property models be developed for EWG3.0, as a large amount of DFHLW glass property data (> 14 × existing data) is expected to be collected in the compositional spaces where no data was previously available. To enable near-term calculations and formulations for designing DFHLW glasses and processing rate estimations, a formulation algorithm with minor modifications will be developed, EWG2.6.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Chemical and Radiological Compatibility Testing of 3D Printed Materials

To dramatically increase the adaptability, performance, and safety of processes in support of the Defense Waste Processing Facility (DWPF), Savannah River National Lab (SRNL) plans to perform chemical and radiological compatibility testing on a wide variety of 3D printed materials of interest. The 3D printing process provides numerous strategic operational benefits such as rapid prototyping of complex designs and geometry specific to the needs of the nuclear waste disposition process, as well as on-demand rapid prototyping and iteration with materials that aren’t as accessible through traditional manufacturing methods. Reaction chemistry in simulated waste batches can be matched closely to its radioactive counterpart, but glass reactor vessels have limitations. Vessel geometry can play a big factor in mixing transport limitations, process chemistry, and degradation reaction kinetics. In addition, additive manufacturing allows for much more detailed vessel design than traditional alternatives. Waste processing techniques in DWPF also encounter extreme chemical environments including high pH, strong acids, abrasive slurries, and significant irradiation. To meet these challenges, a matrix of various polymer, ceramic, and metal additive manufacturing materials have been exposed to a suite of chemical environments of interest as well as radioactive dose (such as gamma radiation from 60 Co) to properly test their durability under these conditions. Mass change has been monitored over a period of up to a week in these conditions, as well as added characterization for surface modification through Scanning Electron Microscopy/Electron Dispersive X-ray analysis (SEM/EDX). Further chemical characterization has been monitored through Fourier-Transform InfraRed Spectroscopy (FTIR), with planned investigation via thermal and tensile strength degradation. While the direct product of this research is identification of material(s) that can withstand specific hazardous environments encountered by the mercury water wash tank in DWPF process simulation experiments, the reference base of materials will be used for many other nuclear processes in the pursuit of rapidly developed, cost-efficient, and highly specific devices for environmental remediation and much more.

Wilson, Nathan W. [Savannah River National Laborat↗
DOE OSTI2024

RESOLVING THE ELECTROCHEMICAL EQUATIONS OF A SOLID OXIDE FUEL CELL FOR USE IN TRANSIENT SIMULATION AND INTEGRATION INTO CYBER-PHYSICAL SYSTEMS

A major challenge with complex cyber-physical systems stems from long model computational time that creates a mismatch between the model system and the physical system. The numerical modeling of solid oxide fuel cells (SOFCs) presents particular challenges due to the highly coupled nature of the underlying equations and the multiphysics needed to fully resolve their behavior during a transient event. To this end current approaches revolve around splitting the computational efforts into resolving temperature effects and resolving electrochemical effects. Current methods employed for the transient simulation of an SOFC for implementation in the Hybrid Performance (HyPer) facility cyberphysical plant at the National Energy Technology Laboratory reveal a distinct need for accelerated results with a high degree of stability. To this aim, an investigation into the computational time for the code reveals that the underlying electrochemical algorithm takes an order of magnitude more time than its thermal counterpart and has a tendency to vary in terms of iteration time and as such a rework of the underlying system is proposed. The primary method for accelerated electrochemical algorithm solutions is to employ higher order root finding recipes for the resolution of the highly coupled electrochemical equations. This is done with the intention to reduce the overall number of subiterations necessary for resolving voltage, current density, and species concentration, properties of the fuel cell that are all directly coupled and require nested iterative approaches. The overall objective of this approach is an order of magnitude reduction in calculation time without sacrificing stability and increasing accuracy. Specific approaches involve using both bounded and unbounded techniques, such as the False Position method and the Secant method (or if applicable Newton-Raphson) respectively, the drawbacks being slower convergence for False Position and instability for the Secant or Newton-Raphson methods. Current preliminary results on simplified versions of the parent functions involved for electrochemical calculations indicate a reduction in computational steps by a factor of two for the secant method and a factor of three for Newton-Raphson. When implemented into new modified electrochemical algorithms, the results indicate a possible order of magnitude reduction in calculation time.

Arias, Jesus↗
DOE OSTI2024

A Fast Algebraic Multigrid Solver and Accurate Discretization for Highly Anisotropic Heat Flux I: Open Field Lines

We present a novel solver technique for the anisotropic heat flux equation, aimed at the high level of anisotropy seen in magnetic confinement fusion plasmas. Such problems pose two major challenges: (i) discretization accuracy and (ii) efficient implicit linear solvers. We simultaneously address each of these challenges by constructing a new finite element discretization with excellent accuracy properties, tailored to a novel solver approach based on algebraic multigrid (AMG) methods designed for advective operators. We pose the problem in a mixed formulation, introducing the directional temperature gradient as an auxiliary variable. The temperature and auxiliary fields are discretized in a scalar discontinuous Galerkin space with upwinding principles used for discretizations of advection. We demonstrate the proposed discretization’s superior accuracy over other discretizations of anisotropic heat flux, achieving error 1000x smaller for anisotropy ratio of 10 9 , for closed field lines. The block matrix system is reordered and solved in an approach where the two advection operators are inverted using AMG solvers based on approximate ideal restriction, which is particularly efficient for upwind discontinuous Galerkin discretizations of advection. To ensure that the advection operators are nonsingular, in this paper we restrict ourselves to considering open (acyclic) magnetic field lines for the linear solvers. We demonstrate fast convergence of the proposed iterative solver in highly anisotropic regimes where other diffusion-based AMG methods fail.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Calculation of permanent magnet arrangements for stellarators: a linear least-squares method

A problem arising in several engineering areas is to design magnets outside a volume that produce a desired magnetic field inside it. One instance of this problem is stellarator design, where it has recently been shown that permanent magnets can provide the required shaping of the magnetic field. Here we demonstrate a robust and efficient algorithm REGCOIL_PM to calculate the spatial distribution of these permanent magnets. The procedure involves a small number of fixed-point iterations, with a linear least-squares problem solved at each step. Furthermore, the method exploits the Biot–Savart Law's exact linearity in magnetization density and approximate linearity in magnet size, for magnets far from the target region. No constraint is placed on the direction of magnetization, so Halbach solutions are found naturally, and the magnitude of the magnetization can be made uniformly equal to a target value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Fourier Analysis Method for Wireless Power Transfer Coil Design

Wireless power transfer (WPT) systems for Electric Vehicles (EVs) are designed to meet specifications such as stray field, power transfer, efficiency, and ground clearance. Typical design approaches include iterative analysis of predetermined coil geometries to identify candidates that meet these constraints. This work instead directly generates WPT coil shapes and magnetic fields to meet specifications and constraints through the optimization of Fourier basis function coefficients. The proposed Fourier Analysis Method (FAM) applies to arbitrary planar coil geometries and does not rely on iterative finite-element analysis (FEA) simulations. This flexibility allows for rapid design evaluation across a larger range of coil geometries and design specifications. A prototype coil is built to compare FAM outputs to experimental measurements and FEA simulations. The FAM is then used to illustrate the tradeoff of coil current and stray field for a given power level showing that the method is capable of generating optimized coil shapes to meet arbitrary field constraints.

Foote, Andrew↗