Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Calculus of Variation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

The exit-time problem for a Markov jump process

The purpose of this paper is to consider the exit-time problem for a finite-range Markov jump process, i.e, the distance the particle can jump is bounded independent of its location. Additionally, such jump diffusions are expedient models for anomalous transport exhibiting super-diffusion or nonstandard normal diffusion. We refer to the associated deterministic equation as a volume-constrained nonlocal diffusion equation. The volume constraint is the nonlocal analogue of a boundary condition necessary to demonstrate that the nonlocal diffusion equation is well-posed and is consistent with the jump process. A critical aspect of the analysis is a variational formulation and a recently developed nonlocal vector calculus. Finally, this calculus allows us to pose nonlocal backward and forward Kolmogorov equations, the former equation granting the various moments of the exit-time distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Identification of the Diffusion Parameter in Nonlocal Steady Diffusion Problems

The problem of identifying the diffusion parameter appearing in a nonlocal steady diffusion equation is considered. The identification problem is formulated as an optimal control problem having a matching functional as the objective of the control and the parameter function as the control variable. The analysis makes use of a nonlocal vector calculus that allows one to define a variational formulation of the nonlocal problem. In a manner analogous to the local partial differential equations counterpart, we demonstrate, for certain kernel functions, the existence of at least one optimal solution in the space of admissible parameters. We introduce a Galerkin finite element discretization of the optimal control problem and derive a priori error estimates for the approximate state and control variables. Using one-dimensional numerical experiments, we illustrate the theoretical results and show that by using nonlocal models it is possible to estimate non-smooth and discontinuous diffusion parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Fluctuations, Saturation, and Diffractive Excitation in High Energy Collisions

Diffractive excitation is usually described by the Good-Walker formalism for low masses, and by the triple-Regge formalism for high masses. In the Good-Walker formalism the cross section is determined by the fluctuations in the interaction. By taking the fluctuations in the BFKL ladder into account, it is possible to describe both low and high mass excitation in the Good-Walker formalism. In high energy pp collisions the fluctuations are strongly suppressed by saturation, which implies that pomeron exchange does not factorise between DIS and pp collisions. The Dipole Cascade Model reproduces the expected triple-Regge form for the bare pomeron, and the triple-pomeron coupling is estimated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A Unified Theory of Fractional Nonlocal and Weighted Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green's identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green's identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. The proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Self-consistent field theory of collisions: Orbital equations with asymptotic sources and self-averaged potentials

The self-consistent field theory of collisions is formulated, incorporating the unique dynamics generated by the self-averaged potentials. The bound state Hartree–Fock approach is extended for the first time to scattering states, by properly resolving the principal difficulties of non-integrable continuum orbitals and imposing complex asymptotic conditions. The recently developed asymptotic source theory provides the natural theoretical basis, as the asymptotic conditions are completely transferred to the source terms and the new scattering function is made fullyintegrable. The scattering solutions can then be directly expressed in terms of bound state HF configurations, establishing the relationship between the bound and scattering state solutions. Alternatively, the integrable spin orbitals are generated by constructing the individual orbital equations that contain asymptotic sources and self-averaged potentials. However, the orbital energies are not determined by the equations, and a special channel energy fixing procedure is developed to secure the solutions. It is also shown that the variational construction of the orbital equations has intrinsic ambiguities that are generally associated with the self-consistent approach. On the other hand, when a small subset of open channels is included in the source term, the solutions are only partiallyintegrable, but the individual open channels can then be treated more simply by properly selecting the orbital energies. The configuration mixing and channel coupling are then necessary to complete the solution. The new theory improves the earlier continuum HF model. - Highlights: • First extension of HF to scattering states, with proper asymptotic conditions. • Orbital equations with asymptotic sources and integrable orbital solutions. • Construction of self-averaged potentials, and orbital energy fixing. • Channel coupling and configuration mixing, involving the new orbitals. • Critical evaluation of the continuum HF and its improvements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Towards a Unified theory of Fractional and Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green’s identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green’s identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. Here, the proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Ponderomotive dynamics of waves in quasiperiodically modulated media

