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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Solar Simulation

Simulation of the surface region and interior of the Sun's convection zone. The objective is to provide understanding of the tachocline at the base of the convection zone. of the differential rotation generated by global convective motions, and of supergranules that are observed on the solar surface. Two large-scale simulation codes are used, one for the global dynamics in spherical coordinates and another for the local surface events in Cartesian geometry. We work in close collaboration with the Helioseismic and Magnetic Imager instrument developers for the Solar Dynamics Observatory satellite. NASA Program: Living with a Star - Sun-Earth Connection.

Wray, Alan↗
NASA NTRSDate not supplied

Implementation of Combinatorial Optimization Techniques for Automated Fiber Placement Through Thickness Defect Stack-Up Minimization

The Computer Aided Process Planning (CAPP) module was developed to facilitate and accelerate the process planning workflow for Automated Fiber Placement (AFP). CAPP assists process planners in identifying optimal starting point locations and layup strategies for each ply of a laminate. Ply optimization operates on measurement and scoring of geometry-based defects such as gaps, overlaps, angle deviation, and steering. This paper expands on the established framework for analyzing defect stack-up through thickness of a laminate. Four different combinatorial optimization algorithms are implemented and evaluated: (1) genetic algorithm, (2) differential evolution, (3) particle swarm, and (4) greedy search. The algorithms identify the optimal combination of ply-level layup strategies, by scoring potential laminates on defect stacking, using two different objective functions. A final optimization approach is also presented which trades some performance for a large gain in efficiency. These approaches are compared to a randomized combination using a complex tool surface in a virtual case study. The result is a streamlined methodology for comparing different laminate-level manufacturing strategies and minimizing the through thickness defect stack up.

CAPP↗
DOE OSTI2017

Multiscale models and stochastic simulation methods for computing rare but key binding events in cell biology

The main difficulty in simulating diffusion processes at a molecular level in cell microdomains is due to the multiple scales involving nano- to micrometers. Few to many particles have to be simulated and simultaneously tracked while there are exploring a large portion of the space for binding small targets, such as buffers or active sites. Bridging the small and large spatial scales is achieved by rare events representing Brownian particles finding small targets and characterized by long-time distribution. These rare events are the bottleneck of numerical simulations. A naive stochastic simulation requires running many Brownian particles together, which is computationally greedy and inefficient. Solving the associated partial differential equations is also difficult due to the time dependent boundary conditions, narrow passages and mixed boundary conditions at small windows. We present here two reduced modeling approaches for a fast computation of diffusing fluxes in microdomains. The first approach is based on a Markov mass-action law equations coupled to a Markov chain. The second is a Gillespie's method based on the narrow escape theory for coarse-graining the geometry of the domain into Poissonian rates. The main application concerns diffusion in cellular biology, where we compute as an example the distribution of arrival times of calcium ions to small hidden targets to trigger vesicular release.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Differential die-away instrument: Report on comparison of fuel assembly experiments and simulations

Experimental results of the assay of mock-up (fresh) fuel with the differential die-away (DDA) instrument were compared to the Monte Carlo N-Particle eXtended (MCNPX) simulation results. Most principal experimental observables, the die-away time and the in tegral of the DDA signal in several time domains, have been found in good agreement with the MCNPX simulation results. The remaining discrepancies between the simulation and experimental results are likely due to small differences between the actual experimental setup and the simulated geometry, including uncertainty in the DT neutron generator yield. Within this report we also present a sensitivity study of the DDA instrument which is a complex and sensitive system and demonstrate to what degree it can be impacted by geometry, material composition, and electronics performance.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2016

State of knowledge and challenges of U-Si compounds for use in light water reactor accident tolerant fuel designs

