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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

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↗
DOE OSTI2025

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↗
DOE OSTI2025

Inferring fracture dilation and shear slip from surface deformation utilising trained surrogate models

An important task in energy and CO 2 storage (sequestration) in the subsurface is to verify that the surrounding fractures and faults are not activated, acting as leakage pathways. This is achievable through effective and efficient Measurement, Monitoring and Verification (MMV) plans. In this work, two surrogate models are trained to captures dilation (opening) and shear deformation of fractures, and the associated surface deformation. The trained surrogate model, based on conditional Generative-Adversarial Networks (cGAN) receives fracture apertures from dilational fractures together with fracture slips from shear fractures and predicts the combined surface deformation. An inversion algorithm based on Bayesian framework is proposed to identify the geometry of both types of fractures, as well as volume of dilational fractures and deformation moment induced by shear fractures, all from the measured surface deformation data. The inversion algorithm utilises the Differential Evolution (DE) optimisation technique that has the superior performance in finding the global minimum of cost function. The proposed surrogate-assisted inversion successfully inferred the unknown dip, dip direction and the volume of the dilational fractures as well as the induced deformation moment in shear fractures. The model was further tested for the inversion of a field hydraulic fracturing tilt dataset applying different scenarios with varying unknowns to show the model's performance, as well as incorporating shear deformation for better match with the observed data.

Dilation and shear↗
NASA NTRS1974

Development of self-acting seals for helicopter engines

An experimental evaluation of a NASA-designed self-acting face seal for use in advanced gas turbine main shaft positions was conducted. The seal incorporated Rayleigh step pads (self-acting geometry) for lift augmentation. Satisfactory performance of the gas film seal was demonstrated in a 500-hour endurance test at speeds to 183 m/s (600 ft/sec, 54,000 rpm) and air pressure differential of 137 newtons per square centimeter (198.7 psi). Carbon wear was minor. Tests were also conducted with seal seat runout greater than that expected in engine operation and in a severe sand and dust environment. Seal operation was satisfactory in both these detrimental modes of operation.

Lynwander, P.↗
NASA NTRS1991

NGC 1068 - Resolution of nuclear structure in the optical continuum

The HST Planetary Camera has been used to obtain an optical continuum image of the Seyfert galaxy NGC 1068. The image reveals a bright nucleus embedded in an irregular cloudlike structure which is well-differentiated against the background of the galaxy. The nucleus is resolved, with an FWHM of about 0.15 arcsec, or 11 pc. There is no evidence for any unresolved nuclear component. The precise geometry of the scattering region cannot yet be inferred. The cloud surrounding the nucleus is elongated in the NE-SSW direction and has extreme dimensions of 3.5 arcsec x 1.7 arcsec. The brightness centroid is situated 0.4 arcsec SW of the nucleus. It is concluded that the light from the cloud is contributed largely by stars, but that the appearance of the cloud is distinctly dissimilar to what is expected from young or old stellar systems and bears to simple relationship to the numerous features that have heretofore been resolved at other wavelengths.

Lynds, Roger↗
DOE OSTI2024

The Computational Analysis of a Generic Hypersonic Store Separation

Many current hypersonic vehicles involve a main vehicle from which components separate during flight and presents simulation challenges such as shock-shock interactions, highly nonlinear interactions, and shock wave boundary layer interactions. These aerodynamic phenomena influence the attitude and trajectory of each object, which needs to be understood before large-scale experiments can be run. To gain a generic understanding of the process, high-fidelity Reynolds-averaged Navier-Stokes solutions are performed on a conical vehicle geometry passing through an oblique shock wave that is representative of different configurations and scenarios. Trajectory and applied forces are tracked and show the vehicle’s dynamics are predominantly the result of differential flow incidence angles causing a strong shock; this leads to a large pressure increase over a fraction of the vehicle, which influences pitch. The vehicle appears to follow conventional stability theory with detached eddy simulation and small disturbances in initial attitude shown to have minimal influence on the scenario.

Savery, Ryan↗
DOE OSTI2023

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗
NASA NTRS1987

Aerodynamic sensitivities from subsonic, sonic and supersonic unsteady, nonplanar lifting-surface theory

The technique of implicit differentiation has been used in combination with linearized lifting-surface theory to derive analytical expressions for aerodynamic sensitivities (i.e., rates of change of lifting pressures with respect to general changes in aircraft geometry, including planform variations) for steady or oscillating planar or nonplanar lifting surfaces in subsonic, sonic, or supersonic flow. The geometric perturbation is defined in terms of a single variable, and the user need only provide simple expressions or similar means for defining the continuous or discontinuous global or local perturbation of interest. Example expressions are given for perturbations of the sweep, taper, and aspect ratio of a wing with trapezoidal semispan planform. In addition to direct computational use, the analytical method presented here should provide benchmark criteria for assessing the accuracy of aerodynamic sensitivities obtained by approximate methods such as finite geometry perturbation and differencing. The present process appears to be readily adaptable to more general surface-panel methods.

