Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Differential Geometry”

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 379 records · Page 21

Modeling of SiC Lateral Resonant Devices Over a Broad Temperature Range

Finite-element analysis (FEA) modal results of 3C-SiC lateral resonant devices anchored to a Si substrate are presented as resonant frequency versus temperature. The suspended elements are etched from a 2 micron, 3C-SiC film grown at 1600 K on a 500 micron-thick, Si substrate. The analysis includes, temperature-dependent properties, shape change due to volume expansion with temperature, and thermal stress caused by differential thermal expansion of different materials. Two designs are considered: type I has anchor locations close to the geometric centroid and a small shuttle; type 11 has a large shuttle with anchors far from the centroid, The resonant frequency decreases approximately 3.5% over a 1000 K temperature increase for the type-I device, and behaves according to theory. The resonant frequency of the type-11 device decreases by 2% over the first 400 K, then rises slightly over the remaining 600 K. This device deviates from theory because of the high thermal stress induced in the beams. The thermal stress is caused by the differential thermal expansion of the suspended element relative to the substrate. The results show that the device geometry must be properly chosen if the resonant frequency of that device will be used to calculate the temperature coefficient of Young's modulus. These results apply only to resonators of one material on a substrate of a different material.

DeAnna, Russell G.↗

Development of a Convolutional Neural Network Classifier for Data Starved Spectra - 20199

The Institute for Clean Energy Technology (ICET) at Mississippi State University is exploring the utility of machine learning in augmenting its mobile radiation surveying platforms, which are currently being developed as means to survey depleted uranium contaminated areas in support of remediation and decommissioning efforts. Mobile survey platforms provide a means to efficiently scan large areas of interest while reducing human exposure to radiation and other hazards. The survey platforms can also be used for scanning for any gamma emitting isotope in addition to depleted uranium. The spectral data that the platforms collect may be data starved with relatively low counts and poorly defined spectral features depending on the speed of the platforms and scintillation detector selection. Such data-starved spectra are difficult to use for isotope identification, requiring advanced knowledge of the possible radionuclides that could be present and environmental factors that could attenuate signals or introduce background noise. These factors in combination with the volume of survey data increases the time it takes to perform analysis of survey data when the source type is unknown. There are a number of algorithms in the field of machine learning that can be used to classify data that would be challenging and time-consuming for a human to identify. Supervised machine learning algorithms train models based on extensive amounts of human-labeled training data. Once sufficiently trained, these models can be used to quickly make high-fidelity predictions on new data. Convolutional neural networks are machine learning algorithms that excel in learning representations of 'shapes'. They do this by taking numerical input data and convolving them with spatial feature detectors referred to as filters. These filters are incrementally adjusted to reduce the prediction error on the data during the backpropagation step of training. Discussed in this paper is the development of a convolutional neural network classifier (CNNC) that can utilize spectral survey data for source discrimination and isotope identification. Bench-top laboratory experiments data using LaBr{sub 3}(Ce) scintillation detectors were used to train and evaluate the performance of the developed CNNC. The CNNC is capable of discriminating a variety of gamma emitting source types, differentiating different forms of uranium (depleted vs. natural), and estimating the amount of uranium for a known geometry. The discussed CNNC may be useful in scenarios where survey systems are deployed in situations where hazardous radioactive material maybe present, but the type is unknown. When used in remediation applications the CNNC can be used to screen-out false positives, helping reduce remediation costs. (authors)

07 ISOTOPE AND RADIATION SOURCES↗

Differential electron yield imaging with STXM

In this report total electron yield (TEY) imaging is an established scanning transmission X-ray microscopy (STXM) technique that gives varying contrast based on a sample's geometry, elemental composition, and electrical conductivity. However, the TEY-STXM signal is determined solely by the electrons that the beam ejects from the sample. A related technique, X-ray beam-induced current (XBIC) imaging, is sensitive to electrons and holes independently, but requires electric fields in the sample. Here we report that multi-electrode devices can be wired to produce differential electron yield (DEY) contrast, which is also independently sensitive to electrons and holes, but does not require an electric field. Depending on whether the region illuminated by the focused STXM beam is better connected to one electrode or another, the DEY-TXM contrast changes sign. DEY-STXM images thus provide a vivid map of a device's connectivity landscape, which can be key to understanding device function and failure. To demonstrate an application in the area of failure analysis, we image a 100 nm, lithographically-defined aluminum nanowire that has failed after being stressed with a large current density.

