Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Arbitrary functions”

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 775 records · Page 43

Electronic density response of warm dense matter

Matter at extreme temperatures and pressures—commonly known as warm dense matter (WDM)—is ubiquitous throughout our Universe and occurs in astrophysical objects such as giant planet interiors and brown dwarfs. Moreover, WDM is very important for technological applications such as inertial confinement fusion and is realized in the laboratory using different techniques. A particularly important property for the understanding of WDM is given by its electronic density response to an external perturbation. Such response properties are probed in x-ray Thomson scattering (XRTS) experiments and are central for the theoretical description of WDM. In this work, we give an overview of a number of recent developments in this field. To this end, we summarize the relevant theoretical background, covering the regime of linear response theory and nonlinear effects, the fully dynamic response and its static, time-independent limit, and the connection between density response properties and imaginary-time correlation functions (ITCF). In addition, we introduce the most important numerical simulation techniques, including path-integral Monte Carlo simulations and different thermal density functional theory (DFT) approaches. From a practical perspective, we present a variety of simulation results for different density response properties, covering the archetypal model of the uniform electron gas and realistic WDM systems such as hydrogen. Moreover, we show how the concept of ITCFs can be used to infer the temperature from XRTS measurements of arbitrary complex systems without the need for any models or approximations. Finally, we outline a strategy for future developments based on the close interplay between simulations and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Chern-Simons invariants from ensemble averages

We discuss ensemble averages of two-dimensional conformal field theories associated with an arbitrary indefinite lattice with integral quadratic form Q. We provide evidence that the holographic dual after the ensemble average is the three-dimensional Abelian Chern-Simons theory with kinetic term determined by Q. The resulting partition function can be written as a modular form, expressed as a sum over the partition functions of Chern-Simons theories on lens spaces. For odd lattices, the dual bulk theory is a spin Chern-Simons theory, and we identify several novel phenomena in this case. We also discuss the holographic duality prior to averaging in terms of Maxwell-Chern-Simons theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Verified Optimizer for Quantum Circuits

We present VOQC, the first verified optimizer for quantum circuits, written using the Coq proof assistant. Quantum circuits are expressed as programs in a simple, low-level language called SQIR, a small quantum intermediate representation, which is deeply embedded in Coq. Optimizations and other transformations are expressed as Coq functions, which are proved correct with respect to a semantics of SQIR programs. SQIR programs denote complex-valued matrices, as is standard in quantum computation, but we treat matrices symbolically to reason about programs that use an arbitrary number of quantum bits. SQIR’s careful design and our provided automation make it possible to write and verify a broad range of optimizations in VOQC, including full-circuit transformations from cutting-edge optimizers.

97 MATHEMATICS AND COMPUTING↗

Alfven waves in spiral interplanetary field

This paper presents a theoretical study of the Alfven waves in the spiral interplanetary magnetic field. The Alfven waves under consideration are arbitrary large-amplitude nonmonochromatic microscale waves of any polarization. They superpose on a mesoscale background flow of thermally anisotropic plasma. When the WKB approximation is used, an analytical solution for the amplitude vectors is obtained as a function of the background flow properties: density, velocity, Alfven speed, thermal anisotropy, and the spiral angle. The relative intensity of fluctuations compared with the magnitude of the background field has its maximum in the region near 1 AU. Thus outside of this region the solar wind is less turbulent. Owing to attenuation of microscale Alfven waves, fluctuation energy is converted into the kinetic energy of the solar wind.

Whang, Y. C.↗

Numerical solution of the Navier-Stokes equations for arbitrary two-dimensional multi-element airfoils

Abstracts are presented on a method of numerical solution of the Navier-Stokes equation for the flow about arbitrary airfoils, using a numerically generated curvilinear coordinate system having a coordinate line coincident with the body contour. Results of continuing research are reported and include: application of the Navier-Stokes solution in the vorticity-stream function formulation to a number of single airfoils at Reynolds numbers up to 2000; programming of the Navier-Stokes solution for multiple airfoils in the primitive variable formulation; testing of the potential flow solution of multiple bodies; and development of a generalized coordinate system program.

Thompson, J. F.↗

A theoretical study of heterojunction and graded band gap type solar cells

A computer program was designed for the analysis of variable composition solar cells and applied to several proposed solar cell structures using appropriate semiconductor materials. The program simulates solar cells made of a ternary alloy of two binary semiconductors with an arbitrary composition profile, and an abrupt or Gaussian doping profile of polarity n-on-p or p-on-n with arbitrary doping levels. Once the device structure is specified, the program numerically solves a complete set of differential equations and calculates electrostatic potential, quasi-Fermi levels, carrier concentrations and current densities, total current density and efficiency as functions of terminal voltage and position within the cell. These results are then recorded by computer in tabulated or plotted form for interpretation by the user.

Sutherland, J. E.↗

Wing and leading-edge thrust

