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 505 records · Page 28

Development of a High Fidelity Dynamic Module of the Advanced Resistive Exercise Device (ARED) Using Adams

NASA's Digital Astronaut Project (DAP) implements well-vetted computational models to predict and assess spaceflight health and performance risks, and enhance countermeasure development. DAP provides expertise and computation tools to its research customers for model development, integration, or analysis. DAP is currently supporting the NASA Exercise Physiology and Countermeasures (ExPC) project by integrating their biomechanical models of specific exercise movements with dynamic models of the devices on which the exercises were performed. This presentation focuses on the development of a high fidelity dynamic module of the Advanced Resistive Exercise Device (ARED) on board the ISS. The ARED module, illustrated in the figure below, was developed using the Adams (MSC Santa Ana, California) simulation package. The Adams package provides the capabilities to perform multi rigid body, flexible body, and mixed dynamic analyses of complex mechanisms. These capabilities were applied to accurately simulate: Inertial and mass properties of the device such as the vibration isolation system (VIS) effects and other ARED components, Non-linear joint friction effects, The gas law dynamics of the vacuum cylinders and VIS components using custom written differential state equations, The ARED flywheel dynamics, including torque limiting clutch. Design data from the JSC ARED Engineering team was utilized in developing the model. This included solid modeling geometry files, component/system specifications, engineering reports and available data sets. The Adams ARED module is importable into LifeMOD (Life Modeler, Inc., San Clemente, CA) for biomechanical analyses of different resistive exercises such as squat and dead-lift. Using motion capture data from ground test subjects, the ExPC developed biomechanical exercise models in LifeMOD. The Adams ARED device module was then integrated with the exercise subject model into one integrated dynamic model. This presentation will describe the development of the Adams ARED module including its capabilities, limitations, and assumptions. Preliminary results, validation activities, and a practical application of the module to inform the relative effect of the flywheels on exercise will be discussed.

Humphreys, B. T.↗

[Interview]: Alexandre Shvartsburg, Pacific Northwest National Laboratory, Richland, WA, USA

Q1. What are your main research activities in ion mobility mass spectrometry (past or present)? My early efforts focused on the structural characterization of atomic (carbon and semiconductor) clusters. After the production of bulk fullerenes, many hoped that other nanoclusters discovered in the gas phase could also coalesce into new materials. As these studies required accurate and robust mobility calculations for any ion geometry, I strived to build the needed theory and implement it in the Mobcal software widely employed today. Since 2004, I have been developing methods and novel applications of differential IMS (FAIMS) at PNNL. The principal achievement has been raising the resolving power by over tenfold (up to ~400 for multiply-charged peptides) using elevated fields, helium and hydrogen-rich buffers, and extended filtering times. This performance broadly allows previously unthinkable separations of very similar species, for example sequence inversions and post-translational modification localization isomers of peptides (including “middle-down” peptides such as histone tails), lipid regioisomers, and even isotopomers. Another major direction is investigating the dipole alignment of larger proteins, which creates an exceptionally strong FAIMS effect that is a potential tool for structural biology. Q2: What have been the most significant instrumentation or applications developments in the history of ion mobility - mass spectrometry? In 1995 when I started graduate research at Northwestern, only two groups worldwide worked with IMS/MS and “the literature” meant papers by Bowers (UCSB). Well-wishers counseled me to “learn something useful like HPLC, as IMS would never have real utility”. This booklet showcases the scale of change since. First, the practical IMS/ToF platforms for complex biological analyses demonstrated by Clemmer have turned IMS/MS from an esoteric physical chemistry technique into a powerful analytical tool. By commercializing the IMS/ToF technology in Synapt instruments, Waters has greatly increased its impact via expanded number and diversity of applications. Concurrently, Guevremont at Canadian NRC has perfected FAIMS coupled to MS, deployed it for real-world bio and environmental analyses, and widely distributed it in the Ionalytics Selectra system (subsequently installed on Thermo MS platforms). The latest breakthrough is ultra-FAIMS by Owlstone, where extreme fields allow numerous qualitatively new separations and operational modes that we just begin to explore. Q3: Where do you see ion mobility - mass spectrometry making the most impact in the next 5 years? Any predictions for where the field will go? Sciences dealing with perturbations in media (such as optics or acoustics) at some point shift from the linear to nonlinear paradigm, where propagation depends on the magnitude of perturbation or its driving force. While the linear part remains industrially important (e.g., eyewear and architectural glass for optics), frontline research moves to nonlinear phenomena. IMS is undergoing that transition now with the rise of FAIMS, which should continue as the fundamental understanding improves, new modalities and applications emerge, and more instrumentation is introduced by vendors. Modifying and augmenting FAIMS separations through vapor dopants that render ion mobilities less linear is becoming routine. I expect this area to advance, extending to more specific interactions and to complexation with solution additives. Another route to higher separation power is integrating FAIMS with conventional IMS; proliferation of both technologies would make such 2-D platforms common. Along with mass spectrometry and conventional IMS, FAIMS will address increasingly large macromolecules, including proteins and their complexes.

