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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2021

CARDINAL

Cardinal is a code that automates the workflow for iteratively solving neutral particle transport via the OpenMC code and computational fluid dynamics via the NekRS code within the Multiphysics Object Oriented Simulation Environment (MOOSE), a finite element framework for solving general partial differential equations. Cardinal establishes a mapping between the Monte Carlo geometry representation used in OpenMC and the spectral element mesh used in NekRS to other physics applications built upon the MOOSE framework, allowing high-resolution particle transport and fluid dynamics physics feedback to the MOOSE application ecosystem. Heat generation rates as solved by OpenMC are used as heat sources in MOOSE, while the temperature and density fields from NekRS are used as conjugate heat transfer boundary conditions and/or source fields in MOOSE. OpenMC, Nek5000, and MOOSE are open source, community-developed simulation codes.

NOVAK, APRILJEAN RODGERS↗

DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

NASA NTRS2002

Dynamics of Diffusion Flames in von Karman Swirling Flows Studied

Von Karman swirling flow is generated by the viscous pumping action of a solid disk spinning in a quiescent fluid media. When this spinning disk is ignited in an oxidizing environment, a flat diffusion flame is established adjacent to the disk, embedded in the boundary layer (see the preceding illustration). For this geometry, the conservation equations reduce to a system of ordinary differential equations, enabling researchers to carry out detailed theoretical models to study the effects of varying strain on the dynamics of diffusion flames. Experimentally, the spinning disk burner provides an ideal configuration to precisely control the strain rates over a wide range. Our original motivation at the NASA Glenn Research Center to study these flames arose from a need to understand the flammability characteristics of solid fuels in microgravity where slow, subbuoyant flows can exist, producing very small strain rates. In a recent work (ref. 1), we showed that the flammability boundaries are wider and the minimum oxygen index (below which flames cannot be sustained) is lower for the von Karman flow configuration in comparison to a stagnation-point flow. Adding a small forced convection to the swirling flow pushes the flame into regions of higher strain and, thereby, decreases the range of flammable strain rates. Experiments using downward facing, polymethylmethacrylate (PMMA) disks spinning in air revealed that, close to the extinction boundaries, the flat diffusion flame breaks up into rotating spiral flames (refs. 2 and 3). Remarkably, the dynamics of these spiral flame edges exhibit a number of similarities to spirals observed in biological systems, such as the electric pulses in cardiac muscles and the aggregation of slime-mold amoeba. The tail of the spiral rotates rigidly while the tip executes a compound, meandering motion sometimes observed in Belousov-Zhabotinskii reactions.

Nayagam, Vedha↗

NASA NTRSDate not supplied

Aviary: An Open-Source Multidisciplinary Design, Analysis, and Optimization Tool for Modeling Aircraft with Analytic Gradients

In recent years demands on aircraft design methods have begun to require higher degrees of coupling between disciplines and optimization in order to satisfy competing objectives involving large numbers of design parameters that define unconventional configurations. These expanding requirements have amplified a need for new and improved aircraft design, analysis, and optimization codes that are capable of performing coupled design and exploiting analytic gradients to perform gradient-based optimization where possible. To address this need, a new multidisciplinary design optimization and analysis tool called Aviary is presented. This tool, built on OpenMDAO, allows for tightly coupled, simultaneous aircraft and subsystem design using analytic gradients. Aviary includes methods from two NASA developed legacy aircraft analysis tools and provides native analytically differentiated calculations for five different disciplines (weights and sizing, aerodynamics, geometry, propulsion, and mission analysis), while also allowing external discipline analysis tools to be coupled, regardless of whether those tools can provide analytic gradients. Verification and preliminary examples and modeling efforts show Aviary’s ability to effectively model novel concepts and explore large and non-intuitive design spaces. Additionally, a multi-level user interface in Aviary creates an easy entry point for users with any level of multidisciplinary design, analysis, and optimization experience.

optimization↗

DOE OSTI2016

Discretization-Accuracy Convergence for Full Algebraic Multigrid

Full multigrid (FMG) is well known for converging to the level of discretization accuracy in a single cycle for a wide class of partial differential equations when the multigrid hierarchy is derived from problem geometry. When applying an FMG cycle to a hierarchy generated by algebraic multigrid (AMG), however, this scalable convergence to discretization accuracy is usually lost. Furthermore, this paper examines the cause of this loss and explores some improvements to standard AMG interpolation that can restore single-cycle convergence to discretization accuracy.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Cu site differentiation in tetracopper(I) sulfide clusters enables biomimetic N 2 O reduction

