Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convex functions”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Adaptive Power Flow Approximations With Second-Order Sensitivity Insights

The power flow equations are fundamental to power system planning, analysis, and control. However, the inherent non-linearity and non-convexity of these equations present formidable obstacles in problem-solving processes. To mitigate these challenges, recent research has proposed adaptive power flow linearizations that aim to achieve accuracy over wide operating ranges. The accuracy of these approximations inherently depends on the curvature of the power flow equations within these ranges, which necessitates considering second-order sensitivities. In this paper, we leverage second-order sensitivities to both analyze and improve power flow approximations. We evaluate the curvature across broad operational ranges and subsequently utilize this information to inform the computation of various sample-based power flow approximation techniques. Additionally, we leverage second-order sensitivities to guide the development of rational approximations that yield linear constraints in optimization problems. In conclusion, this approach is extended to enhance accuracy beyond the limitations of linear functions across varied operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Systematic analysis of melt pool dynamics in laser processing of mixed powder feedstocks

Functionally graded materials (FGMs) fabricated via additive manufacturing of blended powders offer the potential to spatially tailor properties for new technologies, such as fusion first-wall systems, turbine blades, and spacecraft. However, processing these materials is difficult due to the multiplicity of processing parameters to optimize, all of which must be changed as substrate material, powder feedstock compositions, and melt pool dynamics evolve. Here, this work systematically evaluates the qualitative and quantitative effects of these variables on the melt pool size, shape, composition, and particle distribution in an exemplar Ti-Ta system, and connects the experimental results to Marangoni flow behavior and phenomena observed in other systems. Increasing laser power linearly increases melt pool size and layer thickness, driving engineering considerations such as part/geometrical tolerances. Decreasing laser velocity changes the melt pool shape from lenticular to convex and reduces chemical homogeneity due to extreme thermal and compositional gradients between the melt pool center and boundaries. Thermophysical property differences between the powder feedstock and substrate material, as well as the directionality of the gradient, affect dilution and melt pool dynamics, which in turn affect the melt pool boundary characteristics, shape, and uniformity. Mixed powder feedstocks of intermediate compositions do not behave according to linear interpolations between single-material endpoints, instead building taller and wider melt pools. As such, it is recommended to quantify process maps for at least one intermediate composition in the FGM or multi-material system of interest to ensure optimized processing parameters, predictable melt pool sizes and shapes, and compositional and spatial precision.

Dissimilar↗
NASA NTRS1986

Taylor-Goertler instabilities of Tollmien-Schlichting waves and other flows governed by the interactive boundary-layer equations

The Taylor-Goertler vortex instability equations are formulated for steady and unsteady interacting boundary-layer flows. The effective Goertler number is shown to be a function of the wall shape in the boundary layer and the possibility of both steady and unsteady Taylor-Goertler modes exists. As an example the steady flow in a symmetrically constricted channel is considered and it is shown that unstable Goertler vortices exist before the boundary layers at the wall develop the Goldstein singularity discussed by Smith and Daniels (1981). As an example of an unsteady spatially varying basic state, it is considered the instability of high-frequency large-amplitude two- and three-dimensional Tollmien-Schlichting waves in a curved channel. It is shown that they are unstable in the first 'Stokes-layer stage' of the hierarchy of nonlinear states discussed by Smith and Burggraf (1985). This instability of Tollmien-Schlichting waves in an internal flow can occur in the presence of either convex or concave curvature. Some discussion of this instability in external flows is given.

Hall, Philip↗
NASA NTRS2008

Re-entry vehicle shape for enhanced performance

A convex shell structure for enhanced aerodynamic performance and/or reduced heat transfer requirements for a space vehicle that re-enters an atmosphere. The structure has a fore-body, an aft-body, a longitudinal axis and a transverse cross sectional shape, projected on a plane containing the longitudinal axis, that includes: first and second linear segments, smoothly joined at a first end of each the first and second linear segments to an end of a third linear segment by respective first and second curvilinear segments; and a fourth linear segment, joined to a second end of each of the first and second segments by curvilinear segments, including first and second ellipses having unequal ellipse parameters. The cross sectional shape is non-symmetric about the longitudinal axis. The fourth linear segment can be replaced by a sum of one or more polynomials, trigonometric functions or other functions satisfying certain constraints.

Brown, James L.↗
NASA NTRS1974

An analysis of the radiation from apertures in curved surfaces by the geometrical theory of diffraction

