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 451 records · Page 25

Reversed cowl flap inlet thrust augmentor

An adjustable airfoil is described for varying the geometry of a jet inlet and an ejector inlet in a jet engine for providing thrust augmentation and noise reduction. The airfoil comprises essentially a plurality of segments which are extended radially outward and retracted relative to the longitudinal axis of the engine as a function of a change in the pressure differential between the upstream and downstream surfaces of the airfoil. A servo mechanism responsive to the change in the pressure differential is coupled to the airfoil to extend and retract the airfoil segments to maintain the pressure at a maximum on the downstream side of the airfoil relative to the pressure on the upstream side of the airfoil. At low speeds, such as at take-offs and landings, the airfoil is fully extended while at high speeds it is fully retracted.

Cheng, D. Y.↗

Detection, monitoring, and determination of location of changes in metallic structures using multimode acoustic signals

Acoustic transducers generate and receive acoustic signals at multiple locations along a surface of rigid structure, wherein longitudinal spacing between transducer locations define measurement zones. Acoustic signals with chosen amplitude-time-frequency characteristics excite multiple vibration modes in the structure within each zone. Small mechanical changes in inspection zones lead to scattering and attenuation of broadband acoustic signals, which are detectable as changes in received signal characteristics as part of a through-transmission technique. Additional use of short, narrowband pulse acoustic signals as part of a pulse-echo technique allows determination of the relative location of the mechanical change within each zone based on the differential delay profiles. For accurate acoustic modeling and simulation, the mesh size, time step, time delay, and time-window size are optimized. Frequency normalization of the Short-Time Fourier Transform of acoustic response output improves experiment-simulation cross-validation. Applications of the method to structures with arbitrarily complex geometries are also demonstrated.

Findikoglu, Alp Tugrul↗

BISON Robustness and Performance Improvements

BISON is a modern finite-element based nuclear fuel performance code that has been under development at the Idaho National Laboratory (USA) since 2009 [1]. The code is applicable to both steady and transient fuel behavior and can be used to analyze 1D (spherically symmetric), 2D (axisymmetric and generalized plane strain) or 3D geometries. BISON is the fuel performance code used within CASL for LWR fuel under both normal operating and accident conditions. BISON is built using the INL Multiphysics ObjectOriented Simulation Environment, or MOOSE [2, 3]. MOOSE is a massively parallel, finite element-based framework to solve systems of coupled non-linear partial differential equations using the Jacobian-Free Newton Krylov (JFNK) method [4]. This enables investigation of computationally large problems, for example a full stack of discrete pellets in a LWR fuel rod, or every rod in a full reactor core. MOOSE supports the use of complex two and three-dimensional meshes and uses implicit time integration, important for the widely varied time scale in nuclear fuel simulation. An object-oriented architecture is employed which greatly minimizes the programming effort required to add new material and behavioral models. The flexibility of the implicit and fully coupled multiphysics approach comes with a need for constructing suitable approximations for the Jacobian matrix of the coupled system used for either preconditioning a Krylov solve or in a direct Newton solve. Preconditioning options for Bison problems need to be revisited with new preconditioning methods becoming available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quantum Gravity and Laser Interferometry: Towards Observable Predictions

Understanding quantum gravity remains one of the deepest challenges in modern physics, as direct experimental access to Planck-scale effects is beyond current technological reach. However, recent theoretical advances indicate that quantum fluctuations of spacetime may produce measurable effects in precision experiments, particularly near causal horizons. This opens new avenues for testing quantum gravity phenomena through high-precision measurement techniques. This dissertation develops multiple theoretical models to characterize these effects and examines their potential observational signatures in future gravitational wave interferometers. We begin by investigating the role of quantum fluctuations in near-horizon geometries through the lens of the AdS/CFT correspondence, which provides a powerful framework for understanding the interplay between quantum field theory and general relativity via holographic principles. By modeling stochastic energy-momentum sources in Rindler-AdS spacetime, we demonstrate that vacuum fluctuations transform the Einstein equations into a Langevin-type stochastic differential equation, leading to potentially observable fluctuations in photon traversal times. Extending this approach to Minkowski spacetime, we establish a correspondence between gravitational shockwaves and fluid dynamics, showing that near-horizon perturbations satisfy an equation analogous to that governing incompressible fluids, thereby reinforcing the membrane paradigm and hydrodynamic analogies in the context of the fluid/gravity duality. Furthermore, we construct the covariant phase space of a spherically symmetric causal diamond in Minkowski spacetime, identifying two fundamental charges that govern its evolution. These results provide a foundation for quantizing causal horizons and understanding their microscopic degrees of freedom. Building upon these theoretical developments, we further examine a related stochastic phenomenon: the gravitational wave memory background arising from the cumulative memory steps produced by supermassive black hole mergers. After reviewing the standard stochastic gravitational wave background, gravitational memory effects, and BMS symmetries, we model the stochastic memory background using a Brownian motion framework. We show that while the cumulative memory background initially appears above the sensitivity curve of space-based interferometers like LISA, the realistic subtraction of individually resolvable merger events substantially suppresses the residual signal, making its detection more challenging. This highlights the critical importance of source subtraction when evaluating the detectability of gravitational memory effects. By bridging fundamental theory with experimental prospects, this dissertation contributes to the ongoing effort to uncover the quantum nature of spacetime through precision measurement techniques. Whether through detecting quantum spacetime fluctuations, gravitational memory backgrounds, or probing the symmetries of causal horizons, the pursuit of observable quantum gravity phenomena continues to expand the frontiers of both theory and experiment.

