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

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗
DOE OSTI2025

Long-lived zone-boundary magnons in an antiferromagnet

Antiferromagnetic (AFM) insulators exhibit many desirable features for spintronic applications such as fast dynamics in the THz range and robustness to fluctuating external fields. However, large damping typically associated with THz magnons presents a serious challenge for THz magnonic applications. Here, we report long-lived short-wavelength zone boundary magnons in the honeycomb AFM insulator CoTiO 3 , recently found to host topological magnons. We find that its zone-boundary THz magnons exhibit longer lifetimes than its zone-center magnons. This unusual momentum-dependent long magnon lifetime originates from several factors including the antiferromagnetic order, exchange anisotropy, a finite magnon gap, and magnon band dispersion. Our work suggests that magnon-magnon interaction may not be detrimental to magnon lifetimes and should be included in future searches for topological magnons.

42 ENGINEERING↗
DOE OSTI2020

Axial magnetic flux amplification in Hall-magnetohydrodynamic simulations of externally magnetized z-pinches

It is shown that externally magnetized gas-puff z-pinches under compression can exhibit robust axial flux amplification under suitable axial boundary conditions. This effect relies upon the Hall term in the generalized Ohm's law to generate azimuthal currents in the presence of a driving axial electric field. Under dynamic compression, the total current tends to flow in a mostly force-free boundary layer, separating the predominately azimuthal field outside the layer and the predominately axial field inside the layer. The effect only appears to occur if the axial boundaries allow for outflow or absorption of electromagnetic energy. The effect is mitigated by imposing either periodic or conducting axial boundary conditions. A semi-analytic equilibrium analysis agrees with steady-state solutions of the time-dependent electron-magnetohydrodynamic equations and provides an estimate of the scaling of the boundary layer as well as suggesting a scenario for the formation of the boundary layer. When operative, the effect can significantly impede plasma compression due to the increase in axial flux that diffuses into the pre-compressed plasma or through the presence of conductors that inhibit movement of the generated axial flux. Several facilities have noted unusual implosion behavior in z-pinch experiments with applied axial fields that does not appear to be explainable within the standard magnetohydrodynamic model. In conclusion, it is suggested that these experiments can be explained by the axial flux amplification and concentration effect

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Delamination-Assisted Ultrafast Wrinkle Formation in a Freestanding Film

Freestanding films provide a versatile platform for materials engineering thanks to additional structural motifs not found in films with a substrate. A ubiquitous example is wrinkles, yet little is known about how they can develop over as fast as a few picoseconds due to a lack of experimental probes to visualize their dynamics in real time on the nanoscopic scale. Here, we use time-resolved electron diffraction to directly observe light-activated wrinkling formation in freestanding La 2/3 Ca 1/3 MnO 3 films. Via a “lock-in” analysis of oscillations in the diffraction peak position, intensity, and width, we quantitatively reconstructed how wrinkles develop on the time scale of lattice vibration. Contrary to the common assumption of fixed boundary conditions, we found that wrinkle development is associated with ultrafast delamination at the film boundaries. Finally, our work provides a generic protocol to quantify wrinkling dynamics in freestanding films and highlights the importance of the film–substrate interaction in determining the properties of freestanding structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Spatially Mapping the CO 2 Alkaline Sorbent Diffuse Microenvironment Using Operando Raman Spectroscopy

When designing a chemical process, the local balance of transport and kinetics, collectively referred to as the diffuse microenvironment, plays a critical role in performance but is difficult to directly observe. This work demonstrates a method of two-dimensional spatial chemical mapping of the diffuse microenvironment in the context of alkali metal hydroxide direct air capture of carbon dioxide using a custom operando gas-absorption flow cell along with confocal Raman spectroscopy. Notably, we observe the concentration boundary layer near the gas–liquid interface and elucidate the interplay of carbonate and bicarbonate ions within it while inferring local hydroxide depletion through continuum modeling. These first of their kind observations provide a technique to compare the performance of direct air capture solvents based on diffuse microenvironment dynamics while also providing metrics important for air contactor design such as boundary layer thickness. Overall, this work showcases a new experimental platform to study interfacial diffuse microenvironments in and outside of the field of direct air capture of carbon dioxide.

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

