SEARCH · Engineering Papers

Results for “three phase modeling”

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 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Acceleration induced depriming and capillary rewetting of external artery heat pipes - Comparison with SHARE-II flight experiment

A combined analytical and numerical model for the analysis of the deprime and reprime/rewetting characteristics of two high-capacity external artery heat pipe designs undergoing externally induced accelerations was developed using several previously derived analytical expressions. Three distinct phases of the deprime and reprime/rewetting process were analyzed: (1) the effect of longitudinal accelerations on the depriming, (2) the time required for repriming of the liquid artery once the longitudinal acceleration has been terminated, and (3) the rewetting characteristics of the circumferential wall grooves. Combining the three processes, a technique was developed allowing the prediction of the effect of external acceleration on the characteristics of the external artery heat pipes. The predictions made with this technique agreed well with the microgravity flight results.

Ochterbeck, J. M.↗

NASA NTRS2009

Freezable Radiator Coupon Testing and Full Scale Radiator Design

Freezable radiators offer an attractive solution to the issue of thermal control system scalability. As thermal environments change, a freezable radiator will effectively scale the total heat rejection it is capable of as a function of the thermal environment and flow rate through the radiator. Scalable thermal control systems are a critical technology for spacecraft that will endure missions with widely varying thermal requirements. These changing requirements are a result of the space craft s surroundings and because of different thermal loads during different mission phases. However, freezing and thawing (recovering) a radiator is a process that has historically proven very difficult to predict through modeling, resulting in highly inaccurate predictions of recovery time. This paper summarizes tests on three test articles that were performed to further empirically quantify the behavior of a simple freezable radiator, and the culmination of those tests into a full scale design. Each test article explored the bounds of freezing and recovery behavior, as well as providing thermo-physical data of the working fluid, a 50-50 mixture of DowFrost HD and water. These results were then used as a tool for developing correlated thermal model in Thermal Desktop which could be used for modeling the behavior of a full scale thermal control system for a lunar mission. The final design of a thermal control system for a lunar mission is also documented in this paper.

Lillibridge, Sean T.↗

NASA NTRS2013

Neutron Characterization for Additive Manufacturing

Oak Ridge National Laboratory (ORNL) is leveraging decades of experience in neutron characterization of advanced materials together with resources such as the Spallation Neutron Source (SNS) and the High Flux Isotope Reactor (HFIR) shown in Fig. 1 to solve challenging problems in additive manufacturing (AM). Additive manufacturing, or three-dimensional (3-D) printing, is a rapidly maturing technology wherein components are built by selectively adding feedstock material at locations specified by a computer model. The majority of these technologies use thermally driven phase change mechanisms to convert the feedstock into functioning material. As the molten material cools and solidifies, the component is subjected to significant thermal gradients, generating significant internal stresses throughout the part (Fig. 2). As layers are added, inherent residual stresses cause warping and distortions that lead to geometrical differences between the final part and the original computer generated design. This effect also limits geometries that can be fabricated using AM, such as thin-walled, high-aspect- ratio, and overhanging structures. Distortion may be minimized by intelligent toolpath planning or strategic placement of support structures, but these approaches are not well understood and often "Edisonian" in nature. Residual stresses can also impact component performance during operation. For example, in a thermally cycled environment such as a high-pressure turbine engine, residual stresses can cause components to distort unpredictably. Different thermal treatments on as-fabricated AM components have been used to minimize residual stress, but components still retain a nonhomogeneous stress state and/or demonstrate a relaxation-derived geometric distortion. Industry, federal laboratory, and university collaboration is needed to address these challenges and enable the U.S. to compete in the global market. Work is currently being conducted on AM technologies at the ORNL Manufacturing Demonstration Facility (MDF) sponsored by the DOE's Advanced Manufacturing Office. The MDF is focusing on R&D of both metal and polymer AM pertaining to in-situ process monitoring and closed-loop controls; implementation of advanced materials in AM technologies; and demonstration, characterization, and optimization of next-generation technologies. ORNL is working directly with industry partners to leverage world-leading facilities in fields such as high performance computing, advanced materials characterization, and neutron sciences to solve fundamental challenges in advanced manufacturing. Specifically, MDF is leveraging two of the world's most advanced neutron facilities, the HFIR and SNS, to characterize additive manufactured components.