Characterization of the properties of U-Si compounds first required synthesis of high purity material for each of the compositions. The four compounds targeted for investigation in this work were U{sub 3}Si, U{sub 3}Si{sub 2}, USi, and U{sub 3}Si{sub 5}. Specimens of U{sub 3}Si{sub 2}, U{sub 3}Si{sub 5}, and USi with specific geometries for individual testing procedures were prepared using a powder metallurgical method. Thermophysical properties of the U-Si compounds were determined within 100 K of their respective melt point or decomposition temperature. Thermal expansion was determined with a dilatometer, specific heat capacity by the ratio method in a differential scanning calorimeter using a sapphire standard, while thermal diffusivity was determined with a laser flash analyzer. Oxidation performance of the U-Si compounds was determined in a thermogravimetric analyzer (TGA), which measures the weight change of a sample as a function of gas composition and temperature. The results show that each of the U-Si compounds is electronically conductive leading to higher thermal conductivity values than the reference fuel UO{sub 2}. Higher thermal conductivity values will lead to enhanced safety margins within a reactor core and provide better economics over the lifetime of a fuel cycle. However, no improvements have been realized over the reference UO{sub 2} case with regards to oxidation performance in air although more testing is required to substantiate the differences especially in steam and corrosive aqueous LWR environments. Finally, there exists a lack of knowledge in the literature with regards to low burnup irradiation testing on all of the U-Si binary compounds, which would be necessary to simulate their performance in a reactor setting.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Solid-State Syntheses, Crystallographic Spatial Disorders, Thermal Behaviors, and Bandgaps of Hybrid Organic-Inorganic Manganese Halides: A2Mn(Cl/Br)4 (A = NH4+, C(NH2)3+, & C3H4N2+)

By systematically optimizing solid-state synthesis via the concerted use of differential scanning calorimetry and in-situ variable-temperature X-ray diffraction, we report the discovery of four new hybrid organic-inorganic manganese halides A2Mn(Cl/Br)4 (A = NH4+, C(NH2)3+, C3H4N2+), with emphasis on how the organic fragment geometry, polarity, and electron-donating properties influence crystallographic spatial disorders, thermal behaviors, and optical band gaps. With a nonpolar tetrahedral cation, as in (NH4)2MnCl4 (P21/c, mP30) or (NH4)2MnBr4 (C2/c, mC180), the direct optical bandgaps (4.36–4.28 eV) are wider than those with an organic nonpolar planar geometry, as in (C(NH2)3)2MnBr4 (P21/c, mP300, 4.07 eV), or an organic polar planar geometry, as in (C3H4N2+)2MnBr4 (I41/a, tI384, 4.05 eV). Furthermore, the organic components affect the spatial arrangement and dimensionality of the resulting inorganic frameworks comprised of Mn(Cl/Br)6 octahedra, with varying distortions, or spatially disordered MnBr4 tetrahedra. These resulting trends, coupled with the systematic investigation into synthesis, and the motivation behind elucidating the nuanced role of the organic cation, altogether expand upon the phase-space of hybrid materials beyond perovskites and demonstrate the potential for a rational design of hybrid materials with tailored optoelectronic functionalities for advanced energy applications.

Lee, Shannon J.↗
DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2023

Parallel Algebraic Multigrid for Fusion and Higher-Order PDEs

Multigrid methods play a key role in large-scale scientific simulation because they are among the fastest and most scalable approaches for solving the underlying sparse linear systems of equations that arise from a wide array of Partial Differential Equation (PDE) discretizations. Algebraic multigrid (AMG) is a special type of multigrid method that depends only on the description of the linear system, giving it better portability and broader applicability than geometric multigrid, as it requires no explicit knowledge of the problem geometry. Even though these methods are widely used today, there are still applications where further development is needed. In this report, we focus on PDEs with higher-order terms (e.g., fourth order), concentrating on a PDE that arises in tokamak edge plasma simulations (a tokamak is a machine that confines a plasma using magnetic fields and is believed to be the leading plasma confinement concept for future fusion power plants). General multigrid relaxes a linear system on coarser grids and reverses this process with interpolation, but standard AMG methods struggle with the aforementioned higher-order PDEs. We investigate cyclic coarsening and interpolation heuristics, as well as new iterative approximation methods of refining the solution at each grid to improve the existing multigrid approach. To this end, we ensure that these techniques are transferable to a parallelized setting with LLNL’s supercomputers.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1974