Similarly to how charged particles experience time-averaged ponderomotive forces in high-frequency fields, linear waves also experience time-averaged refraction in modulated media. We propose a covariant variational theory of this ponderomotive effect on waves for a general nondissipative linear medium. Using the Weyl calculus, our formulation accommodates waves with temporal and spatial period comparable to that of the modulation (provided that parametric resonances are avoided). This theory also shows that any wave is, in fact, a polarizable object that contributes to the linear dielectric tensor of the ambient medium. Furthermore, the dynamics of quantum particles is subsumed as a special case. As an illustration, ponderomotive Hamiltonians of quantum particles and photons are calculated within a number of models. We also explain a fundamental connection between these results and the well-known electrostatic dielectric tensor of quantum plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme

A fully variational, unstructured, electromagnetic particle-in-cell integrator is developed for integration of the Vlasov-Maxwell equations. Using the formalism of discrete exterior calculus [Desbrun et al., e-print arXiv:math/0508341 (2005)], the field solver, interpolation scheme, and particle advance algorithm are derived through minimization of a single discrete field theory action. As a consequence of ensuring that the action is invariant under discrete electromagnetic gauge transformations, the integrator exactly conserves Gauss's law.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Variational Integration for Ideal Magnetohydrodynamics and Formation of Current Singularities

Coronal heating has been a long-standing conundrum in solar physics. Parker's conjecture that spontaneous current singularities lead to nanoflares that heat the corona has been controversial. In ideal magnetohydrodynamics (MHD), can genuine current singularities emerge from a smooth 3D line-tied magnetic field? To numerically resolve this issue, the schemes employed must preserve magnetic topology exactly to avoid artificial reconnection in the presence of (nearly) singular current densities. Structure-preserving numerical methods are favorable for mitigating numerical dissipation, and variational integration is a powerful machinery for deriving them. However, successful applications of variational integration to ideal MHD have been scarce. In this thesis, we develop variational integrators for ideal MHD in Lagrangian labeling by discretizing Newcomb's Lagrangian on a moving mesh using discretized exterior calculus. With the built-in frozen-in equation, the schemes are free of artificial reconnection, hence optimal for studying current singularity formation. Using this method, we first study a fundamental prototype problem in 2D, the Hahm-Kulsrud-Taylor (HKT) problem. It considers the effect of boundary perturbations on a 2D plasma magnetized by a sheared field, and its linear solution is singular. We find that with increasing resolution, the nonlinear solution converges to one with a current singularity. The same signature of current singularity is also identified in other 2D cases with more complex magnetic topologies, such as the coalescence instability of magnetic islands. We then extend the HKT problem to 3D line-tied geometry, which models the solar corona by anchoring the field lines in the boundaries. The effect of such geometry is crucial in the controversy over Parker's conjecture. The linear solution, which is singular in 2D, is found to be smooth. However, with finite amplitude, it can become pathological above a critical system length. The nonlinear solution turns out smooth for short systems. Nonetheless, the scaling of peak current density vs. system length suggests that the nonlinear solution may become singular at a finite length. With the results in hand, we cannot confirm or rule out this possibility conclusively, since we cannot obtain solutions with system lengths near the extrapolated critical value.

solar corona↗
DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Classification of singularities in the problem of motion of the Kovalevskaya top in a double force field

The problem of motion of the Kovalevskaya top in a double force field is investigated (the integrable case of A.G. Reyman and M.A. Semenov-Tian-Shansky without a gyrostatic momentum). It is a completely integrable Hamiltonian system with three degrees of freedom not reducible to a family of systems with two degrees of freedom. The critical set of the integral map is studied. The critical subsystems and bifurcation diagrams are described. The classification of all nondegenerate critical points is given. The set of these points consists of equilibria (nondegenerate singularities of rank 0), of singular periodic motions (nondegenerate singularities of rank 1), and also of critical two-frequency motions (nondegenerate singularities of rank 2). Bibliography: 32 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme

A fully variational, unstructured, electromagnetic particle-in-cell integrator is developed for integration of the Vlasov-Maxwell equations. Using the formalism of Discrete Exterior Calculus, the field solver, interpolation scheme and particle advance algorithm are derived through minimization of a single discrete field theory action. As a consequence of ensuring that the action is invariant under discrete electromagnetic gauge transformations, the integrator exactly conserves Gauss's law.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Mineral density volume gradients in normal and diseased human tissues

Clinical computed tomography provides a single mineral density (MD) value for heterogeneous calcified tissues containing early and late stage pathologic formations. The novel aspect of this study is that, it extends current quantitative methods of mapping mineral density gradients to three dimensions, discretizes early and late mineralized stages, identifies elemental distribution in discretized volumes, and correlates measured MD with respective calcium (Ca) to phosphorus (P) and Ca to zinc (Zn) elemental ratios. To accomplish this, MD variations identified using polychromatic radiation from a high resolution micro-computed tomography (micro-CT) benchtop unit were correlated with elemental mapping obtained from a microprobe X-ray fluorescence (XRF) using synchrotron monochromatic radiation. Digital segmentation of tomograms from normal and diseased tissues (N=5 per group; 40-60 year old males) contained significant mineral density variations (enamel: 2820-3095mg/cc, bone: 570-1415mg/cc, cementum: 1240-1340mg/cc, dentin: 1480-1590mg/cc, cementum affected by periodontitis: 1100-1220mg/cc, hypomineralized carious dentin: 345-1450mg/cc, hypermineralized carious dentin: 1815-2740mg/cc, and dental calculus: 1290-1770mg/cc). A plausible linear correlation between segmented MD volumes and elemental ratios within these volumes was established, and Ca/P ratios for dentin (1.49), hypomineralized dentin (0.32-0.46), cementum (1.51), and bone (1.68) were observed. Furthermore, varying Ca/Zn ratios were distinguished in adapted compared to normal tissues, such as in bone (855-2765) and in cementum (595-990), highlighting Zn as an influential element in prompting observed adaptive properties. Hence, results provide insights on mineral density gradients with elemental concentrations and elemental footprints that in turn could aid in elucidating mechanistic processes for pathologic formations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Computer-aided detection of renal calculi from noncontrast CT images using TV-flow and MSER features

Purpose: Renal calculi are common extracolonic incidental findings on computed tomographic colonography (CTC). This work aims to develop a fully automated computer-aided diagnosis system to accurately detect renal calculi on CTC images. Methods: The authors developed a total variation (TV) flow method to reduce image noise within the kidneys while maintaining the characteristic appearance of renal calculi. Maximally stable extremal region (MSER) features were then calculated to robustly identify calculi candidates. Finally, the authors computed texture and shape features that were imported to support vector machines for calculus classification. The method was validated on a dataset of 192 patients and compared to a baseline approach that detects calculi by thresholding. The authors also compared their method with the detection approaches using anisotropic diffusion and nonsmoothing. Results: At a false positive rate of 8 per patient, the sensitivities of the new method and the baseline thresholding approach were 69% and 35% (p < 1e − 3) on all calculi from 1 to 433 mm{sup 3} in the testing dataset. The sensitivities of the detection methods using anisotropic diffusion and nonsmoothing were 36% and 0%, respectively. The sensitivity of the new method increased to 90% if only larger and more clinically relevant calculi were considered. Conclusions: Experimental results demonstrated that TV-flow and MSER features are efficient means to robustly and accurately detect renal calculi on low-dose, high noise CTC images. Thus, the proposed method can potentially improve diagnosis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

Matter wave switching in Bose-Einstein condensates via intensity redistribution soliton interactions

