Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SERIES EXPANSION”

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
NASA NTRS1977

Vibrations of incompressible liquid-plastic tank combinations using a series representation of the liquid

A series representation is presented of the dynamic behavior of incompressible non-viscous liquids contained in partially filled elastic shells of revolution. Using a series in which each term satisfies exactly Laplace's equation, the liquid portion of the problem becomes an integral equation on the wetted wall of the shell. Physical insight required in selecting the terms of the series is such that selection may be readily programed on a computer. Liquid-shell vibration interaction is investigated by modifying an existing shell of revolution vibration program with the series expansion of the liquid pressure loading on the wetted shell wall. Application of the modified program is demonstrated on a hemispherical tank and a realistic scaled model of a launch vehicle tank having a more complex geometry. Comparison of the present analysis with a NASTRAN finite element analysis for a hemispherical tank, and comparison with test results for a launch vehicle tank shows excellent agreement and rapid convergence of the series representation of the liquid to the true solution.

Housner, J. M.↗
NASA NTRS1990

On the far-field stream function condition for two-dimensional incompressible flows

The present demonstration of the usefulness of the integral series expansion of the stream function as a far-field computational boundary condition shows the method to require only a 10-percent/time-step increase in computational effort over alternative boundary conditions, in the case of implementation of unsteady problems using a direct elliptic solver. So long as the vorticity was encompassed within the computational domain, the method proved sufficiently accurate to yield virtually identical results for two widely different domains. While the integral-series condition yielded the best results for periodic flow, the Neumann condition gave comparable accuracy with less computation time for the steady-flow case despite its inability to treat periodic flow with vortex shedding.

Sa, Jong-Youb↗
DOE OSTI2014

Extended space expectation values in quantum dynamical system evolutions

The time variant power series expansion for the expectation value of a given quantum dynamical operator is well-known and well-investigated issue in quantum dynamics. However, depending on the operator and Hamiltonian singularities this expansion either may not exist or may not converge for all time instances except the beginning of the evolution. This work focuses on this issue and seeks certain cures for the negativities. We work in the extended space obtained by adding all images of the initial wave function under the system Hamiltonian’s positive integer powers. This requires the introduction of certain appropriately defined weight operators. The resulting better convergence in the temporal power series urges us to call the new defined entities “extended space expectation values” even though they are constructed over certain weight operators and are somehow pseudo expectation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

On the convergence of Neumann series for electrostatic fracture response

Here, the feasibility of Neumann series expansion of Maxwell’s equations in the electrostatic limit is investigated for potentially rapid and approximate subsurface imaging of geologic features proximal to metallic infrastructure in an oilfield environment. While generally useful for efficient modeling of mild conductivity perturbations in uncluttered settings, we raise the question of its suitability for situations, such as oilfield, where metallic artifacts are pervasive, and in some cases, in direct electrical contact with the conductivity perturbation on which the Neumann series is computed. Convergence of the Neumann series and its residual error are computed using the hierarchical finite element framework for a canonical oilfield model consisting of an “L” shaped, steel-cased well, energized by a steady state electrode, and penetrating a small set of mildly conducting fractures near the heel of the well. For a given node spacing h in the finite element mesh, we find that the Neumann series is ultimately convergent if the conductivity is small enough - a result consistent with previous presumptions on the necessity of small conductivity perturbations. However, we also demonstrate that the spectral radius of the Neumann series operator grows as ~ 1/ h , thus suggesting that in the limit of the continuous problem h → 0, the Neumann series is intrinsically divergent for all conductivity perturbation, regardless of their smallness. The hierarchical finite element methodology itself is critically analyzed and shown to possess the h 2 error convergence of traditional linear finite elements, thereby supporting the conclusion of an inescapably divergent Neumann series for this benchmark example. Application of the Neumann series to oilfield problems with metallic clutter should therefore be done with careful consideration to the coupling between infrastructure and geology. Here, the methods used here are demonstrably useful in such circumstances.

58 GEOSCIENCES↗
DOE OSTI2021

An expansion of well tempered gravity