47 OTHER INSTRUMENTATION↗

Generating functionals and Lagrangian partial differential equations

The main goal of this paper is to derive an alternative characterization of the multisymplectic form formula for classical field theories using the geometry of the space of boundary values. We review the concept of Type-I/II generating functionals defined on the space of boundary data of a Lagrangian field theory. On the Lagrangian side, we define an analogue of Jacobi's solution to the Hamilton–Jacobi equation for field theories, and we show that by taking variational derivatives of this functional, we obtain an isotropic submanifold of the space of Cauchy data, described by the so-called multisymplectic form formula. As an example of the latter, we show that Lorentz's reciprocity principle in electromagnetism is a particular instance of the multisymplectic form formula. We also define a Hamiltonian analogue of Jacobi's solution, and we show that this functional is a Type-II generating functional. We finish the paper by defining a similar framework of generating functions for discrete field theories, and we show that for the linear wave equation, we recover the multisymplectic conservation law of Bridges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Numerical calculations of high-altitude differential charging: Preliminary results

A two dimensional simulation program was constructed in order to obtain theoretical predictions of floating potential distributions on geostationary spacecraft. The geometry was infinite-cylindrical with angle dependence. Effects of finite spacecraft length on sheath potential profiles can be included in an approximate way. The program can treat either steady-state conditions or slowly time-varying situations, involving external time scales much larger than particle transit times. Approximate, locally dependent expressions were used to provide space charge, density profiles, but numerical orbit-following is used to calculate surface currents. Ambient velocity distributions were assumed to be isotropic, beam-like, or some superposition of these.

Laframboise, J. G.↗

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗

Kinematic flow from the flow of cuts

The wavefunction coefficients of conformally coupled scalars in power-law FRW cosmologies satisfy differential equations governed by a set of simple combinatorial rules known as the kinematic flow. In this paper we derive the kinematic flow, expressed using a set of differential forms referred to as the cut basis, from a geometric perspective, relying solely on the cosmological hyperplane arrangement and without invoking bulk physics. Each element of the cut basis corresponds to the positive geometry associated to an independent cut of the physical FRW-form and can be labeled by decorating (minors of) the truncated Feynman graph with an acyclic orientation. We provide a straightforward prescription to associate a logarithmic differential form to each element of the cut basis by considering its corresponding decorated graph. Moreover, we show that the residues of the physical FRW-form are canonical forms of certain graphical zonotopes labeled by the same set of decorated graphs. These zonotopes control the cut combinatorics -- flow of cuts -- of the physical FRW-form and the cut basis (by construction). Using the theory of relative twisted cohomology and intersection theory, we derive a closed form formula for the differential equations of the cut basis. We also introduce combinatorial rules that compute the kinematic differential of any basis element without explicit calculation. The combinatorics of our differential equations is a natural consequence of the flow of cuts and is equivalent (up to rescaling) to the kinematic flow for the recently studied time integral basis. In particular, our differential equations decouple into exponentially many sectors, one for each way of cutting a subset of edges of the graph.

General Relativity and Quantum Cosmology↗

A physical zero-knowledge object-comparison system for nuclear warhead verification

Zero-knowledge proofs are mathematical cryptographic methods to demonstrate the validity of a claim while providing no further information beyond the claim itself. The possibility of using such proofs to process classified and other sensitive physical data has attracted attention, especially in the field of nuclear arms control. Here we demonstrate a non-electronic fast neutron differential radiography technique using superheated emulsion detectors that can confirm that two objects are identical without revealing their geometry or composition. Such a technique could form the basis of a verification system that could confirm the authenticity of nuclear weapons without sharing any secret design information. More broadly, by demonstrating a physical zero-knowledge proof that can compare physical properties of objects, this experiment opens the door to developing other such secure proof-systems for other applications.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