Watkins, Thomas↗

NASA NTRS1987

A model El Nino-Southern Oscillation

A coupled atmosphere-ocean model is developed and used to study the ENSO (El Nino/Southern Oscillation) phenomenon. With no anomalous external forcing, the coupled model reproduces certain key features of the observed phenomenon, including the recurrence of warm events at irregular intervals with a preference for three to four years. It is shown that the mean sea surface temperature, wind and ocean current fields determine the characteristic spatial structure of ENSO anomalies. The tendency for phase-locking of anomalies is explained in terms of a variation in coupling strength associated with the annual cycle in the mean fields. Sensitivity studies reveal that both the amplitude and the time scale of the oscillation are sensitive to several parameters that affect the strength of the atmosphere-ocean coupling. Stronger coupling implies larger oscillations with a longer timescale. A critical element of the model oscillation is the variability in the equatorial heat content of the upper ocean. Equatorial heat content increases prior to warm events and decreases sharply during the events. A theory for this variability and the associated transitions between non-El Nino and El Nino states is presented. Implications of the model results for the prediction of El Nino events are discussed.

Zebiak, Stephen E.↗

NASA NTRS1994

Transverse shear stresses and their sensitivity coefficients in multilayered composite panels

A computational procedure is presented for the accurate determination of transverse shear stresses and their sensitivity coefficients in flat multilayered composite panels subjected to mechanical and thermal loads. The sensitivity coefficients measure the sensitivity of the transverse shear stresses to variations in the different lamination and material parameters of the panel. The panel is discretized by using either a three-field mixed finite element model based on a two-dimensional first- order shear deformation plate theory or a two-field degenerate solid element with each of the displacement components having a linear variation throughout the thickness of the laminate. The evaluation of transverse shear stresses can be conveniently divided into two phases. The first phase consists of using a superconvergent recovery technique for evaluating the in-plane stresses in the different layers. In the second phase, the transverse shear stresses are evaluated by using piecewise integration, in the thickness direction, of the three-dimensional equilibrium equations. The same procedure is used for evaluating the sensitivity coefficients of the transverse shear stresses. The effectiveness of the computational procedure is demonstrated by means of numerical examples of multilayered cross-ply panels subjected to transverse loading, uniform temperature change, and uniform temperature gradient through the thickness of the panel. In each case the standard of the comparison is taken to be the exact solution of the three dimensional thermoelasticity equations of the panel.

Noor, Ahmed K.↗

DOE OSTI2022

Coexistence of non-Abelian chiral spin liquid and magnetic order in a spin-1 antiferromagnet

We study the ground-state properties of a spin-1 Heisenberg model on a square lattice with the first- and second-nearest-neighbor antiferromagnetic couplings $J_1$ and $J_2$ and a three-spin scalar chirality term $J_\chi$. Using the density matrix renormalization group calculation, we map out a global phase diagram including various magnetic order phases and an emergent quantum spin liquid phase. The nature of the spin liquid is identified as a bosonic non-Abelian Moore-Read state from the fingerprint of the entanglement spectra and identification of a full set of topological sectors. We further unveil a stripe magnetic order coexisting with this spin liquid. Furthermore, our results not only establish a rare example of non-Abelian spin liquids in simple spin systems but also demonstrate the coexistence of fractionalized excitations and magnetic order beyond mean-field descriptions.

2-dimensional systems↗

NASA NTRS2011

Saturn Ring Data Analysis and Thermal Modeling

