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 739 records · Page 41

The MERiT Onboard the CeREs: A Novel Instrument to Study Energetic Particles in the Earth's Radiation Belts

The Miniaturized Electron pRoton Telescope, MERiT, is a low‐mass, low‐power, compact instrument using an innovative combination of particle detectors, sensor electronics, and onboard processing. MERiT is flying on the Compact Radiation belt Explorer, CeREs, a 3U CubeSat launched into a low earth orbit of 500‐km altitude and inclination of 85° on 16 December 2018. The primary and secondary science goals of CeREs are to investigate electron microbursts and to study solar particles. MERiT comprises a stack of solid state detectors (SSD) behind space facing avalanche photo diodes (APDs) surrounded by W‐Al shielding to reduce side‐penetrating particle background. The APD‐SSD combination enables measurement of electrons from 5 to 200 keV and 1 to 8 MeV; protons from 200–400 keV and 7–100 MeV in differential channels with energy resolution ΔE/E≈30% for both electrons and protons. MERiT measures microbursts with a high time resolution ranging from 4 to 16 ms and solar particles with a cadence of 1 s. MERiT energy channels and cadences are software configurable via algorithms and lookup tables residing on a field‐programmable gate array. The lookup tables can be changed via ground commands. MERiT geometry factor is 31 sq.cm‐sr and optimized to measure microbursts with the instrument viewing the local zenith in orbit. MERiT enables investigation of dynamical processes of radiation belt electron energization and loss, solar electron and proton transport, and their access to the Earth's polar caps. We describe the MERiT sensor design, calibration, operational modes, data products, and science goals.

Kanekal, S. G.↗

Combined UHV/high-pressure catalysis setup for depth-resolved near-surface spectroscopic characterization and catalytic testing of model catalysts

An ultra-high vacuum (UHV) setup for “real” and “inverse” model catalyst preparation, depth-resolved near-surface spectroscopic characterization, and quantification of catalytic activity and selectivity under technologically relevant conditions is described. Due to the all-quartz reactor attached directly to the UHV-chamber, transfer of the catalyst for in situ testing without intermediate contact to the ambient is possible. The design of the UHV-compatible re-circulating batch reactor setup allows the study of reaction kinetics under close to technically relevant catalytic conditions up to 1273 K without contact to metallic surfaces except those of the catalyst itself. With the attached differentially pumped exchangeable evaporators and the quartz-microbalance thickness monitoring equipment, a reproducible, versatile, and standardised sample preparation is possible. For three-dimensional near-surface sample characterization, the system is equipped with a hemispherical analyser for X-ray photoelectron spectroscopy (XPS), electron-beam or X-ray-excited Auger-electron spectroscopy, and low-energy ion scattering measurements. Due the dedicated geometry of the X-ray gun (54.7°, “magic angle”) and the rotatable sample holder, depth analysis by angle-resolved XPS measurements can be performed. Thus, by the combination of characterisation methods with different information depths, a detailed three-dimensional picture of the electronic and geometric structure of the model catalyst can be obtained. To demonstrate the capability of the described system, comparative results for depth-resolved sample characterization and catalytic testing in methanol steam reforming on PdGa and PdZn near-surface intermetallic phases are shown.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

WE-EF-207-09: Single-Scan Dual-Energy CT Using Primary Modulation

Purpose: Compared with conventional CT, dual energy CT (DECT) provides better material differentiation but requires projection data with two different effective x-ray spectra. Current DECT scanners use either a two-scan setting or costly imaging components, which are not feasible or available on open-gantry cone-beam CT systems. We propose a hardware-based method which utilizes primary modulation to enable single-scan DECT on a conventional CT scanner. The CT imaging geometry of primary modulation is identical to that used in our previous method for scatter removal, making it possible for future combination with effective scatter correction on the same CT scanner. Methods: We insert an attenuation sheet with a spatially-varying pattern - primary modulator-between the x-ray source and the imaged object. During the CT scan, the modulator selectively hardens the x-ray beam at specific detector locations. Thus, the proposed method simultaneously acquires high and low energy data. High and low energy CT images are then reconstructed from projections with missing data via an iterative CT reconstruction algorithm with gradient weighting. Proof-of-concept studies are performed using a copper modulator on a cone-beam CT system. Results: Our preliminary results on the Catphan(c) 600 phantom indicate that the proposed method for single-scan DECT is able to successfully generate high-quality high and low energy CT images and distinguish different materials through basis material decomposition. By applying correction algorithms and using all of the acquired projection data, we can reconstruct a single CT image of comparable image quality to conventional CT images, i.e., without primary modulation. Conclusion: This work shows great promise in using a primary modulator to perform high-quality single-scan DECT imaging. Future studies will test method performance on anthropomorphic phantoms and perform quantitative analyses on image qualities and DECT decomposition accuracy. We will use simulations to optimize the modulator material and geometry parameters.