Shvartsburg, Alexandre A.↗

M3 milestone: Advanced contact 2021

This report describes advanced contact approaches and implementations in FY21 for use in emerging reactors with geometries that differ from LWR/EBR-II cylindrical fuel. Contact remains a challenging and crucially important component of our modeling toolkit. While significant improvements have been made through mortar and automatic differentiation, we explored advanced approaches (i.e., dual mortar) that yields better conditioned systems and will improve the robustness and performance of BISON contact problems. This effort is relevant to emerging reactor types, such as micro reactors and the corresponding complex and numerous contact surfaces. In this report, improvements of the mortar-based contact formulation are described. Benchmark problems and BISON assessment cases are included to demonstrate the improvements and advantages that are brought by the recent capability. Remaining issues are discussed and corresponding development plans are outlined.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The Holometer: An Instrument to Probe Planckian Quantum Geometry

This paper describes the Fermilab Holometer, an instrument for measuring correlations of position variations over a four-dimensional volume of space-time. The apparatus consists of two co-located, but independent and isolated, 40 m power-recycled Michelson interferometers, whose outputs are cross-correlated to 25 MHz. The data are sensitive to correlations of differential position across the apparatus over a broad band of frequencies up to and exceeding the inverse light crossing time, 7.6 MHz. A noise model constrained by diagnostic and environmental data distinguishes among physical origins of measured correlations, and is used to verify shot-noise-limited performance. These features allow searches for exotic quantum correlations that depart from classical trajectories at spacelike separations, with a strain noise power spectral density sensitivity smaller than the Planck time. The Holometer in current and future configurations is projected to provide precision tests of a wide class of models of quantum geometry at the Planck scale, beyond those already constrained by currently operating gravitational wave observatories.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Reverse-Mode Algorithmic Differentiation of an OpenMP-Parallel Compressible Flow Solver

y Reverse-mode algorithmic differentiation (AD) is an established method for obtaining adjoint derivatives of computer simulation applications. In computational fluid dynamics (CFD), adjoint derivatives of a cost function output such as drag or lift with respect to design parameters such as surface coordinates or geometry control points are a key ingredient for shape optimization, uncertainty quantification and flow control. The computational cost of CFD applications and their derivatives makes it essential to use high-performance computing hardware efficiently, including multi- and many-core architectures. Nevertheless, OpenMP is not supported in most AD tools, and previously shown methods achieve poor scalability of the derivative code. We present the AD of an OpenMP-parallelized finite volume compressible flow solver for unstructured meshes. Our approach enables us to reuse the parallelization of the original code in the computation of adjoint derivatives. The method works by identifying code segments that can be differentiated in reverse-mode without changing their memory access pattern. The OpenMP parallelization is integrated into the derivative code during the build process in a way that is robust to modifications of the original code and independent of the OpenMP support of the differentiation tool. We show the scalability of our adjoint CFD solver on test cases ranging from thousands to millions of finite volume mesh cells on CPUs with up to 16 threads as well as on an Intel XeonPhi card with 236 threads. We demonstrate that our approach is more practical to implement for production-sized CFD codes and produces more efficient adjoint derivative code than previously shown AD methods.