When faced with two nigh intractable problems in cosmology — how to remove the original cosmological constant problem and how to parametrize modified gravity to explain current cosmic acceleration — we can make progress by counterposing them. The well tempered solution to the cosmological constant through degenerate scalar field dynamics also relates disparate Horndeski gravity terms, making them contrapuntal. We derive the connection between the kinetic term K and braiding term G 3 for shift symmetric theories (including the running Planck mass G 4 ), extending previous work on monomial or binomial dependence to polynomials of arbitrary finite degree. We also exhibit an example for an infinite series expansion. Finally, this contrapuntal condition greatly reduces the number of parameters needed to test modified gravity against cosmological observations, for these "golden" theories of gravity.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

A hybrid neurocomputing/numerical strategy for nonlinear structural analysis

A hybrid neurocomputing/numerical strategy is presented for geometrically nonlinear analysis of structures. The strategy combines model-free data processing capabilities of computational neural networks with a Pade approximants-based perturbation technique to predict partial information about the nonlinear response of structures. In the hybrid strategy, multilayer feedforward neural networks are used to extend the validity of solutions by using training samples produced by Pade approximations to the Taylor series expansion of the response function. The range of validity of the training samples is taken to be the radius of convergence of Pade approximants and is estimated by setting a tolerance on the diverging approximants. The norm of residual vector of unbalanced forces in a given element is used as a measure to assess the quality of network predictions. To further increase the accuracy and the range of network predictions, additional training data are generated by either applying linear regression to weight matrices or expanding the training data by using predicted coefficients in a Taylor series. The effectiveness of the hybrid strategy is assessed by performing large-deflection analysis of a doubly-curved composite panel with a circular cutout, and postbuckling analyses of stiffened composite panels subjected to an in-plane edge shear load. In all the problems considered, the hybrid strategy is used to predict selective information about the structural response, namely the total strain energy and the maximum displacement components only.

Szewczyk, Z. Peter↗
DOE OSTI2018

Incorporating Brueckner-Hartree-Fock correlations in energy density functionals

Recently, a microscopically motivated nuclear energy density functional was derived by applying the density-matrix expansion to the Hartree-Fock (HF) energy obtained from long-range chiral effective-field theory two- and three-nucleon interactions. However, the HF approach cannot account for all many-body correlations. One class of correlations is included by Brueckner-Hartree-Fock (BHF) theory, which gives an improved definition of the one-body HF potential by replacing the interaction by a reaction matrix G. In this paper, we find that the difference between the G matrix and the similarity renormalization group evolved nucleon-nucleon potential V SRG can be well accounted for by a truncated series of contact terms. This is consistent with renormalization-group decoupling generating a series of counterterms as short-distance physics is integrated out. Here, the coefficients C n of the power series expansion Σ C n q n for the counterterms are examined for two potentials at different renormalization-group resolutions and at a range of densities. The success of this expansion for G–V SRG means we can apply the density-matrix expansion at the HF level with low-momentum interactions and density-dependent zero-range interactions to model BHF correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Low-shear three-dimensional equilibria in a periodic cylinder

We carry out expansions of non-symmetric toroidal ideal magnetohydrodynamic (MHD) equilibria with nested flux surfaces about a periodic cylinder, in which physical quantities are periodic of period $2\unicode[STIX]{x03C0}$ in the cylindrical angle $\unicode[STIX]{x1D703}$ and $z$ . The cross-section of a flux surface at a constant toroidal angle is assumed to be approximately circular, and data are given on the cylindrical flux surface $r=1$ . Furthermore, we assume that the magnetic field lines are closed on the lowest-order flux surface, and the magnetic shear is relatively small. We extend earlier work in a flat torus by Weitzner ( Phys. Plasmas , vol. 23, 2016, 062512) and demonstrate that a power series expansion can be carried out to all orders using magnetic flux as an expansion parameter. The cylindrical metric introduces certain new features to the expansions compared to the flat torus. However, the basic methodology of dealing with resonance singularities remains the same. The results, even though lacking convergence proofs, once again support the possibility of smooth, low-shear non-symmetric toroidal MHD equilibria.

