Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic boundary”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Evaluation of the anisotropic grain boundaries and surfaces of α -U via molecular dynamics

Alloys based on uranium-zirconium are gaining renewed interest as fuels for the Versatile Test Reactor and a number of microreactor designs. Implementing metallic fuel in reactors creates the need for robust descriptive and predictive fuel performance modeling. Here, the current state of metallic fuel performance modeling relies on empirical equations derived from historical experiments, which may be unreliable when applied outside of their temperature, power, and composition phase space. One area where such data is lacking is the irradiation behavior of alpha -U, specifically tearing and porosity formation at the early stages of irradiation. While grain boundaries likely play a key role in this fuel behavior, relatively little is known about grain boundaries in alpha-U. Thus, we evaluate the grain boundary, surface energy, and work of adhesion of alpha -U utilizing molecular dynamics. Symmetric tilt grain boundaries (STGBs) are analyzed with the tilt plane oriented along each major crystallographic axis, for a total of eighty unique grain boundaries. The effect of temperature, tilt plane, and misorientation angle on interfacial energies are analyzed. The interfacial energies typically increase with temperature and there is significant variance as a function of misorientation angle, irrespective of the tilt plane. At 500 K, the average surface energy (1.23 J/m 2 ) is approximately 1.5 times the grain boundary energy (0.79 J/m 2 ), and the work of adhesion is approximately twice the grain boundary energy (1.68 J/m 2 ). Orientations for the likely formation of twins and likely failure planes are identified.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Atomistic Origin of Microsecond Carrier Lifetimes at Perovskite Grain Boundaries: Machine Learning-Assisted Nonadiabatic Molecular Dynamics

The polycrystalline nature of perovskites, stemming from their facile solution-based fabrication, leads to a high density of grain boundaries (GBs) and point defects. However, the impact of GBs on perovskite performance remains uncertain, with contradictory statements found in the literature. We developed a machine learning force field, sampled GB structures on a nanosecond time scale, and performed nonadiabatic (NA) molecular dynamics simulations of charge carrier trapping and recombination in stoichiometric and doped GBs. The simulations reveal long, microsecond carrier lifetimes, approaching experimental data, stemming from charge separation at the GBs and small NA coupling, 0.01–0.1 meV. Stoichiometric GBs exhibit transient trap states, which, however, are not particularly detrimental to the carrier lifetime. Halide dopants form interstitial defects in the bulk, but have a stabilizing influence on the GB structure by passivating undersaturated Pb atoms and reducing the transient trap state formation. On the contrary, excess Pb destabilizes GBs, allowing formation of persistent midgap states that trap charges. Still, the charge carrier lifetime reduces relatively little, because the midgap states decouple from the bands, and charges are more likely to escape back into bands upon a GB structural fluctuation. The atomistic study into the structural dynamics of perovskite GBs and its influence on charge carrier trapping and recombination provides valuable insights into the complex properties of perovskites and the intricate role of GBs in the material performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Atomistic study of silicon alloying in the spallation behavior of nanocrystalline aluminum systems

Here, the effect of concentration of the alloying element silicon (Si) on the shock and spall response of aluminum (Al) is studied using classical molecular dynamics (MD) simulations. Silicon is distributed as individual grains and as grain boundaries and the dynamic response under shock is studied for an impact velocity of 1000 m/s. A general trend observed here is that the spallation is almost exclusively intergranular and limited to Al-Si interfaces. The spall strengths are observed to a) decrease with an increase in alloying of Si as grains and b) show little increase for up to 50% grain boundary segregation. An increase in grain size of the whole system showed negligible differences in spall strength, highlighting the significance of concentration of the alloying species rather than its grain size. These results provide an atomistic understanding of light-weight metallic systems under shock compression and pave the way for designing multiphase metal matrix alloys and composites for defense/armor applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hybrid RANS-LES of the Atmospheric Boundary Layer for Wind Farm Simulations: Preprint

Wind farm simulations often do not accurately represent wake-atmospheric boundary layer (ABL) interactions, blade boundary layer (BL) dynamics, and turbine-turbine interactions. In this work, we use Active Model Split (AMS), a new hybrid Reynolds-Averaged Navier Stokes (RANS)-large eddy simulation (LES) model, which is well suited to capture these effects because the model can (i) accurately simulate the ABL with the Coriolis effect, (ii) is accurate in adverse pressure gradients such as those near wind turbine blades, and (iii) has sufficiently low computational cost to simulate multiple turbines while resolving the blade BL. For simplicity and consistency we develop AMS to be used throughout the domain rather than in a zonal method. We implement our work in the massively parallel flow solver, Nalu-Wind, so that our model can access the compute resources needed for blade-resolved simulations of multiple wind turbines. To accomplish these aims, we modify the baseline AMS by changing the RANS contribution to SST k - omega with a length scale limiter, adding the Coriolis effect, and developing an appropriate wall treatment. We show that AMS of the ABL with the Coriolis effect matches LES reference results better than those obtained with RANS. We describe our plans to add buoyancy effects and wind turbines to our AMS simulations.

