Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “integrator”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A machine learning approach for efficient multi-dimensional integration

Many physics problems involve integration in multi-dimensional space whose analytic solution is not available. The integrals can be evaluated using numerical integration methods, but it requires a large computational cost in some cases, so an efficient algorithm plays an important role in solving the physics problems. We propose a novel numerical multi-dimensional integration algorithm using machine learning (ML). After training a ML regression model to mimic a target integrand, the regression model is used to evaluate an approximation of the integral. Then, the difference between the approximation and the true answer is calculated to correct the bias in the approximation of the integral induced by ML prediction errors. Because of the bias correction, the final estimate of the integral is unbiased and has a statistically correct error estimation. Three ML models of multi-layer perceptron, gradient boosting decision tree, and Gaussian process regression algorithms are investigated. The performance of the proposed algorithm is demonstrated on six different families of integrands that typically appear in physics problems at various dimensions and integrand difficulties. The results show that, for the same total number of integrand evaluations, the new algorithm provides integral estimates with more than an order of magnitude smaller uncertainties than those of the VEGAS algorithm in most of the test cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Toward a Fully Integrated Multiphysics Simulation Framework for Fusion Blanket Design

Fusion is an attractive clean-energy solution, thanks to its various advantages, such as reduced radioactivity, little high-level nuclear waste, ample fuel supplies, and increased safety. However, the harsh operating environment introduced by a complex fusion plasma system makes design and integration of fusion blankets incredibly challenging and time-consuming. This work focuses on developing a fully integrated multiphysics simulation framework based on an advanced open-source platform—the Multiphysics Object-Oriented Simulation Environment (MOOSE)—to alleviate the difficulties in fusion blanket design and integration. MOOSE is a massively parallel finite element/volume multiphysics simulation platform that has been widely adopted within the nuclear fission community. Even though fission and fusion are fundamentally different, they involve similar multiphysics phenomena. A fully integrated open-source multiphysics simulation framework tailored for the fusion blanket design will be implemented by leveraging the well-established multiphysics capabilities in MOOSE. Once successfully developed, this fully integrated framework will rapidly evaluate a blanket design concept and offer insights for subsequent iterations. As the first step, we will mainly aim to integrate neutronics analysis, system thermal hydraulics simulation, and full 3-D heat transfer calculations. The efficacy of the integrated framework will be verified using an innovative solid ceramic blanket design. While the project’s final goal is to enable a fully integrated multiphysics simulation platform for various fusion blanket concepts, here this work, as a preliminary step, will mainly focus on a solid ceramic breeder helium-cooled blanket.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Filon-like integration strategy for calculating exact exchange in periodic boundary conditions: a plane-wave DFT implementation

An efficient and accurate approach for calculating exact exchange and other two-electron integrals has been developed for periodic electronic structure methods. Traditional approaches used for integrating over the Brillouin zone in band structure calculations, e.g. trapezoidal or Monkhorst-Pack, are not accurate enough for two-electron integrals. This is because their integrands contain multiple singularities over the double integration of the Brillouin zone, which with simple integration methods lead to very inaccurate results. A common approach to this problem has been to replace the Coulomb interaction with a screened Coulomb interaction that removes singularities from the integrands in the two-electron integrals, albeit at the inelegance of having to introduce a screening factor which must precomputed or guessed. Instead of introducing screened Coulomb interactions in an ad hoc way, the method developed in this work derives an effective screened potential using a Filon-like integration approach that is based only on the lattice parameters. This approach overcomes the limitations of traditionally defined screened Coulomb interactions for calculating two-electron integrals, and makes chemistry many-body calculations tractable in periodic boundary conditions. This method has been applied to several systems for which conventional DFT methods do not work well, including the reaction pathways for the addition of H 2 to phenol and Au$^{–}_{20}$ nanoparticle, and the electron transfer of a charge trapped state in the Fe(II) containing mica, annite.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ion Traps and Packaging for Heterogenous Integration - Chimera