Physics↗
NASA NTRS2011

Validation of Real-Time Data Processing for the Ground and Air-MSPI Systems

JPL is currently developing the multi-angle spectro-polarimetric imager (MSPI), targeted for the Aerosol-Cloud-Ecosystems (ACE) mission, as defined in the National Academies 2007 Decadal Survey. In preparation for the space instrument, the MSPI team has built two incremental camera systems (Ground- and Air-MSPI) to improve understanding of the proposed architecture. Ground-MSPI is a gimballed instrument used primarily for stationary observation and characterization of the imager and optics. The ER-2 based Air-MSPI operates in a step-and-stare mode, providing multi-angle imaging of a static target. This mode-of-operation simulates the observation scenario of the space instrument. Physically, MSPI is a pushbroom camera with a specialized frontend. Before imaging, light entering the camera passes through a pair of photoelastic modulators and a set of pattern polarizers. These optical elements act on the light to make polarimetric extraction computationally feasible. Calculating polarimetric parameters from the imager's data stream requires a real-time least-squares computation that produces coefficients of a truncated time-series expansion of the image. As reported in, the data processing algorithm can operate in real-time on a Xilinx Virtex-5 FPGA. Moving beyond verification with an onboard data source, the algorithm has been validated on a commercial development board interfaced with the ground camera. In addition, the algorithm has been instantiated within the Air-MSPI electronics board's FPGA, and in situ first-light has been achieved.

real-time data processing↗
DOE OSTI2020

Ultrasonic waves from radial mode excitation of a piezoelectric disc on the surface of an elastic solid

In this work, we discuss the ultrasonic wave propagation characteristics in an isotropic elastic solid due to the radial mode excitation of a piezoelectric disc actuator bonded to its surface. Finite element simulations using coupled electromechanical modeling is employed to investigate the wave propagation behavior. We find that the radial mode vibrations on the surface of the elastic solid generates all the three types of ultrasonic waves: longitudinal, shear, and surface waves. The waves in the solid are comprised of a central lobe and multiple side lobes based on the frequency of the radial mode excitation. The central lobe is predominantly composed of longitudinal waves and the side lobes are composed of shear waves. While the longitudinal waves have a strong central lobe and weak side lobes that are not fully developed, shear waves consists of fully developed side lobes. In addition, we observe that the longitudinal waves have fewer side lobes inside the solid compared to the shear waves. A semi-analytical approach is presented to explain the above observations. The interface traction boundary condition between the piezoelectric disc and the elastic solid is approximated as a truncated 'Bessel' excitation over the area of the piezoelectric disc in contact with the solid. A Fourier–Bessel series expansion technique is used to investigate the wave propagation behavior from the above traction boundary condition. The results obtained explain the observations from the finite element simulations with regard to the number and distribution of side lobes pertaining to the longitudinal and shear waves generated from the radial modes of the piezoelectric disc. The proposed semi-analytical approach is general and can be applied to any arbitrary axisymmetric excitations on the surface of an elastic solid.

36 MATERIALS SCIENCE↗
DOE OSTI2011

POLARIZED LINE FORMATION IN MULTI-DIMENSIONAL MEDIA. IV. A FOURIER DECOMPOSITION TECHNIQUE TO FORMULATE THE TRANSFER EQUATION WITH ANGLE-DEPENDENT PARTIAL FREQUENCY REDISTRIBUTION

