Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quadrature”

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

Phase Demodulation by Frequency Chirping in Coherence Microwave Photonic Interferometry

This paper presents a signal processing method to demodulate the optical interference phase of cascaded individual optical fiber intrinsic Fabry-Perot interferometric (IFPI) sensors in a coherent microwave-photonic interferometry (CMPI) distributed sensing system. The new method utilizes the chirp effect of electro-optic modulator (EOM) to create a quasi-quadrature optical interference phase shift between two adjacent pulses which correspond to two adjacent reflection points in the time domain. The phase shift can be controlled by adjusting the bias voltage that is applied to the EOM. The interference phase is calculated by elliptically fitting the phase shift. The interference phase change is proportional to the optical path difference (OPD) change of the interferometer, and the sign can be used to differentiate the increase or decrease of the OPD. The method is demonstrated for distributed strain sensing, showing good linearity, high resolution and large dynamic range.

42 ENGINEERING↗
DOE OSTI2020

Highly Reliable Multi-Port Smart Inverter Modules for PV-Based Energy Systems

Generally, large electrolytic capacitors (E-caps) are required to decouple the double-line frequency power fluctuation inherent in all single-phase dc-ac or ac-dc converters. These E-caps pose the weak link in the photovoltaic (PV) inverter because they have poor reliability as a result of short life expectancy and higher equivalent series resistance (ESR) which leads to lower system efficiency. In this paper, multiport smart dual-inverter modules are proposed for residential PV inverter systems with balanced outputs to eliminate the requirement of large decoupling capacitors, thus leading to more reliable and highly efficient inverter modules. The balanced dual outputs are controlled in phase quadrature such that the double frequency current is completely eliminated. One phase supports the standalone residential loads, and the other phase interfaces with the grid, which is synchronized to the utility grid. Also, the proposed modular configuration is scalable to accommodate higher power rating by stacking and paralleling multiple modules. Simulation results for the dual-inverter modules rated for 1 kW each are presented to validate the proposed concept. Experimental results from a controller hardware-in-the-loop (CHIL) platform also verify its feasibility.

14 SOLAR ENERGY↗
DOE OSTI2021

Detecting Atmospheric Molecules of Nontransiting Temperate Terrestrial Exoplanets Using High-resolution Spectroscopy in the Mid-infrared Domain

Motivated by the development of high-dispersion spectrographs in the mid-infrared (MIR) range, we study their application to the atmospheric characterization of nearby nontransiting temperate terrestrial planets around M-type stars. We examine the detectability of CO{sub 2}, H{sub 2}O, N{sub 2}O, and O{sub 3 }features in high-resolution planetary thermal emission spectra at 12–18 μm assuming an Earth-like profile and a simplified thermal structure. The molecular line width of such planets can be comparable to or broader than the Doppler shift due to the planetary orbital motion. Given the likely difficulty in knowing the high-resolution MIR spectrum of the host star with sufficient accuracy, we propose observing the target system at two quadrature phases and extracting the differential spectra as the planetary signal. In this case, the signals can be substantially suppressed compared with the case where the host star spectrum is perfectly known, as some parts of the spectral features do not remain in the differential spectra. Despite this self-subtraction, the CO{sub 2} and H{sub 2}O features of nearby (≲5 pc) systems with mid-/late-M host stars would be feasible with a 6.5 m class cryogenic space telescope, and orbital inclination could also be constrained for some of them. For CO{sub 2} and N{sub 2}O in a 1 bar Earth-like atmosphere, this method would be sensitive when the mixing ratio is 1–10{sup 3} ppm. The detectability of molecules except O{sub 3} is not significantly improved when the spectral resolution is higher than R≳10,000, although the constraint on the orbital inclination is improved. This study provides some benchmark cases useful for assessing the value of MIR high-resolution spectroscopy in terms of characterization of potentially habitable planets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗
DOE OSTI2025

Dual-Phase Malicious User Detection Scheme for IM-OFDMA Systems Using IQ Imbalance