Numerical case study of the aerosol–cloud interactions in warm boundary layer clouds over the eastern North Atlantic with an interactive chemistry module

The presence of warm boundary layer stratiform clouds over the eastern North Atlantic (ENA) region is commonly influenced by the Azores High, especially during the summer season. To investigate comprehensive aerosol–cloud interactions, this study employs the Weather Research and Forecasting model coupled with a chemistry component (WRF-Chem), incorporating aerosol chemical components that are relevant to the formation of cloud condensation nuclei (CCN) and accounting for aerosol spatiotemporal variation. This study focuses on aerosol indirect effects, particularly the long-range transport of aerosols, in the ENA region under three different weather regimes: a ridge with a surface high-pressure system, a post-trough with a surface high-pressure system, and a weak trough. The WRF-Chem simulations conducted at a near-large-eddy scale offer valuable insights into the model's performance, especially in terms of its ability to use high spatial resolution to capture mesoscale cloud features across various weather regimes. Our result shows that introducing 5 times more aerosols to either non-precipitating or precipitating clouds significantly increases ambient CCN numbers, resulting in, to varying degrees, higher liquid water path (LWP) values. The substantial aerosol–cloud interaction especially occurs in the precipitating clouds and demonstrates the susceptibility of the LWP to changes in CCN under different regimes. Conversely, thin, non-rain clouds at the edges of a cloud system are prone to evaporation, exhibiting an aerosol drying effect. The aerosols released during this process transition back to the accumulation mode, facilitating future activation. This dynamic behavior is not adequately represented in prescribed-aerosol simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Method of information entropy for convergence assessment of molecular dynamics simulations

The lack of a reliable method to evaluate the convergence of molecular dynamics simulations has contributed to discrepancies in different areas of molecular dynamics. Here, the method of information entropy is introduced to molecular dynamics for stationarity assessment. The Shannon information entropy formalism is used to monitor the convergence of the atom motion to a steady state in a continuous spatial domain and is also used to assess the stationarity of calculated multidimensional fields such as the temperature field in a discrete spatial domain. It is demonstrated in this work that monitoring the information entropy of the atom position matrix provides a clear indicator of reaching steady state in radiation damage simulations, non-equilibrium molecular dynamics thermal conductivity computations, and simulations of Poiseuille and Couette flow in nanochannels. A main advantage of the present technique is that it is non-local and relies on fundamental quantities available in all molecular dynamics simulations. Unlike monitoring average temperature, the technique is applicable to simulations that conserve total energy such as reverse non-equilibrium molecular dynamics thermal conductivity computations and to simulations where energy dissipates through a boundary as in radiation damage simulations. The method is applied to simulations of iron using the Tersoff/ZBL splined potential, silicon using the Stillinger–Weber potential, and to Lennard–Jones fluid. Its applicability to both solids and fluids shows that the technique has potential for generalization to other areas in molecular dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Dynamic domain kinematic modelling for predicting interflow over leaky impeding layers

Traditional Boussinesq or kinematic simulations of interflow (i.e., lateral subsurface flow) assume no leakage through the impeding layer and require a no-flow boundary condition at the ridge top. However, recent analyses of many interflow-producing landscapes indicate that leaky impeding layers are common, that most interflow percolates well before reaching the toe slope, and therefore, the downslope contributing length is shorter than the hillslope length. In watersheds characterised by perched interflow over a low conductivity layer through permeable topsoil, interflow with percolation may be modelled with a kinematic wave model using a mobile upslope boundary condition defining the hillslope portion contributing interflow to valleys. Here, we developed and applied a dynamic interflow model to simulate interflow using a downslope travel distance concept such that only the active contributing length is modelled at any time. The model defines a variable active area based on the depth of the perched layer, the topographic slope and the ratio of the hydraulic conductivity of topsoil to that of the impeding layer. It incorporates a two-layer soil moisture accounting water balance analysis, a pedo-transfer function, and percolation and evaporation routines to predict interflow rates in continuous and event-based scenarios. We tested the modelling concept on two sets of data (2-year dataset of rainfall observations for the continuous simulation and a multi-day irrigation experiment for the event simulation) from a 121-m-long open interflow collection trench on an experimental hillslope at the Savannah River Site, South Carolina. The continuous model simulation partially represented the observed interflow hydrograph and perched water depth in the experimental hillslope with correlation coefficients of 0.85 and 0.35, respectively. Model performance improved significantly at event-scale analysis. The modelling approach realistically represents interflow dynamics in hillslopes with leaky impeding layers and can be integrated into catchment-scale hydrology models for more detailed hillslope process modelling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The influence of boundary and edge-plasma modeling in computations of axisymmetric vertical displacement