Copper clusters feature prominently in both metalloenzymes and synthetic nanoclusters that mediate catalytic redox transformations of gaseous small molecules. Such reactions are critical to biological energy conversion and are expected to be crucial parts of renewable energy economies. However, the precise roles of individual metal atoms within clusters are difficult to elucidate, particularly for cluster systems that are dynamic under operating conditions. Here, we present a metal site-specific analysis of synthetic Cu 4 (μ 4 -S) clusters that mimic the Cu Z active site of the nitrous oxide reductase enzyme. Leveraging the ability to obtain structural snapshots of both inactive and active forms of the synthetic model system, we analyzed both states using resonant X-ray diffraction anomalous fine structure (DAFS), a technique that enables X-ray absorption profiles of individual metal sites within a cluster to be extracted independently. Using DAFS, we found that a change in cluster geometry between the inactive and active states is correlated to Cu site differentiation that is presumably required for efficient activation of N 2 O gas. More precisely, we hypothesize that the Cu δ+ ∙∙∙Cu δ- pairs produced upon site differentiation are poised for N 2 O activation, as supported by computational modeling. These results provide an unprecedented level of detail on the roles of individual metal sites within the synthetic cluster system and how those roles interplay with cluster geometry to impact the reactivity function. We expect this fundamental knowledge to inform understanding of metal clusters in settings ranging from (bio)molecular to nanocluster to extended solid systems involved in energy conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Toward Improved Charge Separation through Conformational Control in Copper Coordination Complexes

The continued development of solar energy as a renewable resource necessitates the design and study of new approaches to sustaining photodriven charge separation (CS). To this end, we present a bioinspired approach in which triggered conformational changes are used to control electron transfer (ET) events. We report, photo-induced conformational rearrangements of a ligand are translated into changes in the coordination geometry and environment about a bound metal ion. Taking advantage of the differential coordination properties of Cu I and Cu II , these dynamics facilitate intramolecular ET from Cu I to the ligand to create a CS state. The synthesis and photophysical characterization of CuCl(dpaa R ) (dpaa = dipicolylaminoacetophenone, with R = H and OMe) is presented. These ligands incorporate a fluorophore into their framework that gives rise to a twisted intramolecular charge transfer (TICT) excited state. Excited state ligand twisting provides a tetragonal coordination geometry capable of capturing Cu II in the CS state when an internal ortho-OMe binding site is present (as in dpaaOMe). We employ NMR, IR, EPR, and optical spectroscopies, X-ray diffraction, and electrochemical methods to establish the ground state properties of the Cu I and Cu II complexes. We then investigate the photophysical dynamics of these Cu I complexes via time-resolved photoluminescence (TR-PL), and optical transient absorption (OTA) spectroscopies. We show that relative to controls lacking a TICT-active ligand, the lifetime of the CS state is enhanced ~1000-fold. Further, the presence of the ortho-OMe substituent greatly enhances the lifetime of the TICT* state and biases the coordination environment toward Cu II . The presence of Cu I decreases photoinduced degradation from 14 to <2% but does not result in significant quenching via ET. Factors affecting CS in these systems are discussed, laying the groundwork for our strategy toward solar energy conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Remarks on Hamiltonian structures in G{sub 2}-geometry

In this article, we treat G{sub 2}-geometry as a special case of multisymplectic geometry and make a number of remarks regarding Hamiltonian multivector fields and Hamiltonian differential forms on manifolds with an integrable G{sub 2}-structure; in particular, we discuss existence and make a number of identifications of the spaces of Hamiltonian structures associated to the two multisymplectic structures associated to an integrable G{sub 2}-structure. Along the way, we prove some results in multisymplectic geometry that are generalizations of results from symplectic geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Assessment of MOOSE-Based Tools for Calculating Radial Core Expansion

Radial core expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. An inherently safe reactor design takes advantage of negative reactivity feedback in accident scenarios by utilizing a core restraint system which produces a bowed shape that allows for radial expansion of the fuel regions. Detailed modeling and simulation of core radial expansion itself as well as subsequent reactivity feedback is a challenging task involving contact of many fuel assembly elements and physics feedback from neutronics, thermal hydraulics, and thermal mechanical response. A variety of physics codes have been developed to model aspects of radial core expansion but in general invoke geometrical or physics approximations. No code system currently exists which tightly and robustly couples these physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A high-fidelity code will also be used to benchmark existing lower fidelity, faster running models to understand their benefits, limitations, and range of applications. A code development path using MOOSE-based tools is proposed in order to leverage the detailed geometry capabilities and natural tight coupling and robustness of MOOSE-based applications for modeling this complex phenomena. While simulation of the full phenomenon involves several physics, an assessment has been initiated on the capabilities and readiness of the currently available Tensor Mechanics module within MOOSE for calculation of the structural mechanical responses which occur within the reactor core. This report focuses on modeling the force-deformation response which mimics the physics of a duct contact deformation, as well as differential thermal expansion which produces a thermal bowed shaped for the fuel assemblies in a core. Simple examples were initially performed such as simple supported beam bending under load. The complexity of examples was progressively increased to better mimic the duct behavior by including a differential thermal example and inclusion of hexagonal cross-sections in the geometry. Further assessment of the structural mechanical response simulation capability is still required for modeling duct contact interactions and irradiation creep and swelling. Companion thermal hydraulic and neutronics assessments will also be required; these activities are planned for future years. Finally, integration of the multiple physics components through MOOSE is required to predict the core radial expansion and subsequent feedback effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2014