On current aspects of finite element computational fluid mechanics for turbulent flows

A set of nonlinear partial differential equations suitable for the description of a class of turbulent three-dimensional flow fields in select geometries is identified. On the basis of the concept of enforcing a penalty constraint to ensure accurate accounting of ordering effects, a finite element numerical solution algorithm is established for the equation set and the theoretical aspects of accuracy, convergence and stability are identified and quantized. Hypermatrix constructions are used to formulate the reduction of the computational aspects of the theory to practice. The robustness of the algorithm, and the computer program embodiment, have been verified for pertinent flow configurations.

Baker, A. J.↗

Domain decomposition: A bridge between nature and parallel computers

Domain decomposition is an intuitive organizing principle for a partial differential equation (PDE) computation, both physically and architecturally. However, its significance extends beyond the readily apparent issues of geometry and discretization, on one hand, and of modular software and distributed hardware, on the other. Engineering and computer science aspects are bridged by an old but recently enriched mathematical theory that offers the subject not only unity, but also tools for analysis and generalization. Domain decomposition induces function-space and operator decompositions with valuable properties. Function-space bases and operator splittings that are not derived from domain decompositions generally lack one or more of these properties. The evolution of domain decomposition methods for elliptically dominated problems has linked two major algorithmic developments of the last 15 years: multilevel and Krylov methods. Domain decomposition methods may be considered descendants of both classes with an inheritance from each: they are nearly optimal and at the same time efficiently parallelizable. Many computationally driven application areas are ripe for these developments. A progression is made from a mathematically informal motivation for domain decomposition methods to a specific focus on fluid dynamics applications. To be introductory rather than comprehensive, simple examples are provided while convergence proofs and algorithmic details are left to the original references; however, an attempt is made to convey their most salient features, especially where this leads to algorithmic insight.

Keyes, David E.↗

Features and applications of the Groove Analysis Program (GAP)

An IBM Personal Computer (PC) version of the Groove Analysis program (GAP) was developed to predict the steady state heat transport capability of an axially grooved heat pipe for a specified groove geometry and working fluid. In the model, the capillary limit is determined by the numerical solution of the differential equation for momentum conservation with the appropriate boundary conditions. This governing equation accounts for the hydrodynamic losses due to friction in liquid and vapor flows and due to liquid/vapor shear interaction. Back-pumping in both 0-g and 1-g is accounted for in the boundary condition at the condenser end. Slug formation in 0-g and puddle flow in 1-g are also considered in the model. At the user's discretion, the code will perform the analysis for various fluid inventories (undercharge, nominal charge, overcharge, or a fixed fluid charge) and heat pipe elevations. GAP will also calculate the minimum required heat pipe wall thickness for pressure containment at design temperatures that are greater than or lower than the critical temperature of the working fluid. This paper discusses the theory behind the development of the GAP model. It also presents the many useful and powerful capabilities of the model. Furthermore, a correlation of flight test performance data and the predictions using GAP are presented and discussed.

Ku, Jentung↗

Effect of Molding and Machining on Neoflon CTFE M400H Polychlorotrifluoroethylene Rod Stock and Valve Seat Properties