60 APPLIED LIFE SCIENCES↗

Slow growth of cracks in a rate sensitive Tresca solid.

Extension of the classical theory of fracture to viscoelastic and elastic-plastic materials in which the plasticity effects are confined to a narrow band encompassing the crack front. It is suggested that the Griffith-Irwin criterion of fracture, which requires that the energy release rate computed for a given boundary value problem equal the critical threshold, ought to be replaced by a differential equation governing the slow growth of a crack prior to the onset of rapid propagation. A new term which enters the equation of motion in the dissipative media is proportional to the energy lost within the end sections of the crack, and thus it reflects the extent of inelastic behavior of a solid. A concept of 'apparent' surface energy is introduced to account for the geometry-dependent and the rate-dependent phenomena which influence toughness of an inelastic solid. Three hypotheses regarding the condition for fracture in the subcritical range of load are compared.

Wnuk, M. P.↗

Simulation of confined magnetohydrodynamic flows with Dirichlet boundary conditions using a pseudo-spectral method with volume penalization

A volume penalization approach to simulate magnetohydrodynamic (MHD) flows in confined domains is presented. Here the incompressible visco-resistive MHD equations are solved using parallel pseudo-spectral solvers in Cartesian geometries. The volume penalization technique is an immersed boundary method which is characterized by a high flexibility for the geometry of the considered flow. In the present case, it allows to use other than periodic boundary conditions in a Fourier pseudo-spectral approach. The numerical method is validated and its convergence is assessed for two- and three-dimensional hydrodynamic (HD) and MHD flows, by comparing the numerical results with results from literature and analytical solutions. The test cases considered are two-dimensional Taylor–Couette flow, the z-pinch configuration, three dimensional Orszag–Tang flow, Ohmic-decay in a periodic cylinder, three-dimensional Taylor–Couette flow with and without axial magnetic field and three-dimensional Hartmann-instabilities in a cylinder with an imposed helical magnetic field. Finally, we present a magnetohydrodynamic flow simulation in toroidal geometry with non-symmetric cross section and imposing a helical magnetic field to illustrate the potential of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Probability of Initiation in Neutron Transport

We discuss the numerical solution of the nonlinear integro-differential equation for the probability of a divergent neutron chain in a stationary system (i.e., the probability of initiation (POI)). We follow the development described in Bell’s classic paper on the stochastic theory of neutron transport. As noted by Bell, the linearized form of this equation resembles the linear adjoint neutron transport equation. A matrix formalism for the discretized steady state (or forward) neutron equation in slab geometry is first developed and is then used to derive the discrete adjoint equation. A main advantage of this discrete development is that the resulting discrete adjoint equation does not depend upon how the multigroup cross sections for the forward problem are obtained. That is, we derive the discrete adjoint directly from the discrete forward equations rather than discretizing directly the adjoint equation. This also guarantees that the discrete adjoint operator is consistent with the inner product used to define the adjoint operator. We discuss three approaches for the numerical solution of the POI equations, and present numerical results on several test problems. The three solution methods are a simple fixed-point iteration, a second approach that is akin to a nonlinear Power iteration, and a third approach which uses a Newton-Krylov nonlinear solver. We also give sufficient conditions to guarantee the existence and uniqueness of nontrivial solutions to our discrete POI equations when the discrete system is supercritical, and that only the trivial solution exists when the discrete system is subcritical. Our approach is modeled after the analysis presented for the continuous POI equations by Mokhtar-Kharroubi and Jarmouni-Idrissi, and by Pazy and Rabinowitz.

42 ENGINEERING↗

Multi-scale gravity and cosmology