A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

New insights into the mechanism of ultrasonic atomization for the production of metal powders in additive manufacturing

Ultrasonic atomization is one of the promising technologies for producing metal powders for additive manufacturing, where precise control of particle size and morphology is essential. In this study, we coupled an ultrasonic transducer with a carbon fiber plate and atomized liquid droplets and films under different vibration amplitudes. Water, glycerol, and pure aluminum melt were used to study the atomization mechanism and the resulting droplet/powder characteristics, respectively. High-speed optical and ultrafast synchrotron X-ray imaging were used to study in situ the ultrasonic atomization dynamics, including pulsation and clustering of cavities inside the liquid layer/films, development of capillary waves, and formation of liquid droplets. For the first time, we observed and captured the occurrence of cavitation during the atomization of resting drops, films and impact droplets. The inertial cavitation events interfered with the capillary waves across the interphase boundary, puncturing and breaking the boundary to produce atomized mist. The in situ observation revealed the intricate dynamics of ultrasonic atomization and underscored the pivotal role of cavitation events throughout the entire atomization process. We also conducted experiments on ultrasonic atomization of liquid aluminum, producing particles of perfectly spherical shape. The particle size tended to decrease with reduced vibration amplitude Our work has demonstrated the important processing strategies on how to tailor the particle size while ensuring consistent particle shape and morphology, which is the key processing capability for producing high quality powders for additive manufacturing applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strain-Induced Growth of Twisted Bilayers during the Coalescence of Monolayer MoS 2 Crystals

Tailoring the grain boundaries (GBs) and twist angles between two-dimensional (2D) crystals are two crucial synthetic challenges to deterministically enable envisioned applications such as moiré excitons, emerging magnetism, or single-photon emission. Here, we reveal how twisted 2D bilayers can be synthesized from the collision and coalescence of two growing monolayer MoS 2 crystals during chemical vapor deposition. The twisted bilayer (TB) moiré angles are found to preserve the misorientation angle (θ) of the colliding crystals. The shapes of the TB regions are rationalized by a kink propagation model that predicts the GB formed by the coalescing crystals. Optical spectroscopy measurements reveal a θ-dependent long-range strain in crystals with stitched grain boundaries and a sharp (θ > 20°) threshold for the appearance of TBs, which relieves this strain. Reactive molecular dynamics simulations explain this strain from the continued growth of the crystals during coalescence due to the insertion of atoms at unsaturated defects along the GB, a process that self-terminates when the defects become saturated. The simulations also reproduce atomic-resolution electron microscopy observations of faceting along the GB, which is shown to arise from the growth-induced long-range strain. These facets align with the axes of the colliding crystals to provide favorable nucleation sites for second-layer growth of a TB with twist angles that preserve the misorientation angle θ. This interplay between strain generation and aligned nucleation provides a synthetic pathway for the growth of TBs with deterministic angles.

2D materials↗
DOE OSTI2022

Deep Electric Field Predictions by Drift-Reduced Braginskii Theory with Plasma-Neutral Interactions Based on Experimental Images of Boundary Turbulence

We present 2-dimensional turbulent electric field calculations via physics-informed deep learning consistent with (i) drift-reduced Braginskii theory under the framework of an axisymmetric fusion plasma with purely toroidal field and (ii) experimental estimates of the fluctuating electron density and temperature on open field lines obtained from analysis of gas puff imaging of a discharge on the Alcator C-Mod tokamak. The inclusion of effects from the locally puffed atomic helium on particle and energy sources within the reduced plasma turbulence model are found to strengthen correlations between the electric field and electron pressure. The neutrals are also directly associated with broadening the distribution of turbulent field amplitudes and increasing ${\bf E \times B}$ shearing rates. As a result, this demonstrates a novel approach in plasma experiments by solving for nonlinear dynamics consistent with partial differential equations and data without encoding explicit boundary nor initial conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Ultrafast valley polarization of graphene nanorings