Microfabricated surface ion traps and silicon-based photonics are critical technologies for scaling quantum systems. Current ion trap architectures face scalability and integration challenges due to limitations in optical access, fabrication techniques, and material compatibility. State-of-the-art quantum computers and atomic clocks are investigating monolithic integration, which necessitates custom traps for each ion species and has not overcome the integration hurdles presented by merging these technologies. The Chimera (Ion Traps and Packaging for Heterogeneous Integration) project proposes a novel approach utilizing heterogeneous integration (HI) of ion traps and photonic circuits. This separation of components allows for flexibility in ion trap design and reduces fabrication compromises. The Chimera project specifically designed an ion trap to interface vertically with a separately fabricated waveguide chip and demonstrates the first steps to integrating them at the packaging level. The ion trap features a large area of removed silicon, allowing the photonics chip outputs closer to the ion trap, improving alignment and packaging processes. The alignment must be accurate to < 1 µm to ensure that the light from the waveguide can overlap with the trapping region. This fine alignment must also be maintained through an ultra-high vacuum bake, a critical step in preparing an ion trap experiment. By combining separate chips, we demonstrate a new path for scaling trapped ion technology that is less reliant on monolithic integration. We successfully fabricated a trap with a large area of oxide removed, resulting in a region thinned to about 40 µm, a key milestone toward successful integration.

42 ENGINEERING↗
DOE OSTI2022

m-CUBES An efficient and portable implementation of multi-dimensional integration for gpus

The task of multi-dimensional numerical integration is frequently encountered in physics and other scientific fields, e.g., in modeling the effects of systematic uncertainties in physical systems and in Bayesian parameter estimation. Multi-dimensional integration is often time-prohibitive on CPUs. Efficient implementation on many-core architectures is challenging as the workload across the integration space cannot be predicted a priori. We propose m-Cubes, a novel implementation of the well-known Vegas algorithm for execution on GPUs. Vegas transforms integration variables followed by calculation of a Monte Carlo integral estimate using adaptive partitioning of the resulting space. m-Cubes improves performance on GPUs by maintaining relatively uniform workload across the processors. As a result, our optimized Cuda implementation for Nvidia GPUs outperforms parallelization approaches proposed in past literature. We further demonstrate the efficiency of m-Cubes by evaluating a six-dimensional integral from a cosmology application, achieving significant speedup and greater precision than the CUBA library's CPU implementation of VEGAS. We also evaluate m-Cubes on a standard integrand test suite. m-Cubes outperforms the serial implementations of the Cuba and GSL libraries by orders of magnitude speedup while maintaining comparable accuracy. Our approach yields a speedup of at least 10 when compared against publicly available Monte Carlo based GPU implementations. In summary, m-Cubes can solve integrals that are prohibitively expensive using standard libraries and custom implementations. A modern C++ interface header-only implementation makes m-Cubes portable, allowing its utilization in complicated pipelines with easy to define stateful integrals. Compatibility with non-Nvidia GPUs is achieved with our initial implementation of m-Cubes using the Kokkos framework.

Sakiotis, Ioannis↗
DOE OSTI2024

On-PIC Light Source Integration & Micro-dispensing of Solder Paste for Flip Chip Application

To reach the next level benefits of photonic integrated circuits (PICs), the Integrated Photonic Systems Roadmap-International describes the necessity of either heterogeneous or hybrid integration of light sources [1]. A new approach for hybrid integration is Photonic Wire Bonding where 3D nanolithography is used to pattern a polymer waveguide that connects light sources to PIC waveguides. The waveguides resemble electrical wire bonds and essentially do the same as their electrical counterpart for packaging photonic chips together. In order to reap the full benefits of a photonic wire bonds, the light source must be carefully packaged. In this work I seek to expand RIT’s photonic packaging capability by establishing a packaging process to integrate light sources, specifically an Indium Phosphide distributed feedback lasers and a reflective semiconductor optical amplifiers, directly onto a photonic integrated circuits. Another necessary requirement for PIC packaging is densely integrated electrical connectivity. In this thesis I developed a fine pitch flip chip interconnect technique demonstrated using gold stud bumps in conjunction with micro-dispensed solder paste. This work opens the door to high density photonic flip chip applications. The micro-dispensed solder paste dots having diameters 50-125 µm are currently being tested at a pitch of 150 µm . The gold stud bumps are formed with 1mil gold wire creating bumps with a diameter of 40-60 µm depending on specific parameter values. These capabilities will allow RIT to assemble and test novel photonic integrated devices, cutting down on the time and cost associated with third party assembly and tests facilities. Thereby keeping RIT at the forefront of photonic research.

Wongk, Nicole↗
DOE OSTI2021

A Review of Well Integrity Based on Field Experience at Carbon Utilization and Storage Sites