The gravitational dynamics and cosmological implications of three classes of recently introduced multi-scale spacetimes (with, respectively, ordinary, weighted and q-derivatives) are discussed. These spacetimes are non-Riemannian: the metric structure is accompanied by an independent measure-differential structure with the characteristics of a multi-fractal, namely, different dimensionality at different scales and, at ultra-short distances, a discrete symmetry known as discrete scale invariance. Under this minimal paradigm, five general features arise: (a) the big-bang singularity can be replaced by a finite bounce, (b) the cosmological constant problem is reinterpreted, since accelerating phases can be mimicked by the change of geometry with the time scale, without invoking a slowly rolling scalar field, (c) the discreteness of geometry at Planckian scales can leave an observable imprint of logarithmic oscillations in cosmological spectra and (d) give rise to an alternative mechanism to inflation or (e) to a fully analytic model of cyclic mild inflation, where near scale invariance of the perturbation spectrum can be produced without strong acceleration. Various properties of the models and exact dynamical solutions are discussed. In particular, the multi-scale geometry with weighted derivatives is shown to be a Weyl integrable spacetime.

79 ASTRONOMY AND ASTROPHYSICS↗

Progress Report on Model Development for the Transport of Aerosol through Microchannels

This report summarizes the current progress in the development of a phenomenological model of aerosol transport, deposition, and plugging through microchannels. The purpose is to introduce a generic, reliable numerical model for the prediction of aerosol transport, deposition, and plugging in leak paths while accounting for potential plugging formation, to a user community involving researchers, regulators, and industry. In that regard, a Graphical User Interface (GUI) was generated by integrating the individual MATLAB scripts that make up the model and adding additional features to aid the general understanding of the problem to a user. This report focuses on the model development and features included. Further, predictions from the GUI are compared with experimental data for validation. The strength of the GUI is the ability of a user to plug in basic parameters such as the initial pressure conditions and canister model specifications in order to obtain a first principal approximation of vital information such as the blowdown pressure differential, aerosol penetration, and deposition as a function of time. The user can obtain this without the know-how of the underlying MATLAB scripts and thus enables the model to be readily applied by regulators, industry, and shareholders to reduce the uncertainty in the off-site radiological consequences. The report also lists future tasks under this work scope. This includes ongoing efforts to include additional crack geometries (divergent and divergent-convergent slots) to the model; GUI development and improvement based on feedback from users, and integration of aerosol source term data from Sibling Pin tests into the model as we approach realistic canister stress corrosion cracking–induced aerosol release scenarios.

54 ENVIRONMENTAL SCIENCES↗

A new form of the Saint-Venant equations for variable topography

The solution stability of river models using the one-dimensional (1D) Saint-Venant equations can be easily undermined when source terms in the discrete equations do not satisfy the Lipschitz smoothness condition for partial differential equations. Although instability issues have been previously noted, they are typically treated as model implementation issues rather than as underlying problems associated with the form of the governing equations. This study proposes a new reference slope form of the Saint-Venant equations to ensure smooth slope source terms and eliminate one source of potential numerical oscillations. It is shown that a simple algebraic transformation of channel geometry provides a smooth reference slope while preserving the correct cross-section flow area and the total Piezometric pressure gradient that drives the flow. The reference slope method ensures the slope source term in the governing equations is Lipschitz continuous while maintaining all the underlying complexity of the real-world geometry. The validity of the mathematical concept is demonstrated with the open-source Simulation Program for River Networks (SPRNT) model in a series of artificial test cases and a simulation of a small urban creek. Validation comparisons are made with analytical solutions and the Hydrologic Engineering Center's River Analysis System (HEC-RAS) model. The new method reduces numerical oscillations and instabilities without requiring ad hoc smoothing algorithms.

54 ENVIRONMENTAL SCIENCES↗

Simplified calculation of solar flux distribution on the side wall of cylindrical cavity solar receivers

A simple solar flux calculation algorithm for a cylindrical cavity type solar receiver has been developed and implemented on an IBM PC-AT. Using cone optics, the contour error method is utilized to handle the slope error of a paraboloidal concentrator. The flux distribution on the side wall is calculated by integration of the energy incident from cones emanating from all the differential elements on the concentrator. The calculations are done for any set of dimensions and properties of the receiver and the concentrator, and account for any spillover on the aperture plate. The results of this algorithm compared excellently with those predicted by more complicated programs. Because of the utilization of axial symmetry and overall simplification, it is extremely fast. It can be esily extended to other axisymmetric receiver geometries.