In this paper the geometrical theory of diffraction is extended to treat the radiation from apertures of slots in convex perfectly conducting surfaces. It is assumed that the tangential electric field in the aperture is known so that an equivalent infinitesimal source can be defined at each point in the aperture. Surface rays emanate from this source which is a caustic of the ray system. A launching coefficient is introduced to describe the excitation of the surface ray modes. If the field radiated from the surface is desired, the ordinary diffraction coefficients are used to determine the field of the rays shed tangentially from the surface rays. The field of the surface ray modes is not the field on the surface; hence if the mutual coupling between slots is of interest, a second coefficient related to the launching coefficient must be employed. In the region adjacent to the shadow boundary, the component of the field directly radiated from the source is represented by Fock-type functions. In the illuminated region the incident radiation from the source (this does not include the diffracted field components) is treated by geometrical optics. This extension of the geometrical theory of diffraction is applied to calculate the radiation from slots on elliptic cylinders, spheres, and spheroids.

Pathak, P. H.↗
NASA NTRS1973

The radiation from apertures in curved surfaces

The geometrical theory of diffraction is extended to treat the radiation from apertures or slots in convex, perfectly-conducting surfaces. It is assumed that the tangential electric field in the aperture is known so that an equivalent, infinitesimal source can be defined at each point in the aperture. Surface rays emanate from this source which is a caustic of the ray system. A launching coefficient is introduced to describe the excitation of the surface ray modes. If the field radiated from the surface is desired, the ordinary diffraction coefficients are used to determine the field of the rays shed tangentially from the surface rays. The field of the surface ray modes is not the field on the surface; hence if the mutual coupling between slots is of interest, a second coefficient related to the launching coefficient must be employed. In the region adjacent to the shadow boundary, the component of the field directly radiated from the source is presented by Fock-type functions. In the illuminated region the incident radiation from the source (this does not include the diffracted field components) is treated by geometrical optics. This extension of the geometrical theory of diffraction is applied to calculate the radiation from slots on elliptic cylinders, spheres and spheroids.

Pathak, P. H.↗
DOE OSTI2024

Machine learning assisted search for Fe–Co–C ternary compounds with high magnetic anisotropy

We employ a machine learning (ML)-guided framework to explore rare earth free magnetic materials, specifically focusing on Fe–Co–C ternary compounds for potential use in permanent magnets. Utilizing a specifically trained crystal graph convolutional neural network model, we efficiently screen a vast space of nearly a million substitutional structures to select 620 promising structures for further investigation by first-principles calculation. We predict five low-energy metastable Fe–Co–C compounds with formation energy less than 150 meV/atom above the convex hull. These compounds exhibit high magnetization (Js > 1.0 T) and significant magnetic anisotropy (K1 > 1.0 MJ/m3), making them promising candidates for permanent magnet applications. The phonon calculations indicate these compounds are dynamically stable. Our ML-guided framework demonstrates the utility of rapidly identifying novel materials with tailored magnetic properties.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Organization of P, S, and Fe Inclusions in a Freshwater Magnetococcus

Magnetotactic bacteria are a heterogeneous group of motile, mainly aquatic procaryotes that align and swim along geomagnetic field lines. They are of interest to astrobiologists because of the magnetite crystals found in the Mars meteorite ALH84001 which share many characteristics with the magnetite produced intracellularly by magnetotactic bacteria. These bacteria are diverse morphologically, physiologically and phylogenetically, sharing a few key characteristics: 1) the presence of intracellular membrane-bound magnetic crystals (magnetosomes), usually but not always arranged in chains; 2) motility by means of flagella; and 3) microaerophillic or anaerobic physiology. The bilophotrichous (having two flagella bundles) magnetotactic cocci (MC) are ubiquitous in aquatic habitats but have proven extremely difficult to cultivate. Because only several marine strains have been isolated and grown in axenic culture, little is known about the physiology and the biogeochemical roles of the MC. We studied the composition and distribution of intracellular structures in an uncultured MC, designated ARB-1. To do this, a combination of light microscopy, environmental scanning electron microscopy (ESEM), transmission electron microscopy (TEM) and energy dispersive X-ray spectroscopy (EDS) were used. Cells of ARB-1 were separated from sediments collected from Baldwin Lake (Los Angeles Arboretum, Arcadia, CA). They are large spherical to oblate spheroidal Gram-negative cells, ranging from 1 to 4 micrometers along the maximum dimension, which is perpendicular to the direction of swimming. Cells have two large phosphorus-containing inclusions that comprise a large percentage of the cell volume. Many smaller sulfur inclusions are located at the convex end of the cell. Most of the cellular Fe is present in the magnetosomes. These may be arranged as a clump at the concave end of the cell, near the two flagella bundles, or as chains, or as both a clump and chains. The magnetosomes were identified as magnetite (Fe3O4) by selected area electron diffraction (SAED) and high resolution TEM. We saw a trend between cell size and organization of the magnetosomes. Smaller, more spherical cells were more likely to have chains than were larger, more oblate cells. This may indicate different populations of cells, or it may be attributed to variations in cell growth cycle. The size distribution (length) of magnetosomes in chains was similar to that of magnetosomes in clusters, except that there was a larger size range for clustered magnetosomes. Magnetosomes from ARB-1 cells average 82 nm in length. If plotted on a graph of length as a function of aspect ratio, they fall within the single domain region of the plot. If compared with the size distributions of magnetite from ALH84001 and magnetosomes from the cultured magnetotactic vibrio MV-1, the magnetites produced by ARB-1 cells are, on average, larger and have a wider range of aspect ratio. ARB-1 cells have a specific organization of the P, S, and Fe inclusions. The P inclusions always occupy the majority of the cell volume and separate the S inclusions from the disorganized clumps of magnetosomes and the flagella bundles. The P inclusions may contain polyphosphate, which could play several roles in motility, adaptation to stress, growth and division, buoyancy, and energy. The S inclusions might be a way to store S, a potential energy source, when the cells move from sulfide to oxygen zones. The consistency of P, S, and Fe organization in ARB-1 cells suggests that these inclusions have some specific and interactive functions.

