SEARCH · Engineering Papers

Results for “interface capturing”

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 883 records · Page 49

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1974

Six degree of freedom FORTRAN program, ASTP docking dynamics, users guide

The digital program ASTP Docking Dynamics as outlined is intended to aid the engineer using the program to determine the docking system loads and attendant vehicular motion resulting from docking two vehicles that have an androgynous, six-hydraulic-attenuator, guide ring, docking interface similar to that designed for the Apollo/Soyuz Test Project (ASTP). This program is set up to analyze two different vehicle combinations: the Apollo CSM docking to Soyuz and the shuttle orbiter docking to another orbiter. The subroutine modifies the vehicle control systems to describe one or the other vehicle combinations; the rest of the vehicle characteristics are changed by input data. To date, the program has been used to predict and correlate ASTP docking loads and performance with docking test program results from dynamic testing. The program modified for use on IBM 360 computers. Parts of the original docking system equations in the areas of hydraulic damping and capture latches are modified to better describe the detail design of the ASTP docking system.

Mount, G. O., Jr.↗

DOE OSTI2021

Atomic layer deposition of sodium fluoride thin films

The need for advanced energy conversion and storage devices remains a critical challenge amid the growing worldwide demand for renewable energy. Metal fluoride thin films are of great interest for applications in lithium-ion and emerging rechargeable battery technologies, particularly for enhancing the stability of the electrode-electrolyte interface and thereby extending battery cyclability and lifetime. Reported within, sodium fluoride (NaF) thin films were synthesized via atomic layer deposition. NaF growth experiments were carried out at reactor temperatures between 175 and 250 °C using sodium tert-butoxide and HF-pyridine solution. Here, the optimal deposition temperature range was 175–200 °C, and the resulting NaF films exhibited low roughness (R q ≈ 1.6 nm for films of ~8.5 nm), nearly stoichiometric composition (Na:F = 1:1.05) and a growth per cycle value of 0.85 Å/cycle on SiO 2 substrates. These results are encouraging for future applications of NaF thin films in the development of improved energy capture and storage technologies.

25 ENERGY STORAGE↗

NASA NTRS2022

Design of Flexures for Deployable Reflectarrays Using High Strain Composites

This paper presents a novel flexure design for small scale refelectarrays, specifically CubeSats, that requires high level of deployment precision. Advantages of using a HSC flexure includes additional thermal stability, lesser part count and passive deployment mechanism using the stored strain energy. The proposed design is to use two flat HSC strips of thickness ℎ, length 2𝐿ℎ, width 𝑊ℎ that attach to the core-reflector panel interface. The flexure satisfies two main requirements: ability to fold to a specific compaction 𝑑 without material failure and sufficient stiffness in the deployed configuration. The elastica analysis was implemented to identify the maximum strains in the flexure when folded and provide relations that combine Δ, ℎ, 𝐿ℎ with the maximum allowable strain. A simplified analysis was performed of the vibration of a panel-flexure combination using a two degree of freedom system. The vibration analysis reveals a dimensionless length parameter 𝛽 = (𝐿ℎ/𝐿𝑃) (Δ/ℎ) that determines the transition between bending and shear modes. A non-dimensional relationship was obtained between the natural frequency and 𝛽 that captures the impact of each parameter on the response of the structure. Our current work involves fabricating and testing of the proposed flexures using High Strain Composite materials, with plans to build a prototype for a deployable reflectarray.

López Jiménez, Francisco↗

DOE OSTI2020

FFT-based solver for higher-order and multi-phase-field fracture models applied to strongly anisotropic brittle materials