Protonospheric columnar electron content determination. II - Observations

The results of over 300 hours of measurements of the quantity I sub w - the difference between total columnar electron content, derived from the differential phase-path data, and Faraday electron content - is presented. The interpretation of the measurements in terms of protonospheric content is discussed in the light of a previous analysis which demonstrates that good protonospheric content measurements can only be obtained under certain geometries of observations unless information on the shape and height of the ionospheric layer is available. A strong correlation between I sub w and the ionospheric content is observed in most of the Arecibo data, because they were collected under unfavorable conditions. The Stanford data, on the other hand, yield good measurements of the protonospheric content. The May 1969 observations show a strong depletion of protonospheric content during a magnetic storm and its subsequent recovery. This is in general agreement with flux tube content behavior observed by the whistler technique during similar conditions.

Almeida, O. G.↗
DOE OSTI2020

Estimation of Respirable Aerosol Release Fractions through Stress Corrosion Crack-Like Geometries

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the significant backfill pressurization of the canister, up to approximately 800 kPa. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the carrier gas flow rates and aerosol transmission properties were evaluated for an engineered microchannel with characteristic dimensions similar to those of an SCC. The microchannel was formed by mating two gage blocks with a slot orifice measuring 28.9 μm (0.0011 in.) tall by 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, Ce0 2 , were seeded and mixed inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to demonstrate a new capability to characterize SCCs under well-controlled boundary conditions. Separate modeling efforts are also underway that will be validated using these data. The test apparatus and procedures developed in this study can be easily modified for the evaluation of more complex SCC-like geometries including laboratory-grown SCC samples.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Fractional modeling of viscoelasticity in 3D cerebral arteries and aneurysms

We develop efficient numerical methods for fractional order PDEs, and employ them to investigate viscoelastic constitutive laws for arterial wall mechanics. Recent simulations using one-dimensional models [1] have indicated that fractional order models may offer a more powerful alternative for modeling the arterial wall response, exhibiting reduced sensitivity to parametric uncertainties compared with the integer-calculus-based models. Here, we study three-dimensional (3D) fractional PDEs that naturally model the continuous relaxation properties of soft tissue, and for the first time employ them to simulate flow structure interactions for patient-specific brain aneurysms. To deal with the high memory requirements and in order to accelerate the numerical evaluation of hereditary integrals, we employ a fast convolution method [2] that reduces the memory cost to O(log(N)) and the computational complexity to O(N log(N)). Furthermore, we combine the fast convolution with high-order backward differentiation to achieve third-order time integration accuracy. We confirm that in 3D viscoelastic simulations, the integer order models strongly depends on the relaxation parameters, while the fractional order models are less sensitive. As an application to long-time simulations in complex geometries, we also apply the method to modeling fluid-structure interaction of a 3D patient-specific compliant cerebral artery with an aneurysm. Taken together, our findings demonstrate that fractional calculus can be employed effectively in modeling complex behavior of materials in realistic 3D time-dependent problems if properly designed efficient algorithms are employed to overcome the extra memory requirements and computational complexity associated with the non-local character of fractional derivatives.

Computer Science↗
DOE OSTI2015

Differential rotation in magnetized and non-magnetized stars

