SEARCH · Engineering Papers
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.
In this paper, we construct a generalization of the Jacobi iterative method adapted to a system of linear bounded operators. This new method is used with product integration method to solve a linear Fredholm integral equation of the second kind with a weakly singular kernel. The convergence analysis is proved. The numerical tests developed show its effectiveness.
This paper presents a methodology for optimizing pre-calculated collision-free paths of differential-drive wheeled robots. The main advantage of this methodology is that optimization is done by considering the kinematics and mechanical constraints of the mobile robot. In accordance to this proposal, the optimized path is achieved by applying recursively a local smoothing on an initial path which is originally modeled as a one-dimensional piecewise linear function. By this recursive smoothing, it can be ensured that the original piecewise linear function can be transformed into a smooth one that fulfill the constraints established by the kinematic equations of the wheeled mobile in terms of a minimum radius of curvature. As a result of this, a trajectory which guarantees lower power consumption and lower mechanical wear, is obtained. To show the better performance of the proposed approach, numerical simulation results are contrasted to those obtained from other reported methods with regards to path length, minimum radius of curvature, cross track error, continuity and resulting acceleration.
In this paper, the potential application of the high-level canonical piecewise-linear model in the field of three-dimensional visualization is proposed. This idea emerges from the premise that a three-dimensional surface can be modeled by a continuous two-dimensional function. In this regard, the main contribution consists in providing a methodology for the construction of such function. For this purpose, the well-known algorithm used in the high-level canonical piecewise-linear model has been modified to achieve functions with greater graphing capability. One of the most notorious changes made with respect to the traditional methodology is the use of a new scheme of function domain partition. Contrary to the simplicial partitioning, in our proposal, an equally distributed angular partitioning is considered. With this, the resulting approximate function is expressed in cylindrical coordinates which overcome the typical limitation of graphs in Cartesian coordinates of only representing single-valued surfaces. Another remarkable change can be found in the range of evaluation of the function. Instead of using an open interval, in our approach, a closed interval from 0 to 2π is adopted. To support this proposal, illustrative case studies are reported. Numerical simulations demonstrate the effectiveness of this type of representation model not only in visual perception of three-dimensional surfaces, but also in more complex tasks like volume estimation.
The long-standing resolution of the Abraham–Minkowski electromagnetic momentum controversy is predicated on a decomposition of the total momentum of a closed continuum electrodynamic system into separate field and matter components. Using a microscopic model of a simple linear dielectric, we derive Lagrangian equations of motion for the electric dipoles and show that the dielectric can be treated as a collection of stationary simple harmonic oscillators that are driven by the electric field and produce a polarization field in response. The macroscopic energy and momentum are defined in terms of the electric, magnetic, and polarization fields that travel through the dielectric together as a pulse of electromagnetic radiation. We conclude that both the macroscopic total energy and the macroscopic total momentum are entirely electromagnetic in nature for a simple linear dielectric in the absence of significant reflections. -- Highlights: •The total momentum in a dielectric is identified by conservation principles. •The total momentum in a dielectric cannot be decomposed into field and matter parts. •A component of momentum in a dielectric is due to motion of the polarization field.
The Schrödinger–Langevin equation with linear dissipation is integrated by propagating an ensemble of Bohmian trajectories for the ground state of quantum systems. Substituting the wave function expressed in terms of the complex action into the Schrödinger–Langevin equation yields the complex quantum Hamilton–Jacobi equation with linear dissipation. We transform this equation into the arbitrary Lagrangian–Eulerian version with the grid velocity matching the flow velocity of the probability fluid. The resulting equation is simultaneously integrated with the trajectory guidance equation. Then, the computational method is applied to the harmonic oscillator, the double well potential, and the ground vibrational state of methyl iodide. The excellent agreement between the computational and the exact results for the ground state energies and wave functions shows that this study provides a synthetic trajectory approach to the ground state of quantum systems.
In the context of massless quantum electrodynamics (QED) with a linear covariant gauge fixing, the connection between the counterterm and the Hopf-algebraic approach to renormalization is examined. The coproduct formula of Green’s functions contains two invariant charges, which give rise to different renormalization group functions. All formulas are tested by explicit computations to third loop order. The possibility of a finite electron self-energy by fixing a generalized linear covariant gauge is discussed. An analysis of subdivergences leads to the conclusion that such a gauge only exists in quenched QED.
Based on the Standard Model Extension, we investigate relativistic quantum effects on a scalar particle in backgrounds of the Lorentz symmetry violation defined by a tensor field. We show that harmonic-type and linear-type confining potentials can stem from Lorentz symmetry breaking effects, and thus, relativistic bound state solutions can be achieved. We first analyse a possible scenario of the violation of the Lorentz symmetry that gives rise to a harmonic-type potential. In the following, we analyse another possible scenario of the breaking of the Lorentz symmetry that induces both harmonic-type and linear-type confining potentials. In this second case, we also show that not all values of the parameter associated with the intensity of the electric field are permitted in the search for polynomial solutions to the radial equation, where the possible values of this parameter are determined by the quantum numbers of the system and the parameters associated with the violation of the Lorentz symmetry.
We review the geometry of the space of quantum states S(H) of a finite-level quantum system with Hilbert space H from a group-theoretical point of view. This space carries two stratifications generated by the action of two different Lie groups, namely, the special unitary group SU(H) and its complexification SL(H), the complex special linear group. A stratum of the stratification generated by SU(H) is composed of isospectral states, that is, density operators with the same spectrum. A stratum of the stratification generated by SL(H) is composed of quantum states with the same rank. We prove that on every submanifold of isospectral quantum states there is also a canonical left action of SL(H) which is related with the canonical Kähler structure on isospectral quantum states. The fundamental vector fields of this SL(H)-action are divided into Hamiltonian and gradient vector fields. The former give rise to invertible maps on S that preserve the von Neumann entropy and the convex structure of S(H), while the latter give rise to invertible maps on S(H) that preserve the von Neumann entropy but not the convex structure of S(H). A similar decomposition is given for the fundamental vector fields of the SL(H)-action generating the stratification of S(H) into manifolds of quantum states with the same rank. However, in this case, the gradient vector fields preserve the rank but do not preserve entropy. Finally, some comments on multipartite quantum systems are made, and it is proved that the sets of product states of a multipartite quantum system are homogeneous manifolds for the local action of the complex special linear group associated with the partition.
Optical lattice systems offer the possibility of creating and tuning Dirac points which are present in the tight-binding lattice dispersions. For example, such a behavior can be achieved in the staggered flux lattice or honeycomb type of lattices. Here we focus on the strongly correlated bosonic dynamics in the vicinity of Dirac points. In particular, we investigate bosonic Mott insulator phase in which quasiparticle excitations have a simple particle–hole interpretation. We show that linear response dynamics around Dirac points, can be significantly engineered at least in two ways: by the type of external perturbation or by changing the lattice properties. The key role is played by the interband transitions. Moreover, we explain that the behavior of these transitions is directly connected to different energy scales of the effective hopping amplitudes for particles and holes. Presented in this work theoretical study about tunability of linear response dynamics near the Dirac points can be directly simulated in the optical lattice systems.
Highlights: •The relativistic wave equations with scalar potential are investigated in (1+2)-dimensional rotating space–time. •The bound state solutions are obtained. •The relativistic energy eigenvalues associated with different radial mode are discussed. -- Abstract: In this work, we investigate the linear confinement of spin- 0 relativistic quantum particle in the backgrounds of (1+2)-dimensions Gürses space–time (M. Gürses, Class. Quantum Grav. 11, 2585 (1994)). We solve the Klein–Gordon equation with a linear scalar potential in the considered framework, and evaluate the energy eigenvalues and corresponding eigenfunctions, in detail.
Modern astronomical surveys detect asteroids by linking together their appearances across multiple images taken over time. This approach faces limitations in detecting faint asteroids and handling the computational complexity of trajectory linking. Here, we present a novel method that adapts “digital tracking” – traditionally used for short-term linear asteroid motion across images – to work with large-scale synoptic surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time (Rubin/LSST). Our approach combines hundreds of sparse observations of individual asteroids across their non-linear orbital paths to enhance detection sensitivity by several magnitudes. To address the computational challenges of processing massive data sets and dense orbital phase spaces, we developed a specialized high-performance computing architecture. We demonstrate the effectiveness of our method through experiments that take advantage of the extensive computational resources at Lawrence Livermore National Laboratory. This work enables the detection of significantly fainter asteroids in existing and future survey data, potentially increasing the observable asteroid population by orders of magnitude across different orbital families, from near-Earth objects (NEOs) to Kuiper belt objects (KBOs).
This paper sets forth a method for reachability analysis of linear dynamic systems in continuous time that can be used to compute time-domain bounds of states and outputs with floating-point precision. The focus is on the particular initial conditions and inputs that cause state or output trajectories to attain their extreme values in time. Inputs can be constant, arbitrary, or Lipschitz continuous waveforms. Uncertainties in initial conditions and inputs are modeled using zonotopes. The calculations exploit the definition of zonotopes and the modal information of the linear system dynamics.
Interaction between long DNA molecules and activated cisplatin is believed to be crucial to anticancer activity. However, the exact structural changes of long DNA molecules induced by cisplatin are still not very clear. In this study, structural changes of long linear double-stranded DNA (dsDNA) and short single-stranded DNA (ssDNA) induced by activated cisplatin have been investigated by atomic force microscopy (AFM). The results indicated that long DNA molecules gradually formed network structures, beads-on-string structures and their large aggregates. Electrostatic and coordination interactions were considered as the main driving forces producing these novel structures. An interesting finding in this study is the beads-on-string structures. Moreover, it is worth noting that the beads-on-string structures were linked into the networks, which can be ascribed to the strong DNA–DNA interactions. This study expands our knowledge of the interactions between DNA molecules and cisplatin. - Highlights: • We investigate structural changes of dsDNA and ssDNA induced by cisplatin. • AFM results indicated long dsDNA formed network, beads-on-string and aggregates. • ssDNA can form very similar structures as those of long linear dsDNA. • A possible formation process of theses novel structure is proposed.
In this work, we present a novel method for calculating interface curvature on 3D unstructured meshes from piecewise-linear interface reconstructions typically generated in the volume of fluid method. Interface curvature is a necessary quantity to calculate in order to model surface tension driven flow. Curvature needs only to be computed in cells containing an interface. The approach requires a stencil containing only neighbors sharing a node with a target cell, and calculates curvature from a least-squares paraboloid fit to the interface reconstructions. This involves solving a 6 × 6 symmetric linear system in each mixed cell. We present verification tests where we calculate the curvature of a sphere, an ellipsoid, and a sinusoid in a 3D domain on regular Cartesian meshes, distorted hex meshes, and tetrahedral meshes. Lastly, for both regular and unstructured meshes, we find in all cases the paraboloid fitting method for curvature to converge between first and second order with grid refinement.
In this work, we introduce a new residual distribution (RD) framework for the design of bound-preserving high-resolution finite element schemes. The continuous and discontinuous Galerkin discretizations of the linear advection equation are modified to construct local extremum diminishing (LED) approximations. To that end, we perform mass lumping and redistribute the element residuals in a manner which guarantees the LED property. The hierarchical correction procedure for high-order Bernstein finite element discretizations involves localization to subcells and definition of bound-preserving weights for subcell contributions. Using strong stability preserving (SSP) Runge–Kutta methods for time integration, we prove the validity of discrete maximum principles under CFL-like time step restrictions. The low-order version of our method has roughly the same accuracy as the one derived from a piecewise (multi)-linear approximation on a submesh with the same nodal points. In high-order extensions, we use an element-based flux-corrected transport (FCT) algorithm which can be interpreted as a nonlinear RD scheme. The proposed LED corrections are tailor-made for matrix-free implementations which avoid the rapidly growing cost of matrix assembly for high-order Bernstein elements. The results for 1D, 2D, and 3D test problems compare favorably to those obtained with the best matrix-based approaches.
Here, we formulate a new projection-based reduced-order modeling technique for non-linear dynamical systems. The proposed technique, which we refer to as the Adjoint Petrov–Galerkin (APG) method, is derived by decomposing the generalized coordinates of a dynamical system into a resolved coarse-scale set and an unresolved fine-scale set. A Markovian finite memory assumption within the Mori–Zwanzig formalism is then used to develop a reduced-order representation of the coarse scales. This procedure leads to a closed reduced-order model that displays commonalities with the adjoint stabilization method used in finite elements. The formulation is shown to be equivalent to a Petrov–Galerkin method with a non-linear, time-varying test basis, thus sharing some similarities with the Least-Squares Petrov–Galerkin method. Theoretical analysis examining a priori error bounds and computational cost is presented. Numerical experiments on the compressible Navier–Stokes equations demonstrate that the proposed method can lead to improvements in numerical accuracy, robustness, and computational efficiency over the Galerkin method on problems of practical interest. Improvements in numerical accuracy and computational efficiency over the Least-Squares Petrov–Galerkin method are observed in most cases.
In a previous work, embedded ensemble propagation was proposed to improve the efficiency of sampling-based uncertainty quantification methods of computational models on emerging computational architectures. It consists of simultaneously evaluating the model for a subset of samples together, instead of evaluating them individually. A first approach introduced to solve parametric linear systems with ensemble propagation is ensemble reduction. In Krylov methods for example, this reduction consists in coupling the samples together using an inner product that sums the sample contributions. Ensemble reduction has the advantages of being able to use optimized implementations of BLAS functions and having a stopping criterion which involves only one scalar. However, the reduction potentially decreases the rate of convergence due to the gathering of the spectra of the samples. In this paper, we investigate a second approach: ensemble propagation without ensemble reduction in the case of GMRES. This second approach solves each sample simultaneously but independently to improve the convergence compared to ensemble reduction. This raises two new issues which are solved in this paper: the fact that optimized implementations of BLAS functions cannot be used anymore and that ensemble divergence, whereby individual samples within an ensemble must follow different code execution paths, can occur. We tackle those issues by implementing a high-performing ensemble GEMV and by using masks. The proposed ensemble GEMV leads to a similar cost per GMRES iteration for both approaches, i.e. with and without reduction. For illustration, we study the performances of the new linear solver in the context of a mesh tying problem. Furthermore, this example demonstrates improved ensemble propagation speed-up without reduction.
In this paper, we present a first-order Stress-Hybrid Virtual Element Method (SH-VEM) on six-noded triangular meshes for linear plane elasticity. Here, we adopt the Hellinger–Reissner variational principle to construct a weak equilibrium condition and a stress based projection operator. In each element, the stress projection operator is expressed in terms of the nodal displacements, which leads to a displacement based formulation. This stress-hybrid approach assumes a globally continuous displacement field while the stress field is discontinuous across each element. The stress field is initially represented by divergence-free tensor polynomials based on Airy stress functions, but we also present a formulation that uses a penalty term to enforce the element equilibrium conditions, referred to as the Penalty Stress-Hybrid Virtual Element Method (PSH-VEM). Numerical results are presented for PSH-VEM and SH-VEM, and we compare their convergence to the composite triangle FEM and B-bar VEM on benchmark problems in linear elasticity. The SH-VEM converges optimally in the L 2 norm of the displacement, energy seminorm, and the L 2 norm of hydrostatic stress. Furthermore, the results reveal that PSH-VEM converges in most cases at a faster rate than the expected optimal rate, but it requires the selection of a suitably chosen penalty parameter.