Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “non convex”

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 163 records · Page 9

Dynamical Complexity of Non-Gaussian Many-Body Systems with Dissipation

We characterize the dynamical state of many-body bosonic and fermionic many-body models with intersite Gaussian couplings, on-site non-Gaussian interactions, and local dissipation comprising incoherent particle loss, particle gain, and dephasing. We first establish that, for fermionic systems, if the dephasing noise is larger than the non-Gaussian interactions, irrespective of the Gaussian coupling strength, the system state is a convex combination of Gaussian states at all times. Furthermore, for bosonic systems, we show that if the particle loss and particle gain rates are larger than the Gaussian intersite couplings, the system remains in a separable state at all times. Building on this characterization, we establish that at noise rates above a threshold, there exists a classical algorithm that can efficiently sample from the system state of both the fermionic and bosonic models. Finally, we show that, unlike fermionic systems, bosonic systems can evolve into states that are not convex Gaussian even when the dissipation is much higher than the on-site non-Gaussianity. Similarly, unlike bosonic systems, fermionic systems can generate entanglement even with noise rates much larger than the intersite couplings.

Computational complexity↗

Minimal-Area Metrics on the Swiss Cross and Punctured Torus

The closed string theory minimal-area problem asks for the conformal metric of least area on a Riemann surface with the condition that all non-contractible closed curves have length greater than or equal to an arbitrary constant that can be set to one. Through every point in such a metric there is a geodesic that saturates the length condition, and the saturating geodesics in a given homotopy class form a band. The extremal metric is unknown when bands of geodesics cross, as it happens for surfaces of non-zero genus. Here, we use recently proposed convex programs to numerically find the minimal-area metric on the square torus with a square boundary, for various sizes of the boundary. For large enough boundary the problem is equivalent to the “Swiss cross” challenge posed by Strebel. We find that the metric is positively curved in the two-band region and flat in the single-band regions. For small boundary the metric develops a third band of geodesics wrapping around it, and has both regions of positive and of negative curvature. This surface can be completed to provide the minimal-area metric on a once-punctured torus, representing a closed-string tadpole diagram.

97 MATHEMATICS AND COMPUTING↗

Covariant bit threads

We derive several new reformulations of the Hubeny-Rangamani-Takayanagi covariant holographic entanglement entropy formula. These include: (1) a minimax formula, which involves finding a maximal-area achronal surface on a timelike hypersurface homologous to D(A) (the boundary causal domain of the region A whose entropy we are calculating) and minimizing over the hypersurface; (2) a max V-flow formula, in which we maximize the flux through D(A) of a divergenceless bulk 1-form V subject to an upper bound on its norm that is non-local in time; and (3) a min U-flow formula, in which we minimize the flux over a bulk Cauchy slice of a divergenceless timelike 1-form U subject to a lower bound on its norm that is non-local in space. The two flow formulas define convex programs and are related to each other by Lagrange duality. For each program, the optimal configurations dynamically find the HRT surface and the entanglement wedges of A and its complement. The V-flow formula is the covariant version of the Freedman-Headrick bit thread reformulation of the Ryu-Takayanagi formula. We also introduce a measure-theoretic concept of a “thread distribution”, and explain how Riemannian flows, V-flows, and U-flows can be expressed in terms of thread distributions.

79 ASTRONOMY AND ASTROPHYSICS↗

Running decompactification, sliding towers, and the distance conjecture

We study towers of light particles that appear in infinite-distance limits of moduli spaces of 9-dimensional $\mathcal{N}$ = 1 string theories, some of which notably feature decompactification limits with running string coupling. The lightest tower in such decompactification limits consists of the non-BPS Kaluza-Klein modes of Type I' string theory, whose masses depend nontrivially on the moduli of the theory. We work out the moduli-dependence by explicit computation, finding that despite the running decompactification the Distance Conjecture remains satisfied with an exponential decay rate α ≥ 1 / $\sqrt{d-2}$ in accordance with the sharpened Distance Conjecture. The related sharpened Convex Hull Scalar Weak Gravity Conjecture also passes stringent tests. Our results non-trivially test the Emergent String Conjecture, while highlighting the important subtlety that decompactifcation can lead to a running solution rather than to a higher-dimensional vacuum.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Accelerated Sparse Recovery via Gradient Descent with Nonlinear Conjugate Gradient Momentum

