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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Turbulent properties from Doppler lidar

Permanent and mobile facilities deployed by the ARM program cover a wide range of climatic regions, with a suite of active remote-sensing of cloud and aerosol profiles complemented by comprehensive in-situ measurements of the chemical, physical and optical properties of aerosol measured at the surface. To be able to fully exploit these datasets requires an understanding of the transport from the surface and within the atmosphere, especially within the planetary boundary layer. This requires a method of identifying when the surface and atmosphere are dynamically-linked. Determining the state of the boundary layer, and defining a well-mixed boundary-layer, is therefore suitably described in terms of mixing; here, we will assume that mixing essentially means turbulent mixing. This necessitates accurate measurements of the winds and turbulent mixing in the atmosphere. Doppler lidar is an excellent tool for observing turbulent properties, and ARM has deployed commercially available Doppler lidars at many of the fixed and mobile facilities. Doppler lidars have been shown to provide accurate retrievals of wind and turbulent properties within the boundary layer , however, robust characterization and calculation of the Doppler lidar measurement uncertainties is required before deriving higher level data products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Considerations for testing full-scale wind turbine nacelles with hardware-in-the-loop

In this paper, pole placement control techniques are utilized to compensate for variations in the torsional dynamics resulting from the different boundary conditions experienced by a wind turbine nacelle when mounted on a test bench. Simulation case studies are used to understand the baseline torsional characteristics of the nacelle coupled to the test bench versus that coupled to the rotor. It was found that the frequency shifts by up to 200% and 400% for the first and second torsional modes, respectively. The feasibility of utilizing pole placement technique to match the torsional characteristics of the turbine on the test bench is investigated using PI and PID controllers. Here, the performance of the tuned controller is then verified under two test scenarios: Low Voltage Ride Through (LVRT), and a highly dynamic turbulent wind input. It is demonstrated that tuning the test bench speed controller can effectively shift the poles of the closed-loop system to match the desired first and second modes of the full turbine and replicate transient field events on the test bench.

17 WIND ENERGY↗
DOE OSTI2024

Modeling transient edge plasma transport with dynamic recycling

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Realization of a gas puff imaging system on the Wendelstein 7-X stellarator

A system for studying the spatiotemporal dynamics of fluctuations in the boundary of the W7-X plasma using the “Gas-Puff Imaging” (GPI) technique has been designed, constructed, installed, and operated. This GPI system addresses a number of challenges specific to long-pulse superconducting devices, such as W7-X, including the long distance between the plasma and the vacuum vessel wall, the long distance between the plasma and diagnostic ports, the range of last closed flux surface (LCFS) locations for different magnetic configurations in W7-X, and management of heat loads on the system’s plasma-facing components. The system features a pair of “converging–diverging” nozzles for partially collimating the gas puffed locally ≈135 mm radially outboard of the plasma boundary, a pop-up turning mirror for viewing the gas puff emission from the side (which also acts as a shutter for the re-entrant vacuum window), and a high-throughput optical system that collects visible emission resulting from the interaction between the puffed gas and the plasma and directs it along a water-cooled re-entrant tube directly onto the 8 × 16 pixel detector array of the fast camera. The DEGAS 2 neutral code was used to simulate the Hα (656 nm) and HeI (587 nm) line emission expected from well-characterized gas-puffs of H2 and He and excited within typical edge plasma profiles in W7-X, thereby predicting line brightnesses used to reduce the risks associated with system sensitivity and placement of the field of view. Operation of GPI on W7-X shows excellent signal-to-noise ratios (>100 at 2 Mframes/s) over the field of view for minimally perturbing gas puffs. The GPI system provides detailed measurements of the two-dimensional (radial and poloidal) dynamics of plasma fluctuations in the W7-X edge and scrape-off layer and in and around the magnetic islands outside the LCFS that make up the island divertor configuration employed on W7-X.

Instruments & Instrumentation↗
DOE OSTI2021