This paper presents the application of a fast Fourier transform (FFT) based method to solve two phase field models designed to simulate crack growth of strongly anisotropic materials in the brittle regime. By leveraging the ability of the FFT-based solver to generate solutions with higher-order and global continuities, we design two simple algorithms to capture the complex fracture patterns (e.g. sawtooth, and curved crack growth) common in materials with strongly anisotropic surface energy via the multi-phase-field and high-order phase-field frameworks. A staggered operator-split solver is used where both the balance of linear momentum and the phase field governing equations are formulated in the periodic domain. The unit phase field of the initial failure region is prescribed by the penalty method to alleviate the sharp material contrast between the initial failure region and the base material. The discrete frequency vectors are generalized to estimate the second and fourth order gradients such that the Gibbs effect near shape interfaces or jump conditions can be suppressed. Furthermore, a preconditioner is adopted to improve the convergence rate of the iterative linear solver. Three numerical experiments are used to systematically compare the performance of the FFT-based method in the multi-phase-field and high-order phase-field frameworks.

42 ENGINEERING↗

DOE OSTI2021

Striped electron fluid on (111) KTaO 3

A recent study has revealed that the low carrier density electron gas (2DEG) induced at the interface of EuO and (111) KTaO 3 exhibits a broken symmetry phase with a strong in-plane anisotropy of the resistivity. We present a minimal tight-binding model of this (111) 2DEG, including the large spin-orbit coupling from the Ta ions, which reveals a hexagonal Fermi surface with a highly enhanced 2k F electronic susceptibility. We argue that repulsive electronic interactions, together with a ferromagnetic EuO substrate, favor a magnetic stripe instability leading to a partially gapped Fermi surface. Such a stripe state, or its vestigial nematicity, could explain the observed transport anisotropy. We propose a k · p theory for the low energy j = 3/2 states, which captures the key results from our tight-binding study, and further reveals the intertwined dipolar and octupolar modulations underlying this magnetic stripe order. We conclude by speculating on the relation of this stripe order to the superconductivity seen in this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Standard Library Plant Controller Model Specification for a Grid-Forming Hybrid Control Inverter-Based Resource (REPCGFM_C1)

This document describes a standard library plant controller model to interface with the grid-forming hybrid control inverter-based resource (IBR) model. The initial version of model specification was jointly developed by Pacific Northwest National Laboratory (PNNL), Tesla Energy, and EPRI, and it was revised multiple times later to incorporate suggestions from WECC MVS members. Tesla Energy provided main control algorithms to support the development of this model specification. This standard library model is developed to help the utility industry better understand the GFM technology. The model could be used to represent equipment for long-term planning studies where vendor-specific models are not available. As equipment matures and improves, standard library models will be updated to capture the new functionalities of GFMs. It is not intended that these models will always remain representative of all future GFM technologies.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS2012

Provenance Representation in the Global Change Information System (GCIS)

Global climate change is a topic that has become very controversial despite strong support within the scientific community. It is common for agencies releasing information about climate change to be served with Freedom of Information Act (FOIA) requests for everything that led to that conclusion. Capturing and presenting the provenance, linking to the research papers, data sets, models, analyses, observation instruments and satellites, etc. supporting key findings has the potential to mitigate skepticism in this domain. The U.S. Global Change Research Program (USGCRP) is now coordinating the production of a National Climate Assessment (NCA) that presents our best understanding of global change. We are now developing a Global Change Information System (GCIS) that will present the content of that report and its provenance, including the scientific support for the findings of the assessment. We are using an approach that will present this information both through a human accessible web site as well as a machine readable interface for automated mining of the provenance graph. We plan to use the developing W3C PROV Data Model and Ontology for this system.

Tilmes, Curt↗

DOE OSTI2026

Dynamically reconfigurable topological routing in nonlinear photonic systems