Maintaining well integrity is critical to the success of geologic carbon storage (GCS) and carbon dioxide enhanced oil recovery (CO 2 -EOR) operations. Wells that experience leakage because of integrity issues can potentially become a risk to the environment or human health if they release previously captured CO 2 back into the atmosphere or into freshwater aquifers. There are many GCS and CO 2 -EOR sites in operation around the world. However, well integrity experiences at these sites are not widely documented in the public domain. This study details findings from a survey of well integrity experiences elicited from operators of GCS and CO 2 -EOR sites. The goal of the survey was to obtain information about site characteristics and operator experiences with well integrity, monitoring methods, and risk assessment of legacy wells. Literature relevant to the survey questions was also reviewed and summarized to provide context for survey responses and identify areas where field experiences with well integrity do and do not align with the current state of research. We highlight the current state-of-practice, identify research needs, and provide context for future interactions between researchers, operators, and regulators on issues related to well integrity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A generalized class of strongly stable and dimension-free T-RPMD integrators

Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.

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

Design of Quasi-Integral Experiments at the Flexible Neutron Source at the University of Tennessee

Nuclear reaction data libraries make use of experimental data that are differential in reaction and neutron energy for evaluation and data that are integral in reaction and energy for validation. In many cases, integral measurements can provide lower experimental uncertainties than differential measurements, but the impact of individual reactions and energies cannot generally be determined from the single experimental observable. “Quasi-integral” or “quasi-differential” measurements, which are integral in one characteristic (energy or reaction) and differential in the other, can often provide a middle ground of lower experimental uncertainties combined with some ability to unpack the impact of individual reactions or energies. Furthermore, the Flexible Neutron Source (FNS) at the University of Tennessee, Knoxville, is a new experimental facility that can perform integral and quasi-integral measurements by removing components of the integral configuration to determine the impact of specific reactions. Possible flux-averaged cross-section measurements, flux perturbation measurements, and neutron downscattering measurements are simulated here to show the value of the FNS to evaluation and validation of nuclear data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Twistor Actions for Integrable Systems

Many integrable systems can be reformulated as holomorphic vector bundles on twistor space. This is a powerful organizing principle in the theory of integrable systems. One shortcoming is that it is formulated at the level of the equations of motion. From this perspective, it is mysterious that integrable systems have Lagrangians. In this paper, we study a Chern-Simons action on twistor space and use it to derive the Lagrangians of some integrable sigma models. Our focus is on examples that come from dimensionally reduced gravity and supergravity. The dimensional reduction of general relativity to two spacetime dimensions is an integrable coset sigma model coupled to a dilaton and 2d gravity. The dimensional reduction of supergravity to two spacetime dimensions is an integrable coset sigma model coupled to matter fermions, a dilaton, and 2d supergravity. We derive Lax operators and Lagrangians for these 2d integrable systems using the Chern-Simons theory on twistor space. In the supergravity example, we use an extended setup in which twistor Chern-Simons theory is coupled to a pair of matter fermions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Nyström type exponential integrators for strongly magnetized charged particle dynamics

Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.

general physics↗
DOE OSTI2025

PDB-IHM: A System for Deposition, Curation, Validation, and Dissemination of Integrative Structures

Structures of many large biomolecular assemblies are now being determined using integrative approaches. In these approaches, information derived from multiple experimental and computational methods is combined to compute three-dimensional structures of multi-protein complexes and other macromolecular machines. A standalone prototype data resource for integrative structures called PDB-Dev was built, based on recommendations of the Integrative and Hybrid Methods (IHM) Task Force of the Worldwide Protein Data Bank (wwPDB). This effort included developing data standards and software tools for collecting, curating, validating, visualizing, archiving, and disseminating integrative structures that span diverse spatiotemporal scales and conformational states. Mechanisms have been created to validate integrative structures based on the experimental data underpinning them. Building upon this foundational framework, PDB-Dev has been further expanded to handle large dynamic macromolecular systems and integrative structures that combine, for example, experimental restraints with atomic coordinates computed by machine learning algorithms. Data standards and supporting tools have also been extended to capture information about biomolecular dynamics, such as conformational transitions and related kinetic data derived from biophysical methods. Recently, PDB-Dev was unified with the PDB archive and rebranded as PDB-IHM (pdb-ihm.org), further promoting FAIR (Findable, Accessible, Interoperable, and Reusable) principles of data stewardship for integrative structural biology.

IHMCIF↗
DOE OSTI2022

Integrated control of individual plasma scalars with simultaneous neoclassical tearing-mode suppression