Dynamics of a fractal set of first-order magnetic phase transitions in frustrated Lu 2 CoMnO 6

The axial next-nearest-neighbor Ising model predicts a fractal (infinite) set of phases with incommensurate wave vectors that are separated by first-order phase boundaries. This complexity results from a simple frustration condition between nearest- and next-nearest-neighbor interactions along a chain of Ising spins. Using x-ray photon correlation spectroscopy (XPCS), we investigate the surprising antiferromagnetic dynamics that emerge from such a complex phase diagram over a wide range of temperatures. In this work, we present XPCS measurements of the frustrated magnetic chain compound Lu 2 CoMnO 6 and Monte Carlo simulations. Incommensurate magnetic Bragg peaks slide towards commensurate “up-up-down-down” spin order with decreasing temperature and increasing time. Both simulation and experiment support a counterintuitive “upside-down” temperature dependence of the magnetic dynamics: at higher temperatures in the region of first-order phase boundaries, slower dynamics are observed where the speckle maintains its coherence. At the lowest temperatures, where part of the sample adopts commensurate order, the dynamics speed up and result in fast decoherence.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A vorticity-based criterion to characterise leading edge dynamic stall onset

We propose a more conservative, physically-intuitive criterion, namely, the boundary enstrophy flux ( BEF), to characterise leading-edge-type dynamic stall onset in incompressible flows. Our results are based on wall-resolved large-eddy simulations of pitching aerofoils, with fine spatial and temporal resolution around stall onset. We observe that |BEF| reaches a maximum within the stall onset regime identified. By decomposing the contribution to BEF from the flow field, we find that the dominant contribution arises from the laminar leading edge region, due to the combined effect of large clockwise vorticity and favourable pressure gradient. A relatively small contribution originates from the transitional/turbulent laminar separation bubble (LSB) region, due to LSB-induced counter-clockwise vorticity and adverse pressure gradient. This results in BEF being nearly independent of the integration length as long as the region very close to the leading edge is included. This characteristic of BEF yields a major advantage in that the effect of partial or complete inclusion of the noisy LSB region can be filtered out, without changing the BEF peak location in time significantly. Next, we analytically relate BEF to the net wall shear and show that its critical value (=max(|BEF|) ) corresponds to the instant of maximum net shear prevailing at the wall. Finally, we have also compared BEF with the leading edge suction parameter ( LESP) and find that the former reaches its maximum value between 0.3° and 0.8° of rotation earlier.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Measurement and memory in the periodically driven complex Ginzburg-Landau equation

In the present work we illustrate that classical but nonlinear systems may possess features reminiscent of quantum ones, such as memory, upon suitable external perturbation. As our prototypical example, we use the two-dimensional complex Ginzburg-Landau equation in its vortex glass regime. Here, we impose an external drive as a perturbation mimicking a quantum measurement protocol, with a given “measurement rate” (the rate of repetition of the drive) and “mixing rate” (characterized by the intensity of the drive). Using a variety of measures, we find that the system may or may not retain its coherence, statistically retrieving its original glass state, depending on the strength and periodicity of the perturbing field. The corresponding parametric regimes and the associated energy cascade mechanisms involving the dynamics of vortex waveforms and domain boundaries are discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Inside the hologram: reconstructing the bulk observer’s experience

We develop a holographic framework for describing the experience of bulk observers in AdS/CFT, that allows us to compute the proper time and energy distribution measured along any bulk worldline. Our method is formulated directly in the CFT language and is universal: it does not require knowledge of the bulk geometry as an input. When used to propagate operators along the worldline of an observer falling into an eternal black hole, our proposal resolves a conceptual puzzle raised by Marolf and Wall. Notably, the prescription does not rely on an external dynamical Hamiltonian or the AdS boundary conditions and is, therefore, outlining a general framework for the emergence of time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Simulating hindered grain boundary diffusion using the smoothed boundary method

