{"title":"Multiphysics Numerical Method for Modeling Josephson Traveling-Wave Parametric Amplifiers","authors":"Samuel T. Elkin;Michael Haider;Thomas E. Roth","doi":"10.1109/JMMCT.2024.3428344","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3428344","url":null,"abstract":"Josephson traveling-wave parametric amplifiers (JTWPAs) are wideband, ultralow-noise amplifiers used to enable the readout of superconducting qubits. While individual JTWPAs have achieved high performance, behavior between devices is inconsistent due to wide manufacturing tolerances. Amplifier designs could be modified to improve resilience towards variations in amplifier components; however, existing device models often rely on analytical techniques that typically fail to incorporate component variations. To begin addressing this issue, a 1D numerical method for modeling JTWPAs is introduced in this work. The method treats the Josephson junctions and transmission lines in an amplifier as coupled subsystems and can easily incorporate arbitrary parameter variations. We discretize the transmission line subsystem with a finite element time domain method and the Josephson junction subsystem with a finite difference method, with leap-frog time marching used to evolve the system in time. We validate our method by comparing the computed gain to an analytical model for a traditional JTWPA architecture and one with resonant phase matching. We then use our method to demonstrate the impact of variations in Josephson junctions and phase-matching resonators on amplification. In future work, the method will be adjusted to incorporate additional amplifier architectures and extended to a 3D full-wave approach.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"247-257"},"PeriodicalIF":1.8,"publicationDate":"2024-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141965524","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jinghan Xu;Shengguo Xia;Lixue Chen;Chengxian Li;Hongdan Yang
{"title":"Sheet Element Approximation for Numerical Study of Current on Armature and Rail Interface","authors":"Jinghan Xu;Shengguo Xia;Lixue Chen;Chengxian Li;Hongdan Yang","doi":"10.1109/JMMCT.2024.3422609","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3422609","url":null,"abstract":"The armature and rail (A/R) interface is an imperfect contact that is made at discrete asperities at the microscale resulting from high contact pressure. The current distribution of the interface differs significantly from the bulk behavior. In this paper, based on the contact layer model (CLM) and the Cooper-Mikic-Yoranovich model (CMYM), sheet element approximation and boundary conditions are proposed to analyze the electromagnetic properties of the A/R interface. Assuming zero gradients of the magnetic vector in the thickness direction, there are two ways for the approximation, which are mathematical approximation (MA) and physical approximation (PA). Results from both methods show high agreement, consistent with results from slit boundary conditions. Current distributions on both stationary and sliding A/R interfaces are numerically investigated. On the stationary interface, current diffuses from the edges to the central part of the real contact area, whereas on the sliding interface, current concentration occurs at the trailing edge due to the velocity skin effect (VSE). Furthermore, the contour of the current distribution aligns with the erosion pattern observed in experiments, validating the accuracy of the computational method.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"228-235"},"PeriodicalIF":1.8,"publicationDate":"2024-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141725682","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Two-Dimensional Coupled Electrothermal Method Based on the Unstructured Transmission-Line Modelling Method for Lightning Protection Simulations","authors":"Kaiqi Yan;Ana Vukovic;Phillip Sewell","doi":"10.1109/JMMCT.2024.3421958","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3421958","url":null,"abstract":"This paper outlines a fully coupled electrothermal time-domain method to model the effects of lightning strikes and the formation of plasma. The plasma material is described by using the Drude model. This method predicts the formation of the discharge channel by solving the electromagnetic field and the temperature before, during and after the air breaks down. The proposed method is applied to analyse the performance of a number of segmented lightning diverter strips used for lightning protection.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"208-217"},"PeriodicalIF":1.8,"publicationDate":"2024-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141602542","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Stephen D. Gedney;Nastaran Hendijani;John C. Young;Robert J. Adams
{"title":"Electrostatic Boundary Integral Method for 3D Structures in a Layered Conducting Medium","authors":"Stephen D. Gedney;Nastaran Hendijani;John C. Young;Robert J. Adams","doi":"10.1109/JMMCT.2024.3416688","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3416688","url":null,"abstract":"An integral equation formulation is presented for the modeling of the electrostatic fields surrounding arbitrary three-dimensional structures situated in a conducting layered medium. The layered Green's function for the electrostatic potential and the tensor Green's function for the gradient potential are derived. Closed forms for the 3D layered Green's functions are generated using a discrete complex image method (DCIM) approximation. Improved accuracy of the DCIM approximation is achieved using optimization for the computation of the DCIM poles and residues. The problem is discretized via a high-order locally corrected Nyström method with curvilinear cells. Several examples are shown that demonstrate the accuracy of the DCIM approximation for layered media with disparate layer spacing and conductivities for arbitrary 3D geometries.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"218-227"},"PeriodicalIF":1.8,"publicationDate":"2024-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141602554","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jie Li;Min Tang;Lin-Sheng Wu;Liguo Jiang;Wenliang Dai;Junfa Mao
{"title":"LB-ADI: An Efficient Method for Transient Thermal Simulation of Integrated Chiplets and Packages","authors":"Jie Li;Min Tang;Lin-Sheng Wu;Liguo Jiang;Wenliang Dai;Junfa Mao","doi":"10.1109/JMMCT.2024.3386842","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3386842","url":null,"abstract":"In this article, an efficient Laguerre-based alternating direction implicit (LB-ADI) approach is proposed for the transient thermal simulation of integrated chiplets and packages. The transient heat conduction equation is transformed into the Laguerre domain by the Laguerre basis functions and the Galerkin's testing method. With spatial discretization, the resulting matrix equation based on a marching-on-in-order scheme is established. In order to improve the computational efficiency, a new ADI approach in the Laguerre domain is developed. Only three tridiagonal matrices need to be solved in each order, which significantly reduces the simulation time and memory requirement. The accuracy and efficiency of the proposed method are validated by numerical results.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"149-156"},"PeriodicalIF":2.3,"publicationDate":"2024-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140619562","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A Hybrid Electromagnetic Optimization Method Based on Physics-Informed Machine Learning","authors":"Yanan Liu;Hongliang Li;Jian-Ming Jin","doi":"10.1109/JMMCT.2024.3385451","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3385451","url":null,"abstract":"In this article, we present an optimization method based on the hybridization of the genetic algorithm (GA) and gradient optimization (grad-opt) and facilitated by a physics-informed machine learning model. In the proposed method, the slow-but-global GA is used as a pre-screening tool to provide good initial values to the fast-but-local grad-opt. We introduce a robust metric to measure the goodness of the designs as starting points and use a set of control parameters to fine tune the optimization dynamics. We utilize the machine learning with analytic extension of eigenvalues (ML w/AEE) model to integrate the two pieces seamlessly and accelerate the optimization process by speeding up forward evaluation in GA and gradient calculation in grad-opt. We employ the divide-and-conquer strategy to further improve modeling efficiency and accelerate the design process and propose the use of a fusion module to allow for end-to-end gradient propagation. Two numerical examples are included to show the robustness and efficiency of the proposed method, compared with traditional approaches.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"157-165"},"PeriodicalIF":2.3,"publicationDate":"2024-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=10493126","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140813892","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Zhiwen Dong;Xinlei Chen;Fan Gao;Changqing Gu;Zhuo Li;Wu Yang;Weibing Lu
{"title":"Efficient Iterative Solution of Combined Source Integral Equation Using Characteristic Basis Function Method With Initial Guess","authors":"Zhiwen Dong;Xinlei Chen;Fan Gao;Changqing Gu;Zhuo Li;Wu Yang;Weibing Lu","doi":"10.1109/JMMCT.2024.3382725","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3382725","url":null,"abstract":"Using only the RWG functions, the combined source integral equation (CSIE) with weak form combined source condition can achieve fine accuracy and fast iterative convergence for conductor objects. However, compared with a conventional integral equation in the method of moments (MoM), the conventional CSIE involves more matrices and more complex numerical processing, and these make the CSIE inefficient, especially for multiple excitation problems. In this article, a characteristic basis function (CBF)-based CSIE with initial guess is proposed to mitigate this problem. The CBF is employed to reduce the number of unknowns as well as the storage consumptions and iteration time. In the meantime, an initial guess especially for CBFs is proposed to reduce iterations when solving multiple excitation problems. Numerical results are given to demonstrate the performance of the proposed method.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"142-148"},"PeriodicalIF":2.3,"publicationDate":"2024-03-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140546589","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Kangshuai Du;Shilie He;Chengzhuo Zhao;Na Liu;Qing Huo Liu
{"title":"A 3-D Spectral Element Time-Domain Method With Perfectly Matched Layers for Transient Schrödinger Equation","authors":"Kangshuai Du;Shilie He;Chengzhuo Zhao;Na Liu;Qing Huo Liu","doi":"10.1109/JMMCT.2024.3399911","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3399911","url":null,"abstract":"A spectral element time-domain (SETD) method with perfectly matched layers (PML) is proposed to simulate the behavior of electron waves, interference effects and tunneling effects, in three-dimensional (3-D) devices by solving Schrödinger equation. The proposed method employs Gauss-Lobatto-Legendre (GLL) polynomials to represent the wave function. Easy construction of higher-order element makes refinement straightforward and spectral accuracy can be obtained from the SETD. Meanwhile, by utilizing the GLL quadrature, a diagonal mass matrix is obtained which is meaningful in the time-stepping process. Numerical experiments confirm that, for open boundary problems, employing PML yields results characterized by high numerical efficiency, remarkable flexibility and ease of implementation. These findings underscore the effectiveness of SETD-PML in addressing the challenges posed by open boundary conditions, making it a reliable choice for numerical simulations. Some illustrative numerical examples are presented to demonstrate the performance of the proposed method.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"188-197"},"PeriodicalIF":2.3,"publicationDate":"2024-03-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141292484","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Coupled Mechanical and Electromagnetic Analysis of Current on Armature and Rail Interface With Dynamic Contact","authors":"Jinghan Xu;Shengguo Xia;Hongdan Yang;Lixue Chen","doi":"10.1109/JMMCT.2024.3397464","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3397464","url":null,"abstract":"The electrical contact between the armature and rail (A/R) in a railgun is acknowledged as a dynamic sliding interface, exhibiting properties distinct from bulk. This paper employs a 3-D finite element method (FEM) for coupled mechanical and electromagnetic analysis and proposes boundary conditions for dynamic sliding contact to investigate current distribution on the A/R interface. Results show that contact pressure and area have similar trends as the driving current, which confines current distributed areas. The current distributions on stationary and sliding interfaces reveal different patterns but the distributed areas both locate within the contact areas. In the case of the stationary scenario, the current concentrates at the trailing edge when the current increases and diffuses to the leading edge when the current declines. However, due to the velocity skin effect (VSE), the current fails to diffuse into the interior during all stages. Besides, comparative calculations with constant contact indicate that forced shifts of current occur when the contact is dynamic, dominating the current distributions of the A/R interface. Moreover, the influence of the VSE on forced shifts of current is notable, with significant current variations observed near the trailing edge, whereas those around the leading edge are less pronounced.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"198-207"},"PeriodicalIF":2.3,"publicationDate":"2024-03-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141292551","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Equivalent Circuit Model Development Accounting for Mutual-Coupling Effects","authors":"Chandan Roy;Ke Wu","doi":"10.1109/JMMCT.2024.3396801","DOIUrl":"https://doi.org/10.1109/JMMCT.2024.3396801","url":null,"abstract":"Mutual-coupling effects are of utmost importance in the development of high-frequency circuits and systems. However, it is a common practice to ignore those couplings when establishing equivalent circuit models. Neglecting these couplings leads to inaccurate circuit modelling. Therefore, it becomes imperative to account for mutual couplings in the development of accurate equivalent circuit models. This work presents a holistic process for synthesizing the equivalent circuit model of an electromagnetic (EM) field structure that incorporates mutual couplings of varying orders. The proposed high-order framework begins by developing equivalent circuit models for each individual transmission line discontinuity within the target circuit. Subsequently, the mutual couplings of different orders are extracted in a step-by-step manner. Throughout this process, full-wave EM simulations are deployed, along with a circuit parameter extraction method that utilizes de-embedded circuit responses. By combining these techniques, a comprehensive and accurate equivalent circuit model is generated, enabling a detailed analysis of the target field model structure, and facilitating a deeper understanding of its electrical and magnetic behavior and performance. This paper utilizes a three-step microstrip discontinuity structure and a third-order parallel coupled microstrip filter as examples for theoretical and experimental demonstration of the proposed technique.","PeriodicalId":52176,"journal":{"name":"IEEE Journal on Multiscale and Multiphysics Computational Techniques","volume":"9 ","pages":"166-178"},"PeriodicalIF":2.3,"publicationDate":"2024-03-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140919115","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}