The Effect of Exogenous Zinc Concentration on the Responsiveness of MC3T3-E1 Pre-Osteoblasts to Surface Microtopography: Part II (Differentiation)

Osseointegration of bone implants is a vital part of the recovery process. Numerous studies have shown that micropatterned geometries can promote cell-substrate associations and strengthen the bond between tissue and the implanted material. As demonstrated previously, exogenous zinc levels can influence the responsiveness of pre-osteoblasts to micropatterns and modify their migratory behavior. In this study, we sought to determine the effect of exogenous zinc on differentiation of osteoblasts cultured on micropatterned vs. planar substrates. Levels of activated metalloproteinase-2 (MMP-2) and transforming growth factor-beta 1 (TGF-β1), as well as early stage differentiation marker alkaline phosphatase, were altered with the addition of zinc. These results suggest that exogenous zinc concentration and micropatterning may interdependently modulate osteoblast differentiation.

metalloproteinase-2↗

DOE OSTI2014

A unified radiative magnetohydrodynamics code for lightning-like discharge simulations

A two-dimensional Eulerian finite difference code is developed for solving the non-ideal magnetohydrodynamic (MHD) equations including the effects of self-consistent magnetic field, thermal conduction, resistivity, gravity, and radiation transfer, which when combined with specified pulse current models and plasma equations of state, can be used as a unified lightning return stroke solver. The differential equations are written in the covariant form in the cylindrical geometry and kept in the conservative form which enables some high-accuracy shock capturing schemes to be equipped in the lightning channel configuration naturally. In this code, the 5-order weighted essentially non-oscillatory scheme combined with Lax-Friedrichs flux splitting method is introduced for computing the convection terms of the MHD equations. The 3-order total variation diminishing Runge-Kutta integral operator is also equipped to keep the time-space accuracy of consistency. The numerical algorithms for non-ideal terms, e.g., artificial viscosity, resistivity, and thermal conduction, are introduced in the code via operator splitting method. This code assumes the radiation is in local thermodynamic equilibrium with plasma components and the flux limited diffusion algorithm with grey opacities is implemented for computing the radiation transfer. The transport coefficients and equation of state in this code are obtained from detailed particle population distribution calculation, which makes the numerical model is self-consistent. This code is systematically validated via the Sedov blast solutions and then used for lightning return stroke simulations with the peak current being 20 kA, 30 kA, and 40 kA, respectively. The results show that this numerical model consistent with observations and previous numerical results. The population distribution evolution and energy conservation problems are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

A physically based connection between fractional calculus and fractal geometry

We show a relation between fractional calculus and fractals, based only on physical and geometrical considerations. The link has been found in the physical origins of the power-laws, ruling the evolution of many natural phenomena, whose long memory and hereditary properties are mathematically modelled by differential operators of non integer order. Dealing with the relevant example of a viscous fluid seeping through a fractal shaped porous medium, we show that, once a physical phenomenon or process takes place on an underlying fractal geometry, then a power-law naturally comes up in ruling its evolution, whose order is related to the anomalous dimension of such geometry, as well as to the model used to describe the physics involved. By linearizing the non linear dependence of the response of the system at hand to a proper forcing action then, exploiting the Boltzmann superposition principle, a fractional differential equation is found, describing the dynamics of the system itself. The order of such equation is again related to the anomalous dimension of the underlying geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Molecular Symmetry-Mixed Dichroism in Double Photoionization of $H 2 $

Dichroism in double photoionization of $H_2$ molecules by elliptically polarized extreme ultraviolet pulses is formulated analytically as a sum of atomiclike dichroism (AD) and molecular symmetry-mixed dichroism (MSMD) terms. The MSMD originates from an interplay of $^1Σ$$^+\atop{u}$ and $^1Π$$^+\atop{u}$ continuum molecular ionization amplitudes. For detection geometries in which the AD vanishes, numerical results for the sixfold differential probabilities for opposite pulse helicities show that the MSMD is significant in the electron momentum and angular distributions and is controllable by the ellipticity.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2024

MatCal Users Guide: Release 1.3.0