Bhandari, P.↗

Differential spacecraft tracking by interferometry

This study estimates measurement system errors for two space vehicles on the surface of Mars, and for two Mars orbiting spacecraft, which are being tracked by differential interferometry. In these examples, signals from all spacecraft lie within the same beamwidth of an earth-based radio antenna. The measurements of all spacecraft signals are made simultaneously; errors that scale with angular source separation or with temporal separation between measurement epochs are practically removed. It is shown that errors due to system thermal noise and to systematic effects within ground receiver electronics dominate, except for geometries when signals pass close to the sun, when solar plasma becomes the dominant error source. The instantaneous relative position of two orbiters may be measured to within ten meters, leading to 50-meter three-dimensional orbital accuracy.

Border, James S.↗

Spectral element methods: Algorithms and architectures

Spectral element methods are high-order weighted residual techniques for partial differential equations that combine the geometric flexibility of finite element methods with the rapid convergence of spectral techniques. Spectral element methods are described for the simulation of incompressible fluid flows, with special emphasis on implementation of spectral element techniques on medium-grained parallel processors. Two parallel architectures are considered: the first, a commercially available message-passing hypercube system; the second, a developmental reconfigurable architecture based on Geometry-Defining Processors. High parallel efficiency is obtained in hypercube spectral element computations, indicating that load balancing and communication issues can be successfully addressed by a high-order technique/medium-grained processor algorithm-architecture coupling.

Fischer, Paul↗

Spectral element methods - Algorithms and architectures

Spectral element methods are high-order weighted residual techniques for partial differential equations that combine the geometric flexibility of finite element methods with the rapid convergence of spectral techniques. Spectral element methods are described for the simulation of incompressible fluid flows, with special emphasis on implementation of spectral element techniques on medium-grained parallel processors. Two parallel architectures are considered; the first, a commercially available message-passing hypercube system; the second, a developmental reconfigurable architecture based on Geometry-Defining Processors. High parallel efficiency is obtained in hypercube spectral element computations, indicating that load balancing and communication issues can be successfully addressed by a high-order technique/medium-grained processor algorithm-architecture coupling.

Fischer, Paul↗

Adhesive-Bonded Composite Joint Analysis with Delaminated Surface Ply Using Strain-Energy Release Rate

This paper presents an analytical model to determine the strain energy release rate due to an interlaminar crack of the surface ply in adhesively bonded composite joints subjected to axial tension. Single-lap shear-joint standard test specimen geometry with thick bondline is followed for model development. The field equations are formulated by using the first-order shear-deformation theory in laminated plates together with kinematics relations and force equilibrium conditions. The stress distributions for the adherends and adhesive are determined after the appropriate boundary and loading conditions are applied and the equations for the field displacements are solved. The system of second-order differential equations is solved to using the symbolic computation tool Maple 9.52 to provide displacements fields. The equivalent forces at the tip of the prescribed interlaminar crack are obtained based on interlaminar stress distributions. The strain energy release rate of the crack is then determined by using the crack closure method. Finite element analyses using the J integral as well as the crack closure method are performed to verify the developed analytical model. It has been shown that the results using the analytical method correlate well with the results from the finite element analyses. An attempt is made to predict the failure loads of the joints based on limited test data from the literature. The effectiveness of the inclusion of bondline thickness is justified when compared with the results obtained from the previous model in which a thin bondline and uniform adhesive stresses through the bondline thickness are assumed.

Chadegani, Alireza↗

Biopolymer foams - Relationship between material characteristics and foaming behavior of cellulose based foams

Biopolymers are becoming increasingly important to both industry and consumers. With regard to waste management, CO{sub 2} balance and the conservation of petrochemical resources, increasing efforts are being made to replace standard plastics with bio-based polymers. Nowadays biopolymers can be built for example of cellulose, lactic acid, starch, lignin or bio mass. The paper will present material properties of selected cellulose based polymers (cellulose propionate [CP], cellulose acetate butyrate [CAB]) and corresponding processing conditions for particle foams as well as characterization of produced parts. Special focus is given to the raw material properties by analyzing thermal behavior (differential scanning calorimetry), melt strength (Rheotens test) and molecular weight distribution (gel-permeation chromatography). These results will be correlated with the foaming behavior in a continuous extrusion process with physical blowing agents and underwater pelletizer. Process set-up regarding particle foam technology, including extrusion foaming and pre-foaming, will be shown. The characteristics of the resulting foam beads will be analyzed regarding part density, cell morphology and geometry. The molded parts will be tested on thermal conductivity as well as compression behavior (E-modulus, compression strength)