Huckelheim, Jan↗

Modulating Transition Metal Reactivity with Force

The reactivity and selectivity of a transition metal catalyst is intimately related to its ligand-sphere geometry, and, in many cases, the ideal ligand geometry for one step of a catalytic cycle is poorly matched to the ideal ligand geometry for another. For this reason, methods for reversibly modulating ligand geometry on the time scale of catalytic turnover or monomer enchainment are highly desirable. Mechanical force represents a heretofore untapped approach to modulate catalyst geometry and/or reactivity, with the potential to do so on the timescale of catalytic turnover or monomer enchainment. Macroscopic mechanical forces are large, directional and localized to an extent that differentiates them from other forms of energy input such as heat or light. Here, in this Concept, we describe our efforts to address the fundamental challenges associated with force-modulated transition metal catalysis by employing molecular force probe ligands comprising a stiff stilbene photoswitch tethered to rotationally flexible biaryl bisphosphine ligand. Our efforts to date include the modulation of catalytic activity through force-mediated ligand perturbations, quantification of the force-coupled ligand effects on the energetics of elementary organometallic transformations, and evaluation of the mechanisms of force transduction in these systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nanostructural and rheological transitions of pH-responsive supramolecular systems involving a zwitterionic amphiphile and a triamine

Supramolecular aqueous complexes of a synthesized zwitterionic surfactant, 2-(dimethyl(octadecyl)ammonio)acetate (stearyl betaine) and diethylenetriamine, have been found to display pH-tunable rheological properties without the need for a co-surfactant. For this supramolecular system, the dynamically changing interplay among the constituents is responsible for the rich phase behavior including hollow nanotubules, bilayers, and vesicles. A further hypothesis is that hollow nanotubular structure is highly desirable for achieving chain entanglement geometries with reduced mesh size. The rheological properties of the suspension were measured using a rotational rheometer. An analysis using zero-shear viscosities at various temperatures as well as differential scanning calorimetry(DSC) helped shed light on thermodynamic and phase transitions taking place. The morphology of the suspension was analyzed by using a combination of small angle x-ray scattering(SAXS) and atomic force microscopy(AFM). Large steady-shear viscosities (up to ~160 Pa.s) and viscoelastic behavior was observed in acidic conditions which diminished upon increasing the pH. A novel analysis consisting of measured thermodynamic quantities such as activation energy(E a ) and enthalpy of transition(ΔH T ) further allude to the rich pH-sensitive nanoarchitecture of the system. Rather than commonly observed wormlike micelles, hollow nanotubules, as long as 1.6 µm were observed at pH 4 which disintegrated to bilayer sheets and vesicles upon increasing pH. Finally, gaining a comprehensive understanding from all the results, a transition mechanism between different suspension architectures was proposed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NEK5000 Assessment Milestone Report: Single-Phase Natural Circulation using Direct Numerical Simulation (DNS) & Large Eddy Simulation (LES) Methods

This study evaluates the capabilities of the NEK5000 open-source code as an alternative to the commercially available software that is currently heavily relied upon. The first three case studies initiated the author with (1) the syntax and flow of the program, (2) varying initial and boundary conditions and (3) geometries of greater complexity & dimension. Case 1 studies transient conjugate heat transfer in the entry region of a pipe. Cases 2 & 3 simulate the natural convection flow in an enclosed differentially heated cavity (DHC) of two and three dimensions, respectively. The final case study, Case 4, produces three-dimensional single-phase natural convective flow using the Boussinesq approximation within a simplified representation of an aqueous homogeneous reactor system. The Boussinesq flow in Case 4 clearly depicts both Rayleigh-Taylor and Kelvin-Helmholtz instabilities qualitatively during post-processing analysis. Further analysis of Case 4, up to the point of a quasi-steady state (i.e. thermal-equilibrium condition), is done using the LES solver method in NEK5000. The volume-averaged solution temperatures, when compared between solvers, showed no significant difference. Due to its minimal compromise in fidelity of solution resolution and inexpensive computing cost, compared to the DNS solver, the LES solver is recommended for the current application. Methodology development for the four case studies are reported and a preliminary optimization study with different polynomial resolutions for Case 4 is addressed as a resource when performing future simulations using NEK5000. Hence, this report can be used as a reference starting point to the end user of NEK5000 for natural circulation phenomena in nuclear reactor safety applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimization using pathwise algorithmic derivatives of electromagnetic shower simulations