atmospheric boundary layer↗
DOE OSTI2023

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Mechanisms of Mobility of Grain Boundaries in Tungsten

The objective of this study is to assess the mechanical robustness of tungsten in first-wall applications for tokamaks. These applications unavoidably involve the challenges of high temperature and radiation damage. Here we report progress on determining the mechanisms governing the mobility of tungsten grain boundaries at high temperature, a key property in the process of recrystallization. The determination of mechanisms enables predictive analytic models of grain boundary response to driving forces like the energy stored in the damaged metal lattice that is released as a grain boundary sweeps through it. In prior work we studied tungsten grain boundary properties such as the energies of boundaries with various orientations and misorientations, including energies of grain boundaries whose high-temperature structure differs from the structure at ambient conditions. We have calculated properties associated with the grain boundary mobility including the response of grain boundaries to applied shear stress. We have also developed methods for direct calculation of the motion of curved grain boundaries at high temperature using molecular dynamics. Here we again utilize molecular dynamics simulation of grain boundary motion and go further to analyze how the defect content in a grain boundary evolves as the boundary moves under the action of a driving force.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Upon the horizon’s verge: Thermal particle creation between and approaching horizons

Abstract Quantum particle creation from spacetime horizons, or accelerating boundaries in the dynamical Casimir effect, can have an equilibrium, or thermal, distribution. Using an accelerating boundary in flat spacetime (moving mirror), we investigate the production of thermal energy flux despite non-equilibrium accelerations, the evolution between equilibrium states, and the ‘interference’ between horizons. In particular, this allows us to give a complete solution to the particle spectrum of the accelerated boundary correspondence with Schwarzschild–de Sitter spacetime.

Astronomy & Astrophysics↗
DOE OSTI2022

Spontaneous self-propulsion and nonequilibrium shape fluctuations of a droplet enclosing active particles

Active particles, such as swimming bacteria or self-propelled colloids, spontaneously assemble into large-scale dynamic structures. Geometric boundaries often enforce different spatio-temporal patterns compared to unconfined environment and thus provide a platform to control the behavior of active matter. Here, we report collective dynamics of active particles enclosed by soft, deformable boundary, that is responsive to the particles’ activity. We reveal that a quasi two-dimensional fluid droplet enclosing motile colloids powered by the Quincke effect (Quincke rollers) exhibits strong shape fluctuations with a power spectrum consistent with active fluctuations driven by particle-interface collisions. A broken detailed balance confirms the nonequilibrium nature of the shape dynamics. We further find that rollers self-organize into a single drop-spanning vortex, which can undergo a spontaneous symmetry breaking and vortex splitting. The droplet acquires motility while the vortex doublet exists. Our findings provide insights into the complex collective behavior of active colloidal suspensions in soft confinement.

42 ENGINEERING↗
DOE OSTI2026

Analytical Dynamic Aperture Calculation Compared to Tracking Results

In studying dynamic aperture in the storage rings, it is becoming increasing beneficial to understand the growth of the fractal boundary of the dynamic aperture without tracking for millions of turns in a hadron ring. This note describes the J. Gao theory applied to the Electron Ion Collider (EIC) Hadron Storage Ring (HSR) and compares simulation results to theoretical predictions to show that the analytical method may prove useful in quick determinations of improved dynamic aperture based on sextupole and octupole configurations.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

From boundary data to bound states. Part II. Scattering angle to dynamical invariants (with twist)