Computer program predicts leading-edge thrust for wings of arbitrary planform at supersonic speeds. Methods used in program are based on linearized wing theory. Program first calculates lifting pressures, lift coefficients, drag coefficients, moment coefficients, and lift distributions by using aft-element sensing technique. Next, empirical function calculates set of adjusted pressure-coefficient locations along leading edge. This information is then used to determine limiting value of singularity parameter and value of local leading-edge thrust coefficient.

Carlson, H. W.↗

Formal specification and mechanical verification of SIFT - A fault-tolerant flight control system

The paper describes the methodology being employed to demonstrate rigorously that the SIFT (software-implemented fault-tolerant) computer meets its requirements. The methodology uses a hierarchy of design specifications, expressed in the mathematical domain of multisorted first-order predicate calculus. The most abstract of these, from which almost all details of mechanization have been removed, represents the requirements on the system for reliability and intended functionality. Successive specifications in the hierarchy add design and implementation detail until the PASCAL programs implementing the SIFT executive are reached. A formal proof that a SIFT system in a 'safe' state operates correctly despite the presence of arbitrary faults has been completed all the way from the most abstract specifications to the PASCAL program.

Melliar-Smith, P. M.↗

Real-Time Adaptive Control of Flow-Induced Cavity Tones

An adaptive generalized predictive control (GPC) algorithm was formulated and applied to the cavity flow-tone problem. The algorithm employs gradient descent to update the GPC coefficients at each time step. The adaptive control algorithm demonstrated multiple Rossiter mode suppression at fixed Mach numbers ranging from 0.275 to 0.38. The algorithm was also able t o maintain suppression of multiple cavity tones as the freestream Mach number was varied over a modest range (0.275 to 0.29). Controller performance was evaluated with a measure of output disturbance rejection and an input sensitivity transfer function. The results suggest that disturbances entering the cavity flow are colocated with the control input at the cavity leading edge. In that case, only tonal components of the cavity wall-pressure fluctuations can be suppressed and arbitrary broadband pressure reduction is not possible. In the control-algorithm development, the cavity dynamics are treated as linear and time invariant (LTI) for a fixed Mach number. The experimental results lend support this treatment.

Kegerise, Michael A.↗

Real-Time Feedback Control of Flow-Induced Cavity Tones: Fixed-Gain Control - Part 1

A generalized predictive control (GPC) algorithm was formulated and applied to the cavity flow-tone problem. The control algorithm demonstrated multiple Rossiter-mode suppression at fixed Mach numbers ranging from 0.275 to 0.38. Controller performance was evaluated with a measure of output disturbance rejection and an input sensitivity transfer function. The results suggest that disturbances entering the cavity flow are collocated with the control input at the cavity leading edge. In that case, only tonal components of the cavity wall-pressure fluctuations can be suppressed and arbitrary broadband pressure reduction is not possible with the present sensor/actuator arrangement. In the control-algorithm development, the cavity dynamics were treated as linear and time invariant (LTI) for a fixed Mach number. The experimental results lend support to that treatment.

Kegerise, M. A.↗

Response of a semi-infinite elastic solid to an arbitrary line load along the axis.

The axisymmetric problem of a line load acting along the axis of a semiinfinite elastic solid is solved using Hankel transforms. In this solution the line load is interpreted as a body force loading and by assuming the line load to be of the form of a Dirac delta function the solution of Mindlin's problem of a point load within the interior of the half space is obtained. Solutions of this problem presented in the literature have been obtained using semiinverse techniques whereas the solution given here is obtained in a systematic step-by-step manner.

Agrawal, G. L.↗

Conformal conserved currents in embedding space

We study conformal conserved currents in arbitrary irreducible representations of the Lorentz group using the embedding space formalism. With the help of the operator product expansion, we first show that conservation conditions can be fully investigated by considering only two- and three-point correlation functions. We then find an explicitly conformally-covariant differential operator in embedding space that implements conservation based on the standard position space operator product expansion differential operator ∂μ, although the latter does not uplift to embedding space covariantly. The differential operator in embedding space that imposes conservation is the same differential operator $D_{ijA}$ used in the operator product expansion in embedding space. We provide several examples including conserved currents in irreducible representations that are not symmetric and traceless. With an eye on four-point conformal bootstrap equations for four conserved vector currents $\langle JJJJ\rangle$ and four energy-momentum tensors $\langle TTTT\rangle$, we mostly focus on conservation conditions for $\langle JJ\mathcal{O}\rangle$ and $\langle TT\mathcal{O}\rangle$. Finally, we reproduce and extend the consequences of conformal Ward identities at coincident points by determining three-point coefficients in terms of charges.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Penrose limit of the Weyl double copy

We embed the Penrose limit into the Weyl classical double copy. Thereby, we provide a lift of the double copy properties of plane wave spacetimes into black hole geometries and we open a novel avenue towards taking the classical double copy beyond statements about algebraically special backgrounds. In particular, the Penrose limit, viewed as the leading order Fermi coordinate expansion around a null geodesic, complements approaches leveraging asymptotic flatness such as the asymptotic Weyl double copy. Along the way, we show how our embedding of the Penrose limit within the Weyl double copy naturally fixes the functional ambiguity in the double copy for Petrov type N spacetimes. Here, we also highlight the utility of a spinorial approach to the Penrose limit. In particular, we use this spinorial approach to derive a simple analytical expression for arbitrary Penrose limits of four-dimensional, vacuum type D spacetimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Some remarks on elastic crack-tip stress fields.