Cox, Lea↗
DOE OSTI2022

Cross Inference of Throughput Profiles Using Micro Kernel Network Method

Dedicated network connections are being increasingly deployed in cloud, centralized and edge computing and data infrastructures, whose throughput profiles are critical indicators of the underlying data transfer performance. Due to the cost and disruptions to physical infrastructures, network emulators, such as Mininet, are often used to generate measurements needed to estimate throughput profiles, typically expressed as a function of the connection round trip time. The profiles estimated using measurements from such emulated networks are usually inaccurate for high bandwidth and high latency connections, since they do not accurately reflect the critical network transport dynamics mainly due to computing and memory constraints of the host. We present a machine learning (ML) method to estimate the throughput profiles using emulation measurements to closely match the testbed and production network profiles. In particular, we propose a micro Kernel Network (mKN) that provides baseline throughput measurements on the host running Mininet emulations, which are used to learn a regression map that converts them to the corresponding testbed measurement estimates. Once initially learned, this map is applied to measurements from subsequent network emulations on the same host. We present experimental measurements to illustrate this approach, and derive generalization equations for the proposed mKN-ML method. Using a four-site scenario emulation, we show the effectiveness of this method in providing accurate concave throughput profiles from inaccurate convex or non-smooth ones indicated by Mininet emulation.

Rao, Nageswara↗
NASA NTRS1999

Practical Aspects of Stabilized FEM Discretizations of Nonlinear Conservation Law Systems with Convex Extension

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the Galerkin least-squares (GLS) and the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the POE system. Central to the development of the simplified GLS and DG methods is the Degenerative Scaling Theorem which characterizes right symmetrizes of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobean matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler, Navier-Stokes, and magnetohydrodynamic (MHD) equations. Linear and nonlinear energy stability is proven for the simplified GLS and DG methods. Spatial convergence properties of the simplified GLS and DO methods are numerical evaluated via the computation of Ringleb flow on a sequence of successively refined triangulations. Finally, we consider a posteriori error estimates for the GLS and DG demoralization assuming error functionals related to the integrated lift and drag of a body. Sample calculations in 20 are shown to validate the theory and implementation.

Barth, Timothy↗
DOE OSTI2023

A Novel Multi-Agent Deep Reinforcement Learning-enabled Distributed Power Allocation Scheme for mmWave Cellular Networks

We consider the power allocation problem over shared spectrum for millimeter-Wave (mmWave) cellular downlink. Existing approaches usually find sub-optimal solutions by solving a non-convex optimization which leads to scalability issues due to centralized control. Therefore, distributed and adaptive approaches are desirable. Recently, model-free Deep Reinforcement Learning (DRL) has achieved success in such wireless resource management tasks. By modeling the radio environment as a Markov Decision Process (MDP) with the base stations (BSs) being the agents, power allocation can be automated at the agent level with comparable throughput performance to conventional centralized schemes. The multi-agent setting presents new challenges as the radio environment is impacted by the joint actions of the agents and is no longer stationary from any individual agent’s perspective. Existing literature bypasses this non-stationarity violation by ignoring it which may cause performance degradation. To tackle this issue, we propose a distributed continuous power allocation scheme based on a modified version of multi-agent Deep Deterministic Policy Gradient (MADDPG) that is tailored for the distributed multiple-agent setting. The proposed scheme employs a centralized-training distributed-execution framework where Q-functions are trained over subsets of BSs while each BS determines its transmit power based only on its own local observation. It admits constant per-BS communication and computation complexity and is thus scalable to large networks. Numerical evaluation shows that the proposed scheme adapts well to a wide range of interference conditions and can achieve comparable or better performance than several state-of-the-art non-learning approaches.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Candidate ferroelectrics via ab initio high-throughput screening of polar materials