Yates, E. Carson, Jr.↗
NASA NTRS1990

Boundary-integral method for calculating aerodynamic sensitivities with illustration for lifting-surface theory

The technique of implicit differentiation has been used in combination with linearized lifting-surface theory to derive analytical expressions for aerodynamic sensitivities (i.e., rates of change of lifting pressures with respect to general changes in aircraft geometry, including planform variations) for steady or oscillating planar or nonplanar lifting surfaces in subsonic, sonic, or supersonic flow. The geometric perturbation is defined in terms of a single variable, and the user need only provide simple expressions or similar means for defining the continuous or discontinuous global or local perturbation of interest. Example expressions are given for perturbations of the sweep, taper, and aspect ratio of a wing with trapezoidal semispan planform. The present process appears to be readily adaptable to more general surface-panel methods.

Yates, E. Carson, Jr.↗
DOE OSTI2023

Solution–Doped Donor–Acceptor Copolymers Based on Diketopyrrolopyrrole and 3, 3'–Bis (2–(2–(2–Methoxyethoxy) Ethoxy) ethoxy)–2, 2'–Bithiophene Exhibiting Outstanding Thermoelectric Power Factors with p –Dopants

The design of polymeric semiconductors exhibiting high electrical conductivity (σ) and thermoelectric power factor (PF) will be vital for flexible large-area electronics. In this work, four polymers based on diketopyrrolopyrrole (DPP), 2,3-dihydrothieno[3,4-b][1,4]dioxine (EDOT), thieno[3,2-b]thiophene (TT), and 3, 3'-bis (2-(2-(2-methoxyethoxy) ethoxy) ethoxy)-2, 2'-bithiophene (MEET) are investigated as side-chains, with the MEET polymers newly synthesized for this study. These polymers are systematically doped with tetrafluorotetracyanoquinodimethane ( F 4 TCNQ), CF3SO3H, and the synthesized dopant Cp(CN) 3 -(COOMe) 3 , differing in geometry and electron affinity. The DPP-EDOT-based polymer containing MEET as side-chains exhibits the highest conductivity (σ) ≈700 S cm–1 in this series with the acidic dopant (CF 3 SO 3 H). This polymer also shows the lowest oxidation potential by cyclic voltammetry (CV), the strongest intermolecular interactions evidenced by differential scanning calorimetry (DSC), and has the most oxygen-based functionality for possible hydrogen bonding and ionic screening. Other polymers exhibit high σ ≈300–500 S cm–1 and power factor up to 300 µW m –1 K –2 . The mechanism of conductivity is predominantly electronic, as validated by time-dependent conductance studies and transient thermo voltage monitoring over time, including for those doped with the acid. Furthermore, these materials maintain significant thermal stability and air stability over ≈6 weeks. Density functional theory calculations reveal molecular geometries and inform about frontier energy levels. Raman spectroscopy, in conjunction with scanning electron microscopy (SEM-EDS) and x-ray diffraction, provides insight into the solid-state microstructure and degree of phase separation of the doped polymer films. Infrared spectroscopy enables this study to further quantify the degree of charge transfer from polymer to dopant.

36 MATERIALS SCIENCE↗
NASA NTRS1977

Flow study in the cross sectional planes of a turbine scroll

A numerical study of the nonviscous flow characteristics in the cross-sectional planes of a radial inflow turbine scroll is presented. The velocity potential is used in the formulation to determine the flow velocity in these planes resulting from the continuous mass discharge. The effect of the through flow velocity is simulated by a continuous distribution of source/sink in the cross-section. A special iterative procedure is devised to handle the solution of the resulting Poisson's differential equation with Neumann boundary conditions in a domain with generally curved boundaries. The analysis is used to determine the effects of the radius of curvature, the location of the scroll section and its geometry on the flow characteristics in the turbine scroll.

Hamed, A.↗
NASA NTRS1995

Grid Sensitivity and Aerodynamic Optimization of Generic Airfoils

An algorithm is developed to obtain the grid sensitivity with respect to design parameters for aerodynamic optimization. The procedure is advocating a novel (geometrical) parameterization using spline functions such as NURBS (Non-Uniform Rational B- Splines) for defining the airfoil geometry. An interactive algebraic grid generation technique is employed to generate C-type grids around airfoils. The grid sensitivity of the domain with respect to geometric design parameters has been obtained by direct differentiation of the grid equations. A hybrid approach is proposed for more geometrically complex configurations such as a wing or fuselage. The aerodynamic sensitivity coefficients are obtained by direct differentiation of the compressible two-dimensional thin-layer Navier-Stokes equations. An optimization package has been introduced into the algorithm in order to optimize the airfoil surface. Results demonstrate a substantially improved design due to maximized lift/drag ratio of the airfoil.

