Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2610.40693

Finite-Element Electromagnetic Analysis and Evaluation of Antenna Performance Metrics Using a Quantum Linear System Solver

Hiroaki Ominato, Takahiro Ohyama, Shiro Kawabata

quant-phcs.CESubmitted Oct 28, 2026• Announced Oct 29, 2026
arXivHTMLPDF

Abstract

Electromagnetic analysis based on the finite element method is widely used in antenna design. At higher frequencies, however, the required mesh refinement makes solving large sparse linear systems the main contributor to computational cost. Quantum linear system solvers offer a potential means of accelerating this step, but their output is a superposition state: reconstructing all components of the classical solution vector would forfeit the quantum speedup. Prior studies have proposed methods for obtaining physical quantities such as scattering cross sections and scattering parameters from the solution state. To support antenna design using fault-tolerant quantum computers, we formulate the total radiated power, the radiated power within a prescribed solid angle, and the dielectric and conductive losses as observables for efficient evaluation from a quantum solution state. These observables avoid the measurement cost of reconstructing all components of the output quantum state. We develop an end-to-end hybrid quantum--classical workflow that combines finite-element electromagnetic analysis, matrix preconditioning, a quantum linear system solver, and the evaluation of antenna design metrics. To validate the workflow, we performed numerical experiments evaluating these observables for an axisymmetric coaxially fed monopole antenna. Under idealized conditions, the relative errors in the total radiated power and the dissipated power with respect to the classical finite-element solution were on the order of $10^{-3}$, confirming the consistency of the workflow. The proposed method is expected to provide a basis for using quantum linear system solvers to evaluate antenna performance and, in the future, to design high-frequency wireless devices.

AI Summary

The paper asks how antenna performance metrics can be extracted from a quantum linear-system solver’s normalized finite-element solution state without reconstructing the full electromagnetic field. It expresses total and solid-angle-restricted radiated power, Huygens-surface power flux, and dielectric and conductive losses as quadratic forms $x^\dagger Qx$. Radiation observables are constructed from far-field linear functionals; loss observables follow from volume and surface integrals of finite-element basis fields. The workflow combines Ruiz diagonal equilibration, Hermitian embedding of the complex indefinite Maxwell matrix, and HHL-type solution-state preparation. Observables are transformed consistently with preconditioning, and physical power scales are recovered using the solution norm inferred from postselection probability. Validation uses idealized classical simulation for an axisymmetric coaxially fed monopole with 519 complex degrees of freedom. Equilibration reduces the condition number from $6.16\times10^5$ to $9.34\times10^3$; the simulation uses 14 phase qubits and has postselection probability $1.29\times10^{-5}$ at the reference frequency. Relative errors against the classical discrete solution are $0.769\%$ for total radiated power and $0.662\%$ for angular-band power. Across a 21-frequency sweep, radiated-power errors remain below $0.693\%$, and total-loss errors below $0.313\%$. The classical far-field and Huygens-surface powers differ by $6.64\%$, requiring electromagnetic convergence studies. Finite-shot estimation, gate-synthesis precision, fault-tolerant resources, and end-to-end complexity are not evaluated. Observable formulation alone does not establish efficient measurement or quantum advantage; matrix access, state preparation, preconditioning, spectral bounds, and readout costs remain to be assessed.

Little sign of AI-generated textIncrementalSmall numerics onlyStrong assumptions

Central claim may not hold: The abstract claims efficient evaluation of antenna observables, but the paper establishes quadratic-form representations and ideal-simulator consistency without an efficient quantum measurement construction or end-to-end complexity bound.

Written by gpt-6.1-sol from the paper's full text. Automated: it can be wrong.

Comments
10 pages, 5 figures

Related papers

Solution of the Electric Field Integral Equation Using a Hybrid Quantum-Classical Scheme: Investigation of Accuracy and Efficiency

Rui Chen, Teng-Yang Ma, Meng-Han Dou +1 more