Physical-layer security techniques have contributed to the achievement of various security objectives in an efficient and lightweight manner. Thus, these techniques have been widely considered for limited-resource networks such as Internet of Things networks. Among the different security objectives, malicious user detection by exploiting physical-layer parameters has demonstrated efficient performance. In this work, malicious user detection in the recently proposed index modulation-based orthogonal frequency division multiple access (IM-OFDMA) is addressed. The proposed malicious user detection scheme exploits the hardware impairments, especially the in-phase and quadrature imbalance parameters, for both legitimate and malicious users to design a dual-phase efficient detection scheme. The proposed scheme accounts for the special characteristics of IM-OFDMA transmission that are different from other multiple-access techniques. The performance of the proposed scheme was evaluated considering detection probability and false alarm probability performance metrics. Moreover, closed-form expressions of these metrics were derived for both phases and were validated by Monte Carlo simulation results under different configurations of IM-OFDMA systems.

Alaca, Ozgur [ORNL] (ORCID:0000000153713758)↗
DOE OSTI2021

Extending a 10‐node composite tetrahedral finite element for solid mechanics

Abstract We propose to extend the composite tetrahedral finite element first introduced by Thoutireddy et al. (2002) and recently reformulated by Ostien et al. (2016). We generalize the gradient operator and mass matrix to curved domains through analytical expressions weighted by subtetrahedra Jacobians. Optimal integration weights are constructed to increase the accuracy of Gaussian quadrature. We preserve the variational structure of the formulation through a new five‐field functional with additional, independent fields for the Jacobian and the pressure. This approach not only obviates volumetric locking but also yields symmetry. A deleterious soft mode, common to both the quadratic and composite tetrahedral element with constant pressure formulations is effectively stabilized through a novel convex energy penalty function. Numerous numerical examples spanning a patch test to the impact of a Taylor bar demonstrate the accuracy, robustness, and convergence of the extended composite tetrahedral element for application to structural metals.

Foulk III, James W.↗
DOE OSTI2021

Approximation of periodic Green's operator in real space using numerical integration and its use in fast Fourier transform-based micromechanical models

In this work, we propose an expression for the periodic first derivative of Green's function in real space. The proposed expression allows an alternative way of computing the periodic Green's operator based on periodically summing the free-space Green's operator in terms of an appropriate quadrature rule. We provide computational examples, which show the accuracy of the proposed approach, together with reduced spurious oscillations in the solution fields.

42 ENGINEERING↗
DOE OSTI2022

The strip method for shape derivatives

Abstract A major challenge in shape optimization is the coupling of finite element method (FEM) codes in a way that facilitates efficient computation of shape derivatives. This is particularly difficult with multiphysics problems involving legacy codes, where the costs of implementing and maintaining shape derivative capabilities are prohibitive. The volume and boundary methods are two approaches to computing shape derivatives. Each has a major drawback: the boundary method is less accurate, while the volume method is more invasive to the FEM code. We introduce the strip method , which computes shape derivatives on a strip adjacent to the boundary. The strip method makes code coupling simple. Like the boundary method, it queries the state and adjoint solutions at quadrature nodes, but requires no knowledge of the FEM code implementations. At the same time, it exhibits the higher accuracy of the volume method. As an added benefit, its computational complexity is comparable to that of the boundary method, that is, it is faster than the volume method. We illustrate the benefits of the strip method with numerical examples.

Hardesty, Sean↗
DOE OSTI2024

FLOWERS AEP: An Analytical Model for Wind Farm Layout Optimization

Annual energy production (AEP) is commonly used in objective functions for wind farm layout optimization. AEP is proportional to wind farm power production integrated over an annual distribution of free-stream wind conditions. Physics-based estimates of wind farm power production typically rely on low-fidelity engineering wake models that approximate the steady-state wind farm flow field. AEP estimates are then obtained by performing independent simulations for discrete wind conditions and using rectangular quadrature to account for each condition's expected frequency of occurrence. Depending on the number of simulated discrete wind conditions, this numerical integral could be hampered by poor accuracy or high computational costs. The FLOWERS AEP model instead poses an analytical integral of the engineering wake model over the variable wind conditions, yielding a closed-form, analytical function for wind farm AEP. This paper derives the analytical functions for FLOWERS AEP and its derivatives with respect to turbine position, which are useful for gradient-based wind farm layout optimization, in nondimensional form. We then analyze the benefits of the FLOWERS AEP model over conventional reference models, focusing on its low cost, adequate wake loss predictions, and smooth design space. Although the FLOWERS approach is found to predict the exact value of AEP with some error relative to the reference model (within 14% on average), it dramatically reduces computation time by an order of magnitude, produces a qualitatively similar design space at relatively low resolution, and yields comparable optimal layouts. This significant speed improvement is critical in layout optimization applications, where determining an optimal layout in an efficient manner is more important than precise AEP prediction.