We recently introduced in [9] a boundary-to-bound dictionary between gravitational scattering data and observables for bound states of non-spinning bodies. In this paper, we elaborate further on this holographic map. We start by deriving the following — remarkably simple — formula relating the periastron advance to the scattering angle: $ΔΦ(J,ε)=χ(J,ε)+χ(-J,ε)$, via analytic continuation in angular momentum and binding energy. Using explicit expressions from [9], we confirm its validity to all orders in the Post-Minkowskian (PM) expansion. Furthermore, we reconstruct the radial action for the bound state directly from the knowledge of the scattering angle. The radial action enables us to write compact expressions for dynamical invariants in terms of the deflection angle to all PM orders, which can also be written as a function of the PM-expanded amplitude. As an example, we reproduce our result in [9] for the periastron advance, and compute the radial and azimuthal frequencies and redshift variable to two-loops. Agreement is found in the overlap between PM and Post-Newtonian (PN) schemes. Last but not least, we initiate the study of our dictionary including spin. We demonstrate that the same relation between deflection angle and periastron advance applies for aligned-spin contributions, with J the (canonical) total angular momentum. Explicit checks are performed to display perfect agreement using state-of-the-art PN results in the literature. Using the map between test- and two-body dynamics, we also compute the periastron advance up to quadratic order in spin, to one-loop and to all orders in velocity. We conclude with a discussion on the generalized ‘impetus formula’ for spinning bodies and black holes as ‘elementary particles’. Our findings here and in [9] imply that the deflection angle already encodes vast amount of physical information for bound orbits, encouraging independent derivations using numerical and/or self-force methodologies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Evaluation of Atmospheric Boundary Layer Height From Wind Profiling Radar and Slab Models and Its Responses to Seasonality of Land Cover, Subsidence, and Advection

In this study, we evaluated the effect of land cover, atmospheric subsidence, and advection on the annual dynamics of atmospheric boundary layer (ABL) height from two contrasting sites. The first site is the Walker Branch forest, a deciduous forest of temperate climate, complex topography, and cloudy summers. The second site is the Sacramento-San Joaquin River Delta, a site of Mediterranean climate, flat terrain on a local scale, and clear summers. After testing a new algorithm to calculate ABL heights from 915 MHz radar wind profilers, we evaluated a hierarchy of three slab models to recreate the diurnal and annual patterns of ABL growth. We found that the lower ABL heights in the Delta, particularly during late summer, are driven by the combined effects of increased atmospheric subsidence and marine air advection. In both sites, the annual pattern of ABL height was strongly correlated to total daily incoming radiation, and in the Delta, the annual pattern of ABL height closely followed the seasonal patterns of sensible heat flux from a mosaic of different land covers. A land composite of latent and sensible heat fluxes obtained through a meso-network of eddy covariance measurements and the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) mission resulted in higher model skill, thus showing that land cover heterogeneity is an important driver of ABL growth. To conclude, model simulations show that in the Delta, restoring agricultural land to wetlands with large open water areas could result in a reduction of ABL height during those months with low subsidence and advection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluating fine-resolution, regional outputs of a variable resolution global climate model

Climate models have been used to study water resources and regional hydrologic responses to climate change, but climate model outputs must be downscaled to provide relevant regional data. However, the accuracy of this regional data is limited by uncertainties across and within downscaling methods, uncertainty across global outputs, and discontinuities at downscaled boundaries. A new alternative to traditional downscaling is a variable resolution model that incorporates fine-resolution regions directly into a coarse-resolution, global climate simulation in order to capture contiguous dynamics across resolution boundaries. In this study, we used the Variable-Resolution Community Earth System Model (VR-CESM) to generate one-eighth degree (14 km) fine-resolution outputs for the western U.S. and eastern China from 1970-2006. We focus our evaluation on precipitaiton, temperature, snow pack, solar radiation, and wind. We compare the model outputs with remote-sensing-based precipitation data, and both reanalysis and gridded weather station data for precipitation and temperature. VR-CESM generally has a cold bias in winter and a warm bias in summer in the western U.S., which compensate each other to reduce the annual bias. In eastern China, however, the sign of temperature biases are more consistent throughout the year with cold biases in the higher mountains and warm biases throughout most of the rest of the region. Precipitation biases are dependent upon reference data, and show slight overestimation in high mountain regions in both the U.S. and China with respect to gridded weather station data. Simulated snow cover in the western U.S. is reasonable compared to remote sensing data, but snow cover and snow water equivalent have larger biases when compared to reanalysis data. In eastern China there are widespread snow cover biases compared to remote sensing data. VR-CESM underestimates downward shortwave radiation to a greater degree in summer than in winter, and underestimates surface layer windspeed over mountains to a greater degree than in other areas. Comparison between VR-CESM and a coarser simulation (1-degree Beijing Climate Center model) shows reduced precipitation biases in the mountainous regions with finer resolution, indicating the value of variable-resolution modeling for reigonal studies.

Di Vittorio, Alan↗
DOE OSTI2021

A Poisson equation method for prescribing fully developed non-Newtonian inlet conditions for computational fluid dynamics simulations in models of arbitrary cross-section