The effects of magnetic field on stellar differential rotation (DR) are studied by comparing magnetohydrodynamic (MHD) models and their hydrodynamic (HD) counterparts in the broad range of rotation rates and across varying initial rotation profiles. Fully compressible MHD simulations of rotating penetrative convection are performed in a full-spherical shell geometry. Critical conditions for the transition of the DR between a faster equator (solar type) and a slower equator (anti-solar type) are explored by focusing on the “Rossby number (Ro)” and the “convective Rossby number (Ro{sub conv}).” It is confirmed that the transition is more gradual and the critical value for it is higher in the MHD model than in the HD model in view of the Ro{sub conv} dependence. As observed in earlier studies, the rotation profile shows a bistability near the transition in the HD model, whereas it disappears when allowing the growth of magnetic fields except in the model taking the anti-solar-type solution as the initial condition. We find that the transition occurs at Ro≃1 both in the MHD and HD models independently of the hysteresis. Not only the critical value but also the sharpness of the transition is similar between the two models in view of the Ro dependence. The influences of the dynamo-generated magnetic field and/or the hysteresis on convective motion are reflected by Ro. This is why the transition is unified in view of the Ro dependence. We also discuss the Ro dependence of magnetic dynamo activities with emphasis on their possible relation to the kinetic helicity profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Optimization of a cyclone using MFIX and Nodeworks

Video depicting the optimization process of a cyclone on NETL's chemical looping reactor (CLR) using MFIX and Nodeworks. MFIX is used to model the cyclone using PIC. Nodeworks is then used to generate proposed geometry changes using a Latin hypercube. Each design is simulated, with an objective value being computed based on the cyclone efficiency and pressure drop. A Gaussian Process surrogate model is then constructed from the objective values. This surrogate model is then used by a differential evolution optimization algorithm to identify the optimal cyclone design. Details published here: Weber, J., Fullmer, W., Gel, A., and Musser, J. (February 4, 2020). "Optimization of a Cyclone Using Multiphase Flow Computational Fluid Dynamics." ASME. J. Fluids Eng. March 2020; 142(3): 031111. https://doi.org/10.1115/1.4045952 OSTI: https://www.osti.gov/pages/servlets/purl/1763893

cyclone↗
DOE OSTI2020

EBS Task Force: Task 9/FEBEX Modeling Final Report: Thermo-Hydrological Modeling with PFLOTRAN

This report outlines Sandia National Laboratories modeling studies applied to Stage 1 and Stage 2 of the Full-scale Engineered Barriers Experiment in Crystalline Host Rock (FEBEX) in situ test for the SKB EBS Task Force Task 9. The FEBEX test was a full-scale test conducted over ~18 years at the Grimsel, Switzerland Underground Research Laboratory (URL) managed by NAGRA. It involved emplacing simulated waste packages, in the form of welded cylindrical heaters, inside a tunnel in crystalline granitic rock and surrounded by a bentonite barrier and cement plug. Sensors emplaced within the bentonite monitored the wetting-up, heating, and drying out of the bentonite barrier, and the large resulting data set provides an excellent opportunity for validation of multiphysics Thermal-Hydrological (TH), Thermal-Hydrologic-Chemical (THC), and Thermal-Hydrological-Mechanical (THM) modeling approaches for underground nuclear waste storage and the performance of engineered bentonite barriers. The present status of the EBS Task Force is finalizing Task 9, which follows years of modeling studies of the FEBEX test, by many notable modeling teams (Gens et al., 2009; Sanchez et al. 2010; 2012; Samper et al., 2018). These modeling studies generally use two-dimensional axisymmetric meshes, ignoring threedimensional effects, gravity and asymmetric wetting and dry out of the bentonite engineered barrier. This study investigates these effects with use of the PFLOTRAN THC code with massively parallel computational methods in modeling FEBEX Stage 1 and Stage 2 results. The PFLOTRAN numerical code is an open source, state-of-the-art, massively parallel subsurface flow and reactive transport code operating in a high-performance computing environment (Hammond et al., 2014). Section 2 describes the applied partial differential equations describing mass, momentum and energy balance used in this study, considerations derived by assuming phase equilibrium between gas and liquid phases, constitutive equations for granite, cement plug, and bentonite domains, and specific approaches for use inthe PFLOTRAN code. Section 3 describes the geometry, meshing, and model set-up. Section 4 describes modeling results, Section 5 compares modeling results to field testing data, and Section 6 gives conclusions. The Appendix provides detailed information required by the EBSTask Force for final reporting.