Among the well-known methods to approximate derivatives of expectancies computed by Monte-Carlo simulations, averages of pathwise derivatives are often the easiest one to apply. Computing them via algorithmic differentiation typically does not require major manual analysis and rewriting of the code, even for very complex programs like simulations of particle-detector interactions in high-energy physics. However, the pathwise derivative estimator can be biased if there are discontinuities in the program, which may diminish its value for applications. This work integrates algorithmic differentiation into the electromagnetic shower simulation code HepEmShow based on G4HepEm, allowing us to study how well pathwise derivatives approximate derivatives of energy depositions in a sampling calorimeter with respect to parameters of the beam and geometry. We found that when multiple scattering is disabled in the simulation, means of pathwise derivatives converge quickly to their expected values, and these are close to the actual derivatives of the energy deposition. Additionally, we demonstrate the applicability of this novel gradient estimator for stochastic gradient-based optimization in a model example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling of Orbital and Attitude Dynamics of Nanosatellites Controlled via Active Electrostatic Charging

The large-scale exploration of airless bodies, such as asteroids and moons, is gaining interest, however it is limited by mobility issues: the lack of atmosphere, low gravity, and unknown soil properties pose difficult challenges for many forms of traditional locomotion. The environment in proximity of these bodies is also electrically charged due to interactions with solar wind and UV radiation. The EGlider (Electrostatic Glider) concept would be able to overcome these mobility issues by leveraging the natural environment, allowing operations in close proximity of the surface, while enabling long duration missions by minimizing propellant consumption. The EGlider is an advanced concept for small satellite mobility and propulsion, which relies on the electric fields naturally present around airless bodies in order to generate forces and torques useful for maneuvering. It does so by extending electrically charged appendages, which enable it to electrostatically soar above the surface. By differentially charging its electrodes it can also produce torques to control its attitude. The charges are maintained by continuous active ion or electron emission from the spacecraft, in order to cancel out the neutralizing influx of charges from surrounding plasma. An investigation of the spacecraft-plasma interaction was carried out. This included studying the effect of electrode geometry and calculating the charge-to-mass ratios required to enable several mission scenarios. Long, thin wire electrodes were identified to be the most power-efficient and would allow power-to-weight ratios achievable with current nanosatellite technologies. High electrode potential represents the main limiting factor for the system design. In order to test the feasibility of active control by means of differential charging, a simple 2D interaction model was developed, and a feedback controller to stabilize the vehicle was tested in a simulation environment. The results confirmed that good performance could be obtained for both position and attitude control. Finally, a dedicated software was developed for future simulation and testing of control strategies for the EGlider. This software allows to study the trajectory and attitude of an arbitrarily configured spacecraft in the proximity of an arbitrarily defined main airless body. The spacecraft can be assembled from basic parts, each with specified electrical, mass and optical properties. Efficient models allow to calculate gravitational and electrical interactions with the rotating main body and the local plasma field, as well as solar radiation pressure effects. Control models can be implemented as simple plug-in functions and easily tested. The preliminary validation campaign showed good matching with the reference cases that have been analyzed.

Corpino, Sabrina↗

Self-acting geometry for noncontact seals