A novel integrated-control architecture has been tested in nonlinear, one-dimensional simulations using the control-oriented transport simulator (COTSIM©) and in DIII-D experiments. Integrated architectures that can perform continuous-mission control while also handling off-normal events will be vital in future reactor-grade tokamaks. Continuous-mission controllers for individual magnetic and kinetic scalars (thermal stored-energy (W), volume-average toroidal rotation (Ω Φ ), and safety factor profile (q) at different spatial locations) have been integrated in this work with event-triggered neoclassical tearing-mode (NTM) suppression controllers by combining them into an architecture augmented by a supervisory and exception handling (S&EH) system and an actuator management (AM) system. Here, the AM system, which enables the integration of competing controllers, solves in real time a nonlinear optimization problem that takes into account the high-level control priorities dictated by the S&EH system. The resulting architecture offers a high level of integration and some of the functionalities that will be required to fulfill the advanced-control requirements anticipated for ITER. Initial simulations using COTSIM suggest that the plasma performance and its MHD stability may be improved under integrated feedback control. In addition, the integrated-control architecture has been implemented in the DIII-D plasma control system and tested experimentally for the first time ever in DIII-D in a high-q min scenario, which is a candidate for steady-state operation in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Technical and Economic Assessment and Gap Analysis of Advanced Nuclear Reactor Integration with a Reference Oil Refinery

Efforts to identify the most-economic methods to decarbonize several sectors of the U.S. economy are underway. Industrial processes such as crude-oil refining rely heavily on energy-dense and easily stored and transported fossil fuels for powering their operations. Refineries use large amounts of energy, primarily derived from fossil sources to separate crude-oil components, break down heavier hydrocarbons into lighter compounds, remove impurities, reform hydrocarbon molecules, and generate steam and electricity for pumps and compressors and other various auxiliary systems. Crude-oil refining operations such as distillation, cracking, desulfurization, reforming, utilities systems and some offsite facilities collectively account for most of the energy consumption. Other operations such as hydrocracking or hydrotreating also require hydrogen for developing hydrogenation reactions which involve substantial heating to keep the reactors at high-temperature and pressure levels. All heat and energy demands are typically provided by natural gas (NG), oil, or other fuels, which makes refinery industry one of the most-difficult sectors to decarbonize. Nuclear power is a viable and energy-dense source of clean electricity, heat, and hydrogen to provide the large, sustainable energy supply that the refining industry demands. The U.S. Department of Energy’s (DOE’s) Integrated Energy Systems (IES) program is working to perform research and development, design, economic siting, and risk analysis. This state-of-the-art work will enable the first on-site demonstrations and commercial deployments of advanced small modular nuclear reactors (SMNRs) integrated with industries such as chemical production, refining, iron and steel making, and more. IES seeks to demonstrate the ability of advanced nuclear reactors to meet the heat and power demands of these industries while reducing carbon emissions in a sustainable and cost-competitive way. The primary objective of this research effort is to analyze industrial-scale SMNR integration intended to decarbonize refining facilities. The foreseen outcome is the provision of reliable, cost-competitive, and sustainable clean energy, alongside a reduction of carbon emissions. Specifically, the focus of this work lies on meeting the reference facilities’ heat and electricity demands with nuclear power while also supplying clean hydrogen via integrated high-temperature steam electrolysis (HTSE). This report presents a comprehensive technical and economic assessment of the integration of advanced nuclear reactors into a reference refinery, leveraging financial incentives from the Inflation Reduction Act (IRA). The evaluation aims to explore the potential economic benefits and challenges associated with incorporating advanced nuclear reactors into refinery operations, particularly in terms of energy efficiency, economic implications and environmental impact. By examining both the technical feasibility and economic viability, this analysis seeks to identify existing gaps and propose solutions for successful nuclear integration implementation. The findings are intended to provide valuable insights for stakeholders considering the adoption of advanced nuclear reactors in the refining sector. A refinery reference-plant was developed, using an open-source refinery model, Petroleum Refinery Lifecycle Inventory Model (PRELIM) and expert assessment, as a base case for comparison with various nuclear integration options. The capacity of 100 kbd/day (KBD) of heavy crude-oil feed was selected to represent a general coking-type refinery with deep conversion capabilities (incorporating heavy-oil upgrading with FCC, coking, and associated hydrotreating process units), using a heavy crude-oil feed, which represents about 70% of U.S. refineries configurations. A summary of all cases considered in this study is shown in Table 1.

13 HYDRO ENERGY↗
DOE OSTI2020