Since 1997 numerous fires have been reported to the Food and Drug Administration involving cylinder valves installed on medical use oxygen cylinders sold and operated within the United States. All of the cylinder valves in question had polychlorotrifluoroethylene (PCTFE) valve seats. Subsequent failure analysis showed that the main seat was the primary source of ignition. A review of the incidents involving cylinder valve fires indicated three possible ignition mechanisms: contaminant promotion, flow friction, and resonance. However, gas purity analysis showed that uncombusted, residual oxygen was within specification. Infrared and energy dispersive spectroscopy further showed that no contaminants or organic compounds were present in the remaining, uncombusted valve seat material or on seat plug surfaces. Therefore, contaminant-promoted ignition did not appear to be responsible for the failures. Observations of extruded material along the outer edge of the coined or loaded seat area produced by cylinder overuse or poppet overload led to concerns that accelerated gas flow across a deformed seat surface could generate enough localized heating to ignite the polymeric seat. Low molecular weight or highly amorphous quick-quenched PCTFE grades might be expected to be especially prone to this type of deformation. Such a failure mechanism has been described as flow friction; however, the corresponding mechanistic parameters are poorly understood. Subsequent revelation of low-temperature dimensional instability by thermomechanical analysis (TMA) in a variety of PCTFE sheet and rod stock samples led to new concerns that PCTFE valve seats could undergo excessive expansion or contraction during service. During expansion, additional extrusion and accompanying flow friction could occur. During contraction, a gap between the seal and adjacent metal surfaces could form. Gas flowing past the gap could, in turn, lead to resonance heating and subsequent ignition as described in ASTM Guide for Evaluation Nonmetallic Materials for Oxygen Service (G 63). Attempts to uncover the origins of the observed dimensional instability were hindered by uncertainties about resin grade, process history, and post-process heat history introduced by machining, annealing, and sample preparation. An approach was therefore taken to monitor property changes before and after processing and machining using a single, well-characterized lot of Neoflon CTFE.1 M400H resin. A task group consisting of the current PCTFE resin supplier, two molders, and four valve seat manufacturers was formed, and phased testing on raw resin, intermediate rod stock, and finished valve seats initiated. The effect of processing and machining on the properties of PCTFE rod stock and oxygen gas cylinder valve seats was then determined. Testing focused on two types of extruded rod stock and one type of compression-molded rod stock. To accommodate valve seat manufacturer preferences for certain rod stock diameters, two representative diameters were used (4.8 mm (0.1875 in.) and 19.1 mm (0.75 in.)). To encompass a variety of possible sealing configurations, seven different valve seat types with unique geometries or machining histories were tested. The properties investigated were dimensional stability as determined by TMA, specific gravity, differential scanning calorimetry (DSC), compressive strength, zero strength time, and intrinsic viscosity. Findings are discussed in the context of polymer structure-process-property relationships whenever possible.

Waller, Jess M.↗

Microfluidic Mixing Technology for a Universal Health Sensor

A highly efficient means of microfluidic mixing has been created for use with the rHEALTH sensor an elliptical mixer and passive curvilinear mixing patterns. The rHEALTH sensor provides rapid, handheld, complete blood count, cell differential counts, electrolyte measurements, and other lab tests based on a reusable, flow-based microfluidic platform. These geometries allow for cleaning in a reusable manner, and also allow for complete mixing of fluid streams. The microfluidic mixing is performed by flowing two streams of fluid into an elliptical or curvilinear design that allows the combination of the flows into one channel. The mixing is accomplished by either chaotic advection around micro - fluidic loops. All components of the microfluidic chip are flow-through, meaning that cleaning solution can be introduced into the chip to flush out cells, plasma proteins, and dye. Tests were performed on multiple chip geometries to show that cleaning is efficient in any flowthrough design. The conclusion from these experiments is that the chip can indeed be flushed out with microliter volumes of solution and biological samples are cleaned readily from the chip with minimal effort. The technology can be applied in real-time health monitoring at patient s bedside or in a doctor s office, and real-time clinical intervention in acute situations. It also can be used for daily measurement of hematocrit for patients on anticoagulant drugs, or to detect acute myocardial damage outside a hospital.

Chan, Eugene Y.↗

Thermal and hydraulic performance of volumetrically heated triply periodic minimal surface heaters