Two hydrodynamic self-acting seal designs for a LOX turbopump were analyzed in order to predict performance. A radial face seal-to-seal LOX at 310 N/sq cm and 32,000 rpm (130 m/s) was analyzed for pressure differentials of 172 to 448 N/sq cm and speeds from 98 to 147 m/s. A segmented circumferential seal-to-seal helium at 34.5 or 69 N/sq cm and 157 m/s was analyzed for pressures of 35 to 86 N/sq cm (10 N/sq cm ambient) and speeds from 94 to 189 m/s. Test results confirmed noncontact operation near the design speed and pressure, and relatively good qualitative agreement between test and analysis. The face seal was found to operate with mostly liquid in the pads and mostly gas across the dam.

Allen, G. P.↗

A Consistent Spatially Adaptive Smoothed Particle Hydrodynamics Method for fluid-structure interactions

A new consistent, spatially adaptive, smoothed particle hydrodynamics (SPH) method for Fluid-Structure Interactions (FSI) is presented. The method combines several attributes that have not been simultaneously satisfied by other SPH methods. Specifically, it is second-order convergent; it allows for resolutions spatially adapted with moving (translating and rotating) boundaries of arbitrary geometries; and, it accelerates the FSI solution as the adaptive approach leads to fewer degrees of freedom without sacrificing accuracy. The key ingredients in the method are a consistent discretization of differential operators, a posteriori error estimator/distance-based criterion of adaptivity, and a particle-shifting technique. The method is applied in simulating six different flows or FSI problems. The new method's convergence, accuracy, and efficiency attributes are assessed by comparing the results it produces with analytical, finite element, and consistent SPH uniform high-resolution solutions as well as experimental data.

Hu, Wei↗

Picosecond magneto-optic thermometry measurements of nanoscale thermal transport in AlN thin films

The thermal conductivity Λ of wide bandgap semiconductor thin films, such as AlN, affects the performance of high-frequency devices, power devices, and optoelectronics. However, accurate measurements of Λ in thin films with sub-micrometer thicknesses and Λ > 100 W m -1 K -1 is challenging. Widely used pump/probe metrologies, such as time–domain thermoreflectance (TDTR) and frequency–domain thermoreflectance, lack the spatiotemporal resolution necessary to accurately quantify thermal properties of sub-micrometer thin films with high Λ. In this work, we use a combination of magneto-optic thermometry and TiN interfacial layers to significantly enhance the spatiotemporal resolution of pump/probe thermal transport measurements. We use our approach to measure Λ of 100, 400, and 1000 nm AlN thin films. We coat AlN thin films with a ferromagnetic thin-film transducer with the geometry of (1 nm-Pt/0.4 nm-Co) x3 /(2 nm-TiN). This PtCo/TiN transducer has a fast thermal response time of <50 ps, which allows us to differentiate between the thermal response of the transducer, AlN thin film, and substrate. For the 100, 400, and 1000 nm thick AlN films, we determine Λ to be 200 ± 80, 165 ± 35, and 300 ± 70 W m -1 K -1 , respectively. We conclude with an uncertainty analysis that quantifies the errors associated with pump/probe measurements of thermal conductivity, as a function of transducer type, thin-film thermal conductivity, and thin-film thickness. Time resolved magneto-optic Kerr effect experiments can measure films that are three to five times thinner than is possible with standard pump/probe metrologies, such as TDTR. This advance in metrology will enable better characterization of nanoscale heat transfer in high thermal conductivity material systems like wide bandgap semiconductor heterostructures and devices.

36 MATERIALS SCIENCE↗

The analysis of high spatial resolution UV and X-ray images by computational modeling

Very high resolution stereoscopic images of high temperature loop structures observed at UV and X-ray wavelengths in the solar corona can be used to understand physical processes in the corona. An existing computational model is described and sample results are given to demonstrate that computational modeling of coronal structures can indeed take advantage of very high resolution images. The sample results include the run of temperature and number density along a typical loop and the variation of the differential emission measure with temperature. The integration of the differential emission measure with temperature along a column commensurate with an instrument's spatial resolution is the relevant parameter obtained from UV and X-ray observations. The effects of loop geometry and energy input are examined.

Vesecky, J. F.↗

