SEARCH · Engineering Papers
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.
This is the final report on the project, and describes the main accomplishments in the development of the Shifted Boundary Method (SBM), a new computational framework for embedded/immersed computational mechanics applications.
We provide corrections to the slot capacitance and inverse inductance per unit length for slot gasket groove geometries using an approximate conformal mapping approach. We also provide corrections for abrupt step changes in slot width along with boundary discontinuity conditions for implementation in the various slot models.
Accurate measurement of atmospheric turbulent fluctuations is critical for understanding environmental dynamics and improving models in applications such as wind energy. Advanced remote sensing technologies are essential for capturing instantaneous velocity and temperature fluctuations. Acoustic tomography (AT) offers a promising approach that utilizes sound travel times between an array of transducers to reconstruct turbulence fields. This study presents a systematic evaluation of the time-dependent stochastic inversion (TDSI) algorithm for AT using synthetic travel-time measurements derived from large-eddy simulation (LES) fields under both neutral and convective atmospheric boundary-layer conditions. Unlike prior work that relied on field observations or idealized fields, the LES framework provides a ground-truth atmospheric state, enabling quantitative assessment of TDSI retrieval reliability, sensitivity to travel-time measurement noise, and dependence on covariance model parameters and temporal data integration. A detailed sensitivity analysis was conducted to determine the best-fit model parameters, identify the tolerance thresholds for parameter mismatch, and establish a maximum spatial resolution. The TDSI algorithm successfully reconstructed large-scale velocity and temperature fluctuations with root mean square errors ( RMSE s) below 0.35 m/s and 0.12 K, respectively. Spectral analysis established a maximum spatial resolution of approximately 1.4 m, and reconstructions remained robust for travel-time measurement uncertainties up to 0.002 s. These findings provide critical insights into the operational limits of TDSI and inform future applications of AT for atmospheric turbulence characterization and system design.
Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.
The role of turbulence in setting boundary plasma conditions is presently a key uncertainty in projecting to fusion energy reactors. To robustly diagnose edge turbulence, we develop and demonstrate a technique to translate brightness measurements of HeI line radiation into local plasma fluctuations via a novel integrated deep learning framework that combines neutral transport physics and collisional radiative theory for the $3^3 D - 2^3 P$ transition in atomic helium. The tenets for experimental validity are reviewed, illustrating that this turbulence analysis for ionized gases is transferable to both magnetized and unmagnetized environments with arbitrary geometries. Based upon fast camera data on the Alcator C-Mod tokamak, we present the first 2-dimensional time-dependent experimental measurements of the turbulent electron density, electron temperature, and neutral density revealing shadowing effects in a fusion plasma using a single spectral line.
Abstract While numerous modeling studies have focused on the interaction of ocean surface waves with the atmospheric boundary layer, most employ idealized waves that are either monochromatic or synthetically generated from a theoretical wave spectrum, and the atmospheric solvers are typically incompressible. To study wind–wave coupling in real-world scenarios, a model that can simulate both realistic meteorological and wave conditions is necessary. In this paper we describe the implementation of a moving bottom boundary condition into the Weather Research and Forecasting Model for large-eddy simulation applications. We first describe the moving bottom boundary conditions within WRF’s pressure-based vertical coordinate system. We then validate our code with idealized test cases that have analytical solutions, including flow over a monochromatic wave with and without viscosity. Finally, we present results from turbulent flows over a moving monochromatic wave with different wave ages, and demonstrate satisfactory agreement of the wave growth rate with results from the literature. We also compare atmospheric stress and wind parameters from two physically equivalent cases. The first specifies a wind moving in the same direction as a propagating wave, while the second involves a stationary wave with the wind adjusted such that the wind relative to the wave is the same as in the first case. Results indicate that the velocity and Reynolds stress profiles for the two cases match, further validating the moving bottom implementation.
The mixed boundary value problem of a tubular conductor is solved using an approximate representation of its Fourier coefficients. A two term solution is derived, which represents the solution over an extremely broad range of aspect ratios. This representation is used to find the electrostatic solution and capacitance of a charged tube as well as the solution of a tube in a uniform field and its dipole moment. This second case is directly useful as a model for a monopole electric field probe. This approximation is a special case of a representation using a combination of Chebyshev and Legendre polynomials. Combining the charged tube and tube in a uniform field allows the solution of voltage driven tubular antennas. Comparisons are made with numerical solutions using piecewise sinusoidal representations of the current. The results are also generalized to the dynamic case and up to and beyond the first resonance. Simple corrections for finite gap and delta gap drives to magnetic frill drives are examined using infinite tube integral transform representations. Corrections between magnetic frill drives and coaxial drives are also given. Approximate drive corrections using conformal mapping and an effective radius are also discussed. Finally, this efficient current representation is applied to the magnetic problem involving simple tubular solenoids.
Understanding the influence of radiation damage on the mechanical properties of HT-9 and other tempered martensitic alloys is part of the mission in developing radiation-tolerant materials for the next generation of nuclear reactors. Although there has been extensive data on the macroscopic mechanical properties of the irradiated HT-9 alloys and microstructural changes, it is not well-understood how radiation damage and the resulting microstructural changes influence the local mechanical properties. In this study, we utilized in situ SEM microtensile testing to investigate the deformation behavior of specific martensitic boundaries in 1 dpa proton-irradiated HT-9 and answered the question of what is the weakest link. Additionally, we provided the direct observation of the failure modes of microtensiles containing high angle and low angle martensitic boundaries. In the unirradiated condition, the deformation is ductile and no martensitic boundary failure is observed. In the irradiated condition, the high angle martensitic boundaries are more susceptible to radiation-induced boundary failure as compared to the low angle martensitic boundaries.
The classical Neumann boundary value problem of an isotropic, homogeneous elastic half-plane under plane strain conditions is readdressed as the limiting case of the fully three-dimensional problem. Analytical solutions of the stress and strain tensors are obtained by taking the limit from known three-dimensional solutions. It is shown that the displacement fields for the plane strain problem are not well defined. A small number of simple expressions are developed, which provide a general solution for linearly-varying traction over arbitrary regions on the boundary. A simple, efficient, and rapidly convergent algorithm is developed which uses these solutions as analytic elements and provides a solution approach to the general boundary value problem. The method is verified against known solutions for Hertzian contact between parallel cylinders. Two numerical examples are presented for the analysis of shallow foundation systems. In the first, the boundary conditions are informed by analytical elastoplastic calculations and a strain influence analysis is performed and compared with the Schmertmann method. Subsequently, empirical laboratory contact traction distributions measured by Bauer et al., in both the normal and tangential directions are employed as boundary conditions for an analysis of the underlying stress field.
Accurate knowledge of the electrical and thermal conductivities and structural properties of hydrogen–helium mixtures under thermodynamic conditions within and beyond the immiscibility range is very important to predict the thermal evolution and internal structure of gas giant planets like Jupiter and Saturn. Here, we propose a novel method to determine the immiscibility boundary accurately without the need for free energy calculations, while providing consistent insights into structural and transport properties of mixtures. We show with direct large-scale ab initio simulations that the insulator–metal transition (IMT) of the hydrogen subsystem is strongly affected by an admixture with a small fraction of helium and occurs at temperatures significantly higher than those of pure hydrogen. At pressures below 150 GPa, the IMT boundary is not related anymore to the H 2 subsystem dissociation, the system remains insulating even after the full dissociation of H 2 molecules and its transition to an H–He mixture. The offset of the IMT in the H–He mixture relative to the dissociation region in the hydrogen subsystem and the significant reduction of static electrical and thermal conductivity by a factor between two and a few thousand relative to pure hydrogen found in mixtures have consequences for Jupiter and Saturn’s thermal evolution, internal structure, and dynamo action, affecting a large fraction of the interior of both planets.
Enhanced fluctuations, steep gradients, and intensified heat transfer are characteristics of wall-bounded turbulence at transcritical conditions. Although such conditions are prevalent in numerous technical applications, the structure of the thermal boundary layer under realistic density gradients and heating conditions remains poorly understood. Specifically, statistical descriptions of the temperature field in such flows are provided inconsistently using existing models. To address this issue, direct numerical simulations are performed by considering fully developed transcritical turbulent channel flow at pressure and temperature conditions that cause density changes of a factor of up to $O(20)$ between the hot and cold walls. As a consequence of the proximity of the Widom line to the hot wall, significant asymmetries are observed when comparing regions near the cold wall and near the hot wall. Previous transformations that attempt to collapse the near-wall mean temperature profiles among different cases to a single curve are examined. By addressing model deficiencies of these transformations, a formulation for an improved mean temperature transformation is proposed, with appropriate considerations for real fluid effects that involve strong variations in thermodynamic quantities. Our proposed transformation is shown to perform well in collapsing the slope of the logarithmic region to a single universal value with reduced uncertainty. Coupled with a predictive framework to estimate the non-universal shift parameter of the logarithmic region using a priori information, our transformation provides an analytic profile to model the near-wall mean temperature. These results thus provide a framework to guide the development of models for wall-bounded transcritical turbulence.
Grain boundaries (GBs) accommodate the misorientation between adjacent grains in a polycrystalline material. GBs are geometrically described by the macroscopic and microscopic degrees of freedom. Besides, at the atomistic level, GBs exhibit complicated behaviors under varying thermodynamic conditions. The complexity of atomistic GB structures demands stochastic searching for possible states. The effectiveness of stochastic search methods relies on techniques to recreate and select atomistic structures. In this work, we developed a new mutation operator that can induce direct and collective atomistic structure changes to boost the search efficiency of exploring GB structures with evolutionary algorithms (EA). We implemented the mutation methods along with innovative selection, crossover, boundary condition preprocessing methods to form an EA-based package to explore GB structures in grand canonical ensembles with atomistic simulations. We used this package to study the [001] symmetric tilt grain boundaries (STGBs) in FCC copper (Cu), the [110] STGBs in BCC tungsten (W), and the $[1\bar{2}10]$ STGBs in HCP magnesium (Mg). The results show that our design and implementation based on new mutation procedures, selection, and boundary conditions provide a high-quality search of atomistic GB structures in the grand canonical ensemble for different crystal lattices.
It is proposed that a salt-saturated, water-based mud might improve borehole stability for UT-GOM2-2 relative to a fresh-water-based mud. However, the primary objective for UT-GOM2-2 is to safely and successfully acquire uncompromised hydrate-bearing sediment cores. This can only be done if conditions stay within the hydrate stability zone, and well away from the hydrate stability boundary. The addition of salt to the drilling mud, shifts the hydrate stability boundary closer to estimated conditions. This shift, by some estimates of pressure and temperature, where conditions now fall outside of the hydrate stability zone, may result in borehole enlargement and the release of free gas into the borehole. This shift by all estimates, even when coring conditions stay inside the hydrate stability zone, shrinks the window between the estimated conditions and the hydrate stability boundary. This shift is likely to compromise the hydrate-bearing cores while they are being recovered from the bottom-hole to the rig floor. In this report we present two examples: 1) Assuming temperature and pressure from in-situ estimates: Hydrate within the target reservoirs (Orange and Blue) at Walker Ridge Block 313 (WR 313) will be outside the hydrate stability zone in the presence of a drilling mud with 9.5 wt.% salinity (the salinity of the proposed 10.5 ppg salt-based mud) at in situ pressure and temperature. Thus, the hydrate will be unstable. A 10.5 ppg salt-based mud may result in dissociation of the hydrate into its components: water and gas. 10.5 ppg salt-based mud may enlarge the borehole, release free gas into the borehole, and compromise the cores. 2) Assuming temperatures equal to the measured LWD borehole temperatures at this location: Hydrate within the target reservoirs will be stable with a salinity of 9.5 wt.% (the salinity of the proposed 10.5 ppg salt-based mud), but possibly unstable (just at the methane hydrate stability boundary) for a salinity of 17.2 wt.% (the salinity of the proposed 13.5 ppg salt-based kill mud). if the borehole has the same temperature as recorded during previous LWD drilling at this location. A 13.5 ppg salt-based kill mud may destabilize the borehole. More importantly, a 10.5 ppg salt-based mud does not provide a large enough window between the estimated conditions and the hydrate stability boundary. Core temperatures and pressures fluctuate from the bottom-hole conditions as the core is brought up from the bottom-hole to the rig floor. These fluctuations are likely to cause the hydrate in the cores to touch or cross the stability boundary. Therefore, cores captured using a 10.5 ppg salt-based mud are likely to be compromised during core recovery to the rig floor.
We present theoretical and experimental studies of superconductivity and low temperature structural phase boundaries in lithium. We mapped the structural phase diagram of 6Li and 7Li under hydrostatic conditions between 5 – 55 GPa and within the temperature range of 15 – 75 K, observing the FCC → hR1 → cI16 phase transitions. 6Li and 7Li show some differences at the structural boundaries, with a potential shift of the phase boundaries of 6Li to lower pressures. Density functional theory calculations and topological analysis of the electron density elucidates the superconducting properties and interatomic interactions within these phases of lithium.
This large-eddy simulation (LES) study examines how wind-farm-induced atmospheric gravity waves (AGWs) and wind farm performance depend on non-dimensional atmospheric parameters and simulation configuration. A hypothetical aligned wind farm of actuator disks is simulated under neutral surface conditions, with a stable capping inversion and a mildly stable free atmosphere, to assess the effects of stratification beyond the atmospheric boundary layer (ABL) on ABL flow. Simulation set-ups fully resolving AGWs are validated to minimize spurious wave generation and reflection from the domain boundaries. The validated set-up is then used to analyze AGW types and characteristics, as well as stratification impacts under conventionally neutral boundary layer (CNBL) conditions. These conditions are governed by four non-dimensional parameters: the Froude numbers of the free atmosphere and capping inversion (Fr, Fr i ), and the aspect ratios of the ABL and wind farm (H̃ i , S h ). Simulation configurations that fully resolve AGWs – capturing at least one wavelength both horizontally and vertically – yield the most realistic stratification effects on ABL flow, whereas partial or unresolved configurations produce non-physical, channel-like behavior. A coherent description of the AGW phenomena is provided, highlighting the central role of capping inversion displacement in linking ABL fluctuations with AGWs. Trapped waves are confined within the capping inversion, while interfacial and internal waves aloft are identified as the AGW types most relevant to wind farm performance. The wavy inversion, analogous to an interfacial wave, forms converging and diverging zones that drive power fluctuations across the farm. The interfacial wavelength, measured over the wind farm, corresponds to one diverging, one converging, and one mildly diverging zone. As the interfacial wavelength decreases with Fr i , multiple convergence–divergence zones develop under sub-critical conditions (Fr i <1.0), while for super-critical conditions (Fr i > 1.0), the wavelength approaches the farm length. Wave amplitude increases with decreasing H̃ i (i.e., shallower capping inversions). Wind farm performance is most sensitive to H̃ i : shallow boundary layers increase blockage and reduce efficiency, while deeper layers enhance efficiency. Increasing Fr and Fr i mitigates blockage, and increasing S h mainly improves wake recovery. Although local power fluctuations arise from AGWs, overall wind farm efficiency remains nearly constant with Fr and Fr i , improving primarily with larger H̃ i and S h .
FeO is a crucial component of the Earth’s core, and its thermodynamic properties are essential to developing more accurate core models. It is also a notorious correlated insulator in the NaCl-type (B1) phase at ambient conditions. It undergoes two polymorphic transitions at 300 K before it becomes metallic in the NiAs-type (B8) structure at ~100 GPa. Although its phase diagram is not fully mapped, it is well established that the B8 phase transforms to the CsCl-type (B2) phase at core pressures and temperatures. Here, we report a successful ab initio calculation of the B8↔B2 phase boundary in FeO at Earth’s core pressures. We show that fully anharmonic free energies computed with the Perdew–Burke–Ernzerhof-generalized gradient approximation coupled with thermal electronic excitations reproduce the experimental phase boundary within uncertainties at P > 255 GPa, including the largely negative Clapeyron slope of –52 MPa/K. This study validates the applicability of a standard density functional theory functional to FeO under Earth’s core conditions and demonstrates the theoretical framework that enables complex predictive studies of this region.
One distinguished property of the topological insulator (TI) is its robust quantized edge conductance against edge defect. However, this robustness, underlined by the topological principle of bulk-boundary correspondence, is conditioned by assuming a perfect bulk. Here, we investigate the robustness of the TI phase against bulk defects, including vacancy (VA), vacancy cluster (VC), and grain boundary (GB), instead of edge defect. In this work, based on a tight-binding model analysis, we show that a two-dimensional (2D) TI phase, as characterized by a nonzero spin Bott index, will vanish beyond a critical VA concentration ($n^\text{c}_\text{v}$). Generally, $n^\text{c}_\text{v}$ decreases monotonically with the decreasing topological gap induced by spin-orbit coupling. Interestingly, the $n^\text{c}_\text{v}$ to destroy the topological order, namely, the robustness of the TI phase, is shown to be increased by the presence of VCs but decreased by GBs. As a specific example of a large-gap 2D TI, we further show that the surface-supported monolayer Bi can sustain a nontrivial topology up to $n^\text{c}_\text{v}$ ~ 17%, based on a density-functional theory–Wannier-function calculation. Our findings should provide useful guidance for future experimental studies of effects of defects on TIs.
Here, a high-order accurate implicit-mesh discontinuous Galerkin framework for wave propagation in single-phase and bi-phase solids is presented. The framework belongs to the embedded-boundary techniques and its novelty regards the spatial discretization, which enables boundary and interface conditions to be enforced with high-order accuracy on curved embedded geometries. High-order accuracy is achieved via high-order quadrature rules for implicitly-defined domains and boundaries, whilst a cell-merging strategy addresses the presence of small cut cells. The framework is used to discretize the governing equations of elastodynamics, written using a first-order hyperbolic momentum-strain formulation, and an exact Riemann solver is employed to compute the numerical flux at the interface between dissimilar materials with general anisotropic properties. The space-discretized equations are then advanced in time using explicit high-order Runge–Kutta algorithms. Several two- and three-dimensional numerical tests including dynamic adaptive mesh refinement are presented to demonstrate the high-order accuracy and the capability of the method in the elastodynamic analysis of single- and bi-phases solids containing complex geometries.