Abstract Grain boundaries can greatly affect the transport properties of polycrystalline materials, particularly when the grain size approaches the nanoscale. While grain boundaries often enhance diffusion by providing a fast pathway for chemical transport, some material systems, such as those of solid oxide fuel cells and battery cathode particles, exhibit the opposite behavior, where grain boundaries act to hinder diffusion. To facilitate the study of systems with hindered grain boundary diffusion, we propose a model that utilizes the smoothed boundary method to simulate the dynamic concentration evolution in polycrystalline systems. The model employs domain parameters with diffuse interfaces to describe the grains, thereby enabling solutions with explicit consideration of their complex geometries. The intrinsic error arising from the diffuse interface approach employed in our proposed model is explored by comparing the results against a sharp interface model for a variety of parameter sets. Finally, two case studies are considered to demonstrate potential applications of the model. First, a nanocrystalline yttria-stabilized zirconia solid oxide fuel cell system is investigated, and the effective diffusivities are extracted from the simulation results and are compared to the values obtained through mean-field approximations. Second, the concentration evolution during lithiation of a polycrystalline battery cathode particle is simulated to demonstrate the method’s capability.

Materials Science↗
DOE OSTI2020

Coordinate transformation methodology for simulating quasistatic elastoplastic solids

Molecular dynamics simulations frequently employ periodic boundary conditions where the positions of the periodic images are manipulated in order to apply deformation to the material sample. For example, Lees-Edwards conditions use moving periodic images to apply simple shear. Here, we examine the problem of precisely comparing this type of simulation to continuum solid mechanics. We employ a hypo-elastoplastic mechanical model, and develop a projection method to enforce quasi-static equilibrium. We introduce a simulation framework that uses a fixed Cartesian computational grid on a reference domain, and imposes deformation via a time-dependent coordinate transformation to the physical domain. As a test case for our method, we consider the evolution of shear bands in a bulk metallic glass using the shear transformation zone theory of amorphous plasticity. We examine the growth of shear bands in simple shear and pure shear conditions as a function of the initial preparation of the bulk metallic glass.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Biophysical model of eelgrass and water quality in Coos Bay, OR shows greater mitigation potential for ocean acidification than hypoxia

Seagrass beds provide important ecosystem services and are valued, in part, for their potential to mediate stressors such as ocean acidification and hypoxia (OAH) for sensitive species. However, the susceptibility of seagrasses to anthropogenic impacts and recent declines motivate the need to better understand the drivers of seagrass and the water quality consequences that occur with variation in seagrass abundance. To meet this need, we leveraged existing monitoring data (water quality and seagrass), hydrodynamic circulation model, and biogeochemical model framework with seagrass submodel, to produce a biophysical model of Coos Bay estuary, Oregon, U.S. The model includes biogeochemical processes involving water quality, plankton, seagrass, and sediment-water interactions. Ecosystem models like this are useful for evaluating complex estuarine systems because they allow us to extend our understanding of system dynamics beyond existing observations and perform experiments to identify the processes driving observed patterns. We used the biophysical model of Coos Bay to evaluate the dynamics of water quality and native eelgrass (Zostera marina) under three eelgrass abundance scenarios (zero eelgrass, current extent, and maximum observed extent) to elucidate the relationship between eelgrass and OAH. Including eelgrass in the Coos Bay model produced results that more closely resembled water quality observations - dissolved oxygen (DO) and pH were more dynamic in simulations with eelgrass, often having both higher highs and lower lows. While there were some areas of the estuary where DO improved with the addition of eelgrass to the model there was overall a small net increase in harmful DO conditions (based on a salmon physiological threshold). In contrast, ocean acidification conditions, pH and calcium carbonate saturation state for aragonite (Ω), were improved (based on oyster requirements) with the addition of eelgrass - although the magnitude of improvement differed seasonally and spatially. Our new model represents a useful tool - one which accounts for and controls the relevant physical and biogeochemical processes - to evaluate conditions that confer resilience or enhance vulnerability to OAH in an important Pacific Northwest coastal estuary and results can inform the OAH-related dynamics occurring in other eastern boundary current estuaries.