Here, we study theoretically electron dynamics of a graphene nanoring placed in the field of an ultrashort optical pulse. We describe the graphene nanoring within an effective model with infinite mass boundary conditions. For an optical pulse with a duration of just a few femtoseconds, the electron dynamics is coherent and is described by a time-dependent Schrödinger equation. If the optical pulse is circularly polarized, then two valleys of graphene are populated differently, resulting in a finite valley polarization of the system after the pulse. Such a valley polarization is a unique property of graphene nanoscale systems, while for a graphene monolayer, a circularly polarized pulse does not produce any valley polarization. The valley polarization of the graphene nanoring depends on parameters of the system, such as inner and outer radii. With the system's size increasing, the valley polarization monotonically decreases, converging to its zero value for the infinite graphene monolayer.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Kinetics of HCP-BCC Phase Transition Boundary in Magnesium at High Pressure

Under high pressures, many crystalline metals undergo solid–solid phase transformations. In order to accurately model the behavior of materials under extreme loading conditions, it is essential to understand the kinetics of phase transition. Using molecular dynamics simulations, this work demonstrates the feasibility of characterizing the speeds of a moving phase boundary using atomistic simulations employing a suitable empirical potential for single-crystal magnesium. The model can provide temperature- and tensorial stress-dependent velocity of a moving phase boundary as a rate-limiting contribution to the kinetics of phase transformation in continuum codes. Results demonstrate that a nonlinear interaction exists between plasticity and phase transition, facilitating a jump in the velocity of a moving phase boundary, facilitated by activated plastic deformation mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dynamic experiments to study the α - ε phase transition in cerium

The ability to understand and predict the response of matter at extreme conditions requires knowledge of a material’s equation-of-state including the location of phase boundaries, transition kinetics, and the evolution of material strength. Cerium is a material with a complex phase diagram that continues to attract significant scientific interest. Recent dynamic experiments have provided information on the low-pressure γ–α phase transition, sound speed, and Hugoniot data for the higher-pressure α phase, as well as the incipient shock melt transition. Despite these efforts, there are still regions of the phase diagram that are largely unexplored dynamically, including the high-pressure region below the melt boundary. Along a room temperature isotherm, diamond anvil cell data report a transition to the ϵ phase between 13 and 17 GPa. At higher temperatures, similar diamond anvil cell data show significant disagreement regarding the existence, location, and slope of the ϵ-phase boundary. In this work, double-shock loading was used to access the α–ϵ region of the phase diagram to obtain equation-of-state information and to determine the location of the ϵ-phase boundary for shock loading.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Insights into elevated temperature tensile deformation mechanisms and kink banding in additively manufactured tungsten

Additive manufacturing (AM) potentially enables fabrication and repair of plasma facing components (PFCs) of tungsten. However, deformation mechanisms in AM-fabricated tungsten under service-relevant thermomechanical conditions remain unknown majorly due to challenges in achieving crack-free W. This study reveals the deformation mechanisms operative under tension at 800 °C and 1200 °C in the electron beam melting powder bed fusion (EBM-PBF) fabricated W. Textured columnar grains with mixed <001 >/<111 >|| build-direction were observed. Strong anisotropy in tensile properties persisted across test temperatures, and all specimens exhibited extensive deformation-induced banding phenomena. Grain boundary character analysis of the deformed specimens indicated that these deformation bands were kink bands with tilt grain boundaries prominently present at the band-matrix interface. Interestingly, the material within the kink bands rotated toward higher resolved shear stress, providing insights into the origins of kink band formation in body-centered cubic refractory metals. Intragranular misorientation axis analysis revealed that plasticity was dominated by {110} <111 > slip at 800 °C, whereas additional (213)[11⁢̄1] and (112)[11⁢̄1] slip systems activated at 1200 °C. A dense low-angle boundary networks observed near fractured surface indicated the onset of dynamic recrystallization. Results reveal plasticity governing mechanisms in AM-fabricated W under power plant-relevant conditions, and provide insights into kink band attributes in refractory metals.

Mayes, Riley [ORNL] (ORCID:000900089307010X)↗