To explain the linear polarization observed in spatially resolved structures in the solar atmosphere, the solution of polarized radiative transfer (RT) equation in multi-dimensional (multi-D) geometries is essential. For strong resonance lines, partial frequency redistribution (PRD) effects also become important. In a series of papers, we have been investigating the nature of Stokes profiles formed in multi-D media including PRD in line scattering. For numerical simplicity, so far we have restricted our attention to the particular case of PRD functions which are averaged over all the incident and scattered directions. In this paper, we formulate the polarized RT equation in multi-D media that takes into account the Hanle effect with angle-dependent PRD functions. We generalize here to the multi-D case the method for Fourier series expansion of angle-dependent PRD functions originally developed for RT in one-dimensional geometry. We show that the Stokes source vector S = (S{sub I} , S{sub Q} , S{sub U} ){sup T} and the Stokes vector I = (I, Q, U){sup T} can be expanded in terms of infinite sets of components S-tilde{sup (k)}, I-tilde{sup (k)}, respectively, k in [0, +{infinity}). We show that the components S-tilde{sup (k)} become independent of the azimuthal angle ({psi}) of the scattered ray, whereas the components I-tilde{sup (k)} remain dependent on {psi} due to the nature of RT in multi-D geometry. We also establish that S-tilde{sup (k)} and I-tilde{sup (k)} satisfy a simple transfer equation, which can be solved by any iterative method such as an approximate Lambda iteration or a Bi-Conjugate Gradient-type projection method provided we truncate the Fourier series to have a finite number of terms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

POLARIZED LINE TRANSFER WITH F-STATE INTERFERENCE IN A NON-MAGNETIC MEDIUM: PARTIAL FREQUENCY REDISTRIBUTION EFFECTS IN THE COLLISIONLESS REGIME

Quantum interference phenomena manifest themselves in several ways in the polarized solar spectrum formed due to coherent scattering processes. One such effect arises due to interference between the fine structure (J) states giving rise to multiplets. Another effect is that which arises due to interference between the hyperfine structure (F) states. We extend the redistribution matrix derived for the J-state interference to the case of F-state interference. We then incorporate it into the polarized radiative transfer equation and solve it for isothermal constant property slab atmospheres. The relevant transfer equation is solved using a polarized approximate lambda iteration (PALI) technique based on operator perturbation. An alternative method derived from the Neumann series expansion is also proposed and is found to be relatively more efficient than the PALI method. The effects of partial frequency redistribution and the F-state interference on the shapes of the linearly polarized Stokes profiles are discussed. The emergent Stokes profiles are computed for hypothetical line transitions arising due to hyperfine structure splitting of the upper J = 3/2 and lower J = 1/2 levels of a two-level atom model with nuclear spin I{sub s} = 3/2. We confine our attention to the non-magnetic scattering in the collisionless regime.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Sensitivity analysis of a layered piezoelectric system using ZFEM

The complex variable finite element method (ZFEM) is a numerical technique which aims to find the partial derivatives of the independent variables with respect to variation in dependent parameters declared in the physics. This is done by combining the complex Taylor series expansion within the weak formulation of the governing equation in a coupled system of linear equations forming a complex valued block matrix given by the Cauchy–Riemann matrix representation. In this work, two-dimensional linear first-order elements have been implemented in ZFEM to predict the design derivatives of the mechanical displacement field and the voltage potential field for a layered piezoelectric system in a steady-state study with Dirichlet boundary condition applied at the top and bottom edges of the geometry. This approach allows the standard FEM solution to quantify the sensitivity of the mechanical displacement and voltage potential fields with respect to small variations in the material properties through the information obtained from the computation of the derivatives. The domain is formed by a layered body with PZT-4 and PZT-5 stacked together. For result verification, the numerical solution obtained with ZFEM was compared to results from a commercial FEM package and the solution from the imaginary part was compared to the exact solution of a well-known benchmark problem. In conclusion, comparison of the results showed good agreement for both the real and imaginary parts of the solution and the largest sensitivities were found in PZT-5 specifically in C 13 , C 33 , and ε 33 .

42 ENGINEERING↗
DOE OSTI2019

A Framework for Distributed Approximation of Moments with Higher-Order Derivatives through Automatic Differentiation

We present a framework for the distributed approximation of moments, enabling the evaluation of the uncertainty in a dynamical system. The first and second moment, mean, and variance are computed with up to third-order Taylor series expansion. The required derivatives for the expansion are generated automatically by automatic differentiation and propagated through an implicit time stepper. The computational kernels are the accumulation of the derivatives (Jacobian, Hessian, tensor) and the covariance matrix. We apply distributed parallelism to the Hessian or third-order tensor, and the user merely has to provide a function for the differential equation, thus achieving similar ease of use as Monte Carlo-based methods. We demonstrate our approach using with benchmarks on Theta, a KNL-based system at the Argonne Leadership Computing Facility.