This paper applies an idea of adaptive momentum for the nonlinear conjugate gradient to accelerate optimization problems in sparse recovery. Specifically, we consider two types of minimization problems: a (single) differentiable function and the sum of a non-smooth function and a differentiable function. In the first case, we adopt a fixed step size to avoid the traditional line search and establish the convergence analysis of the proposed algorithm for a quadratic problem. This acceleration is further incorporated with an operator splitting technique to deal with the non-smooth function in the second case. As a result, we use the convex ι 1 and the nonconvex ι 1 – ι 2 functionals as two case studies to demonstrate the efficiency of the proposed approaches over traditional methods.

97 MATHEMATICS AND COMPUTING↗

Improved image decompression for reduced transform coding artifacts

The perceived quality of images reconstructed from low bit rate compression is severely degraded by the appearance of transform coding artifacts. This paper proposes a method for producing higher quality reconstructed images based on a stochastic model for the image data. Quantization (scalar or vector) partitions the transform coefficient space and maps all points in a partition cell to a representative reconstruction point, usually taken as the centroid of the cell. The proposed image estimation technique selects the reconstruction point within the quantization partition cell which results in a reconstructed image which best fits a non-Gaussian Markov random field (MRF) image model. This approach results in a convex constrained optimization problem which can be solved iteratively. At each iteration, the gradient projection method is used to update the estimate based on the image model. In the transform domain, the resulting coefficient reconstruction points are projected to the particular quantization partition cells defined by the compressed image. Experimental results will be shown for images compressed using scalar quantization of block DCT and using vector quantization of subband wavelet transform. The proposed image decompression provides a reconstructed image with reduced visibility of transform coding artifacts and superior perceived quality.

Orourke, Thomas P.↗

A positivity-preserving and conservative intersection-distribution-based remapping algorithm for staggered ALE hydrodynamics on arbitrary meshes

In this study, we introduce new intersection-distribution-based remapping tools for indirect staggered arbitrary Lagrangian-Eulerian (ALE) simulations of multi-material shock hydrodynamics on arbitrary meshes. In addition to conserving momentum and total energy, the three-stage remapper proposed in this work preserves non-negativity of the internal energy. At the first stage, we construct slope-limited piecewise-linear reconstructions of all conserved quantities on zones of the source mesh and perform intersection-based remap to obtain bound-preserving zonal quantities on the target mesh. At the second stage, we define bound-preserving nodal quantities of the staggered ALE discretization as convex combinations of corner quantities. The nodal internal energy is corrected in a way which keeps it non-negative, while providing exact conservation of total energy. At the final stage, we distribute the non-negative nodal internal energy to corners, zones and materials using non-negative weights. Proofs of positivity preservation are provided for each stage. This work is a natural extension of our paper [14] in which a similar intersection-distribution-based remapping procedure was employed. The original version used a nodal kinetic energy fix which did not provably ensure positivity preservation for the zonal internal energy after the final distribution stage. The new algorithm cures this potential drawback by using ‘coordinated’ limiters for piecewise-linear reconstructions, remapping the internal energy to nodes and correcting it before redistribution. The effectiveness of the new nodal fix is illustrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗

Transmission interference fringe ( TIF ) technique for the dynamic visualization of evaporating droplet

Here, the transmission interference fringe (TIF) technique was developed to visualize the dynamics of evaporating droplets based on the Reflection Interference Fringe (RIF) technique for micro-sized droplets. The geometric formulation was conducted to determine the contact angle (CA) and height of macro-sized droplets without the need for the prism used in RIF. The TIF characteristics were analyzed through experiments and simulations to demonstrate a wider range of contact angles from 0 to 90°, in contrast to RIF's limited range of 0–30°. TIF was utilized to visualize the dynamic evaporation of droplets in the constant contact radius (CCR) mode, observing the droplet profile change from convex-only to convex-concave at the end of dry-out from the interference fringe formation. The TIF also observed the contact angle increase from the fringe radius increase. This observation is uniquely reported as the interference fringe (IF) technique can detect the formation of interference fringe between the reflection from the center convex profile and the reflection from the edge concave profile on the far-field screen. Unlike general microscopy techniques, TIF can detect far-field interference fringes as it focuses beyond the droplet-substrate interface. The formation of the convex-concave profile during CCR evaporation is believed to be influenced by the non-uniform evaporative flux along the droplet surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nastic curvatures of wheat coleoptiles that develop in true microgravity

