Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 145 records · Page 8

Enhancing the cooling performance of thermocouples: a power-constrained topology optimization procedure

Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

Gutiérrez, G. Reales↗

McCormick envelopes in mixed-integer PDE-constrained optimization

McCormick envelopes are a standard tool for deriving convex relaxations of optimization problems that involve polynomial terms. Such McCormick relaxations provide lower bounds, for example, in branch-and-bound procedures for mixed-integer nonlinear programs but have not gained much attention in PDE-constrained optimization so far. This lack of attention may be due to the distributed nature of such problems, which on the one hand leads to infinitely many linear constraints (generally state constraints that may be difficult to handle) in addition to the state equation for a pointwise formulation of the McCormick envelopes and renders bound-tightening procedures that successively improve the resulting convex relaxations computationally intractable. We analyze McCormick envelopes for a model problem class that is governed by a semilinear PDE involving a bilinearity and integrality constraints. We approximate the nonlinearity and in turn the McCormick envelopes by averaging the involved terms over the cells of a partition of the computational domain on which the PDE is defined. This yields convex relaxations that underestimate the original problem up to an a priori error estimate that depends on the mesh size of the discretization. These approximate McCormick relaxations can be improved by means of an optimization-based bound-tightening procedure. We show that their minimizers converge to minimizers to a limit problem with a pointwise formulation of the McCormick envelopes when driving the mesh size to zero. We provide a computational example, for which we certify all of our imposed assumptions. The results point to both the potential of the methodology and the gaps in the research that need to be closed. Our methodology provides a framework first for obtaining pointwise underestimators for nonconvexities and second for approximating them with finitely many linear inequalities in an infinite-dimensional setting.

Approximations and Expansions↗

Convergence analysis for a nonlocal gradient descent method via directional Gaussian smoothing

We analyze the convergence of a nonlocal gradient descent method for minimizing a class of high-dimensional non-convex functions, where a directional Gaussian smoothing (DGS) is proposed to define the nonlocal gradient (also referred to as the DGS gradient). The method was first proposed in [Zhang et al., Enabling long-range exploration in minimization of multimodal functions, UAI 2021], in which multiple numerical experiments showed that replacing the traditional local gradient with the DGS gradient can help the optimizers escape local minima more easily and significantly improve their performance. However, a rigorous theory for the efficiency of the method on nonconvex landscape is lacking. In this work, we investigate the scenario where the objective function is composed of a convex function, perturbed by deterministic oscillating noise. We provide a convergence theory under which the iterates exponentially converge to a tightened neighborhood of the solution, whose size is characterized by the noise wavelength. Here, we also establish a correlation between the optimal values of the Gaussian smoothing radius and the noise wavelength, thus justifying the advantage of using moderate or large smoothing radii with the method. Furthermore, if the noise level decays to zero when approaching the global minimum, we prove that DGS-based optimization converges to the exact global minimum with linear rates, similarly to standard gradient-based methods in optimizing convex functions. Several numerical experiments are provided to confirm our theory and illustrate the superiority of the approach over those based on the local gradient.