36 MATERIALS SCIENCE↗

Deformation T-Cup: A new multi-anvil apparatus for controlled strain-rate deformation experiments at pressures above 18 GPa

A new multi-anvil deformation apparatus, based on the widely used 6-8 split-cylinder, geometry, has been developed which is capable of deformation experiments at pressures in excess of 18 GPa at room temperature. In 6-8 (Kawai-type) devices eight cubic anvils are used to compress the sample assembly. In our new apparatus two of the eight cubes which sit along the split-cylinder axis have been replaced by hexagonal cross section anvils. Combining these anvils hexagonal-anvils with secondary differential actuators incorporated into the load frame, for the first time, enables the 6-8 multi-anvil apparatus to be used for controlled strain-rate deformation experiments to high strains. Testing of the design, both with and without synchrotron-X-rays, has demonstrated the Deformation T-Cup (DT-Cup) is capable of deforming 1–2 mm long samples to over 55% strain at high temperatures and pressures. To date the apparatus has been calibrated to, and deformed at, 18.8 GPa and deformation experiments performed in conjunction with synchrotron X-rays at confining pressures up to 10 GPa at 800 °C.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Machine Learning to Select Experiments Driven by Fundamental Science and Applications for Targeted Nuclear Data Improvement

This work describes a blueprint for a process that accelerates progress in science by quantitatively answering the following question: What is the optimal combination of fundamental-science and application-driven experiments to maximally reduce pertinent data uncertainties? Answering this question entails solving a high-dimensional and complex optimization problem that is best solved with advanced statistic techniques often classified as machine learning. We apply this process within the framework of nuclear data with the aim to select an experiment combination that will reduce uncertainties in 239 Pu nuclear data for neutron energies between 1 and 600 keV. In this field, fundamental-physics driven data, called differential, look at one nuclear physics observable at a time. They are contrasted to application-driven, integral, data where one or few resulting values inform a broad set of nuclear data across several nuclides and energies. The candidates for integral experiments are criticality measurements that were refined by a genetic algorithm to be maximally sensitive to 239 Pu fission cross sections in the desired energy range. Twenty-three candidate differential experiments were investigated and span multiple nuclear physics observables (e.g., total, capture cross sections) for isotopes appearing in the integral experiments. The optimal combination among these candidate experiments was investigated via generalized least squares fitting, augmented with Gaussian processes to ameliorate statistical irregularities in data, and the D-optimality criterion. The latter evaluates for each pair of candidates the joint reduction in uncertainties of all 12200 nuclear data appearing in the integral experiments compared to the knowledge we have from 168 past experiments, theory, and nuclear data. We chose as differential measurements those that investigate 63 Cu and 239 Pu total cross sections, based on D-optimality rank and feasibility constraints. Two integral (criticality) experiments were selected: An experiment with Al 2 ⁢O 3 and graphite interleaved with Pu and a thick Cu reflector explores 1–30 keV, while we target the 30–600 keV range with an experiment that swaps boron in place of graphite with a different geometry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On covariant Poisson brackets in classical field theory

How to give a natural geometric definition of a covariant Poisson bracket in classical field theory has for a long time been an open problem—as testified by the extensive literature on “multisymplectic Poisson brackets,” together with the fact that all these proposals suffer from serious defects. On the other hand, the functional approach does provide a good candidate which has come to be known as the Peierls–De Witt bracket and whose construction in a geometrical setting is now well understood. Here, we show how the basic “multisymplectic Poisson bracket” already proposed in the 1970s can be derived from the Peierls–De Witt bracket, applied to a special class of functionals. This relation allows to trace back most (if not all) of the problems encountered in the past to ambiguities (the relation between differential forms on multiphase space and the functionals they define is not one-to-one) and also to the fact that this class of functionals does not form a Poisson subalgebra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