Zhang, Yiwen [Caltech] (ORCID:0000000323559416)↗

The discrete Green's function paradigm for two-way coupled Euler-Lagrange simulation

We outline a methodology for the simulation of particle-laden flows whereby the dispersed and fluid phases are two-way coupled. The drag force which couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. The method is general in that it can be extended to other partial differential equations which may be associated with particle-laden flows, such as the thermal energy equation or Maxwell's equations. In this work, we demonstrate the method of discrete Green's functions by obtaining these functions for the Navier-Stokes equations at low particle Reynolds number in a two-plane channel geometry. We perform verification at low and finite particle Reynolds number for the case of a point-particle settling under gravity parallel to a plane wall, for different wall normal separations. In comparing to other point-particle schemes the discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations and is the most accurate in the near wall region at finite Reynolds number. We discuss how the accuracy away from the wall at finite Reynolds number could be improved by appealing to Oseen-like discrete Green's functions. Finally we demonstrate that the discrete Green's function approach can have important implications on statistics of particle-laden turbulent channel flow.

42 ENGINEERING↗

The Method of Finite Averages

The Method of Finite Averages (MoFA) is a rigorous multiscale modeling methodology for efficiently modeling multi-physical phenomena in heterogeneous porous media. The code developed in this project aims to perform the numerical calculations required to formulate, implement, and verify MoFA models for Earth and Energy systems (i.e., model verification refers to performing fully-resolved simulations of the systems and comparing their results to those of the models). In general, MoFA transforms partial differential equations (PDEs) describing the fine-scale physics of a system into coupled ordinary differential equations (ODEs)---in time---that describe the coarse-scale---or "average"---physical behaviors of the system. This transformation significantly expedites system simulation, as the coarse-scale ODEs involve vastly fewer degrees of freedom than the fine-scale PDEs. The code developed under this project will allow users to 1.) generate system geometries and numerical meshes, 2.) solve the PDE and ODE systems required for MoFA model formulation and implementation, 3.) solve the PDE systems required to obtain fully-resolved simulation results for model verification, and 4.) compare and plot results (e.g., the model and fully-resolved simulation solutions, the error between the solutions, etc.).