CIRS, VIMS, UVIS, and ISS (Cassini's Composite Infrared Specrtometer, Visual and Infrared Mapping Spectrometer, Ultra Violet Imaging Spectrometer and Imaging Science Subsystem, respectively), have each operated in a multidimensional observation space and have acquired scans of the lit and unlit rings at multiple phase angles. To better understand physical and dynamical ring particle parametric dependence, we co-registered profiles from these three instruments, taken at a wide range of wavelengths, from ultraviolet through the thermal infrared, to associate changes in ring particle temperature with changes in observed brightness, specifically with albedos inferred by ISS, UVIS and VIMS. We work in a parameter space where the solar elevation range is constrained to 12 deg - 14 deg and the chosen radial region is the B3 region of the B ring; this region is the most optically thick region in Saturn's rings. From this compilation of multiple wavelength data, we construct and fit phase curves and color ratios using independent dynamical thermal models for ring structure and overplot Saturn, Saturn ring, and Solar spectra. Analysis of phase curve construction and color ratios reveals thermal emission to fall within the extrema of the ISS bandwidth and a geometrical dependence of reddening on phase angle, respectively. Analysis of spectra reveals Cassini CIRS Saturn spectra dominate Cassini CIRS B3 Ring Spectra from 19 to 1000 microns, while Earth-based B Ring Spectrum dominates Earth-based Saturn Spectrum from 0.4 to 4 microns. From our fits we test out dynamical thermal models; from the phase curves we derive ring albedos and non-lambertian properties of the ring particle surfaces; and from the color ratios we examine multiple scattering within the regolith of ring particles.

Saturn's rings↗

DOE OSTI2023

Multiband Fermi surface in $1T–$$\mathrm{VSe_2}$ and its implication for the charge density wave phase

Angle-resolved photoemission spectroscopy experiments reveal a surprisingly richer surface electronic structure in $1T–$$\mathrm{VSe_2}$ than previously predicted or probed. Earlier claims supporting a charge density wave phase in this material are reexamined in terms of these findings and are found to be untenable. Here, the Fermi surface is found to be gapless, while band warping effects, currently attributed to three-dimensional lattice distortion, result from the simultaneous dispersion of the closely lying multiple bands. Based on these findings, a charge density wave scenario in $1T–$$\mathrm{VSe_2}$ is unlikely. On the other side, the presence of multiple states crossing the Fermi level should constitute relevant constraints for any viable microscopic model of the structural phase transition of $\mathrm{VSe_2}$.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Elucidating Mass Transport Regimes in Gas Diffusion Electrodes for CO 2 Electroreduction

Gas diffusion electrodes (GDEs) have shown promising performance for the electrochemical reduction of CO 2 (CO 2 R). In this study, a resolved, pore scale model of electrochemical reduction of CO 2 within a liquid-filled catalyst layer is developed. Three CO 2 mass transport regimes are identified in which the CO 2 penetration depth is controlled by CO 2 consumption in the electrolyte, CO 2 conversion along the solid-electrolyte double-phase boundaries (DPBs), and CO 2 conversion concentrated around the gas–solid–electrolyte triple-phase boundaries (TPBs). While it is possible for CO 2 R to be localized around the TPBs, in systems with submicron pore radii operating at <1 A cm –2 CO 2 R will be distributed across the DPBs within the catalyst layer. This validates the assumption of pore-scale uniformity implicit in popular, volume-averaged GDE models. Overall, the CO 2 conversion efficiency depends strongly on the governing mass transport regime, and operational-phase diagrams are constructed to guide the catalyst layer design.

10 SYNTHETIC FUELS↗

NASA NTRSDate not supplied

Multi-Scale Thermo-Mechanical Modeling of Porous 3D Woven TPS Materials

This work summarizes the process to compute and analyze the thermal conductivity and mechanical properties of 3D woven TPS materials such as 3MDCP (3D Mid-Density Carbon-Phenolic). The PuMA software, developed at NASA Ames, was used to characterize different 3MDCP samples from their constituents' data, averaging their thermal conductivity and elastic properties in the three main directions to obtain their effective orthotropic thermo-mechanical properties. This was performed in multiple steps: firstly, TPS samples were digitally reconstructed using micro computed tomography (µCT) and their constituents were segmented; then, the porous matrix phase was analyzed at the micro-scale and these results were used, along with the fibers’ constituents information, to model the tows at the meso-scale; finally, results from the constituents were homogenized and used to model the thermo-mechanical behavior of the unit cell at the macro-scale. Additionally, to gain a better understanding of the tows’ morphology and distribution of the carbon-phenolic blended fibers, microscopy images of the tows’ cross-section were segmented using deep learning techniques and analyzed with PuMA. This workflow can be applied to any TPS material to computationally obtain its thermo-mechanical properties, enabling more informed TPS design and manufacturing choices.

Multi-Scale↗

DOE OSTI2024

Emergence of the isotropic Kitaev honeycomb lattice $α$–RuCl 3 and its magnetic properties

Here, we present a comprehensive investigation of the crystal and magnetic structures of the van der Waals antiferromagnet $α$–RuCl 3 using single crystal x-ray and neutron diffraction. The crystal structure at room temperature is a monoclinic (C2/m). However, with decreasing temperature, a remarkable first-order structural phase transition is observed, leading to the emergence of a rhombohedral ($R\bar3$) structure characterized by three-fold rotational symmetry forming an isotropic honeycomb lattice. On further cooling, a zigzag-type antiferromagnetic order develops below $T_{\textrm{N}} = 6\sim6.6$ K. The critical exponent of the magnetic order parameter was determined to be $\beta = 0.11(1)$, which is close to the two-dimensional Ising model. Additionally, the angular dependence of the magnetic critical field of the zigzag antiferromagnetic order for the polarized ferromagnetic phase reveals a six-fold rotational symmetry within the ab–plane. These findings reflect the symmetry associated with the Ising-like bond-dependent Kitaev spin interactions and underscore the universality of the Kitaev interaction-dominated antiferromagnetic system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

A space-time adaptive finite element method with exponential time integrator for the phase field model of pitting corrosion

In this paper we propose a space-time adaptive finite element method for the phase field model of pitting corrosion, which is a parabolic partial differential equation system consisting of a phase variable and a concentration variable. A major challenge in solving this phase field model is that the problem is very stiff, which makes the time step size extremely small for standard temporal discretizations. Another difficulty is that a high spatial resolution is required to capture the steep gradients within the diffused interface, which results in very large number of degrees of freedom for uniform meshes. To overcome the stiffness of this model, we combine the Rosenbrock–Euler exponential integrator with Crank–Nicolson scheme for the temporal discretization. Moreover, by exploiting the fact that the speed of the corroding interface decreases with time, we derive an adaptive time stepping formula. For the spatial approximation, we propose a simple and efficient strategy to generate adaptive meshes that reduces the computational cost significantly. Thus, the proposed method utilizes local adaptivity and mesh refinement for efficient simulation of the corrosive dissolution over long times in heterogeneous media with complex microstructures. We also present an extensive set of numerical experiments in both two and three dimensional spaces to demonstrate efficiency and robustness of the proposed method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗

NASA NTRSDate not supplied

Multi-Scale Thermo-Mechanical Modeling of Porous 3D Woven TPS Materials

This work summarizes the process to compute and analyze the thermal conductivity and mechanical properties of 3D woven TPS materials such as 3MDCP (3D Mid-Density Carbon-Phenolic). The PuMA [1,2] software, developed at NASA Ames, was used to characterize different 3MDCP samples from their constituents' data, averaging their thermal conductivity and elastic properties in the three main directions to obtain their effective orthotropic thermo-mechanical properties. This was performed in multiple steps: firstly, TPS samples were digitally reconstructed using micro computed tomography (µCT) and their constituents were segmented; then, the porous matrix phase was analyzed at the micro-scale and these results were used, along with the fibers’ constituents information, to model the tows at the meso-scale; finally, results from the constituents were homogenized and used to model the thermo-mechanical behavior of the unit cell at the macro-scale. Additionally, to gain a better understanding of the tows’ morphology and distribution of the carbon-phenolic blended fibers, microscopy images of the tows’ cross-section were segmented using deep learning techniques and analyzed with PuMA. This workflow can be applied to any TPS material to computationally obtain its thermo-mechanical properties, enabling more informed TPS design and manufacturing choices.

Multi-Scale↗

NASA NTRSDate not supplied

A Methodology for Actuating RCS Jets in a Continuous, Time-Accurate CFD Simulation

Computational fluid dynamics simulations of Mars entry vehicle reaction control systems have been used to gain a better understanding of the wide range of flow phenomena encountered during various phases of flight. One aspect of this is simulating the system during controlled flight. Current practices model the system in a continuously actuated state, but this methodology is unable to accurately model the dynamic behavior of the system during EDL. One method to improve this is to create a coupled CFD-RBD-Control simulation. The first step to this is establishing a method of accurately actuating an RCS jet during a time accurate CFD simulation. This paper presents a method to achieve this actuation. There are three phases to the methodology: establishing the off state, creating a ramping function that will initialize the jet flow, and establishing the on state. A test case of a jet in a supersonic crossflow is implemented to establish the success of the methodology in being able to actuate a jet in the midst of an ongoing CFD simulation. It shows significant improvement in boundary condition responsiveness for linear and exponential ramping functions in comparison to the step function. A preliminary demonstration of the methodology is presented for the MSL vehicle which shows significant advantages for the implementation of a linear ramping function over a step function.

Alexandra Hickey↗

NASA NTRSDate not supplied

A Methodology for Actuating RCS Jets in a Continuous, Time-Accurate CFD Simulation

Alexandra Hickey↗

DOE OSTI2021

Robust d-Wave Superconductivity in the Square-Lattice t–J Model

Unravelling competing orders emergent in doped Mott insulators and their interplay with unconventional superconductivity is one of the major challenges in condensed matter physics. Here, to explore the possible superconducting state in a doped Mott insulator, we study the square-lattice t-J model with both the nearest-neighbor and next-nearest-neighbor electron hoppings and spin interactions. By using the state-of-the-art density matrix renormalization group calculation with imposing charge U(1) and spin SU(2) symmetries on the six-leg cylinders, we establish a quantum phase diagram including three phases: a stripe charge density wave phase, a superconducting phase without static charge order, and a superconducting phase coexistent with a weak charge stripe order. Crucially, we demonstrate that the superconducting phase has a power-law pairing correlation that decays much slower than the charge density and spin correlations, which is a quasi-1D descendant of the uniform d-wave superconductor in two dimensions. These findings reveal that enhanced charge and spin fluctuations with optimal doping is able to produce robust d-wave superconductivity in doped Mott insulators, providing a foundation for connecting theories of superconductivity to models of strongly correlated systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Projected Changes in Temperature and Precipitation Over the United States, Central America, and the Caribbean in CMIP6 GCMs

The Coupled Model Intercomparison Project Phase 6 (CMIP6) dataset is used to examine projected changes in temperature and precipitation over the United States (U.S.), Central America and the Caribbean. The changes are computed using an ensemble of 31 models for three future time slices (2021–2040, 2041–2060, and 2080–2099) relative to the reference period (1995–2014) under three Shared Socioeconomic Pathways (SSPs; SSP1-2.6, SSP2-4.5, and SSP5-8.5). The CMIP6 ensemble reproduces the observed annual cycle and distribution of mean annual temperature and precipitation with biases between – 0.93 and 1.27 °C and – 37.90 to 58.45%, respectively, for most of the region. However, modeled precipitation is too large over the western and Midwestern U.S. during winter and spring and over the North American monsoon region in summer, while too small over southern Central America. Temperature is projected to increase over the entire domain under all three SSPs, by as much as 6 °C under SSP5-8.5, and with more pronounced increases in the northern latitudes over the regions that receive snow in the present climate. Annual precipitation projections for the end of the twenty-first century have more uncertainty, as expected, and exhibit a meridional dipole-like pattern, with precipitation increasing by 10–30% over much of the U.S. and decreasing by 10–40% over Central America and the Caribbean, especially over the monsoon region. Seasonally, precipitation over the eastern and central subregions is projected to increase during winter and spring and decrease during summer and autumn. Over the monsoon region and Central America, precipitation is projected to decrease in all seasons except autumn. The analysis was repeated on a subset of 9 models with the best performance in the reference period; however, no significant difference was found, suggesting that model bias is not strongly influencing the projections.

54 ENVIRONMENTAL SCIENCES↗