FVCOM-ICM↗
DOE OSTI2022

Development of A Multi-Spectral Pyrometry Sensor for High-Speed Transient Surface-Temperature Measurements in Combustion-Relevant Harsh Environments

Accurate and high-speed transient surface-temperature measurements of combustion devices including internal combustion (IC) engines, gas turbines, etc., provide validation targets and boundary conditions for computational fluid dynamics models, and are broadly relevant to technology advancements such as performance improvement and emissions reduction. Development and demonstration of a multi-infrared-channel pyrometry-based optical instrument for high-speed surface-temperature measurement is described. The measurement principle is based on multi-spectral radiation thermometry (MRT) and uses surface thermal radiation at four discrete spectral regions and a corresponding emissivity model to obtain surface temperature via non-linear least squares (NLLS) optimization. Rules of thumb for specifying the spectral regions and considerations to avoid interference with common combustion products are developed; the impact of these along with linear and non-linear MRT analysis are assessed as a function of temperature and signal-to-noise ratio. A multi-start method to determine the MRT-solution global optimum is described and demonstrated. The resulting multi-channel transient pyrometry instrument is described along with practical considerations including optical-alignment drift, matching intra-channel transient response, and solution-confidence indicators. The instrument demonstrated excellent >97% accuracy and >99% 2-sigma precision over the 400–800 °C range, with ~20 µs (50 kHz, equivalent to 0.2 cad at 2000 RPM IC-engine operation) transient response in the bench validation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Parametric Decay of Alfvénic Wave Packets in Nonperiodic Low-beta Plasmas

The parametric decay of finite-size Alfvén waves in nonperiodic low-beta plasmas is investigated using one-dimensional (1D) hybrid simulations. Compared with the usual small periodic system, a wave packet in a large system under the absorption boundary condition shows different decay dynamics, including reduced energy transfer, localized density cavitation, and ion heating. The resulting Alfvén wave dynamics are influenced by several factors relating to this instability, including the growth rate, central wave frequency, and unstable bandwidth. A final steady state of the wave packet may be achieved when the instability does not have enough time to develop within the residual packet, and the packet size shows well-defined scaling dependencies on the growth rate, wave amplitude, and plasma beta. Under the proper conditions, enhanced secondary decay can also be excited in the form of a narrow, amplified wave packet. These results may help to interpret laboratory and spacecraft observations of Alfvén waves, and to refine our understanding of the associated energy transport and ion heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

RELAP5-3D in the OECD-NEA HTGR Thermal Hydraulics Benchmark

The OECD-NEA HTGR Thermal Hydraulics benchmark uses data from the High Temperature Test Facility for a series of code-to-code and code-to-data exercises aimed at improving the state of knowledge on existing thermal hydraulics modeling and simulation tools for prismatic HTGR applications. This benchmark includes problems representing hot gas mixing in the lower plenum, the depressurized conduction cooldown accident and the pressurized conduction cooldown accident. This presentation discusses the benchmark and RELAP5-3D's role in that benchmark, including as a tool for predicting behavior in the core and as a tool for providing boundary conditions to computational fluid dynamics analysis in the lower plenum.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Probing chiral Kitaev spin liquids via dangling boundary fermions

Identifying experimental probes capable of diagnosing extreme quantum behavior is widely regarded as one of the foremost challenges in modern condensed matter physics. Here, we propose a novel approach for detecting chiral Kitaev spin-liquid states through measurements of the local dynamical spin structure factor on the boundary using scanning tunneling microscopy (STM). We specifically focus on unpaired (“dangling”) Majorana fermions, which naturally emerge along boundaries of Kitaev spin liquids, and can serve as indicators of chiral boundary modes under broad conditions, thereby offering a clear signature of these exotic quantum states.

Zhang, Shang-Shun [University of Tennessee, Knoxvi↗