On systems of orthogonal functions
Hilberts method of solving linear integral equations using system of orthogonal functions
LINEAR EQUATION↗
SEARCH · Search NASA
Results for “Functional equations”
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.
A vortex wake capturing method for potential flow calculations
A method is presented for modifying finite difference solutions of the potential equation to include the calculation of non-planar vortex wake features. The approach is an adaptation of Baker's 'cloud in cell' algorithm developed for the stream function-vorticity equations. The vortex wake is tracked in a Lagrangian frame of reference as a group of discrete vortex filaments. These are distributed to the Eulerian mesh system on which the velocity is calculated by a finite difference solution of the potential equation. An artificial viscosity introduced by the finite difference equations removes the singular nature of the vortex filaments. Computed examples are given for the two-dimensional time dependent roll-up of vortex wakes generated by wings with different spanwise loading distributions.
Murman, E. M.↗
Method of controlling a resin curing process
An analytical technique for controlling the curing process of fiber-reinforced composite materials that are formed using thermosetting resins. The technique is the percent gel method and involves development of a time-to-gel equation as s function of temperature. From this equation a rate-of-gel equation is then determined, and a percent gel is calculated which is the product of the rate-of-gel times time. Percent gel accounting is used to control the proper pressure application point in an autoclave cure process to achieve desired properties in a production composite part.
Charles Neal Webster↗
A Probabilistic Framework for the Validation and Certification of Computer Simulations
The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.
Ghanem, Roger↗
Spherical balloon response to three-dimensional time-dependent flows
The concept of the Lagrangian displacement of a balloon is introduced. It is shown that the general balloon response problem is extremely complicated because the wind-forcing functions in the balloon equations of motion are functions of the wind velocity vector and its Eulerian first derivatives evaluated at the location of the balloon. The linear perturbation equations for a spherical balloon are derived by perturbing the components of velocity of the balloon about a terminal velocity state which is in equilibrium with a space-time invariant mean horizontal flow. The atmospheric flow is also perturbed such that the resulting equations can be used to analyze the responses of spherical balloons to three-dimensional time-dependent flows. The wind field is represented in terms of a four-fold Fourier integral that involves three orthogonal wave numbers and a frequency, while the balloon components of velocity are represented as Fourier integrals involving a frequency which, in turn, is a function of the wind field wave numbers and frequency and the unperturbed flow components of velocity.
Fichtl, G. H.↗
Electromagnetic scattering by coated convex surfaces and wedges simulated by approximate boundary conditions
Asymptotic/high-frequency solutions are developed for analyzing the non-specular scattering mechanisms associated with coated convex surfaces and edges simulated by approximate boundary conditions. In particular, the standard impedance boundary conditions (SIBC's) and the second order generalized impedance boundary conditions (GIBC's) are employed for a characterization of the edge diffraction, creeping wave, and surface diffracted wave contributions. To study the creeping wave and surface diffracted wave mechanisms, rigorous UTD (uniform geometrical theory of diffraction) diffraction coefficients are developed for a convex coated cylinder simulated with SIBC's and GIBC's. The ray solutions obtained remain valid in the transition region and reduce uniformly to those in the deep lit and shadow regions. A uniform asymptotic solution is also presented for observations in the close vicinity of the cylinder. The diffraction coefficient for a convex cylinder are obtained via a generalization of the corresponding ones of the circular cylinder. To validate the asymptotic/high-frequency solution, integral equations are derived for both E and H-polarization and solved numerically using the method of moments. Results are presented for a single and three layered coated convex cylinder. Some insights are also provided on the accuracy of the employed GIBC's versus SIBC's for application to curved surfaces. To characterize the scattering by impedance wedges illuminated at skew incidence, diffraction coefficients are derived from an approximate solution of the governing functional difference equations. This solution exactly recovers the known ones for an impedance half plane or an arbitrary wedge at normal incidence, and to validate it for other wedge angles, a moment method code was used. Finally, to test the usefulness of the approximate skew incidence impedance wedge diffraction coefficient for three dimensional structures, equivalent currents are derived in the context of the physical theory of diffraction (PTD) for a finite length impedance wedge of arbitrary internal angles. These are incorporated in a standard general purpose PTD code and results are presented for a number of different impedance structures.
Syed, H. H.↗
State elimination for mixed‐integer optimal control of partial differential equations by semigroup theory
Abstract Mixed‐integer optimal control problems governed by partial differential equations (MIPDECOs) are powerful modeling tools but also challenging in terms of theory and computation. We propose a highly efficient state elimination approach for MIPDECOs that are governed by partial differential equations that have the structure of an abstract ordinary differential equation in function space. This allows us to avoid repeated calculations of the states for all time steps, and our approach is applied only once before starting the optimization. The presentation of theoretical results is complemented by numerical experiments.
97 MATHEMATICS AND COMPUTING↗
Accurate parameterization of the kinetic energy functional
The absence of a reliable formulation of kinetic energy density functional has hindered the development of orbital free density functional theory. Using the data-aided learning paradigm, we propose a simple prescription to accurately model the kinetic energy density of any system. Our method relies on a dictionary of functional forms for local and nonlocal contributions which have been proposed in the literature and the appropriate coefficients are calculated via a linear regression framework. To model the nonlocal contributions, we explore two new nonlocal functionals - a functional that captures fluctuations in electronic density and a functional that incorporates gradient information. Since, the analytical functional forms of the kernels present in these nonlocal terms are not known from theory, we propose a basis function expansion to model these seemingly difficult nonlocal quantities. This allows us to easily reconstruct kernels for any system using only a few structures. The proposed method is able to learn kinetic energy densities and total kinetic energies of molecular and periodic systems, such as H 2 , LiH, LiF and a one-dimensional chain of 8 hydrogens using data from Kohn-Sham density functional theory calculations for only a few structures. For the ease of reproduction, codes used to generate the models are provided in the supporting materials.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Theoretical prediction of magnetic exchange coupling constants from broken-symmetry coupled cluster calculations
Exchange coupling constants ( J ) are fundamental to the understanding of spin spectra of magnetic systems. In this study, we investigate the broken-symmetry (BS) approaches of Noodleman and Yamaguchi in conjunction with coupled cluster (CC) methods to obtain exchange couplings. J values calculated from CC in this fashion converge smoothly toward the full configuration interaction result with increasing level of CC excitation. We compare this BS-CC scheme to the complementary equation-of-motion CC approach on a selection of bridged molecular cases and give results from a few other methodologies for context.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Foundations of variational discrete action theory
Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Coupling solidification kinetics with phase-behavior computations in hydrodynamic simulations of high-pressure, dynamic-compression processes
In this study, we report a numerical scheme to integrate models for the kinetics of solidification processes together with phase-behavior computations in the context of continuum-scale hydrodynamic simulations. The objective of the phase-behavior computations is to determine the pressure and temperature, given the following three sets of inputs: (1) an appropriate equation of state to describe our system, (2) the phase fraction(s) produced by the kinetic models, (3) and the volume and internal energy obtained by solving the conservation equations that govern the hydrodynamic behavior. The kinetics are assumed to be governed by the Kolmogorov–Johnson–Mehl–Avrami equation, and the nucleation and growth rates that enter into that equation are functions of the pressure and temperature produced by the phase-behavior computations. Our formulation allows for the fluid and solid phases to be at different temperatures (thermal nonequilibrium) and pressures (arising from surface-tension-induced Laplace contributions). The formulation is presented in a fairly general setting that is independent of any particular material, although we demonstrate it in some examples related to high-energy-density science applications where materials are rapidly compressed to pressures exceeding several gigapascals in less than a microsecond. We conclude with a critical evaluation of our approach and provide suggestions for future work to improve the predictive capabilities and generality of the models.
42 ENGINEERING↗
Stochastic renormalization group and gradient flow
A non-perturbative and continuous definition of RG transformations as stochastic processes is proposed, inspired by the observation that the functional RG equations for effective Boltzmann factors may be interpreted as Fokker-Planck equations. The result implies a new approach to Monte Carlo RG that is amenable to lattice simulation. Long-distance correlations of the effective theory are shown to approach gradient-flowed correlations, which are simpler to measure. The Markov property of the stochastic RG transformation implies an RG scaling formula which allows for the measurement of anomalous dimensions when transcribed into gradient flow expectation values.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
On the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers: Selective mode-type excitation and specific mode enhancement
Here, we investigate the in-plane vibrations and electromechanical resonance characteristics of non-uniformly polarized rectangular piezoelectric wafers. Non-uniform polarization is represented as a non-uniform electromechanical coupling coefficient using a polarization function. Governing equations are derived for the forced in-plane vibrations of a thin wafer under the assumption of generalized plane stress. The effect of non-uniform polarization is explicitly obtained in the forcing terms of the governing equations. These equations are then recast into a variational weak form that is then solved using the finite element method to obtain the displacement fields for different modes. The electromechanical response of a non-uniformly polarized piezoelectric wafer is derived in terms of the out-of-plane displacement profile on the surface of the wafer. Using the derived analytical expression, a necessary and sufficient condition for the presence/absence of a vibrational mode in the electromechanical impedance spectrum is obtained. Based on this condition, criteria for selective mode-type excitation and specific mode enhancement of vibrational modes in the electromechanical impedance spectrum are postulated. Selective mode-type excitation of in-plane extensional and shear modes is demonstrated for a square wafer and that of in-plane bending modes is demonstrated for a rectangular wafer. Specific mode enhancement is demonstrated for both square and rectangular wafers. In addition, it is also demonstrated how the criteria can be used to suppress specific vibrational modes in the electromechanical impedance spectrum. The proposed methodology of using non-uniformly polarized piezoelectric wafers finds application in the design of single element transducers with multi-frequency operation, frequency-tuned receivers/sensors, acoustic holograms, designing acoustic beams of prescribed shape/lobes, and other non-traditional applications such as information storage.
36 MATERIALS SCIENCE↗
Sudden collapse of a pressure profile generated by off-axis heating in a linear magnetized plasma
Here, the features of an unexpected, large event that arises spontaneously during a basic heat transport experiment are presented. It consists of the sudden collapse of the radial plasma pressure profile, akin to disruption events observed in toroidal magnetic confinement devices. The experiment is performed on the Large Plasma Device at the University of California, Los Angeles (UCLA). It uses a LaB6 thermionic emitter of annular shape to induce off-axis heating of a cold, afterglow plasma, in a linear magnetic geometry. The temporal evolution consists of three regimes. During an early, quiescent period, classical heat transport along and across the magnetic field arises from Coulomb collisions. After significant pressure gradients develop, drift-Alfvén waves become unstable. Upon reaching large amplitude, they trigger avalanche events that flatten the outer part of the heated region, which, in turn, quenches the instability. Due to the sustained heating, the pressure profile rebuilds and the process repeats, leading to a relatively long, second regime that displays multiple avalanches, but suddenly, the annular pressure profile is observed to collapse. After this collapse, the system enters a third regime with large fluctuations. Before the collapse, a rapid, runaway heating environment arises whose time evolution exhibits a self-similar dependence on the applied voltage. The time evolution, morphology, and scaling of the collapse event are presented, and an examination is made of the underlying mechanisms.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Relativistic hybrid stars in light of the NICER PSR J0740 + 6620 radius measurement
We explore the implications of the recent radius determination of PSR J0740 + 6620 by the NICER experiment, combined with the neutron skin measurement by the PREX-II experiment and the associated inference of the slope of symmetry energy, for the structure of hybrid stars with a strong first-order phase transition from nucleonic to quark matter. Furthermore, we combine a covariant density-functional nucleonic equation of state (EOS) with a constant-speed-of-sound EOS for quark matter. We show that the radius and tidal deformability ranges obtained from GW170817 can be reconciled with the implication of the PREX-II experiment if there is a phase transition to quark matter in the low-mass compact star. In the high-mass segment, the EOS needs to be stiff to comply with the large-radius inference for PSR J0740 + 6620 and J0030 + 0451 with masses M ≃ 2 M ⊙ and M≃1.4 M⊙. We show that twin stars are not excluded, but the mass and radius ranges (with M ≥ M ⊙ ) are restricted to narrow domains ΔM twin ≲ 0.05 M ⊙ and ΔR twin ~ 1.0 km. We also show that the existence of twin configurations is compatible with the light companion in the GW190814 event being a hybrid star in the case of values of the sound-speed square s = 0.6 and s = 1/3.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Dense plasma opacity via the multiple-scattering method
The calculation of the optical properties of hot dense plasmas with a model that has self-consistent plasma physics is a grand challenge for high energy density science. Here we exploit a recently developed electronic structure model that uses multiple scattering theory to solve the Kohn-Sham density functional theory equations for dense plasmas. Here we calculate opacities in this regime, validate the method, and apply it to recent experimental measurements of opacity for Cr, Ni, and Fe. Good agreement is found in the quasicontinuum region for Cr and Ni, while the self-consistent plasma physics of the approach cannot explain the observed difference between models and the experiment for Fe.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
pnnl/neuromancer
Dynamics-based deep learning methods to modernize current scientific computing methods. Neuromancer is currently capable of solving inverse problems for a system of ordinary differential equations. The functionality includes system identification and constrained optimal control of unknown or partially known ODEs.
Skomski, Elliott↗
Connecting Functions of Nonlinear Ordinary Differential Equations.
Abstract not provided.
Steyer, Andrew↗