The propagation path of topologically protected states is bound to the interface between regions with different topology, and as such, the functionality of linear photonic devices leveraging these states is fixed during fabrication. Here, we propose a mechanism for dynamic control over a driven dissipative system’s local topology, yielding reconfigurable topological interfaces and thus tunable paths for protected routing. We illustrate our approach in non-resonantly pumped polariton lattices, where the nonlinear interaction between the polaritons and the exciton reservoir due to non-resonant pumping can yield picosecond-scale changes in the propagation paths of the chiral edge states. To analytically confirm the numerically observed topological dynamics, we generalize the spectral localizer framework to non-linear non-Hermitian Chern materials and apply this framework to a continuous model of the polariton system based on a driven-dissipative Gross-Pitaevskii equation. In doing so, we show that the local changes in the polariton lattice’s topology are captured by a local Chern marker. Looking forward, we anticipate such reconfigurable topological routing will enable the realization of novel classes of topological photonic devices.

Wong, Stephan [Sandia National Laboratories (SNL-N↗

DOE OSTI2021

Kinetic mix at gas-shell interface in inverted corona fusion targets

Gas-filled, laser-driven “inverted corona” fusion targets have attracted interest as a low-convergence neutron source and platform for studying kinetic physics. At the fill pressures under investigation, ejected particles from the shell can penetrate deeply into the gas before colliding, leading to significant mixing across the gas–shell interface. Here, we use kinetic-ion, fluid-electron hybrid particle-in-cell simulations to explore the nature of that mix. Simulations of the system demonstrate characteristics of a weakly collisional electrostatic shock, whereby a strong electric field accelerates shell ions into the rarefied gas and reflects upstream gas ions. This interpenetration is mediated by collisional processes: At higher initial gas pressure, fewer shell particles pass into the mix region and reach the hotspot. This effect is detectable through neutron yield scaling vs gas pressure. Predictions of neutron yield scaling show excellent agreement with experimental data recorded at the OMEGA laser facility, suggesting that 1D kinetic mechanisms are sufficient to capture the mix process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Redistribution of subsurface neutrons caused by ground ice on Mars

The purpose of this study is to build on that of Squyres and Evans to approach the specific problem of the buried ice deposits explored by Mellon and Jakosky. Therefore, we calculated the depth dependence of equilibrium neutron energy spectra and related neutron-capture gamma ray production rates for a semi-infinite layer of regolith having a Chryse-fines composition containing 5, 10, and 18 wt percent of water. This high water-content layer was buried beneath the same regolith containing 1 percent water at depths of 20, 40, and 80 g/ sq. cm. A 16-g/sq. cm atmosphere having the abundances measured by Owen et al. was included throughout. We also calculated the effects of a 30-g/sq. cm-thick layer containing 18 percent water embedded within a semi-finite regolith containing 1 percent water, at depths of 20, 40, and 80 g/sq. cm below the atmosphere-surface interface.

Feldman, W. C.↗

NASA NTRS2016

Biomechanical Modeling of Split-leg Squat and Heel Raise on the Hybrid Ultimate Lifting Kit (HULK)

Long duration space travel will expose astronauts to extended periods of reduced gravity. Since gravity is not present to aid loading, astronauts will use resistive and aerobic exercise regimes for the duration of the space flight to minimize the loss of bone density, muscle mass and aerobic capacity that occurs during exposure to a reduced gravity environment. Unlike the International Space Station (ISS), the area available for an exercise device in the next generation of spacecraft is limited and therefore compact resistance exercise device prototypes are being developed. The Advanced Resistive Exercise Device (ARED) currently on the ISS is being used as a benchmark for the functional performance of these new devices. Biomechanical data collection and computational modeling aid the device design process by quantifying the joint torques and musculoskeletal forces that occur during exercises performed on the prototype devices. Computational models currently use OpenSim software, an open source code for musculoskeletal modeling, with biomechanical input data from subjects for estimation of muscle and joint loads. Subjects are instrumented with reflective markers for motion capture data collection while exercising on the Hybrid Ultimate Lifting Kit (HULK) prototype device. Ground reaction force data is collected with force plates under the feet and device loading is recorded through load cells internal to the HULK. This data is input into the OpenSim biomechanical model, which has been scaled to match the anthropometrics of the test subject, to calculate the loads on the body. Multiple exercises are performed and evaluated during a test session such as a full squat, single leg squat, heel raise and dead lift. Variables for these exercises include applied device load, narrow or wide foot stance, slow or fast cadence and the harness or long bar interface between the test subject and the device. Data from free weights are compared to the resistively loaded exercise device. The focus of this presentation is to summarize the results from the single-leg squat and heel raise exercises performed during three sessions occurring in 2015. Differences in loading configuration, cadence and stance produce differences in kinematics, joint toques and force and muscle forces.

computational modeling↗

NASA NTRS2016

Structural Design Considerations for a 50 kW-Class Solar Array for NASA's Asteroid Redirect Mission

NASA is planning an Asteroid Redirect Mission (ARM) to take place in the 2020s. To enable this multi-year mission, a 40 kW class solar electric propulsion (SEP) system powered by an advanced 50 kW class solar array will be required. Powered by the SEP module (SEPM), the ARM vehicle will travel to a large near-Earth asteroid, descend to its surface, capture a multi-metric ton (t) asteroid boulder, ascend from the surface and return to the Earth-moon system to ultimately place the ARM vehicle and its captured asteroid boulder into a stable distant orbit. During the years that follow, astronauts flying in the Orion multipurpose crew vehicle (MPCV) will dock with the ARM vehicle and conduct extra-vehicular activity (EVA) operations to explore and sample the asteroid boulder. This paper will review the top structural design considerations to successfully implement this 50 kW class solar array that must meet unprecedented performance levels. These considerations include beyond state-of-the-art metrics for specific mass, specific volume, deployed area, deployed solar array wing (SAW) keep in zone (KIZ), deployed strength and deployed frequency. Analytical and design results are presented that support definition of stowed KIZ and launch restraint interface definition. An offset boom is defined to meet the deployed SAW KIZ. The resulting parametric impact of the offset boom length on spacecraft moment of inertias and deployed SAW quasistatic and dynamic load cases are also presented. Load cases include ARM spacecraft thruster plume impingement, asteroid surface operations and Orion docking operations which drive the required SAW deployed strength and damping. The authors conclude that to support NASA's ARM power needs, an advanced SAW is required with mass performance better than 125 W/kg, stowed volume better than 40 kW/cu m, a deployed area of 200 sq m (100 sq m for each of two SAWs), a deployed SAW offset distance of nominally 3-4 m, a deployed SAW quasistatic strength of nominally 0.1 g in any direction, a deployed loading displacement under 2 m, a deployed fundamental frequency above 0.1 Hz and deployed damping of at least 1%. These parameters must be met on top of challenging mission environments and ground testing requirements unique to the ARM project.

solar arrays↗

DOE OSTI2022

A level set approach for the computational study of a yield stress fluid filling a thin mold

Many important engineering and scientific applications such as cement slurries, foams, crude oil, and granular avalanches involve the concept of yield stress. Therefore, modeling yield stress fluids in different flow configurations, including the accurate prediction of the yield surface, is important. In this paper, we present a computational model based on the finite element method to study the flow of yield stress fluids in a thin mold and compare the results with data from flow visualization experiments. Here, we use the level set method to describe the interface between the filling fluid and air. We use polypropylene glycol as a model Newtonian fluid and Carbopol for the model yield stress fluid, as the Carbopol solution demonstrates yielding without thixotropy. To describe the yielding and shear-thinning behavior, we use a generalized Newtonian constitutive equation with a Bingham–Carreau–Yasuda form. We compare the results obtained from the mold filling experiments with the results from the three-dimensional (3D) model and from a reduced-order Hele-Shaw (HS) model that is two-dimensional, including the effect of shear-thinning along the thin direction only approximately. We show that both the 3D and the HS model can capture the experimental meniscus shape reasonably well for all the fluids considered at three different flow rates. This indicates that the shape evolution is insensitive to the dimensionality of the model. However, the viscosity and yield surfaces predicted by the 3D and HS models are different. The HS model underestimates the high viscosity and unyielded regions compared to the estimation by the 3D model.

42 ENGINEERING↗

DOE OSTI2020

Infrared spectroscopy data- and physics-driven machine learning for characterizing surface microstructure of complex materials

There is a need to characterize complex materials and their dynamics under reaction conditions to accelerate materials design. Adsorbate vibrational excitations are selective to adsorbate/surface interactions and infrared (IR) spectra associated with activating adsorbate vibrational modes are accurate, capture details of most modes, and can be obtained operando. Current interpretation depends on heuristic peak assignments for simple spectra, precluding the possibility of obtaining detailed structural information. Here, we combine data-based approaches with chemistry-dependent problem formulation to develop physics-driven surrogate models that generate synthetic IR spectra from first-principles calculations. Using synthetic IR spectra of carbon monoxide on platinum, we implement multinomial regression via neural network ensembles to learn probability distributions functions (pdfs) that describe adsorption sites and quantify uncertainty. We use these pdfs to infer detailed surface microstructure from experimental spectra and extend this methodology to other systems as a first step towards characterizing complex interfaces and closing the materials gap.

Science & Technology - Other Topics↗

NASA NTRS2016

Modeling Fatigue Damage Onset and Progression in Composites Using an Element-Based Virtual Crack Closure Technique Combined With the Floating Node Method

A new methodology is proposed to model the onset and propagation of matrix cracks and delaminations in carbon-epoxy composites subject to fatigue loading. An extended interface element, based on the Floating Node Method, is developed to represent delaminations and matrix cracks explicitly in a mesh independent fashion. Crack propagation is determined using an element-based Virtual Crack Closure Technique approach to determine mixed-mode energy release rates, and the Paris-Law relationship to obtain crack growth rate. Crack onset is determined using a stressbased onset criterion coupled with a stress vs. cycle curve and Palmgren-Miner rule to account for fatigue damage accumulation. The approach is implemented in Abaqus/Standard® via the user subroutine functionality. Verification exercises are performed to assess the accuracy and correct implementation of the approach. Finally, it was demonstrated that this approach captured the differences in failure morphology in fatigue for two laminates of identical stiffness, but with layups containing θdeg plies that were either stacked in a single group, or distributed through the laminate thickness.

De Carvalho, Nelson V.↗

DOE OSTI2022

A hierarchical stochastic modeling approach for representing point bar geometries and petrophysical property variations

The flow of fluids in point bars is affected by the existence of heterogeneities like shale drapes that are found on the surfaces of inclined heterolithic stratifications. In fact, these shale drapes can act as fluid flow baffles; therefore, developing a framework for modeling point bars and their associated heterogeneities is vital. In this study, a stochastic process-based modeling approach is presented for capturing the main point bar heterogeneities: accretion surfaces (i.e., the aerial heterogeneity) and inclined heterolithic stratifications (i.e., the vertical heterogeneity). The former was modeled using a sine-generation function and the latter, with a sigmoidal function after which they were combined into a 3D point bar model. To ensure proper modeling of petrophysical properties, we developed a more representative gridding scheme which generates curvilinear grids representative of the point bar geometry. This grid was then transformed into a rectilinear grid to allow for geostatistical simulation after which all petrophysical properties were mapped back into the original curvilinear grid. An essential element of this modeling approach is the stochastic representation of shale drapes at the interface between successive accretion surfaces. The workflow was tested using a real field dataset for the Cranfield field, Mississippi. The constructed model was then subjected to a flow simulation study mimicking a CO 2 storage scenario. Various sensitivities were simulated to evaluate the effect of heterogeneities on CO 2 flow within the point-bar. Results demonstrate the importance of representing point-bar related heterogeneity and the spatial distribution of shale drapes on CO 2 plume migration and storage.

58 GEOSCIENCES↗