Conventional classical solvers are commonly used for solving matrix equation systems resulting from the discretization of SIEs in computational electromagnetics (CEM). However, the memory requirement would become a bottleneck for classical computing as the electromagentic problems become much larger. As an alternative, quantum computing has a natural "parallelization" advantage with much lower storage complexity due to the superposition and entanglement in quantum mechanics. Even though several quantum algorithms have been applied for the SIEs-based methods in the literature, the size of the matrix equation systems solvable using them is still limited. In this work, we use a hybrid quantum-classical scheme to solve the EFIE for analyzing electromagentic scattering from three-dimensional (3D) perfect electrically conducting objects with arbitrary shapes in CEM for the first time. Instead of directly solving the original EFIE matrix equation system using the quantum algorithms, the hybrid scheme first designs the preconditioned linear system and then uses a double-layer iterative strategy for its solution, where the external iteration layer builds subspace matrix equation systems with smaller dimension and the internal iteration layer solves the smaller systems using the quantum algorithms. Two representative quantum algorithms, HHL and VQLS, are considered in this work, which are executed on the quantum simulator and quantum computer platforms. We present the theoretical time complexity analysis of the hybrid quantum-classical scheme and perform numerical experiments to investigate the accuracy and efficiency of the hybrid scheme. The results show that the computational complexity of the hybrid VQLS-classical scheme is lower than the conventional fast solvers in classical computing, which indicates the hybrid scheme is more promising for analyzing large-scale electromagnetic problems.

quant-pharXiv:2512.03808•Dec 4, 2025•0 citations

End-to-end performance of quantum-accelerated large-scale linear algebra workflows

Daiwei Zhu, Miguel Angel Lopez-Ruiz, François-Henry Rouet +5 more

Solving large-scale sparse linear systems is a challenging computational task due to the introduction of non-zero elements, or "fill-in". The Graph Partitioning Problem (GPP) arises naturally when minimizing fill-in and accelerating solvers. In this paper, we measure the end-to-end performance of a hybrid quantum-classical framework designed to accelerate Finite Element Analysis (FEA) by integrating a quantum solver for GPP into Synopsys/Ansys' LS-DYNA multiphysics simulation software. The quantum solver we use is based on Iterative-QAOA, a scalable, non-variational quantum approach for optimization. We focus on two specific classes of FEA problems, namely vibrational (eigenmode) analysis and transient simulation. We report numerical simulations on up to 150 qubits done on NVIDIA's CUDA-Q/cuTensorNet and implementation on IonQ's Forte quantum hardware. The potential impact on LS-DYNA workflows is quantified by measuring the wall-clock time-to-solution for complex problem instances, including vibrational analysis of large finite element models of a sedan car and a Rolls-Royce jet engine, as well as transient simulations of a drill and an impeller. We performed end-to-end performance measurements on meshes comprising up to 35 million elements. Measurements were conducted using LS-DYNA in distributed-memory mode via Message Passing Interface (MPI) on AWS and Synopsys compute clusters. Our findings indicate that with a quantum computer in the loop, amortized LS-DYNA wall-clock time can be improved by up to 14.6% for specific cases and by at least 5.9% for all models considered. These results highlight the significant potential of quantum computing to reduce time-to-solution for large-scale FEA simulations within the Noisy Intermediate-Scale Quantum (NISQ) era, offering an approach that is scalable and extendable into the fault-tolerant quantum computing regime.

quant-pharXiv:2603.15515•Mar 17, 2026•0 citations

Quantum Method of Calculating Multiport S-Parameter Magnitudes Based on FEM for Microwave Components: An HHL Framework

Xiaolong Li, Feng Feng, Jiali Zhang +5 more

Quantum computing in the electromagnetic (EM) field represents a transformative frontier with significant potential for advancement. Recent studies have demonstrated the potential of quantum computing for solving waveguide mode problems, solving finite element method (FEM) equations for EM problems, and scaling the FEM solution. Among these applications, the Harrow–Hassidim–Lloyd (HHL) algorithm is an algorithm with large quantum advantages and has attracted attention. However, existing quantum methods remain constrained by the problem that the quantum-based FEM solution only obtains electric field magnitude and lacks phase information. The lack of phase information results in the inability to solve an important characteristic of EM problems: scattering parameters (S-parameters). To address the above issue, this work proposes a novel quantum method of calculating multiport S-parameter magnitudes based on FEM and HHL algorithm for microwave components. The proposed method introduces novel mathematical formulations optimized for quantum architectures. To streamline the quantum computation process and avoid the need for multiple quantum circuit simulations for each S-parameter, we design an integrated quantum circuit capable of computing magnitudes of all multiport S-parameters within a single execution. Furthermore, a specific matrix-scaling method is also developed to simplify the proposed circuits. The specific matrix-scaling method is used to improve the condition of the FEM matrix and reduce the complexity of the designed quantum circuits. The proposed method is validated using two microwave structures, demonstrating how our proposed quantum method computes S-parameter magnitudes in practical EM applications.