42 ENGINEERING↗
NASA NTRS1993

Grid and aerodynamic sensitivity analyses of airplane components

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 wing-section geometry. An interactive algebraic grid generation technique, known as Two-Boundary Grid Generation (TBGG) is employed to generate C-type grids around wing-sections. 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 wing-section surface. Results demonstrate a substantially improved design due to maximized lift/drag ratio of the wing-section.

Sadrehaghighi, Ideen↗
NASA NTRS1995

Two-dimensional mesh embedding for Galerkin B-spline methods

A number of advantages result from using B-splines as basis functions in a Galerkin method for solving partial differential equations. Among them are arbitrary order of accuracy and high resolution similar to that of compact schemes but without the aliasing error. This work develops another property, namely, the ability to treat semi-structured embedded or zonal meshes for two-dimensional geometries. This can drastically reduce the number of grid points in many applications. Both integer and non-integer refinement ratios are allowed. The report begins by developing an algorithm for choosing basis functions that yield the desired mesh resolution. These functions are suitable products of one-dimensional B-splines. Finally, test cases for linear scalar equations such as the Poisson and advection equation are presented. The scheme is conservative and has uniformly high order of accuracy throughout the domain.

Shariff, Karim↗
NASA NTRS1993

Modeling bristle lift-off in idealized brush seal configurations

We attempt in this paper to develop a model for the flow through brush seals and determine their elastic behavior in order to predict the dependence of brush/journal clearance on geometry and operating conditions. Several idealizations regarding brush seal configuration, flow conditions, and elastic behavior are made in the analysis in order to determine closed form parametric dependence. This formulation assumes that there is no initial interference between the bristle tip and the rotor. Also, interbristle, bristle-backing plate, and bristle-rotor friction is neglected. The bristle bundle or the brush seal as it is alternately called is assumed homogeneous and isotropic on a macroscopic scale so that a physical property like permeability is uniform. The fluid is assumed to be homogeneous, incompressible, viscous, and flowing under steady conditions. A schematic of a brush seal is shown. If the nominal bristle-shaft interference is absent then under static conditions the bristles may deflect axially due to the imposed pressure differential. This axial deflection may create a clearance permitting leakage flow in excess of that which occurs through the porous matrix formed by the bristle bundles. Under dynamic conditions the Couette flow created by shaft motion could be strong enough to cause bristle deflection and once again a clearance may develop. The paper proposes a means to determine this clearance (or at least describe its parametric dependence on geometry and operating conditions) under static as well as dynamic conditions.

Modi, Vijay↗
DOE OSTI2019

3D Printing Amorphous Polysiloxane Terpolymers via Vat Photopolymerization

Abstract Photocuring and vat photopolymerization (VP) additive manufacturing (AM) is reported for two families of fully amorphous poly(dimethyl siloxane) (PDMS) terpolymers containing either diphenylsiloxy (DiPhS) or diethylsiloxy (DiEtS) repeating units. A thiol‐functionalized PDMS crosslinker enables rapid crosslinking in air using efficient thiol–ene addition. Differential scanning calorimetry and dynamic mechanical analysis (DMA) confirm the absence of crystallinity for the DiPhS‐containing systems, while DMA shows a rubbery plateau extending to greater than 200 °C for the DiEtS‐containing system. VP‐AM of both photopolymer systems afford well‐defined 3D geometries, including high aspect ratio structures, which demonstrate feasibility of these photopolymers for the 3D printing of unique geometric objects that require elastomeric performance to temperatures as low as −120 °C.

Sirrine, Justin M.↗