Assessment of Turbulence Models for a Single-Injector Cooling Flow

This work assesses the ability of Reynolds-averaged Navier-Stokes turbulence models to predict a fundamental thermal cooling flow configuration. This geometry consists of a plate with a single cooling hole positioned downstream of a heated exhaust nozzle. The linear Shear Stress Transport (SST), non-linear Quadratic Constitutive Relation (QCR), and differential SSG/LRR Reynolds Stress Model (RSM) turbulence models were investigated for flow conditions with a blowing ratio of 1.0 and 2.0. Despite significant differences in turbulence model complexity, the results miss important details observed in the data and large eddy simulation (LES). Injection penetration is over-predicted. Mixing and centerline wake recovery is too slow. Temperature predictions are worse than velocity, suggesting that perhaps the simple gradient diffusion approximation used to model the turbulent heat flux may be to blame. These observations are put into context using other studies to identify shortcomings in the modeling of both the turbulent stress and turbulent heat flux.

Film Cooling↗

Assessment of Turbulence Models for a Single-Injector Cooling Flow

This work assesses the ability of Reynolds-averaged Navier-Stokes turbulence models to predict a fundamental thermal cooling flow configuration. This geometry consists of a plate with a single cooling hole positioned downstream of a heated exhaust nozzle. The linear Shear Stress Transport (SST), non-linear Quadratic Constitutive Relation (QCR), and differential SSG/LRR Reynolds Stress Model (RSM) turbulence models were investigated for flow conditions with a blowing ratio of 1.0 and 2.0. Despite significant differences in turbulence model complexity, the results miss important details observed in the data and large eddy simulation (LES). Injection penetration is over-predicted. Mixing and centerline wake recovery is too slow. Temperature predictions are worse than velocity, suggesting that perhaps the simple gradient diffusion approximation used to model the turbulent heat flux may be to blame. These observations are put into context using other studies to identify shortcomings in the modeling of both the turbulent stress and turbulent heat flux.

Film Cooling↗

Hysteresis in amplitudes of self-excited oscillations for co-axial electrode-geometry DC glow discharge plasma

Hysteresis in amplitudes of the self-excited oscillations of the floating potential and discharge current is observed in an unmagnetized co-axial electrode-geometry DC glow discharge plasma system. The nonlinearities of these oscillations are studied using standard dynamical analysis tools. The characterization revealed the transition of low-amplitude high-frequency period-n oscillations to a large amplitude low frequency period-1 oscillations through a chaotic intermediate route. The transition of the low amplitude, high frequency period-n oscillations to chaotic type is observed to be linked to the dynamical change in the plasma system, i.e., after a negative differential resistance (NDR) region, whereas the transition from chaotic to period-1 is observed to be linked to a discharge current threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Strongly nonlinear wave propagation in elasto-plastic metamaterials: Low-order dynamic modeling

Nonlinear elastic metamaterials are known to support a variety of dynamic phenomena that enhance our capacity to manipulate elastic waves. Since these properties stem from complex, subwavelength geometry, full-scale dynamic simulations are often prohibitively expensive at scales of interest. Prior studies have therefore utilized low-order effective medium models, such as discrete mass-spring lattices, to capture essential properties in the long-wavelength limit. While models of this type have been successfully implemented for a wide variety of nonlinear elastic systems, they have predominantly considered dynamics depending only on the instantaneous kinematics of the lattice, neglecting history-dependent effects, such as wear and plasticity. Here, to address this limitation, the present study develops a lattice-based modeling framework for nonlinear elastic metamaterials undergoing plastic deformation. Due to the history- and rate-dependent nature of plasticity, the framework generally yields a system of differential-algebraic equations whose computational cost is significantly greater than an elastic system of comparable size. We demonstrate the method using several models inspired by classical lattice dynamics and continuum plasticity theory and explore means to obtain empirical plasticity models for general geometries, thereby gaining insight into the influence of microstructural plasticity on effective material performance, which can be used to improve the design of nonlinear mechanical metamaterials.

Dynamic simulation↗