Advanced heat transfer surfaces, mainly Triply Periodic Minimal Surfaces (TPMSs) have great potential to increase heat transfer performance over traditional heat transfer surfaces due to their tortuous flow path and higher surface area-volume ratio. The emergence of additive manufacturing of nuclear fuel motivates experimental investigation of such advanced geometries in the context of a nuclear reactor core. Here, this study employed a novel method to explore the performance of gyroid and diamond TPMS geometries. Volumetrically heated heaters were additively manufactured from conductive polymer filament, and airflow was used to simulate the convective heat transfer from the core within a nuclear reactor. The mass flow and volumetric generation rates were varied, local Nusselt number correlations were obtained using embedded thermocouples to measure the solid and fluid temperature, and friction factor correlations were developed using differential pressure measurements. The gyroid TPMS presented a friction factor 1.5 times higher than the diamond TPMS and 90 times that of the rod bundle. TPMS geometries showed Nusselt numbers 8–10 times higher than that of the rod bundle, also maintaining colder internal temperatures, which suggested that they could sustain 250–275 % higher power density than the rod bundle for a given centerline temperature. Hence, having the potential to greatly increase the safety margins and power output of a nuclear reactor core with similar volume.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Differential equations from unitarity cuts: nonplanar hexa-box integrals

We compute ϵ-factorized differential equations for all dimensionally-regularized integrals of the nonplanar hexa-box topology, which contribute for instance to 2-loop 5-point QCD amplitudes. A full set of pure integrals is presented. For 5-point planar topologies, Gram determinants which vanish in 4 dimensions are used to build compact expressions for pure integrals. Using unitarity cuts and computational algebraic geometry, we obtain a compact IBP system which can be solved in 8 hours on a single CPU core, overcoming a major bottleneck for deriving the differential equations. Alternatively, assuming prior knowledge of the alphabet of the nonplanar hexa-box, we reconstruct analytic differential equations from 30 numerical phase-space points, making the computation almost trivial with current techniques. We solve the differential equations to obtain the values of the master integrals at the symbol level. Full results for the differential equations and solutions are included as supplementary material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analytical solution of the Poisson-Nernst-Planck equations for an electrochemical system close to electroneutrality

Single charge densities and the potential are used to describe models of electrochemical systems. These quantities can be calculated by solving a system of time dependent nonlinear coupled partial differential equations, the Poisson-Nernst-Planck equations. Assuming small deviations from the electroneutral equilibrium, the linearized and decoupled equations are solved for a radial symmetric geometry, which represents the interface between a cell and a sensor device. The densities and the potential are expressed by Fourier-Bessels series. The system considered has a ratio between the Debye-length and its geometric dimension on the order of 10{sup −4} so the Fourier-Bessel series can be approximated by elementary functions. The time development of the system is characterized by two time constants, τ{sub c} and τ{sub g}. The constant τ{sub c} describes the approach to the stationary state of the total charge and the potential. τ{sub c} is several orders of magnitude smaller than the geometry-dependent constant τ{sub g}, which is on the order of 10 ms characterizing the transition to the stationary state of the single ion densities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Kinematical vortices in double photoionization of helium by attosecond pulses

Two-armed helical vortex structures are predicted in the two-electron momentum distributions produced in double photoionization (DPI) of the He atom by a pair of time-delayed elliptically polarized attosecond pulses with opposite helicities. These predictions are based upon both a first-order perturbation theory analysis and numerical solutions of the two-electron, time-dependent Schrödinger equation in six spatial dimensions. The helical vortex structures originate from Ramsey interference of a pair of ionized two-electron wave packets, each having a total angular momentum of unity, and appear in the sixfold differential DPI probability distribution for any energy partitioning between the two electrons. The vortex structures are exquisitely sensitive to the time delay between the two pulses, their relative phase, their ellipticity, and their handedness; moreover, they occur in a variety of electron detection geometries. However, the vortex structures only occur when the angular separation β = cos -1 (ˆ p 1 ∙ ˆ p 2 ) between the electron momenta p 1 and p 2 is held fixed. The vortex structures can also be observed in the fourfold differential DPI probability distribution obtained by averaging the sixfold differential probability over the emission angles of one electron. Such kinematical vortices are a general phenomenon that may occur in any ionization process, initiated by two time-delayed short pulses with opposite ellipticities, for particular detection geometries.

74 ATOMIC AND MOLECULAR PHYSICS↗