Dark-grown wheat coleoptiles developed strong curvatures within 5 h of being transferred in orbit from a 1 g centrifuge to microgravity during an experiment flown on the IML-1 shuttle mission. The curving tendency was strongest in seedlings that were immature, with coleoptiles shorter than 10 mm at the time of transfer. The curvature direction was non-random, and directed away from the caryopsis (the coleptile face adjacent to the caryopsis becoming convex). The curvatures were most marked in the basal third of the coleoptiles, contrasting with phototropic responses, which occur in the apical third. We interpret these curvatures as being nastic, and related to the curvatures commonly reported to occur during clinostat rotation treatments.

manned↗

Symbiotic niche mapping reveals functional specialization by two ectomycorrhizal fungi that expands the host plant niche

Mutualisms are ubiquitous in natural systems, but less is known about how these positive interactions influence species distributions compared with antagonistic interactions, such as competition and predation. The niche concept is one useful approach for thinking about factors that shape species ranges, which we apply here towards understanding how the nature of plant-mycorrhizal symbioses change across large environmental gradients. We used a continuous niche mapping approach to examine how two ectomycorrhizal fungi (Thelephora terrestris, Suillus pungens) impact pine seedling growth across a two-dimensional soil nitrogen (N) and phosphorus (P) gradient. We found that ectomycorrhizal fungi improved plant growth most in nutrient addition treatments with highly imbalanced N:P ratios, demonstrating that mycorrhizal benefits depend on interactions between niche axes. Nevertheless, T. terrestris (highN:lowP) and S. pungens (lowN:highP) benefited plants most at opposite ends of the resource ratio spectrum, consistent with niche partitioning and functional specialization. While ectomycorrhizal fungi are often thought of as being most beneficial for nitrogen uptake, our results suggest that members of the Thelephoraceae may specialize in improving plant P uptake. Ectomycorrhizal colonization by a single fungus increased plant niche volume (calculated as convex hull volumes of plant growth response surfaces across N and P gradients) compared to non-mycorrhizal control plants and shows the overall positive effects of mutualisms on plant niche volume. Despite plant host benefits in S. pungens and T. terrestris single species treatments, the presence of both fungi together decreased plant niche volume. The lack of functional complementarity, despite functional specialization, indicates that in some environments, either fungal competition or the cost of maintaining a suboptimal mycorrhizal partner can limit the benefits of a higher quality partner. The niche mapping approach we present has the potential to answer fundamental questions about the dimensions of functional diversity in ectomycorrhizal fungi and the distributions of mycorrhizal symbioses.

59 BASIC BIOLOGICAL SCIENCES↗

Gratings Fabricated on Flat Surfaces and Reproduced on Non-Flat Substrates

A method has been developed for fabricating gratings on flat substrates, and then reproducing the groove pattern on a curved (concave or convex) substrate and a corresponding grating device. First, surface relief diffraction grating grooves are formed on flat substrates. For example, they may be fabricated using photolithography and reactive ion etching, maskless lithography, holography, or mechanical ruling. Then, an imprint of the grating is made on a deformable substrate, such as plastic, polymer, or other materials using thermoforming, hot or cold embossing, or other methods. Interim stamps using electroforming, or other methods, may be produced for the imprinting process or if the same polarity of the grating image is required. The imprinted, deformable substrate is then attached to a curved, rigid substrate using epoxy or other suitable adhesives. The imprinted surface is facing away from the curved rigid substrate. As an alternative fabrication method, after grating is imprinted on the deformable substrate as described above, the grating may be coated with thin conformal conductive layer (for example, using vacuum deposition of gold). Then the membrane may be mounted over an opening in a pressured vessel in a manner of a membrane on a drum, grating side out. The pressure inside of the vessel may be changed with respect to the ambient pressure to produce concave or convex membrane surface. The shape of the opening may control the type of the surface curvature (for example, a circular opening would create spherical surface, oval opening would create toroidal surface, etc.). After that, well-known electroforming methods may be used to create a replica of the grating on the concave or convex membrane. For example, the pressure vessel assembly may be submerged into an electro-forming solution and negative electric potential applied to the metal coated membrane using an insulated wire. Positive electric potential may be then applied to a nickel or other metal plate submerged into the same solution. Metal ions would transfer from the plate through the solution into the membrane, producing high fidelity metal replica of the grating on the membrane. In one variation, an adhesive may be deposited on the deformable substrate, and then cured without touching the rigid, curved substrate. Edges of the deformable substrate may be attached to the rigid substrate to ensure uniform deformation of the deformable substrate. The assembly may be performed in vacuum, and then taken out to atmospheric pressure conditions to ensure that no air is trapped between the deformable and rigid substrates. Alternatively, a rigid surface with complementary curvature to the rigid substrate may be used to ensure uniform adhesion of the deformable substrate to the rigid substrate. Liquid may be applied to the surface of the deformable substrate to uniformly distribute pressure across its surface during the curing or hardening of the adhesive, or the film may be pressed into the surface using a deformable object or surface. After the attachment is complete, the grooves may be coated with reflective or dielectric layers to improve diffraction efficiency.