IEEE Transactions on Microwave Theory and Techniques•Nov 5, 2025•4 citations

Quantum Computing for Industrial Electromagnetics: Applicability and Case Studies in Solving Maxwell's Equations

Francesco Turro, Marco Maronese, Daniele Dragoni

Computational electromagnetics plays a central role in many industrial applications but often requires substantial computational resources, particularly when fine spatial discretizations are needed. While classical approaches remain the standard, quantum computing offers the potential to accelerate large-scale simulations by encoding them with a limited number of qubits. Here, we investigate the performance and resource scaling of the Harrow-Hassidim-Lloyd (HHL) and Quantum Singular Value Transformation (QSVT) algorithms for solving linear systems generated by the finite-difference time-domain (FDTD) method, a widely adopted numerical scheme for discretizing Maxwell's equations. We benchmark their performance across representative industrial use cases, including radar propagation, lens simulations, and beamforming processes. Our results demonstrate the validity of the approaches, achieving state infidelities smaller than $2\cdot 10^{-3}$ with success probabilities greater than $10^{-3}$, compatible with practical quantum state sampling. Overall, we observe that the QSVT method consistently delivers higher accuracy. We further observe that the condition number of the linear matrix, a key factor governing the performance of quantum solvers, saturates as the number of spatial lattice points increases. This implies that the spatial grid can be scaled to realistic industrial dimensions without increasing the HHL or QSVT circuit depth due to ill-conditioned matrices.

quant-pharXiv:2608.10671•Aug 12, 2026•0 citations

Quantum-Based Methods for Solving Partial Differential Equations in Electromagnetics Engineering

Emanuel Colella, Luca Bastianelli, Valter Mariani Primiani +2 more

Quantum methods for solving partial differential equations (PDEs) generally fall into a few general categories, each leveraging different quantum mechanisms to achieve potential speedups. One well-known approach is the quantum linear system solver, epitomized by the Harrow-Hassidim-Lloyd (HHL) algorithm, which uses phase estimation and eigenvalue inversion to solve large, sparse linear systems arising from discretized linear PDEs [1]. This method can yield exponential improvements over classical solvers under certain conditions (mostly when the matrix is well-conditioned and sparse). More recent developments, such as Quantum Singular Value Transformation (QSVT), provide a framework for matrix inversion in a more flexible way, often improving on the cost and error dependence of earlier solvers [2]. Time-dependent PDEs can be tackled via Hamiltonian simulations. The Trotter-Suzuki method decomposes the exponential of a Hamiltonian into products of exponentials of simpler sub-Hamiltonians, allowing a piecewise approximation of continuous dynamics [3]. Continuous-time quantum walks, which treat a PDE operator as a Hamiltonian that governs the walk, similarly simulate processes like diffusion or wave propagation on discretized grids or graphs [4]. Beyond direct simulation, variational quantum algorithms have emerged as a promising route in the noisy intermediate-scale quantum era. These approaches encode the PDE solution as a parameterized quantum state and optimize the parameters via classical feedback, thereby minimizing a cost function tied to the PDE residual or energy functional [5]. Quantum annealing or adiabatic methods, meanwhile, can solve optimization problems derived from PDE constraints by adiabatically evolving a simple initial Hamiltonian to the final problem Hamiltonian and staying in the system's ground state [3]. In our study, we implement a new parameterized quantum model, specifically designed to solve large linear systems on quantum hardware. This model is less deep and computationally faster than the proposed one in [5], reducing the complexity and speeding up convergence times. Such improvements are particularly valuable in cases where classical algorithms may fail due to the problem's size or complexity. Ongoing research is focusing on generalized quantum models that can be adapted to practical PDE problems for solving large electromagnetic problems where classical computers can fail.

USNC-URSI Radio Science Meeting•Jul 13, 2025•0 citations

A Cost-Efficient Approach towards Computational Fluid Dynamics Simulations on Quantum Devices

Szabolcs Jóczik, Zoltán Zimborás, Tamás Majoros +1 more