Sadrehaghighi, Ideen↗
DOE OSTI2022

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1976

A differential game solution to the Coplanar tail-chase aerial combat problem

Numerical results obtained in a simplified version of the one on one aerial combat problem are presented. The primary aim of the data is to specify the roles of pursuer and evader as functions of the relative geometry and of the significant physical parameters of the problem. Numerical results are given in a case in which the slower aircraft is more maneuverable than the faster aircraft. A third order dynamic model of the relative motion is described, for which the state variables are relative range, bearing, and heading. The ranges at termination are arbitary in the present version of the problem, so the weapon systems of both aircraft can be visualized as forward firing high velocity weapons, which must be aimed at the tail pipe of the evader. It was found that, for the great majority of the ralative geometries, each aircraft can evade the weapon system of the other.

Merz, A. W.↗
DOE OSTI2022

Tunable phononic bandgap materials designed via topology optimization

Topology optimization is used to design phononic bandgap materials that are tunable by mechanical deformation. A periodic media is considered, which due to the assumption of length scale separation, allows the dispersion relations to be obtained by analyzing a single unit cell subjected to Floquet–Bloch boundary conditions. A finite macroscopic deformation is applied to the unit cell to affect its geometry and hence dispersion. We tune the dispersion–deformation relation to our liking by solving a topology optimization problem using nonlinear programming. The adjoint method is employed to compute the sensitivities, and the non-differentiability of degenerate eigenvalues is avoided using symmetric polynomials. Several tunable phononic crystal designs are presented. Also, a verification analysis is performed, wherein the optimized design is interpreted and analyzed using a conforming finite element mesh.

42 ENGINEERING↗
DOE OSTI2022

Yield strength of CeO 2 measured from static compression in a radial diamond anvil cell

Cerium oxide (ceria, CeO 2 ) is frequently used as a standard in applications such as synchrotron and x-ray free electron lasers for calibrating x-ray wavelengths and offers the potential for understanding the high pressure properties and deformation mechanisms in a wide range of similar face centered cubic (fcc) materials. In this study, the pressure dependence of the strength of ceria was investigated up to 38 GPa using angle dispersive x-ray diffraction in a radial geometry in a diamond anvil cell. In this experiment, the difference in the stress along the axis of compression and perpendicular to the direction of compression can be determined, giving a quantity known as the differential stress. It was found that the differential stress (t), a measure of the lower bound for yield strength, initially increases rapidly from 0.35 ± 0.06 GPa to 2.2 ± 0.4 GPa at pressures of 1.8 and 3.8 GPa, respectively. Above 4 GPa, t increases more slowly to 13.8 ± 2.6 GPa at a pressure of 38 GPa. The changes in the preferred orientation (texture) of CeO 2 with pressure were also measured, allowing for the determination of active deformation mechanisms using an elasto-viscoplastic self-consistent model (EVPSC). It was found that as pressure increased, the [001] direction had a slight preferred orientation along the axis of compression. Our EVPSC model of experimental fiber (cylindrically symmetric) textures and lattice strains were most consistent with dominant slip activity along ${111}$ $\langle$$1\bar10$$\rangle$.

36 MATERIALS SCIENCE↗
DOE OSTI2024

CurvilinearGrids.jl: A Julia package for curvilinear coordinate transformations

Finite-difference discretizations of partial differential equations are widespread throughout the scientific community. Oftentimes finite-differences are used to compute spatial gradients of fields on a discrete grid, which is typically a uniform or rectilinear Cartesian mesh. Arbitrary multidimensional geometry is difficult to discretize directly with finite differences, however, due to non-uniform grid spacing and non-orthogonality. Curvilinear coordinate transformations can be used as a strategy to enable arbitrary geometry. While these curvilinear transformations are straightforward, the governing PDEs require additional terms (metrics) and must adhere to strict conservation laws; these criteria complicate the application of the transformation and require careful implementation.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

Periodic motion in a class of nth-order autonomous differential equations

Sufficient conditions are obtained for the existence of periodic motion in a class of autonomous nonlinear differential equations of order greater than two. The approach is based on the decomposition of an equation into a linear and a nonlinear part. The analysis relies on some basic ideas from linear analysis and geometry. Sufficient conditions for a periodic solution are derived by means of a general topological principle referred to as the torus principle. The existence of a periodic solution is concluded by an appropriate use of the Brouwer fixed-point theorem.

Williamson, D.↗