Content, David↗

Shape Preserving Spline Interpolation

A rational spline solution to the problem of shape preserving interpolation is discussed. The rational spline is represented in terms of first derivative values at the knots and provides an alternative to the spline-under-tension. The idea of making the shape control parameters dependent on the first derivative unknowns is then explored. The monotonic or convex shape of the interpolation data can then be preserved automatically through the solution of the resulting non-linear consistency equations of the spline.

Gregory, J. A.↗

Machine learning guided search for energetically favorable metal borocarbide ternary compounds

In this work, we employ machine-learning (ML) combined with first principles calculations to efficiently search for the energetically favorable metal borocarbide (M-B-C) ternary compounds with M being the group 1–3 metal elements. Using a crystal graph convolutional neural network (CGCNN) ML approach followed by first-principles calculations, we predicted 47 energetically favorable stable and metastable ternary Na-B-C, Ca-B-C, and La-B-C ternary compounds with their decomposition energy (E d ) below or within 100 meV/atom from the currently known convex hulls. Phonon spectra and electronic structures of the 14 energetically favorable stable structures are also investigated by first-principles calculations. By substituting the metal atoms in the 29 energetically favorable non-equivalent template structures of Na (Ca, La)-B-C with other group 1–3 elements in the periodic table, we further obtain 22 stable structures and 52 metastable structures (E d ≤100 meV/atom with respect to the known convex hulls) for Li-B-C, K-B-C, Rb-B-C, Mg-B-C, Sr-B-C, Ba-B-C, Sc-B-C and Y-B-C ternary compounds. New convex hulls including our newly predicted stable ternary structures and the known stable structures are constructed for the M-B-C systems. The results obtained by our ML guided first-principles calculations enrich our knowledge in the structure and energy landscape of metal borocarbide ternary compounds and provide useful guidance for further experimental synthesis and discovery.

36 MATERIALS SCIENCE↗

Convex Programs for Minimal-Area Problems

The minimal-area problem that defines string diagrams in closed string field theory asks for the metric of least area on a Riemann surface with the condition that all non-contractible closed curves have length at least 2π. This is an extremal length problem in conformal geometry as well as a problem in systolic geometry. Here, we consider the analogous minimal-area problem for homology classes of curves and, with the aid of calibrations and the max flow-min cut theorem, formulate it as a local convex program. We derive an equivalent dual program involving maximization of a concave functional. These two programs give new insights into the form of the minimal-area metric and are amenable to numerical solution. We explain how the homology problem can be modified to provide the solution to the original homotopy problem.

97 MATHEMATICS AND COMPUTING↗

Computational optimal transport for molecular spectra: The fully continuous case

Computational optimal transport is used to analyze the difference between pairs of continuous molecular spectra. It is demonstrated that transport distances which are derived from this approach may be a more appropriate measure of the difference between two continuous spectra than more familiar measures of distance under many common circumstances. Associated with the transport distances is the transport map which provides a detailed analysis of the difference between two molecular spectra and is a key component of our study of quantitative differences between two continuous spectra. The use of optimal transport for comparing molecular spectra is developed in detail here with a set of model spectra, so that the discussion is self-contained. The difference between the transport distance and more common definitions of distance is elucidated for some well-chosen examples and it is shown where transport distances may be very useful alternatives to standard definitions of distance. The transport distance between a theoretical and experimental electronic absorption spectrum for SO 2 is studied and it is shown how the theoretical spectrum can be modified to fit the experimental spectrum better adjusting the theoretical band origin and the resolution of the theoretical spectrum. In conclusion, this analysis includes the calculation of transport maps between the theoretical and experimental spectra suggesting future applications of the methodology.