Automatic Differentiation↗
DOE OSTI2024

One- and two-dimensional higher-point conformal blocks as free-particle wavefunctions in $$ {\textrm{AdS}}_3^{\otimes m} $$

Abstract We establish that all of the one- and two-dimensional global conformal blocks are, up to some choice of prefactor, free-particle wavefunctions in tensor products of AdS 3 or limits thereof. Our first core observation is that the six-point comb-channel conformal blocks correspond to free-particle wavefunctions on an AdS 3 constructed directly in cross-ratio space. This construction generalizes to blocks for a special class of diagrams, which are determined as free-particle wavefunctions in tensor products of AdS 3 . Conformal blocks for all the remaining topologies are obtained as limits of the free wavefunctions mentioned above. Our results show directly that the integrable models associated with all one- and two-dimensional conformal blocks can be seen as limits of free theory, and manifest a relation between AdS and CFT kinematics that lies outside of the standard AdS/CFT dictionary. We complete the discussion by providing explicit Feynman-like rules that can be used to work out blocks for all topologies, as well as a Mathematica notebook that allows simple computation of Casimir equations and series expansions for blocks, by requiring just an OPE diagram as input.

Physics↗
DOE OSTI2017

An algebraic approach to the analytic bootstrap

We develop an algebraic approach to the analytic bootstrap in CFTs. By acting with the Casimir operator on the crossing equation we map the problem of doing large spin sums to any desired order to the problem of solving a set of recursion relations. We compute corrections to the anomalous dimension of large spin operators due to the exchange of a primary and its descendants in the crossed channel and show that this leads to a Borel-summable expansion. Here, we analyse higher order corrections to the microscopic CFT data in the direct channel and its matching to infinite towers of operators in the crossed channel. We apply this method to the critical O(N ) model. At large N we reproduce the first few terms in the large spin expansion of the known two-loop anomalous dimensions of higher spin currents in the traceless symmetric representation of O(N ) and make further predictions. At small N we present the results for the truncated large spin expansion series of anomalous dimensions of higher spin currents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Two-loop master integrals for leading-color $$ pp\to t\overline{t}H $$ amplitudes with a light-quark loop

Abstract We compute the two-loop master integrals for leading-color QCD scattering amplitudes including a closed light-quark loop in$$ t\overline{t}H $$ t t ¯ H production at hadron colliders. Exploiting numerical evaluations in modular arithmetic, we construct a basis of master integrals satisfying a system of differential equations inϵ-factorized form. We present the analytic form of the differential equations in terms of a minimal set of differential one-forms. We explore properties of the function space of analytic solutions to the differential equations in terms of iterative integrals which can be exploited for studying the analytic form of related scattering amplitudes. Finally, we solve the differential equations using generalized series expansions to numerically evaluate the master integrals in physical phase space. As the first computation of a set of two-loop seven-scale master integrals, our results provide valuable input for analytic studies of scattering amplitudes in processes involving massive particles and a large number of kinematic scales.

Physics↗
DOE OSTI2022

Exploring the strong-coupling region of SU( N ) Seiberg-Witten theory

We consider the Seiberg-Witten solution of pure N = 2 gauge theory in four dimensions, with gauge group SU(N). A simple exact series expansion for the dependence of the 2(N – 1) Seiberg-Witten periods a I (u), a DI (u) on the N – 1 Coulomb-branch moduli un is obtained around the Z 2N -symmetric point of the Coulomb branch, where all u n vanish. This generalizes earlier results for N = 2 in terms of hypergeometric functions, and for N = 3 in terms of Appell functions. Using these and other analytical results, combined with numerical computations, we explore the global structure of the Kähler potential K = 1/2Σ I Im(a¯ I a DI ), which is single valued on the Coulomb branch. Evidence is presented that K is a convex function, with a unique minimum at the Z 2N -symmetric point. Finally, we explore candidate walls of marginal stability in the vicinity of this point, and their relation to the surface of vanishing Kähler potential.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