It is shown that if the displacement field and stress intensity factor are known as functions of crack length for any symmetrical load system acting on a linear elastic body in plane strain, then the stress intensity factor for any other symmetrical load system whatsoever on the same body may be directly determined. The result is closely related to Bueckner's (1970) weight function, through which the stress intensity factor is expressed as a sum of work-like products between applied forces and values of the weight function at their points of application. An example of the method is given wherein the solution for a crack in a remotely uniform stress field is used to generate the expression for the stress intensity factor due to an arbitrary traction distribution on the faces of a crack. A corresponding theory is developed in an appendix for three-dimensional crack problems, although this appears to be directly useful chiefly for problems in which there is axial symmetry.

Rice, J. R.↗

Radiation from a localized Langmuir oscillation in a uniformly magnetized plasma

The electromagnetic radiation at the first and second harmonics of the electron plasma frequency from a localized Langmuir perturbation is computed for the case of a uniformly magnetized plasma. It is assumed that the localized perturbations in both the electrostatic field and plasma density have cylindrical symmetry about the direction of the ambient magnetic field. The analysis is initially performed for such an arbitrary localized perturbation, and then applied to treat the case of a quasi-planar Langmuir soliton propagating in the direction of the magnetic field. An extensive numerical study of the angular dependence of the radiation spectrum as a function of the ratio of the electron plasma and cyclotron frequencies is described.

Freund, H. P.↗

Working with Bézier Curves as bases for Functional Expansion Tallies

Functional expansion tallies (FETs) are powerful tools for getting more information per history from Monte Carlo simulations, but in the past they have been constrained to orthogonal bases. Bézier curves are used widely in computer aided design (CAD) geometry kernels and could be well-suited for FETs due to their ability to assume many arbitrary shapes, but they use nonorthogonal bases. Recent developments in 2021 have made nonorthogonal FETs possible. The convergence of Bézier curve FETs in both polynomial order and number of samples is explored in this work. It is shown that these bases are well-suited for representing normal distributions and this opens the door to the possibility of other CAD-derived FET bases.

97 - MATHEMATICS AND COMPUTING↗

Regression Based Approach for Robust Finite Element Analysis on Arbitrary Grids. LDRD Final Report

This report summarizes the work performed under a one-year LDRD project aiming to enable accurate and robust numerical simulation of partial differential equations for meshes that are of poor quality. Traditional finite element methods use the mesh to both discretize the geometric domain and to define the finite element shape functions. The latter creates a dependence between the quality of the mesh and the properties of the finite element basis that may adversely affect the accuracy of the discretized problem. In this project, we propose a new approach for defining finite element shape functions that breaks this dependence and separates mesh quality from the discretization quality. At the core of the approach is a meshless definition of the shape functions, which limits the purpose of the mesh to representing the geometric domain and integrating the basis functions without having any role in their approximation quality. The resulting non-conforming space can be utilized within a standard discontinuous Galerkin framework providing a rigorous foundation for solving partial differential equations on low-quality meshes. We present a collection of numerical experiments demonstrating our approach in a wide range of settings: strongly coercive elliptic problems, linear elasticity in the compressible regime, and the stationary Stokes problem. We demonstrate convergence for all problems and stability for element pairs for problems which usually require inf-sup compatibility for conforming methods, also referring to a minor modification possible through the symmetric interior penalty Galerkin framework for stabilizing element pairs that would otherwise be traditionally unstable. Mesh robustness is particularly critical for elasticity, and we provide an example that our approach provides a greater than 5x improvement in accuracy and allows for taking an 8x larger stable timestep for a highly deformed mesh, compared to the continuous Galerkin finite element method. The report concludes with a brief summary of ongoing projects and collaborations that utilize or extend the products of this work.

97 MATHEMATICS AND COMPUTING↗

Elliptic grid generation with orthogonality and spacing control on an arbitrary number of boundaries

A procedure for the generation of two and quasi-three-dimensional grids with control of orthogonality and spacing with respect to any and/or all boundaries of the domain is described. The elliptic grid generation equations of Thompson are solved implicitly. Control of the grid behavior is achieved through the introduction of forcing functions terms in the manner of Steger and Sorenson or in a modification of the method of Hilgenstock. The forcing function terms are constructed on the boundaries and propagated into the domain using transfinite Lagrangian bivariate interpolation. An anisotropic transfinite stencil is introduced and is shown to produce excellent grid behavior particularly in the vicinity of corner singularities. Emphasis is placed on the generation of viscous grids and the method is shown to be suited for use in the generation of grids for internal as well as external flow geometries. A FORTRAN program named PISCES has been written to implement the algorithm. Examples of grids for internal and external flows are given that highlight the characteristics and behavior of the algorithm.

White, J. A.↗