Pietrzyk, KyleM [Lawrence Livermore National Labor↗

Judging rolling wheels: Dynamic and kinematic aspects of rotation-translation coupling

Four experiments were carried out to investigate observers' abilities to judge rolling motions. The experiments were designed to assess whether two important aspects of such motions are appreciated: the kinematic coupling of rotation and translation, and the dynamic effects of gravity. Different motion contexts of rolling wheels were created using computer-generated displays. The first experiment involved wheels rolling down an inclined plane. Observers spontaneously appreciated the anomaly of wheels that failed to accelerate, but they were not able to differentiate between different acceleration functions. Moreover, their judgements were almost exclusively based on the translation component of the rolling motion, neglecting the rotation component. In a second experiment it was found that observers could accurately estimate the perimeter of various objects. Thus, their inability to consider rotation information is not attributable to misperceptions of the geometry of wheels. In a third experiment the finding that rolling wheels appear to overrotate was replicated; however, findings from this experiment also showed, together with those from a fourth experiment, that observers are able to make very accurate judgments about translation-rotation coupling in rolling wheels when information is provided about the orientation of the wheel and the texture of the surface on which it rolls.

Hecht, Heiko↗

Simulation of Ultra-Small MOSFETs Using a 2-D Quantum-Corrected Drift-Diffusion Model

The continued down-scaling of electronic devices, in particular the commercially dominant MOSFET, will force a fundamental change in the process of new electronics technology development in the next five to ten years. The cost of developing new technology generations is soaring along with the price of new fabrication facilities, even as competitive pressure intensifies to bring this new technology to market faster than ever before. To reduce cost and time to market, device simulation must become a more fundamental, indeed dominant, part of the technology development cycle. In order to produce these benefits, simulation accuracy must improve markedly. At the same time, device physics will become more complex, with the rapid increase in various small-geometry and quantum effects. This work describes both an approach to device simulator development and a physical model which advance the effort to meet the tremendous electronic device simulation challenge described above. The device simulation approach is to specify the physical model at a high level to a general-purpose (but highly efficient) partial differential equation solver (in this case PROPHET, developed by Lucent Technologies), which then simulates the model in 1-D, 2-D, or 3-D for a specified device and test regime. This approach allows for the rapid investigation of a wide range of device models and effects, which is certainly essential for device simulation to catch up with, and then stay ahead of, electronic device technology of the present and future. The physical device model used in this work is the density-gradient (DG) quantum correction to the drift-diffusion model [Ancona, Phys. Rev. B 35(5), 7959 (1987)]. This model adds tunneling and quantum smoothing of carrier density profiles to the drift-diffusion model. We used the DG model in 1-D and 2-D (for the first time) to simulate both bipolar and unipolar devices. Simulations of heavily-doped, short-base diodes indicated that the DG quantum corrections do not have a large effect on the IN characteristics of electronic devices without heteroj unction s. On the other hand, ultra-small MOSFETs certainly exhibit important quantum effects that the DG model will include: quantum repulsion of the inversion and gate charges from the oxide interfaces, and quantum tunneling through thin gate oxides. We present initial results of 2-D DG simulations of ultra-small MOSFETs. Subtle but important issues involving the specification of the model, boundary conditions, and interface constraints for DG simulation of MOSFETs will also be illuminated.

Biegal, Bryan A.↗

Synergistic learning with multi-task DeepONet for efficient PDE problem solving

Multi-task learning (MTL) is an inductive transfer mechanism designed to leverage useful information from multiple tasks to improve generalization performance compared to single-task learning. It has been extensively explored in traditional machine learning to address issues such as data sparsity and overfitting in neural networks. In this work, we apply MTL to problems in science and engineering governed by partial differential equations (PDEs). However, implementing MTL in this context is complex, as it requires task-specific modifications to accommodate various scenarios representing different physical processes. To this end, we present a multi-task deep operator network (MT-DeepONet) to learn solutions across various functional forms of source terms in a PDE and multiple geometries in a single concurrent training session. We introduce modifications in the branch network of the vanilla DeepONet to account for various functional forms of a parameterized coefficient in a PDE. Additionally, we handle parameterized geometries by introducing a binary mask in the branch network and incorporating it into the loss term to improve convergence and generalization to new geometry tasks. Our approach is demonstrated on three benchmark problems: (1) learning different functional forms of the source term in the Fisher equation; (2) learning multiple geometries in a 2D Darcy Flow problem and showcasing better transfer learning capabilities to new geometries; and (3) learning 3D parameterized geometries for a heat transfer problem and demonstrate the ability to predict on new but similar geometries. Finally, our MT-DeepONet framework offers a novel approach to solving PDE problems in engineering and science under a unified umbrella based on synergistic learning that reduces the overall training cost for neural operators.

42 ENGINEERING↗

Comparison of Computed and Measured Performance of a Pulsed Inductive Thruster Operating on Argon Propellant

Pulsed inductive plasma accelerators are electrodeless space propulsion devices where a capacitor is charged to an initial voltage and then discharged through a coil as a high-current pulse that inductively couples energy into the propellant. The field produced by this pulse ionizes the propellant, producing a plasma near the face of the coil. Once a plasma is formed if can be accelerated and expelled at a high exhaust velocity by the Lorentz force arising from the interaction of an induced plasma current and the magnetic field. A recent review of the developmental history of planar-geometry pulsed inductive thrusters, where the coil take the shape of a flat spiral, can be found in Ref. [1]. Two concepts that have employed this geometry are the Pulsed Inductive Thruster (PIT)[2, 3] and the Faraday Accelerator with Radio-frequency Assisted Discharge (FARAD)[4]. There exists a 1-D pulsed inductive acceleration model that employs a set of circuit equations coupled to a one-dimensional momentum equation. The model was originally developed and used by Lovberg and Dailey[2, 3] and has since been nondimensionalized and used by Polzin et al.[5, 6] to define a set of scaling parameters and gain general insight into their effect on thruster performance. The circuit presented in Fig. 1 provides a description of the electrical coupling between the current flowing in the thruster I1 and the plasma current I2. Recently, the model was upgraded to include an equation governing the deposition of energy into various modes present in a pulsed inductive thruster system (acceleration, magnetic flux generation, resistive heating, etc.)[7]. An MHD description of the plasma energy density evolution was tailored to the thruster geometry by assuming only one-dimensional motion and averaging the plasma properties over the spatial dimensions of the current sheet to obtain an equation for the time-evolution of the total energy. The equation set governing the dynamics of the coupled electrodynamic-current sheet system is composed of first-order, coupled ordinary differential equations that can be easily solved numerically without having to resort to much more complex 2-D finite element plasma simulations.

Polzin, Kurt A.↗

The Distribution of Methane in Titan's Stratosphere from Cassini/CIRS Observations

Cassini/CIRS spectra in the far- and mid-infrared region are used to determine the abundance of methanein Titans lower stratosphere and investigate its distribution with latitude. The CIRS spectra include emissionfrom both the CH4 v4 band at 7.7 micron and pure rotational lines longwards of 50 micron, which show differential sensitivities to thermal profile and methane mole fraction. We analyze nadir and limb data taken over the first part of the Cassini mission (August 2005 to June 2010), including a selection of 12 latitudes that provides a reasonably complete and regular sampling of both hemispheres. Unexpectedly, but in a consistent manner for limb and nadir geometries, large variations of the methane mole fraction near 15 mbar (approx. 85 km) are found, with values ranging from approx. 1.0% (at low latitudes and near +/-50-55 deg) to approx. 1.5% (at +/-30-35 deg and polar latitudes). Error bars on the retrieved methane mole fraction are 0.07-0.12% at low latitudes in the Southern hemisphere and 0.14-0.21% northward of 40 deg N. A 1.0% methanemole fraction at low latitudes permits us to reconcile the HASI-measured temperatures below 147 kmaltitude (2.7 mbar) with inferences from CIRS. The roughly hemispherically-symmetric distribution ofmethane gas is reminiscent of that observed or predicted for the tropospheric methane clouds, whichon a yearly-averaged basis, show preferential occurrences at tropical and polar latitudes. We speculatethat convective events at these latitudes result into local stratospheric methane enrichment, whichmay persist year-round due to dynamical mixing times in the lower stratosphere only moderately shorterthan a Titan year.

Methane↗

Generalized fiducial inference on differentiable manifolds

We introduce a novel approach to inference on parameters that take values in a Riemannian manifold embedded in a Euclidean space. Parameter spaces of this form are ubiquitous across many fields, including chemistry, physics, computer graphics, and geology. Here, this new approach uses generalized fiducial inference (GFI) to obtain a posterior-like distribution on the manifold, without needing to know local parameterizations that map to the constrained space from an unconstrained Euclidean space. Using mathematical tools from Riemannian geometry, we construct a constrained generalized fiducial distribution (CGFD). A Bernstein-von Mises-type result for the CGFD, which provides intuition for how the desirable asymptotic qualities of the unconstrained generalized fiducial distribution are inherited by the CGFD, is provided. To illustrate the practical use of the CGFD, we provide a proof-of-concept example in the context of a linear logspline density estimation problem, and demonstrate that CGFD-based confidence sets exhibit desirable coverage properties via simulation. As an application, we fit a CGFD to COVID-19 case count data from North Carolina, USA.

97 MATHEMATICS AND COMPUTING↗

Analysis of a finite composite plate with smooth rigid pin

An elastic continuum method of analysis developed for a smooth rigid pin in a finite composite plate by using complex potentials that satisfy the governing orthotropic differential equation is presented. The plate is subjected to uniaxial loading and the nonlinear problems of partial contact in the pin-plate are studied. Three types of pin fits are examined: interference, push, and clearance. The problems of stresses due to pure interference, load contact behavior, and stresses variations are addressed. The effects of the rectangular geometry and the orthotropic property of the plate are demonstrated through numerical results. It is found that the use of large hole diameter/plate width ratios have deleterious effects and the use of interference has beneficial effects. The solutions obtained will help in the development of special finite elements for fastener joints. Finally, the boundary conditions on the finite outer boundaries are found to be satisfied by a successive integration technique.

Naidu, A. C. B.↗

Resonant torque differential magnetometry with high frequency quartz oscillators

Sensitive magnetometry has been a powerful probe for investigating quantum materials. Extreme conditions, such as sub-kelvin cryogenic temperatures and ultrahigh magnetic fields, demand further durability for sensitive magnetometry. However, significant mechanical vibrations and rapid magnetic field changes give enormous challenges to conventional magnetometry. Here, this article presents a possible solution to this problem by developing a new magnetometry technique using high-frequency quartz oscillators. The technique takes advantage of the symmetry and geometry of mechanical vibration configurations of standard commercially available MHz quartz oscillators, and the setup keeps the high quality factor resonance with the sample mounted on the oscillator. We further demonstrate the sensitivity of the technique using bismuth single crystals and a Fe 0.25 TaS 2 ferromagnetic material. Quantum oscillations are observed in the magnetometry response below 1 T, and the detected oscillation frequency is shown to come from the electron pockets of the bismuth.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Rindler fluids from gravitational shockwaves

We study a correspondence between gravitational shockwave geometry and its fluid description near a Rindler horizon in Minkowski spacetime. Utilizing the Petrov classification that describes algebraic symmetries for Lorentzian spaces, we establish an explicit mapping between a potential fluid and the shockwave metric perturbation, where the Einstein equation for the shockwave geometry is equivalent to the incompressibility condition of the fluid, augmented by a shockwave source. Then we consider an Ansatz of a stochastic quantum source for the potential fluid, which has the physical interpretation of shockwaves created by vacuum energy fluctuations. Under such circumstance, the Einstein equation, or equivalently, the incompressibility condition for the fluid, becomes a stochastic differential equation. By smearing the quantum source on a stretched horizon in a Lorentz invariant manner with a Planckian width (similarly to the membrane paradigm), we integrate fluctuations near the Rindler horizon to find an accumulated effect of the variance in the round-trip time of a photon traversing the horizon of a causal diamond.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Small-scale production of bespoke accelerated aging plutonium alloy

Here, this study demonstrates the 100 g scale manufacture of a plutonium alloy that ages at an accelerated rate. The resulting alloy ages six times faster than typical weapons-grade plutonium due to the addition of 238 Pu. As a major innovation, the process involved using a partial direct oxide reduction technique. This method was achieved by developing a new, complex geometry stirrer using additive manufacturing to reduce the 238Pu oxide and efficiently incorporate it into weapons-grade plutonium metal. The material was then purified by molten salt extraction and electrorefining before being alloyed with gallium. The alloy was then cold-rolled and annealed in a homogenization heat treatment. The resulting disk was characterized by metallography and differential scanning calorimetry, and the impurity content was determined using analytical chemistry techniques. The results show that a homogeneous delta phase plutonium alloy was achieved with expected microstructure and minimal impurities. This study was also successful in changing the plutonium isotopic composition by incorporating additional 238 Pu to accelerate the effects of radiation damage. This enables researchers to study long-term aging phenomena in a reduced time frame, thus avoiding the need for large-scale material production and circumventing the limitations of using naturally aged, archived plutonium.

Aging theory↗

Accelerating FEM-Based Corrosion Predictions Using Machine Learning

Atmospheric corrosion of metallic parts is a widespread materials degradation phenomena that is challenging to predict given its dependence on many factors (e.g. environmental, physiochemical, and part geometry). For materials with long expected service lives, accurately predicting the degree to which corrosion will degrade part performance is especially difficult due to the stochastic nature of corrosion damage spread across years or decades of service. The Finite Element Method (FEM) is a computational technique capable of providing accurate estimates of corrosion rate by numerically solving complex differential Eqs. characterizing this phenomena. Nevertheless, given the iterative nature of FEM and the computational expense required to solve these complex equations, FEM is ill-equipped for an efficient exploration of the design space to identify factors that accelerate or deter corrosion, despite its accuracy. In this work, a machine learning based surrogate model capable of providing accurate predictions of corrosion with significant computational savings is introduced. Specifically, this work leverages AdaBoosted Decision trees to provide an accurate estimate of corrosion current per width given different values of temperature, water layer thickness, molarity of the solution, and the length of the cathode for a galvanic couple of aluminum and stainless steel.

36 MATERIALS SCIENCE↗

The CMC:3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 1: Theoretical

An order-of-magnitude analysis of the subsonic three dimensional steady time averaged Navier-Stokes equations, for semibounded aerodynamic juncture geometries, yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement of first order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. Numerical results document accuracy, convergence, and utility of the developed finite element algorithm, and the CMC:3DPNS computer code applied to an idealized wing-body juncture region. Additional results document accuracy aspects of the algorithm turbulence closure model.

Baker, A. J.↗