Any continuum mechanics model will require three components: (1) a discretized geometry of the boundary value problem being studied, (2) the partial differential equations to be solved, and (3) the initial conditions and boundary conditions for the problem. To describe material behavior in these computational models, material models contribute to (2) the underlying equations and, occasionally, to (3) the initial conditions for the simulation. These material models can exhibit a mathematical form that is empirically based, based on first principles, or developed from both empirical observations and known physics. In general, these models are meant to represent a class of materials with well understood behavior. As a result, material models have parameters that must be tuned or calibrated so that the model response matches characterization data available for the specific material it is intended to represent when used to simulate a specific system. For simple models, such as isotropic, linear elastic materials in solid mechanics, this calibration process can be a simple analytical calculation directly extracting the parameters from experimental measurements. For complex models that have many inputs and require many characterization datasets to adequately identify the material behavior, the model calibration process can require an inverse problem approach where an optimization is performed to tune the model parameters to the available data.

36 MATERIALS SCIENCE↗

DOE OSTI2016

KELT-17B: A HOT-JUPITER TRANSITING AN A-STAR IN A MISALIGNED ORBIT DETECTED WITH DOPPLER TOMOGRAPHY

We present the discovery of a hot Jupiter transiting the V = 9.23 mag main-sequence A-star KELT-17 (BD+14 1881). KELT-17b is a 1.31{sub −0.29}{sup +0.28} M{sub J}, 1.525{sub −0.060}{sup +0.065} R{sub J} hot-Jupiter in a 3.08-day period orbit misaligned at −115.°9 ± 4.°1 to the rotation axis of the star. The planet is confirmed via both the detection of the radial velocity orbit, and the Doppler tomographic detection of the shadow of the planet during two transits. The nature of the spin–orbit misaligned transit geometry allows us to place a constraint on the level of differential rotation in the host star; we find that KELT-17 is consistent with both rigid-body rotation and solar differential rotation rates (α<0.30 at 2σ significance). KELT-17 is only the fourth A-star with a confirmed transiting planet, and with a mass of 1.635{sub −0.061}{sup +0.066} M{sub ⊙}, an effective temperature of 7454 ± 49 K, and a projected rotational velocity of vsinI{sub ∗}=44.2{sub −1.3}{sup +1.5} km s{sup −1}; it is among the most massive, hottest, and most rapidly rotating of known planet hosts.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1992

Three-dimensional co-culture process

The present invention relates to a 3-dimensional co-culture process, more particularly to methods or co-culturing at least two types of cells in a culture environment, either in space or in unit gravity, with minimum shear stress, freedom for 3-dimensional spatial orientation of the suspended particles and localization of particles with differing or similar sedimentation properties in a similar spatial region to form 3-dimensional tissue-like structures. Several examples of multicellular 3-dimensional experiences are included. The protocol and procedure are also set forth. The process allows simultaneous culture of multiple cell types and supporting substrates in a manner which does not disrupt the 3-dimensional spatial orientation of these components. The co-cultured cells cause a mutual induction effect which mimics the natural hormonal signals and cell interactions found in the intact organism. This causes the tissues to differentiate and form higher 3-dimensional structures such as glands, junctional complexes polypoid geometries, and microvilli which represent the corresponding in-vitro structures to a greater degree than when the cell types are cultured individually or by conventional processes. This process was clearly demonstrated for the case of two epithelial derived colon cancer lines, each co-cultured with normal human fibroblasts and with microcarrier bead substrates. The results clearly demonstrate increased 3-dimensional tissue-like structure and biochemical evidence of an increased differentiation state. With the present invention a variety of cells may be co-cultured to produce tissue which has 3-dimensionality and has some of the characteristics of in-vitro tissue. The process provides enhanced 3-dimensional tissue which create a multicellular organoid differentiation model.

Wolf, David A.↗

DOE OSTI2026

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Spatiotemporal dynamics of Gaussian laser pulse in a multi ions plasma

Spatiotemporal evolutions of Gaussian laser pulse propagating through a plasma with multiple charged ions are studied, taking into account the ponderomotive nonlinearity. Coupled differential equations for beam width and pulse length parameters are established and numerically solved using paraxial ray approximation. In one-dimensional geometry, effects of laser and plasma parameters such as laser intensity, plasma density, and temperature on the longitudinal pulse compression and the laser intensity distribution are analyzed for plasmas with singly and doubly charged ions. The results demonstrate that self-compression occurs in a laser intensity range with a turning point intensity in which the self-compression process has its strongest extent. The results also show that the multiply ionized ions have different effect on the pulse compression above and below turning point intensity. Finally, three-dimensional geometry is used to analyze the simultaneous evolution of both self-focusing and self-compression of Gaussian laser pulse in such plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