Integrated Transmission-and-Distribution System Modeling of Power Systems: State-of-the-Art and Future Research Directions

Integrated transmission-and-distribution (T&D) modeling is a new and developing method for simulating power systems. Interest in integrated T&D modeling is driven by the changes taking place in power systems worldwide that are resulting in more decentralized power systems with increasingly high levels of distributed energy resources. Additionally, the increasing role of the hitherto passive energy consumer in the management and operation of power systems requires more capable and detailed integrated T&D modeling to understand the interactions between T&D systems. Although integrated T&D modeling has not yet found widespread commercial application, its potential for changing the decades-old power system modeling approaches has led to several research efforts in the last few years that tried to (i) develop algorithms and software for steady-state and dynamic modeling of power systems and (ii) demonstrate the advantages of this modeling approach compared with traditional, separated T&D system modeling. In this paper, we provide a review of integrated T&D modeling research efforts and the methods employed for steady-state and dynamic modeling of power systems. We also discuss our current research in integrated T&D modeling and the potential directions for future research. This paper should be useful for power systems researchers and industry members because it will provide them with a critical summary of current research efforts and the potential topics where research efforts are needed to further advance and demonstrate the utility of integrated T&D modeling.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Decomposition of Feynman integrals by multivariate intersection numbers

We present a detailed description of the recent idea for a direct decomposition of Feynman integrals onto a basis of master integrals by projections, as well as a direct derivation of the differential equations satisfied by the master integrals, employing multivariate intersection numbers. We discuss a recursive algorithm for the computation of multivariate intersection numbers, and provide three different approaches for a direct decomposition of Feynman integrals, which we dub the straight decomposition, the bottom-up decomposition, and the top-down decomposition. These algorithms exploit the unitarity structure of Feynman integrals by computing intersection numbers supported on cuts, in various orders, thus showing the synthesis of the intersection-theory concepts with unitarity-based methods and integrand decomposition. We perform explicit computations to exemplify all of these approaches applied to Feynman integrals, paving a way towards potential applications to generic multi-loop integrals.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Amplitudes, supersymmetric black hole scattering at $\mathcal{O}\left({G}^5\right)$, and loop integration

We compute the potential-graviton contribution to the scattering amplitude, the radial action, and the scattering angle of two extremal black holes in $\mathcal{N}$ = 8 supergravity at the fifth post-Minkowskian order and to next-to-leading order in a large mass expansion (first self-force order). Properties of classical unitarity cuts allow us to focus on the integration-by-parts reduction of planar integrals, while nonplanar integrals at this order are obtained from the planar ones by straightforward manipulations. We present the solution to the differential equations for all master integrals necessary to evaluate the classical scattering amplitudes of massive scalar particles at this order in all gravitational theories, in particular in $\mathcal{N}$ = 8 supergravity, and in general relativity. Despite the appearance of higher-weight generalized polylogarithms and elliptic functions in the solution to the differential equation for master integrals, the final supergravity answer is remarkably simple and contains only (harmonic) polylogarithmic functions up to weight 2. The systematic analysis of elliptic integrals discussed here, as well as the particular organization of boundary integrals in $\mathcal{N}$ = 8 observables are independent of supersymmetry and may have wider applications, including to aspects of collider physics.

Black Holes↗
DOE OSTI2021

Scaled boundary cubature scheme for numerical integration over planar regions with affine and curved boundaries

This paper introduces the scaled boundary cubature (SBC) scheme for accurate and efficient integration of functions over polygons and two-dimensional regions bounded by parametric curves. Over two-dimensional domains, the SBC method reduces integration over a region bounded by m curves to integration over m regions (referred to as curved triangular regions), where each region is bounded by two line segments and a curve. With proper (counterclockwise) orientation of the boundary curves, the scheme is applicable to convex and nonconvex domains. Additionally, for star-convex domains, a tensor-product cubature rule with positive weights and integration points in the interior of the domain is obtained. If the integrand is homogeneous, we show that this new method reduces to the homogeneous numerical integration scheme; however, the SBC scheme is more versatile since it is equally applicable to both homogeneous and non-homogeneous functions. This paper also introduces several methods for smoothing integrands with point singularities and near-singularities. When these methods are used, highly efficient integration of weakly singular functions is realized. Finally, the SBC method is applied to a number of benchmark problems, which reveal its broad applicability and superior performance (in terms of time to generate a rule and accuracy per cubature point) when compared to existing methods for integration.

42 ENGINEERING↗