Tran, Hoang [Oak Ridge National Laboratory (ORNL),↗

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↗

Investigating the origin of the far-field reflection interference fringe (RIF) of microdroplets

We show that the reflection interference fringe (RIF) is formed on a screen far away from the microdroplets placed on a prism-based substrate, which have low contact angles and thin droplet heights, caused by the dual convex–concave profile of the droplet, not a pure convex profile. The geometric formulation shows that the interference fringes are caused by the optical path difference when the reflected rays from the upper convex profile at the droplet–air interface interfere with reflection from the lower concave profile at oblique angles lower than the critical angle. Analytic solutions are obtained for the droplet height and the contact angle out of the fringe number and the fringe radius in RIF from the geometric formulation. Furthermore, the ray tracing simulation is conducted using the custom-designed code. The geometric formulation and the ray tracing show excellent agreement with the experimental observation in the relation between the droplet height and the fringe number and the relation between the contact angle and the fringe radius. This study is remarkable as the droplet's dual profile cannot be easily observed with the existing techniques. However, the RIF technique can effectively verify the existence of a dual profile of the microdroplets in a simple setup. In this work, the RIF technique is successfully developed as a new optical diagnostic technique to determine the microdroplet features, such as the dual profile, the height, the contact angle, the inflection point, and the precursor film thickness, by simply measuring the RIF patterns on the far-field screen.

42 ENGINEERING↗

Data-Driven Compositional Optimization in Misspecified Regimes

With a manifold growth in the scale and intricacy of systems, the challenges of parametric misspecification become pronounced. These concerns are further exacerbated in compositional settings, which emerge in problems complicated by modeling risk and robustness. In “Data-Driven Compositional Optimization in Misspecified Regimes,” the authors consider the resolution of compositional stochastic optimization problems, plagued by parametric misspecification. In considering settings where such misspecification may be resolved via a parallel learning process, the authors develop schemes that can contend with diverse forms of risk, dynamics, and nonconvexity. They provide asymptotic and rate guarantees for unaccelerated and accelerated schemes for convex, strongly convex, and nonconvex problems in a two-level regime with extensions to the multilevel setting. Surprisingly, the nonasymptotic rate guarantees show no degradation from the rate statements obtained in a correctly specified regime and the schemes achieve optimal (or near-optimal) sample complexities for general T-level strongly convex and nonconvex compositional problems.

Business & Economics↗

Precision Computations in Strongly Coupled Conformal Field Theories (Final Technical Report)

Conformal Field Theories (CFTs) are quantum field theories that are invariant under the conformal symmetry group (which includes translations and rotations, but also local rescalings of spacetime). They are building blocks of general quantum field theories, and appear in many areas of physics, including statistical physics, condensed matter physics, particle physics, and quantum gravity. Because of their extra symmetries, the mathematical structure of CFTs is tightly constrained, and this leads to the idea of the ``conformal bootstrap," which is to use these mathematical structures to constrain, and in some cases determine, CFT observables. A new numerical implementation of the conformal bootstrap idea appeared in 2008 with the work of Rattazzi, Rychkov, Tonni, and Vichi. Their observation was that certain bootstrap constraints (conformal symmetry and unitarity) could be combined to yield a convex optimization problem that constraints CFT data. By solving this convex optimization problem on a computer, one could obtain bounds on observables like critical exponents and operator product expansion (OPE) coefficients. Over the course of this award, the PI has improved numerical bootstrap techniques by optimizing known algorithms and finding new ones for performing the required convex optimization computations. The PI has applied these techniques to compute high-precision observables in several important strongly-coupled systems. The PI has also explored both analytical and numerical bootstrap methods for constraining the space of low energy effective field theories of quantum gravity, and developed new analytical techniques for CFT and QFT more broadly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An experimental investigation of turbulent boundary layers along curved surfaces

A curved wall tunnel was designed, and an equilibrium turbulent boundary layer was set up on the straight section preceding the curved test section. Turbulent boundary layer flows with uniform and adverse pressure distributions along convex and concave walls were investigated. Hot-wire measurements along the convex surface indicated that turbulent mixing between fluid layers was very much reduced. However, the law of the wall held and the skin friction, thus determined, correlated well with other measurements. Hot-wire measurements along the concave test wall revealed a system of longitudinal vortices inside the boundary layer and confirmed that concave curvature enhances mixing. A self-consistent set of turbulent boundary layer equations for flows along curved surfaces was derived together with a modified eddy viscosity. Solution of these equations together with the modified eddy viscosity gave results that correlated well with the present data on flows along the convex surface with arbitrary pressure distribution. However, it could only be used to predict the mean characteristics of the flow along concave walls because of the existence of the system of longitudinal vortices inside the boundary layer.

So, R. M. C.↗

A Natural Bias Approach to Cardinal Spline Curves

The cardinal spline approach to defining interpolatory curves is of interest in the application of computer graphics to modelling and animation problems. The formulation of this interpolation problem (C) a vector valued function V is obtained which forms a C1 map from 0,n into the plane in such a way that V(i) and V prime (i) attain prespecified (vector) values P sub i and T sub i for the integers i=0, 1,...n. The intent is that the planar curves V(s) pass through each point P sub i in the direction corresponding to T sub i. In actual practice it is common for only the interpolation points P sub i to be user specified. Suitable tangents T sub i are then provided by some heuristic to complete the specification of (C). A convex control polygon yields a convex interpolant. Examples exist where the interpolant fails to be convex if the tangents are left unbiased.

Fletcher, G. Y.↗

Mean flow in turbulent boundary layers disturbed to alter skin friction

A comparative evaluation is made of recent developments in methods for the reduction of boundary layer drag, encompassing longitudinal surface riblets, 'outer layer' devices, (OLDs) and longitudinal convex surface curvature. The boundary layer of a surface with a longitudinal concave curvature is also studied, to complement the convex case results. The net drag reductions achievable by both riblets and OLDs are noted to be a rather modest 10 percent. Boundary layers exhibit asymmetric response to streamwise surface curvatures, with the response being slower for the case of a concave than a convex curvature.

Bandyopadhyay, P. R.↗

Space ultra-vacuum facility and method of operation

A wake shield facility providing an ultrahigh vacuum level for space processing is described. The facility is in the shape of a truncated, hollow hemispherical section, one side of the shield convex and the other concave. The shield surface is preferably made of material that has low out-gassing characteristics such as stainless steel. A material sample supporting fixture in the form of a carousel is disposed on the convex side of the shield at its apex. Movable arms, also on the convex side, are connected by the shield in proximity to the carousel, the arms supporting processing fixtures, and providing for movement of the fixtures to predetermined locations required for producing interations with material samples. For MBE processes a vapor jet projects a stream of vaporized material onto a sample surface. The fixtures are oriented to face the surface of the sample being processed when in their extended position, and when not in use they are retractable to a storage position. The concave side of the shield has a support structure including metal struts connected to the shield, extending radially inward. The struts are joined to an end plate disposed parallel to the outer edge of the shield. This system eliminates outgassing contamination.

Naumann, Robert J.↗

Asteroid lightcurve inversion

One of the most fundamental physical properties of any asteroid is its shape. Lightcurves provide the only source of shape information for most asteroids. Unfortunately, the functional form of a lightcurve is determined by the viewing/illumination geometry and the asteroid's light scattering characteristics as well as its shape, and in general it is impossible to determine an asteroid's shape from lightcurves. A technique called convex-profile inversion (CPI) that obtains a convex profile, P, from any lightcurve is introduced. If certain ideal conditions are satisfied, then P is an estimator for the asteroid's mean cross section, C, a convex set defined as the average of all cross sections C(z) cut by planes a distance z above the asteroids's equatorial plane. C is therefore a 2-D average of the asteroid's 3-D shape.

Ostro, Steven J.↗

Curvature effects on the stability of three-dimensional laminar boundary layers

The linear stability equations for compressible, three-dimensional laminar boundary layer flow are derived in an orthogonal curvilinear coordinate system. The system of equations is solved using a finite difference scheme in order to study the effects of streamline and surface curvature and compressibility on the stability of the flow past a swept wing. It is known that convex surface curvature can have a stabilizing effect on the laminar boundary layer. Conversely, concave surface curvature can be destabilizing. The magnitude of these effects for swept wings is determined. Results indicate that amplification rates and hence, N-factors, for the flow over the convex upper surface of a swept wing can be reduced by about 15 to 45 percent when curvature effects are included in the linear stability analysis. The results of the calculations show that concave curvature destabilizes crossflow type disturbances with a significant increase in amplification rate. In addition, comparisons are made with some experimental results on a swept concave-convex surface. Calculated velocity vector plots show good agreement with observed disturbances in the laminar boundary layer over the concave surface.

Collier, F. S., Jr.↗

The distributed virtual windtunnel

An implementation of a distributed virtual environment for the shared interactive visualization of large unsteady three-dimensional flowfields is described. Computation of the visualizations is performed on a Convex C3240 computer, and the visualization data are transferred over a high-speed network to a Silicon Graphics Iris workstation for rendering. A boom-mounted six-degree-of-freedom head-position-sensitive stereo CRT (cathode ray tube) system is used for display. A hand-position-sensitive glove controller is used for controlling various tracers for the visualization of the flow. User commands are sent to the Convex, which interprets these commands and computes the corresponding visualization. With this architecture, several users may share and cooperatively control the visualization generated by the Convex. The distributed architecture is also interesting to those using conventional screen and mouse interfaces.

Bryson, Steve↗

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Elishakoff, Isaac↗

Undulated Nozzle for Enhanced Exit Area Mixing

A nozzle having an undulating surface for enhancing the mixing of a primary flow with a secondary flow or ambient air, without requiring an ejector. The nozzle includes a nozzle structure and design for introducing counter-rotating vorticity into the primary flow either through (i) internal surface corrugations where an axisymmetric line through each corrugation is coincident with an axisymmetric line through the center of the flow passageway or (ii) through one or more sets of alternating convexities and cavities in the internal surface of the nozzle where an axisymmetric line through each convexity and cavity is coincident with an axisymmetric line through the center of the flow passageway, and where the convexities contract from the entrance end towards the exit end. Exit area mixing is also enhanced by one or more chevrons attached to the exit edge of the nozzle. The nozzle is ideally suited for application as a jet engine nozzle. When used as a jet engine nozzle, noise suppression with simultaneous thrust augmentation/minimal thrust loss is achieved.

Seiner, John M.↗

Computer-based video digitizer analysis of surface extension in maize roots: kinetics of growth rate changes during gravitropism

We used a video digitizer system to measure surface extension and curvature in gravistimulated primary roots of maize (Zea mays L.). Downward curvature began about 25 +/- 7 min after gravistimulation and resulted from a combination of enhanced growth along the upper surface and reduced growth along the lower surface relative to growth in vertically oriented controls. The roots curved at a rate of 1.4 +/- 0.5 degrees min-1 but the pattern of curvature varied somewhat. In about 35% of the samples the roots curved steadily downward and the rate of curvature slowed as the root neared 90 degrees. A final angle of about 90 degrees was reached 110 +/- 35 min after the start of gravistimulation. In about 65% of the samples there was a period of backward curvature (partial reversal of curvature) during the response. In some cases (about 15% of those showing a period of reverse bending) this period of backward curvature occurred before the root reached 90 degrees. Following transient backward curvature, downward curvature resumed and the root approached a final angle of about 90 degrees. In about 65% of the roots showing a period of reverse curvature, the roots curved steadily past the vertical, reaching maximum curvature about 205 +/- 65 min after gravistimulation. The direction of curvature then reversed back toward the vertical. After one or two oscillations about the vertical the roots obtained a vertical orientation and the distribution of growth within the root tip became the same as that prior to gravistimulation. The period of transient backward curvature coincided with and was evidently caused by enhancement of growth along the concave and inhibition of growth along the convex side of the curve, a pattern opposite to that prevailing in the earlier stages of downward curvature. There were periods during the gravitropic response when the normally unimodal growth-rate distribution within the elongation zone became bimodal with two peaks of rapid elongation separated by a region of reduced elongation rate. This occurred at different times on the convex and concave sides of the graviresponding root. During the period of steady downward curvature the elongation zone along the convex side extended farther toward the tip than in the vertical control. During the period of reduced rate of curvature, the zone of elongation extended farther toward the tip along the concave side of the root. The data show that the gravitropic response pattern varies with time and involves changes in localized elongation rates as well as changes in the length and position of the elongation zone. Models of root gravitropic curvature based on simple unimodal inhibition of growth along the lower side cannot account for these complex growth patterns.

NASA Discipline Plant Biology↗

Interface Shape Control using Localized Heating during Bridgman Growth

Numerical calculations were performed to assess the effect of localized radial heating on the melt-crystal interface shape during vertical Bridgman growth. System parameters examined include the ampoule, melt and crystal thermal conductivities, the magnitude and width of localized heating, and the latent heat of crystallization. Concave interface shapes, typical of semiconductor systems, could be flattened or made convex with localized heating. Although localized heating caused shallower thermal gradients ahead of the interface, the magnitude of the localized heating required for convexity was less than that which resulted in a thermal inversion ahead ofthe interface. A convex interface shape was most readily achieved with ampoules of lower thermal conductivity. The conditions under which convection in the melt must be considered were determined.

Volz, M. P.↗