17 WIND ENERGY↗
DOE OSTI2024

Tensor network representation of non-abelian gauge theory coupled to reduced staggered fermions

We show how to construct a tensor network representation of the path integral for reduced staggered fermions coupled to a non-abelian gauge field in two dimensions. The resulting formulation is both memory and computation efficient because reduced staggered fermions can be represented in terms of a minimal number of tensor indices while the gauge sector can be approximated using Gaussian quadrature with a truncation. Numerical results obtained using the Grassmann TRG algorithm are shown for the case of SU(2) lattice gauge theory and compared to Monte Carlo results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Fuzzy spheres in stringy matrix models: quantifying chaos in a mixed phase space

We consider a truncation of the BMN matrix model to a configuration of two fuzzy spheres, described by two coupled non-linear oscillators dependent on the mass parameter μ. The classical phase diagram of the system generically (μ ≠ 0) contains three equilibrium points: two centers and a center-saddle; as μ → 0 the system exhibits a pitchfork bifurcation. We demonstrate that the system is exactly integrable in quadratures for μ = 0, while for very large values of μ, it approaches another integrable point characterized by two harmonic oscillators. The classical phase space is mixed, containing both integrable islands and chaotic regions, as evidenced by the classical Lyapunov spectrum. At the quantum level, we explore indicators of early and late time chaos. The eigenvalue spacing is best described by a Brody distribution, which interpolates between Poisson and Wigner distributions; it dovetails, at the quantum level, the classical results and reemphasizes the notion that the quantum system is mixed. We also study the spectral form factor and the quantum Lyapunov exponent, as defined by out-of-time-ordered correlators. These two indicators of quantum chaos exhibit weak correlations with the Brody distribution. We speculate that the behavior of the system as μ → 0 dominates the spectral form factor and the quantum Lyapunov exponent, making these indicators of quantum chaos less effective in the context of a mixed phase space.

AdS-CFT correspondence↗
DOE OSTI2022

Study of ${\mathrm{B}}_{\mathrm{c}}^{+}\to \mathrm{J}/\psi {\mathrm{D}}_{\mathrm{s}}^{+}$ and ${\mathrm{B}}_{\mathrm{c}}^{+}\to \mathrm{J}/\psi {\mathrm{D}}_{\mathrm{s}}^{\ast +}$ decays in pp collisions at $\sqrt{\mathrm{s}}$ = 13 TeV with the ATLAS detector

A study of ${B}_c^{+}\to J/\psi {D}_s^{+}$ and ${B}_c^{+}\to J/\psi {D}_s^{\ast +}$ decays using 139 fb –1 of integrated luminosity collected with the ATLAS detector from $\sqrt{s}$ = 13 TeV pp collisions at the LHC is presented. The ratios of the branching fractions of the two decays to the branching fraction of the ${B}_c^{+}$ → J/ψπ + decay are measured: $\mathcal{B}\left({B}_c^{+}\to J/\psi {D}_s^{+}\right)/\mathcal{B}\left({B}_c^{+}\to J/{\psi \pi}^{+}\right)$ = 2.76 ± 0.47 and $\mathcal{B}\left({B}_c^{+}\to J/\psi {D}_s^{\ast +}\right)/\mathcal{B}\left({B}_c^{+}\to J/{\psi \pi}^{+}\right)$ = 5.33 ± 0.96. The ratio of the branching fractions of the two decays is found to be $\mathcal{B}\left({B}_c^{+}\to J/\psi {D}_s^{\ast +}\right)/\mathcal{B}\left({B}_c^{+}\to J/\psi {D}_s^{\ast +}\right)$ = 1.93 ± 0.26. For the ${B}_c^{+}\to J/\psi {D}_s^{\ast +}$ decay, the transverse polarization fraction, Γ ±± /Γ, is measured to be 0.70 ± 0.11. The reported uncertainties include both the statistical and systematic components added in quadrature. The precision of the measurements exceeds that in all previous studies of these decays. These results supersede those obtained in the earlier ATLAS study of the same decays with $\sqrt{s}$ = 7 and 8 TeV pp collision data. A comparison with available theoretical predictions for the measured quantities is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A variational multiscale immersed meshfree method for heterogeneous materials