Ferroelectrics are a class of polar and switchable functional materials with diverse applications, from microelectronics to energy conversion. Computational searches for new ferroelectric materials have been constrained by accurate prediction of the polarization and switchability with electric field, properties that, in principle, require a comparison with a nonpolar phase whose atomic-scale unit cell is continuously deformable from the polar ground state. For most polar materials, such a higher-symmetry nonpolar phase does not exist or is unknown. Here, we introduce a general high-throughput workflow that screens polar materials as potential ferroelectrics. We demonstrate our workflow on 1978 polar structures in the Materials Project database, for which we automatically generate a nonpolar reference structure using pseudosymmetries, and then compute the polarization difference and energy barrier between polar and nonpolar phases, comparing the predicted values to known ferroelectrics. Focusing on a subset of 182 potential ferroelectrics, we implement a systematic ranking strategy that prioritizes candidates with large polarization and small polar-nonpolar energy differences. To assess stability and synthesizability, we combine information including the computed formation energy above the convex hull, the Inorganic Crystal Structure Database id number, a previously reported machine learning-based synthesizability score, and ab initio phonon band structures. To distinguish between previously reported ferroelectrics, materials known for alternative applications, and lesser-known materials, we combine this ranking with a survey of the existing literature on these candidates through Google Scholar and Scopus databases, revealing ~130 promising materials uninvestigated as ferroelectric. Our workflow and large-scale high-throughput screening lays the groundwork for the discovery of novel ferroelectrics, revealing numerous candidates materials for future experimental and theoretical endeavors.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Statistics of base polytopes in F-theory

We propose a new statistical ensemble of toric bases for elliptic Calabi-Yaus used in F-theory models, by focusing on only the convex hull of the base, i.e., the base polytope. This physically motivated coarse-graining greatly simplifies the combinatorial complexity of the part of the 4d F-theory landscape with toric bases. We develop a Monte Carlo approach that randomly samples the base polytopes within fixed boxes, with proper statistical weights. We first apply the algorithm to the set of 2d base polytopes, generating an enlarged set of toric 2d bases that include certain types of codimension-two (4,6) points, and we validate our approach against exact numbers. We then explore the set of 3d base polytopes which fit in a set of “maximal” 3d boxes, and estimate the total number of inequivalent 3d base polytopes to be 10 85 –10 90 . We provide statistical data such as the distribution of non-Higgsable gauge groups on these bases. Amusingly, a similar method can also be applied to generate reflexive polytopes in various dimensions. In both the reflexive and base polytope cases, the number of relevant polytopes obeys a Gaussian distribution as a function of the number of vertices, which can be understood in terms of other results on random polytopes in the math literature.

Differential and algebraic geometry↗
NASA NTRS1999

Edge-Flames in Von Karman Swirling Flows

Classical understanding of diffusion flames dictates that they, unlike the premixed flames, do not possess a characteristic propagation velocity and are constrained by stoichiometric requirements at the flame surface. However, it has been commonly observed that when local extinction occurs within a diffusion flame sheet, the edges that are formed propagate with distinct speeds. In general, the propagation speed of these edges depend on their geometrical shape (concave, convex, or straight) among other factors. Recently, Buckmaster investigated the dynamics of straight diffusion flame edges separating burning and quenched regions using simplified one-dimensional models. He showed that these flame edges can have positive, negative, or zero velocity depending on the Damkoehler number of the equilibrium diffusion flame that support them. It was also shown that this unsteady flame-edge behavior is intrinsically linked to S-curve behavior of the diffusion flame with varying Damkoehler number. When the system Damkoehler number lies between the extinction and ignition limits, flame edges can propagate as an "ignition wave" or as a "failure wave," and for a critical Damkoehler number remain as a stationary flame-edge. We have extend Buckmaster's 1-d model to more general edge-flame configurations where the edges appear as "flame holes" or as "flame disks". These two configurations along with the straight-edge case cover the entire range of possible edge-flame geometry observable in planar diffusion-flame sheets. A generalized map of edge-flame propagation velocities as a function of the system Damkoehler number and the edge-flame radius is presented. Experimentally we show that edge flames can be created using diffusion flames embedded in von Karman boundary layers. In a von Karman boundary layer, the flow is generated by spinning a solid (fuel) disk in a quiescent ambient gas. Under normal gravity we were able to produce "flame disks" over a range of fuel-disk rotational velocities varying from 0 to 20 revolutions per second, by orienting the burning surface of the fuel disk facing downward.