Numerical simulations of physical systems are found in many industries, as they currently play a crucial role in product development. There are many numerical methods for solving differential equations that describe the underlying physics behind the mathematical models in the simulation, among which, the finite element method (FEM) is one of the most commonly used. Although in many applications the FEM seems to provide an acceptable solution to the problem, there are still many complex real-life processes that can be challenging to simulate numerically due to their complexity and large size. Recently, there has been a shift in research towards efficiently applying quantum algorithms in finite element analysis (FEA), as the potential and speedup that they could offer have been shown, but little to no effort has been made towards the applicability and cost efficiency of these algorithms in real-world quantum devices. In this paper, we propose a cost-efficient method for applying quantum algorithms in FEA for industrial problems post-processed by classical algorithms in order to address the limitations of available quantum hardware and their cost when accessing them through different cloud-based services. We carry this out by approximating the solution of the initially large system with a suitable quantum algorithm and using the obtained solutions to generate a set of reduced-order models (ROMs) that are much smaller in complexity and size than the original model. This allows the simulation of the original model with different parameter sets and excitations to be run efficiently on classical computers without having the need to access quantum subroutines again. This way, we have reduced the usage of quantum hardware (and thus the development cost) while still taking advantage of its quantum speedup.

Applied Sciences•Mar 10, 2022•12 citations
Electromagnetic analysis based on the finite element method is widely used in antenna design. At higher frequencies, however, the required mesh refinement makes solving large sparse linear systems the main contributor to computational cost. Quantum linear system solvers offer a potential means of accelerating this step, but their output is a superposition state: reconstructing all components of the classical solution vector would forfeit the quantum speedup. Prior studies have proposed methods for obtaining physical quantities such as scattering cross sections and scattering parameters from the solution state. To support antenna design using fault-tolerant quantum computers, we formulate the total radiated power, the radiated power within a prescribed solid angle, and the dielectric and conductive losses as observables for efficient evaluation from a quantum solution state. These observables avoid the measurement cost of reconstructing all components of the output quantum state. We develop an end-to-end hybrid quantum--classical workflow that combines finite-element electromagnetic analysis, matrix preconditioning, a quantum linear system solver, and the evaluation of antenna design metrics. To validate the workflow, we performed numerical experiments evaluating these observables for an axisymmetric coaxially fed monopole antenna. Under idealized conditions, the relative errors in the total radiated power and the dissipated power with respect to the classical finite-element solution were on the order of 10−310^{-3}10−3, confirming the consistency of the workflow. The proposed method is expected to provide a basis for using quantum linear system solvers to evaluate antenna performance and, in the future, to design high-frequency wireless devices.
The paper asks how antenna performance metrics can be extracted from a quantum linear-system solver’s normalized finite-element solution state without reconstructing the full electromagnetic field. It expresses total and solid-angle-restricted radiated power, Huygens-surface power flux, and dielectric and conductive losses as quadratic forms x†Qxx^\dagger Qxx†Qx. Radiation observables are constructed from far-field linear functionals; loss observables follow from volume and surface integrals of finite-element basis fields.
The workflow combines Ruiz diagonal equilibration, Hermitian embedding of the complex indefinite Maxwell matrix, and HHL-type solution-state preparation. Observables are transformed consistently with preconditioning, and physical power scales are recovered using the solution norm inferred from postselection probability. Validation uses idealized classical simulation for an axisymmetric coaxially fed monopole with 519 complex degrees of freedom. Equilibration reduces the condition number from 6.16×1056.16\times10^56.16×105 to 9.34×1039.34\times10^39.34×103; the simulation uses 14 phase qubits and has postselection probability 1.29×10−51.29\times10^{-5}1.29×10−5 at the reference frequency. Relative errors against the classical discrete solution are 0.769%0.769\%0.769% for total radiated power and 0.662%0.662\%0.662% for angular-band power. Across a 21-frequency sweep, radiated-power errors remain below 0.693%0.693\%0.693%, and total-loss errors below 0.313%0.313\%0.313%.
The classical far-field and Huygens-surface powers differ by 6.64%6.64\%6.64%, requiring electromagnetic convergence studies. Finite-shot estimation, gate-synthesis precision, fault-tolerant resources, and end-to-end complexity are not evaluated. Observable formulation alone does not establish efficient measurement or quantum advantage; matrix access, state preparation, preconditioning, spectral bounds, and readout costs remain to be assessed.
Conventional classical solvers are commonly used for solving matrix equation systems resulting from the discretization of SIEs in computational electromagnetics (CEM). However, the memory requirement would become a bottleneck for classical computing as the electromagentic problems become much larger. As an alternative, quantum computing has a natural "parallelization" advantage with much lower storage complexity due to the superposition and entanglement in quantum mechanics. Even though several quantum algorithms have been applied for the SIEs-based methods in the literature, the size of the matrix equation systems solvable using them is still limited. In this work, we use a hybrid quantum-classical scheme to solve the EFIE for analyzing electromagentic scattering from three-dimensional (3D) perfect electrically conducting objects with arbitrary shapes in CEM for the first time. Instead of directly solving the original EFIE matrix equation system using the quantum algorithms, the hybrid scheme first designs the preconditioned linear system and then uses a double-layer iterative strategy for its solution, where the external iteration layer builds subspace matrix equation systems with smaller dimension and the internal iteration layer solves the smaller systems using the quantum algorithms. Two representative quantum algorithms, HHL and VQLS, are considered in this work, which are executed on the quantum simulator and quantum computer platforms. We present the theoretical time complexity analysis of the hybrid quantum-classical scheme and perform numerical experiments to investigate the accuracy and efficiency of the hybrid scheme. The results show that the computational complexity of the hybrid VQLS-classical scheme is lower than the conventional fast solvers in classical computing, which indicates the hybrid scheme is more promising for analyzing large-scale electromagnetic problems.
Solving large-scale sparse linear systems is a challenging computational task due to the introduction of non-zero elements, or "fill-in". The Graph Partitioning Problem (GPP) arises naturally when minimizing fill-in and accelerating solvers. In this paper, we measure the end-to-end performance of a hybrid quantum-classical framework designed to accelerate Finite Element Analysis (FEA) by integrating a quantum solver for GPP into Synopsys/Ansys' LS-DYNA multiphysics simulation software. The quantum solver we use is based on Iterative-QAOA, a scalable, non-variational quantum approach for optimization. We focus on two specific classes of FEA problems, namely vibrational (eigenmode) analysis and transient simulation. We report numerical simulations on up to 150 qubits done on NVIDIA's CUDA-Q/cuTensorNet and implementation on IonQ's Forte quantum hardware. The potential impact on LS-DYNA workflows is quantified by measuring the wall-clock time-to-solution for complex problem instances, including vibrational analysis of large finite element models of a sedan car and a Rolls-Royce jet engine, as well as transient simulations of a drill and an impeller. We performed end-to-end performance measurements on meshes comprising up to 35 million elements. Measurements were conducted using LS-DYNA in distributed-memory mode via Message Passing Interface (MPI) on AWS and Synopsys compute clusters. Our findings indicate that with a quantum computer in the loop, amortized LS-DYNA wall-clock time can be improved by up to 14.6% for specific cases and by at least 5.9% for all models considered. These results highlight the significant potential of quantum computing to reduce time-to-solution for large-scale FEA simulations within the Noisy Intermediate-Scale Quantum (NISQ) era, offering an approach that is scalable and extendable into the fault-tolerant quantum computing regime.
Quantum computing in the electromagnetic (EM) field represents a transformative frontier with significant potential for advancement. Recent studies have demonstrated the potential of quantum computing for solving waveguide mode problems, solving finite element method (FEM) equations for EM problems, and scaling the FEM solution. Among these applications, the Harrow–Hassidim–Lloyd (HHL) algorithm is an algorithm with large quantum advantages and has attracted attention. However, existing quantum methods remain constrained by the problem that the quantum-based FEM solution only obtains electric field magnitude and lacks phase information. The lack of phase information results in the inability to solve an important characteristic of EM problems: scattering parameters (S-parameters). To address the above issue, this work proposes a novel quantum method of calculating multiport S-parameter magnitudes based on FEM and HHL algorithm for microwave components. The proposed method introduces novel mathematical formulations optimized for quantum architectures. To streamline the quantum computation process and avoid the need for multiple quantum circuit simulations for each S-parameter, we design an integrated quantum circuit capable of computing magnitudes of all multiport S-parameters within a single execution. Furthermore, a specific matrix-scaling method is also developed to simplify the proposed circuits. The specific matrix-scaling method is used to improve the condition of the FEM matrix and reduce the complexity of the designed quantum circuits. The proposed method is validated using two microwave structures, demonstrating how our proposed quantum method computes S-parameter magnitudes in practical EM applications.
Computational electromagnetics plays a central role in many industrial applications but often requires substantial computational resources, particularly when fine spatial discretizations are needed. While classical approaches remain the standard, quantum computing offers the potential to accelerate large-scale simulations by encoding them with a limited number of qubits. Here, we investigate the performance and resource scaling of the Harrow-Hassidim-Lloyd (HHL) and Quantum Singular Value Transformation (QSVT) algorithms for solving linear systems generated by the finite-difference time-domain (FDTD) method, a widely adopted numerical scheme for discretizing Maxwell's equations. We benchmark their performance across representative industrial use cases, including radar propagation, lens simulations, and beamforming processes. Our results demonstrate the validity of the approaches, achieving state infidelities smaller than 2⋅10−32\cdot 10^{-3}2⋅10−3 with success probabilities greater than 10−310^{-3}10−3, compatible with practical quantum state sampling. Overall, we observe that the QSVT method consistently delivers higher accuracy. We further observe that the condition number of the linear matrix, a key factor governing the performance of quantum solvers, saturates as the number of spatial lattice points increases. This implies that the spatial grid can be scaled to realistic industrial dimensions without increasing the HHL or QSVT circuit depth due to ill-conditioned matrices.
Quantum methods for solving partial differential equations (PDEs) generally fall into a few general categories, each leveraging different quantum mechanisms to achieve potential speedups. One well-known approach is the quantum linear system solver, epitomized by the Harrow-Hassidim-Lloyd (HHL) algorithm, which uses phase estimation and eigenvalue inversion to solve large, sparse linear systems arising from discretized linear PDEs [1]. This method can yield exponential improvements over classical solvers under certain conditions (mostly when the matrix is well-conditioned and sparse). More recent developments, such as Quantum Singular Value Transformation (QSVT), provide a framework for matrix inversion in a more flexible way, often improving on the cost and error dependence of earlier solvers [2]. Time-dependent PDEs can be tackled via Hamiltonian simulations. The Trotter-Suzuki method decomposes the exponential of a Hamiltonian into products of exponentials of simpler sub-Hamiltonians, allowing a piecewise approximation of continuous dynamics [3]. Continuous-time quantum walks, which treat a PDE operator as a Hamiltonian that governs the walk, similarly simulate processes like diffusion or wave propagation on discretized grids or graphs [4]. Beyond direct simulation, variational quantum algorithms have emerged as a promising route in the noisy intermediate-scale quantum era. These approaches encode the PDE solution as a parameterized quantum state and optimize the parameters via classical feedback, thereby minimizing a cost function tied to the PDE residual or energy functional [5]. Quantum annealing or adiabatic methods, meanwhile, can solve optimization problems derived from PDE constraints by adiabatically evolving a simple initial Hamiltonian to the final problem Hamiltonian and staying in the system's ground state [3]. In our study, we implement a new parameterized quantum model, specifically designed to solve large linear systems on quantum hardware. This model is less deep and computationally faster than the proposed one in [5], reducing the complexity and speeding up convergence times. Such improvements are particularly valuable in cases where classical algorithms may fail due to the problem's size or complexity. Ongoing research is focusing on generalized quantum models that can be adapted to practical PDE problems for solving large electromagnetic problems where classical computers can fail.
Numerical simulations of physical systems are found in many industries, as they currently play a crucial role in product development. There are many numerical methods for solving differential equations that describe the underlying physics behind the mathematical models in the simulation, among which, the finite element method (FEM) is one of the most commonly used. Although in many applications the FEM seems to provide an acceptable solution to the problem, there are still many complex real-life processes that can be challenging to simulate numerically due to their complexity and large size. Recently, there has been a shift in research towards efficiently applying quantum algorithms in finite element analysis (FEA), as the potential and speedup that they could offer have been shown, but little to no effort has been made towards the applicability and cost efficiency of these algorithms in real-world quantum devices. In this paper, we propose a cost-efficient method for applying quantum algorithms in FEA for industrial problems post-processed by classical algorithms in order to address the limitations of available quantum hardware and their cost when accessing them through different cloud-based services. We carry this out by approximating the solution of the initially large system with a suitable quantum algorithm and using the obtained solutions to generate a set of reduced-order models (ROMs) that are much smaller in complexity and size than the original model. This allows the simulation of the original model with different parameter sets and excitations to be run efficiently on classical computers without having the need to access quantum subroutines again. This way, we have reduced the usage of quantum hardware (and thus the development cost) while still taking advantage of its quantum speedup.