Prescribing inlet boundary conditions for computational fluid dynamics (CFD) simulations of internal flow in complex geometries such as anatomical vascular models is challenging. In the absence of patient-specific inlet velocity data, a common approach for long blood vessels is to assume that the inlet flow is fully developed. In vessels of irregular cross section, however, prescribing fully developed conditions is complicated due to the lack of a general closed-form analytical solution. In this study, we develop a simple Poisson equation method for prescribing fully developed inlet conditions for the flow of either Newtonian or non-Newtonian fluids in CFD models of arbitrary cross-section. We first derive the generalized Poisson equation for fully developed flow of a non-Newtonian fluid and we then develop and verify a methodology for numerically computing the solution on any planar boundary domain. In addition, we develop a simple extension of the method for prescribing a non-orthogonal inlet velocity that represents fully developed flow from an upstream tube that is connected to the CFD inlet at a non-orthogonal angle. This may be used to investigate a common source of uncertainty in CFD simulations of internal flow that is due to a lack of information concerning the exact streamwise flow direction at the inlets. Comparison to several Newtonian and non-Newtonian benchmark verification solutions shows the method to be extremely accurate. As a practical demonstration case, we use the method to prescribe fully developed conditions on multiple non-circular inlets for the non-Newtonian flow of blood in a patient-specific model of the inferior vena cava (IVC). Finally, we further demonstrate the utility of the method by performing a sensitivity study using the patient-specific IVC model, wherein we investigate the influence of inlet velocity flow direction on the non-Newtonian IVC hemodynamics. Given its simplicity and computational efficiency, the method is shown to be far superior to alternative approaches for prescribing fully developed inlet conditions in such complicated geometries. In conclusion, to facilitate the adoption of our Poisson equation method, we have distributed our OpenFOAM source code and the associated test cases from this study as open-source software.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Performance Analysis and Tank Test Validation of a Hybrid Wave-Current Energy Converter with a Single Power Takeoff

Marine and hydrokinetic (MHK) energy, including ocean waves, tidal current, ocean current and river current, has been recognized as a promising power source due to its full-day availability and high energy potential. At this stage, ocean current energy, tidal energy and ocean wave energy are currently the most competitive sourves among all the categories of MHK. The state of art MHK energy harvesting technology mainly focus on harvesting either ocean wave energy or current energy, but not both. However, a significant amount of ocean waves and tidal/ ocean current coexist in many sites and traditional devices that harvest from a single form of MHK energy, cannot make full use of the coexisting ocean energy. Furthermore, MHK energy harvesting devices need to advance to be cost-effective and competitive with other energy sources. This is difficult to achieve. Ocean wave excitation is irregular, which means that ocean wave height and wave periods are unpredictable and excitation forces on energy harvesting devices can have large variance in amplitude and frequency. Such problems/ restrictions can be possibly addressed by the concept of a hybrid energy converter. In this sense, a hybrid wave-current ocean energy conveter (HWCEC) that simutaneously harvests energy from current and wave with one single power takeoff (PTO) is designed.The wave energy is extracted through relative heaving motion between a floating buoy and a submerged second body, while the current energy is extracted using a marine current turbine (MCT). Energy from both sources are integrated by a hybrid PTO whose concept is based on a mechanical motion rectifier (MMR). In this study, different working modes are investigated together with switching criteria.Simulations were conducted with hydrodynamic coefficients obtained from computational fluid dynamics analysis and boundary element method. Tank tests were conducted for a HWCEC under coexisting wave and current inputs. For comparison, separate baseline tests of a turbine and a two-body point absorber, each acting in isolation, are conducted. Experimental results validate the dynamic modeling and show that a HWCEC can increase the output power with a range between 29-87 percent over either current turbine and wave energy converter acting individually, and it can reduce by up to 70 percent the peak-to-average power ratio compared with the wave energy converter on the tested conditions.Such results demonstrate the potential of the HWCEC as an efficient and cost-effective design.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Measuring dynamics of differentially rotating, unmagnetized, free-boundary plasma produced by soft x-ray irradiation

We present results from a new experiment, fielded on the MAGPIE pulsed power generator, producing a rotating, boundary-free plasma with no pre-imposed magnetic field. Angular momentum is introduced to the rotating system by the oblique collision of multiple plasma flows, which are driven by x-ray ablation of solid targets using the soft x-rays emitted from the implosion of wire array z-pinches. This produces a hydrodynamically stable plasma which undergoes ~2-3 rotations over the duration of the experiment, significantly more than previous pulsed-power platforms. Estimating the angular frequency from the electron density profile in the inner part of the rotating plasma shows that it has a quasi-Keplerian rotation profile. The system also allows for the addition of a controllable magnetic field. This, combined with the sufficiently large Reynolds (~10 5 ) and magnetic Reynolds (~10) numbers, will enable investigation of the effect of magnetic field on the structure and stability of the rotating plasma.

Physics - Plasma physics↗