Nayagam, Vedha↗
NASA NTRS1993

Techniques for estimating Space Station aerodynamic characteristics

A method was devised and calculations were performed to determine the effects of reflected molecules on the aerodynamic force and moment coefficients for a body in free molecule flow. A procedure was developed for determining the velocity and temperature distributions of molecules reflected from a surface of arbitrary momentum and energy accommodation. A system of equations, based on momentum and energy balances for the surface, incident, and reflected molecules, was solved by a numerical optimization technique. The minimization of a 'cost' function, developed from the set of equations, resulted in the determination of the defining properties of the flow reflected from the arbitrary surface. The properties used to define both the incident and reflected flows were: average temperature of the molecules in the flow, angle of the flow with respect to a vector normal to the surface, and the molecular speed ratio. The properties of the reflected flow were used to calculate the contribution of multiply reflected molecules to the force and moments on a test body in the flow. The test configuration consisted of two flat plates joined along one edge at a right angle to each other. When force and moment coefficients of this 90 deg concave wedge were compared to results that did not include multiple reflections, it was found that multiple reflections could nearly double lift and drag coefficients, with nearly a 50 percent increase in pitching moment for cases with specular or nearly specular accommodation. The cases of diffuse or nearly diffuse accommodation often had minor reductions in axial and normal forces when multiple reflections were included. There were several cases of intermediate accommodation where the addition of multiple reflection effects more than tripled the lift coefficient over the convex technique.

Thomas, Richard E.↗
NASA NTRS1993

Large eddy simulation of a boundary layer with concave streamwise curvature

One of the most exciting recent developments in the field of large eddy simulation (LES) is the dynamic subgrid-scale model. The dynamic model concept is a general procedure for evaluating model constants by sampling a band of the smallest scales actually resolved in the simulation. To date, the procedure has been used primarily in conjunction with the Smagorinsky model. The dynamic procedure has the advantage that the value of the model constant need not be specified a priori, but rather is calculated as a function of space and time as the simulation progresses. This feature makes the dynamic model especially attractive for flows in complex geometries where it is difficult or impossible to calibrate model constants. The dynamic model was highly successful in benchmark tests involving homogeneous and channel flows. Having demonstrated the potential of the dynamic model in these simple flows, the overall direction of the LES effort at CTR shifted toward an evaluation of the model in more complex situations. The current test cases are basic engineering-type flows for which Reynolds averaged approaches were unable to model the turbulence to within engineering accuracy. Flows currently under investigation include a backward-facing step, wake behind a circular cylinder, airfoil at high angles of attack, separated flow in a diffuser, and boundary layer over a concave surface. Preliminary results from the backward-facing step and cylinder wake simulations are encouraging. Progress on the LES of a boundary layer on a concave surface is discussed. Although the geometry of a concave wall is not very complex, the boundary layer that develops on its surface is difficult to model due to the presence of streamwise Taylor-Gortler vortices. These vortices arise as a result of a centrifugal instability associated with the convex curvature.

Lund, Thomas S.↗
DOE OSTI2022

Scaling and Benchmarking an Evolutionary Algorithm for Constructing Biophysical Neuronal Models

Single neuron models are fundamental for computational modeling of the brain's neuronal networks, and understanding how ion channel dynamics mediate neural function. A challenge in defining such models is determining biophysically realistic channel distributions. Here, we present an efficient, highly parallel evolutionary algorithm for developing such models, named NeuroGPU-EA. NeuroGPU-EA uses CPUs and GPUs concurrently to simulate and evaluate neuron membrane potentials with respect to multiple stimuli. We demonstrate a logarithmic cost for scaling the stimuli used in the fitting procedure. NeuroGPU-EA outperforms the typically used CPU based evolutionary algorithm by a factor of 10 on a series of scaling benchmarks. We report observed performance bottlenecks and propose mitigation strategies. Finally, we also discuss the potential of this method for efficient simulation and evaluation of electrophysiological waveforms.

59 BASIC BIOLOGICAL SCIENCES↗