74 ATOMIC AND MOLECULAR PHYSICS↗

Refinement of a Method for Identifying Probable Archaeological Sites from Remotely Sensed Data

To facilitate locating archaeological sites before they are compromised or destroyed, we are developing approaches for generating maps of probable archaeological sites, through detecting subtle anomalies in vegetative cover, soil chemistry, and soil moisture by analyzing remotely sensed data from multiple sources. We previously reported some success in this effort with a statistical analysis of slope, radar, and Ikonos data (including tasseled cap and NDVI transforms) with Student's t-test. We report here on new developments in our work, performing an analysis of 8-band multispectral Worldview-2 data. The Worldview-2 analysis begins by computing medians and median absolute deviations for the pixels in various annuli around each site of interest on the 28 band difference ratios. We then use principle components analysis followed by linear discriminant analysis to train a classifier which assigns a posterior probability that a location is an archaeological site. We tested the procedure using leave-one-out cross validation with a second leave-one-out step to choose parameters on a 9,859x23,000 subset of the WorldView-2 data over the western portion of Ft. Irwin, CA, USA. We used 100 known non-sites and trained one classifier for lithic sites (n=33) and one classifier for habitation sites (n=16). We then analyzed convex combinations of scores from the Archaeological Predictive Model (APM) and our scores. We found that that the combined scores had a higher area under the ROC curve than either individual method, indicating that including WorldView-2 data in analysis improved the predictive power of the provided APM.

Tilton, James C.↗

Fragmentation analysis of a bar with the Lip-field approach

The Lip-field approach was introduced in Moës and Chevaugeon (2021) as a new way to regularize softening material models. It was tested in 1D quasistatic in Moës and Chevaugeon (2021) and 2D quasistatic in Chevaugeon and Moës (2021): this paper extends it to 1D dynamics, on the challenging problem of dynamic fragmentation. The Lip-field approach formulates the mechanical problem to be solved as an optimization problem, where the incremental potential to be minimized is the non-regularized one. Spurious localization is prevented by imposing a Lipschitz constraint on the damage field. Here, the displacement and damage field at each time step are obtained by a staggered algorithm, that is the displacement field is computed for a fixed damage field, then the damage field is computed for a fixed displacement field. Indeed, these two problems are convex, which is not the case of the global problem where the displacement and damage fields are sought at the same time. The incremental potential is obtained by equivalence with a cohesive zone model, which makes material parameters calibration simple. A non-regularized local damage equivalent to a cohesive zone model is also proposed. It is used as a reference for the Lip-field approach, without the need to implement displacement jumps. These approaches are applied to the brittle fragmentation of a 1D bar with randomly perturbed material properties to accelerate spatial convergence. Both explicit and implicit dynamic implementations are compared. Favorable comparison to several analytical, numerical and experimental references serves to validate the modeling approach.

36 MATERIALS SCIENCE↗

Moments-based interface reconstruction, remap and advection

Here, we present a new moment-of-fluid (MOF 2 ) interface reconstruction method. It uses the zeroth, first, and second moments of the fragment of material inside a cell of the mesh to reconstruct a convex material polygon or a union of convex polygons that approximate the respective material fragment. The new method requires information about the material moments only for the cell under consideration. The MOF 2 method allows to exactly reproduce several convex shapes: corners, filaments, and some concave shapes: cell-complements to corners and filaments. Interface reconstruction is formulated as a local (for each cell), non-linear, equality constrained optimization problem, which does not require additional communication and allows for an efficient parallel implementation. We present an extensive set of test problems, both for interface reconstruction on a single cell, and for reconstruction of a variety of shapes on a variety of meshes. We describe how to perform two-material advection using the MOF 2 method and present the results for the classical advection tests. We also show the examples of material interface remapping needed in the framework of multi-material arbitrary Lagrangian-Eulerian methods, and give a brief description of a procedure that can be used to update the material moments on the Lagrangian stage of those methods.

97 MATHEMATICS AND COMPUTING↗