Using time dependent nonlinear (s-wave scattering length) coupling between the components of a weakly interacting two component Bose-Einstein condensate (BEC), we show the possibility of matter wave switching (fraction of atoms transfer) between the components via shape changing/intensity redistribution (matter redistribution) soliton interactions. We investigate the exact bright-bright N-soliton solution of an effective one-dimensional (1D) two component BEC by suitably tailoring the trap potential, atomic scattering length, and atom gain or loss. In particular, we show that the effective 1D coupled Gross-Pitaevskii equations with time dependent parameters can be transformed into the well known completely integrable Manakov model described by coupled nonlinear Schroedinger equations by effecting a change of variables of the coordinates and the wave functions under certain conditions related to the time dependent parameters. We obtain the one-soliton solution and demonstrate the shape changing/matter redistribution interactions of two and three-soliton solutions for the time-independent expulsive harmonic trap potential, periodically modulated harmonic trap potential, and kinklike modulated harmonic trap potential. The standard elastic collision of solitons occur only for a specific choice of soliton parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Hydrogen recombination due to collisions with He and Ar

Quantum-mechanical calculations are reported for hydrogen recombination in the presence of a chemically inert spectator. The calculations employ a square integrable Sturmian basis set to provide a discrete representation of the H{sub 2} continuum. Direct three-body recombination is approximated by computing transitions from the nonresonant continuum. Resonant and nonresonant states are handled on equal footing within the sequential two-step energy transfer mechanism. Theoretical rate coefficients are computed within the equilibrium and steady-state approximations for the density of intermediate molecules. The results are compared with existing experimental data for He and Ar. The sensitivity of the calculations to pressure variations and to changes in the potential energy surface are investigated for He. The reliability of these calculations and their relevance for astrophysical models are discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2010

Correlations of RMT characteristic polynomials and integrability: Hermitean matrices

Integrable theory is formulated for correlation functions of characteristic polynomials associated with invariant non-Gaussian ensembles of Hermitean random matrices. By embedding the correlation functions of interest into a more general theory of {tau} functions, we (i) identify a zoo of hierarchical relations satisfied by {tau} functions in an abstract infinite-dimensional space and (ii) present a technology to translate these relations into hierarchically structured nonlinear differential equations describing the correlation functions of characteristic polynomials in the physical, spectral space. Implications of this formalism for fermionic, bosonic, and supersymmetric variations of zero-dimensional replica field theories are discussed at length. A particular emphasis is placed on the phenomenon of fermionic-bosonic factorisation of random-matrix-theory correlation functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Convective mixing in vertically-layered porous media: The linear regime and the onset of convection

In this paper, we study the effect of permeability heterogeneity on the stability of gravitationally unstable, transient, diffusive boundary layers in porous media. Permeability is taken to vary periodically in the horizontal plane normal to the direction of gravity. In contrast to the situation for vertical permeability variation, the horizontal perturbation structures are multimodal. We therefore use a two-dimensional quasi-steady eigenvalue analysis as well as a complementary initial value problem to investigate the stability behavior in the linear regime, until the onset of convection. We find that thick permeability layers enhance instability compared with thin layers when heterogeneity is increased. On the contrary, for thin layers the instability is weakened progressively with increasing heterogeneity to the extent that the corresponding homogeneous case is more unstable. For high levels of heterogeneity, we find that a small change in the permeability field results in large variations in the onset time of convection, similar to the instability event in the linear regime. However, this trend does not persist unconditionally because of the reorientation of vorticity pairs due to the interaction of evolving perturbation structures with heterogeneity. Consequently, an earlier onset of instability does not necessarily imply an earlier onset of convection. A resonant amplification of instability is observed within the linear regime when the dominant perturbation mode is equal to half the wavenumber of permeability variation. On the other hand, a substantial damping occurs when the perturbation mode is equal to the harmonic and sub-harmonic components of the permeability wavenumber. Finally, the phenomenon of such harmonic interactions influences both the onset of instability as well as the onset of convection.

36 MATERIALS SCIENCE↗