Abstract We introduce an immersed meshfree formulation for modeling heterogeneous materials with flexible non-body-fitted discretizations, approximations, and quadrature rules. The interfacial compatibility condition is imposed by a volumetric constraint, which avoids a tedious contour integral for complex material geometry. The proposed immersed approach is formulated under a variational multiscale based formulation, termed the variational multiscale immersed method (VMIM). Under this framework, the solution approximation on either the foreground or the background can be decoupled into coarse-scale and fine-scale in the variational equations, where the fine-scale approximation represents a correction to the residual of the coarse-scale equations. The resulting fine-scale solution leads to a residual-based stabilization in the VMIM discrete equations. The employment of reproducing kernel (RK) approximation for the coarse- and fine-scale variables allows arbitrary order of continuity in the approximation, which is particularly advantageous for modeling heterogeneous materials. The effectiveness of VMIM is demonstrated with several numerical examples, showing accuracy, stability, and discretization efficiency of the proposed method.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An efficient method to integrate polynomials over polytopes and curved solids

Here in this paper, we present an efficient approach to compute the integral of monomials and polynomials over polyhedra and regions defined by parametric curved boundary surfaces. We use Euler's theorem for homogeneous functions in combination with Stokes's theorem to reduce the integration of a monomial over a three-dimensional solid to its boundary. If the solid is a polytope, through a recursive application of these theorems, the integral is further reduced to just the evaluation of the monomial and its derivatives at the vertices of the polytope. The present approach is simpler than existing techniques that rely on repeated use of the divergence theorem, which require the antiderivative of the monomials and the projection of these functions onto hyperplanes. For convex and nonconvex polytopes, our approach does not introduce any approximation for the integration of monomials. For curved solid regions bounded by surfaces that admit a parameterization, the same approach yields simplified formulas to compute the integral of any homogeneous function, including monomials. For surfaces parameterized by polynomial surfaces (such as Bezier surface triangles and B-spline patches), the method yields machine-precision accuracy for the volumetric integration of monomials with an appropriate quadrature rule. Numerical examples over regions bounded by polynomial surfaces and rational surfaces are presented to establish the accuracy and efficiency of the method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗
DOE OSTI2021

Parallel projection—An improved return mapping algorithm for finite element modeling of shape memory alloys

Here, we present a novel finite element analysis of inelastic structures containing Shape Memory Alloys (SMAs). Phenomenological constitutive models for SMAs lead to material nonlinearities, that require substantial computational effort to resolve. Finite element analysis methods, which rely on Gauss quadrature integration schemes, must solve two sets of coupled differential equations: one at the global level and the other at the local, i.e. Gauss point level. In contrast to the conventional return mapping algorithm, which solves these two sets of coupled differential equations separately using a nested Newton procedure, we propose a scheme to solve the local and global differential equations simultaneously. In the process we also derive closed-form expressions used to update the internal/constitutive state variables, and unify the popular closest-point and cutting plane methods with our formulas. Numerical testing indicates that our method allows for larger thermomechanical loading steps and provides increased computational efficiency, over the standard return mapping algorithm.

42 ENGINEERING↗
DOE OSTI2022

Modeling wave propagation in elastic solids via high-order accurate implicit-mesh discontinuous Galerkin methods

Here, a high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.

42 ENGINEERING↗