Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Real-Time Optical Monitoring and Simulations of Gas Phase Kinetics in InN Vapor Phase Epitaxy at High Pressure

Understanding the kinetics of nucleation and coalescence of heteroepitaxial thin films is a crucial step in controlling a chemical vapor deposition process, since it defines the perfection of the heteroepitaxial film both in terms of extended defect formation and chemical integrity of the interface. The initial nucleation process also defines the film quality during the later stages of film growth. The growth of emerging new materials heterostructures such as InN or In-rich Ga(x)In(1-x)N require deposition methods operating at higher vapor densities due to the high thermal decomposition pressure in these materials. High nitrogen pressure has been demonstrated to suppress thermal decomposition of InN, but has not been applied yet in chemical vapor deposition or etching experiments. Because of the difficulty with maintaining stochiometry at elevated temperature, current knowledge regarding thermodynamic data for InN, e.g., its melting point, temperature-dependent heat capacity, heat and entropy of formation are known with far less accuracy than for InP, InAs and InSb. Also, no information exists regarding the partial pressures of nitrogen and phosphorus along the liquidus surfaces of mixed-anion alloys of InN, of which the InN(x)P(1-x) system is the most interesting option. A miscibility gap is expected for InN(x)P(1-x) pseudobinary solidus compositions, but its extent is not established at this point by experimental studies under near equilibrium conditions. The extension of chemical vapor deposition to elevated pressure is also necessary for retaining stoichiometric single phase surface composition for materials that are characterized by large thermal decomposition pressures at optimum processing temperatures.

Dietz, Nikolaus↗
NASA NTRS1994

Generation of radiation in solar corona and interplanetary space by energetic electrons

Emissions of electromagnetic waves with frequencies close to the plasma frequency and/or its second harmonic have been frequently observed in the solar corona and interplanetary space. In the past, a number of theories have been put forward to esplain the generation mechanism of the observed radiation. In this paper, a new model is proposed. The essential point of the present theory is that the Langmuir waves amplified as a result of the usual beam instability can lead to two important effects: first, electrostatic waves with frequencies close to twice the plasma frequency can be excited; and second, a significant modification of the dispersion relation can occur, so that these electrostatic waves can naturally change into electomagnetic waves as they propagate in a plasma in which the plasma density decreases spatially. The latter effect is attributed to a mode couplng process. In addition to the second harmonic emission, emission at the fundamental is also briefly discussed. In this case, as in many other theories, the presence of a very low frequency electostatic wave such as the ion-acoustic wave is assumed. The emission process discussed in the present theory stresses the importance of mode coupling and conversion rather than kinetic processes such as a nonlinear wave-wave or wave-particle scattering.

Wu, C. S.↗
DOE OSTI2023

Twin Boundaries Contribute to The First Cycle Irreversibility of LiNiO 2

Abstract LiNiO 2 remains a target for layered oxide Li‐ion cathode development as it can theoretically deliver the highest energy density of this materials class. In practice, LiNiO 2 suffers from poor capacity retention due to electrochemically‐induced structural changes. While the impact of Ni off‐stoichiometry on the electrochemical performance has been extensively studied, that of planar defects present in the as‐synthesized cathode is not well understood. Advanced ex situ and operando structure probes are used to identify and quantify point and planar defects present in as‐synthesized Li 1‐y Ni 1+y O 2 cathodes and monitor their evolution during the first cycle. Specifically, a 7 Li nuclear magnetic resonance (NMR) signature characteristic of Li environments near twin boundaries is identified; an assignment supported by first‐principles calculations and scanning transmission electron microscopy (STEM) images of twin boundary defects. The NMR results suggest that the concentration of twin boundaries depends on the amount of Ni excess. Moreover, operando magnetometry and ex situ synchrotron X‐ray diffraction and NMR demonstrate that these planar defects impede Li reinsertion into the bulk cathode at reasonable discharge rates and contribute to the first cycle irreversible capacity. These findings provide new design rules for Li 1‐y Ni 1+y O 2 cathodes, whereby a reduced concentration of twin boundaries in the pristine material leads to reduced kinetic limitations and improved cathode utilization.

Nguyen, H.↗
NASA NTRS2006

Raised Relief Mars Globe Brings the Red Planet Closer

The Mars Globe 20x is the first digitally produced raised relief globe to be sold at retail establishments. The rises and dips of the Martian landscape have been multiplied by 20 to make the shapes more noticeable to human touch and sight. They make this globe, with its over 1 million elevation points, a visually stunning introduction to the planet. Spectrum 3D used the NASA digital land elevation data to aim lasers that then shaped and defined the master globe s surface. Subsequent copies were then made by creating a master globe mold. The molded copies are hand finished by workers who remove errant edges or lines that may appear on the raw globes and then paint the surfaces. The result is a globe that measures in at 18 inches in diameter, roughly 1:15,729,473 scale of the actual planet. The exaggerated raised relief is like having a 3-D digital microscope for planetary shapes. This makes the landmarks easier to learn and understand, as it provides easy visuals for orientation. People have a natural propensity for understanding 3-D shapes more easily than numbers or words. The 3-D globes appeal to both the kinetic and visual learning aspects of the brain, making it easy for people to readily memorize the landmarks and to make a mental model that they will remember for a long time.

Source record↗
DOE OSTI2024

Exploring ship track spreading rates with a physics-informed Langevin particle parameterization

Abstract. The rate at which aerosols spread from a point source injection, such as from a ship or other stationary pollution source, is critical for accurately representing subgrid plume spreading in a climate model. Such climate model results will guide future decisions regarding the feasibility and application of large-scale intentional marine cloud brightening (MCB). Prior modeling studies have shown that the rate at which ship plumes spread may be strongly dependent on meteorological conditions, such as precipitating versus non-precipitating boundary layers and shear. In this study, we apply a Lagrangian particle model (PM-ABL v1.0), governed by a Langevin stochastic differential equation, to create a simplified framework for predicting the rate of spreading from a ship-injected aerosol plume in sheared, precipitating, and non-precipitating boundary layers. The velocity and position of each stochastic particle is predicted with the acceleration of each particle being driven by the turbulent kinetic energy, dissipation rate, momentum variance, and mean wind. These inputs to the stochastic particle velocity equation are derived from high-fidelity large-eddy simulations (LES) equipped with a prognostic aerosol–cloud microphysics scheme (UW-SAM) to simulate an aerosol injection from a ship into a cloud-topped marine boundary layer. The resulting spreading rate from the reduced-order stochastic model is then compared to the spreading rate in the LES. The stochastic particle velocity representation is shown to reasonably reproduce spreading rates in sheared, precipitating, and non-precipitating cases using domain-averaged turbulent statistics from the LES.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗
DOE OSTI2024

A Gaussian process autoregressive model capturing microstructure evolution paths in a Ni–Mo–Nb alloy

Additive manufacturing is increasingly being employed to produce components of complex geometries in structural alloys because of the expected energy savings associated with the near-net-shape capability and the ability to build in novel internal features that are not possible with many conventional manufacturing approaches. However, because of the extreme thermal conditions encountered, the non-equilibrium microstructures produced during powder bed-based additive manufacturing processes must be subjected to custom post-heat treatment processes to recover the target mechanical properties. Phase-field models and simulation techniques have matured to a state where the microstructure evolution paths, and the morphologies of the resulting precipitate phases can be predicted reasonably accurately, considering alloy-specific thermodynamic and kinetic aspects of the nucleation and growth processes. However, phase-field simulations are computationally intensive, which precludes the ability to apply the simulations directly to the length scale of the entire component. Therefore, it is highly desirable to develop low-computational-cost surrogate models that effectively capture the physics at the microstructural length scale, while facilitating the design of optimized processing conditions resulting in location-specific targeted microstructures at the component scale. The work presented here demonstrates the application of the materials knowledge system framework to develop a surrogate model that effectively captures the microstructural path during annealing of a Ni–Mo–Nb alloy containing different Mo and Nb compositions known to segregate during solidification under additive manufacturing conditions. Specifically, the surrogate model built in this work is based on a Gaussian process autoregressive model informed by statistical representation of simulated microstructures using two-point correlations and dimensionality reduction through principal component analysis. In conclusion, this surrogate model is shown to capture the bifurcation of the microstructural path during precipitation, which yields a microstructure dominated by the $\gamma^{\prime\prime}$ phase at high Nb concentrations and the $\delta$ phase at low Nb concentrations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Emerging investigator series: kinetics of diopside reactivity for carbon mineralization in mafic–ultramafic rocks

The ongoing use of fossil fuels to supply modern energy demands has necessitated research on combating carbon dioxide (CO 2 ) emissions and climate change. Carbon storage via mineral trapping in basalt and related rocks is a promising strategy. However, mineralization rates depend on the variable minerology that makes up these rock formations. Diopside (CaMgSi 2 O 6 ) is a common pyroxene mineral in ultramafic and mafic rocks including basalt, but relatively little work has been done to understand its carbon mineralization kinetics using hydrated supercritical CO 2 , which induces the formation of reactive nanoscale interfacial water films. Here, in situ XRD experiments at 50–110 °C and 90 bar indicate that diopside transforms into a myriad of Mg/Ca carbonates, including huntite [Mg 3 Ca(CO 3 ) 4 ] and very high magnesium calcite (VHMC, i.e., protodolomite). Through ex situ characterization, we were able to constrain reaction pathways for the dissolution–precipitation transformation process including metastable intermediate precipitates. Experiments performed at variable temperatures enabled Avrami-derived rate constants and an apparent activation energy of 97 ± 16 kJ mol –1 , implying the dissolution of diopside is the rate-limiting step. Density functional theory (DFT) calculations, used to gain molecular insight into the surface stability of the diopside during dissolution, suggest that exposed calcium cations are susceptible to dissolution when put in contact with water given their coordination environment. The collective results point to the high CO 2 mineralization potential of diopside in basalts, which could help guide parameterization of reactive transport models needed to design and permit commercial-scale subsurface carbon storage operations.

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

Demonstration of low-mode shape control in indirect-drive double shell implosions at the NIF

Double shell inertial confinement fusion is a concept for achieving robust thermonuclear burn that uses dense metal shells to compress deuterium-tritium (DT) fuel to fusion conditions. Double shell implosions are typically indirectly driven and involve a target that consists of a low-Z ablator, a foam layer, and a high-Z pusher surrounding the DT fuel. The goal of the campaign is to achieve a volumetric burn as radiation losses from the DT fuel are trapped by the opaque high-Z shell. The overall performance of double shell implosions relies on the efficient collisional transfer of kinetic energy between layers. The efficiency of this transfer (and therefore the overall performance of a given implosion) is degraded by the presence of low-mode asymmetries. P2 asymmetries are often observed in spatially resolved 2D radiographs of nominal double shell implosions. This work discusses three such experiments: one with an oblate P2 asymmetry, one with a prolate P2 asymmetry, and one with an approximate spherical symmetry. After performing a shape analysis of the oblate and prolate implosions to quantify asymmetries, these experimental results were compared with the results of hydrodynamic simulations for the two experiments. Differences between the experiment and simulation were then used to design an approximately spherical implosion by altering the incident laser cone fraction. Radiographs from the experiment that implemented the modified cone fraction show evidence of an implosion that is approximately spherical until bang time. This design is intended to serve as a point design for future studies that will seek to optimize various aspects of the double shell target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Dipole potentials indicate restructuring of the membrane interface induced by gadolinium and beryllium ions

The dipole component of the membrane boundary potential, phi(d), is an integral parameter that may report on the conformational state of the lipid headgroups and their hydration. In this work, we describe an experimental approach to measurements of the dipole potential changes, Deltaphi(d), and apply it in studies of Be(2+) and Gd(3+) interactions with membranes composed of phosphatidylserine (PS), phosphatidylcholine (PC), and their mixtures. Deltaphi(d) is determined as the difference between the changes of the total boundary potential, phi(b), measured by the IFC method in planar lipid membranes and the surface potential, phi(s), determined from the electrophoretic mobility of liposomes. The Gouy-Chapman-Stern formalism, combined with the condition of mass balance, well describes the ion equilibria for these high-affinity cations. For the adsorption of Be(2+) and Gd(3+) to PC membranes, and of Mg(2+) to PS membranes, the values of Deltaphi(b) and Deltaphi(s) are the same, indicative of no change of phi(d). Binding of Gd(3+) to PS-containing membranes induces changes of phi(d) of opposite signs depending on the density of ionized PS headgroups in the bilayer. At low density, the induced Deltaphi(d) is negative (-30 mV), consistent with the effect of dehydration of the surface. At maximal density (pure PS, neutral pH), adsorption of Be(2+) or Gd(3+) induces an increase of phi(d) of 35 or 140 mV, respectively. The onset of the strong positive dipole effect on PS membranes with Gd(3+) is observed near the zero charge point and correlates with a six-fold increase of membrane tension. The observed phenomena may reflect concerted reorientation of dipole moments of PS headgroups as a result of ion adsorption and lipid condensation. Their possible implications to in-vivo effects of these high-affinity ions are discussed.

Non-NASA Center↗
NASA NTRS2002

Solutal Convection Around Growing Protein Crystal and Diffusional Purification in Space

This work theoretically addressed two subjects: 1) onset of convection, 2) distribution of impurities. Onset of convection was considered analytically and numerically. Crystal growth was characterized by slow surface incorporation kinetics, i.e. growth kinetic coefficient beta (cm/s) small as compared to the typical bulk diffusion rate, D(sub 1)/h, where D(sub 1) is diffusivity of major crystallizing protein and h is the crystal size. Scaling type analysis predicted two laws on how the convection rate, v, essentially the Peclet number, Pe exactly equal to vh/D(sub 1), depends on dimensionless kinetic coefficient a exactly equal to beta h/D(sub 1). Namely: Pe = C(sub 2/5)(aRa(sup 2/5)) and Pe = C(sub 1) aRa. Here, Reynolds number Ra = rho(sub 1)(sup 0)gh(sup 3)(rho(sub p) - rho(sub w))/rho(sup p)rho(sub 1)vD(sub 1), v being solution viscosity. The constants C(sub 2/5), exactly equal to 0.28 and C(sub 1) exactly equal to 10(exp -2) found from the full scale computer simulation for a cylindrical crystal inside big cylindrical vessel. The linear boundary conditions connecting protein and impurity concentration at the interface with the flux to/from the interface was applied. No-slip condition for Navier-Shocker equations was employed. With these conditions, flow and concentration distributions were calculated. Validity of the Pe(Ra) dependencies follows for wide range of parameters for which numerical calculations have been accomplished and presented by various points.

Chernov, A. A.↗
NASA NTRS1995

A pseudo-sound constitutive relationship for the dilatational covariances in compressible turbulence: An analytical theory

The mathematical consequences of a few simple scaling assumptions about the effects of compressibility are explored using a simple singular perturbation idea and the methods of statistical fluid mechanics. Representations for the pressure-dilation and dilatational dissipation covariances appearing in single-point moment closures for compressible turbulence are obtained. While the results are expressed in the context of a second-order statistical closure they provide some interesting and very clear physical metaphors for the effects of compressibility that have not been seen using more traditional linear stability methods. In the limit of homogeneous turbulence with quasi-normal large-scales the expressions derived are - in the low turbulent Mach number limit - asymptotically exact. The expressions obtained are functions of the rate of change of the turbulence energy, its correlation length scale, and the relative time scale of the cascade rate. The expressions for the dilatational covariances contain constants which have a precise and definite physical significance; they are related to various integrals of the longitudinal velocity correlation. The pressure-dilation covariance is found to be a nonequilibrium phenomena related to the time rate of change of the internal energy and the kinetic energy of the turbulence. Also of interest is the fact that the representation for the dilatational dissipation in turbulence, with or without shear, features a dependence on the Reynolds number. This article is a documentation of an analytical investigation of the implications of a pseudo-sound theory for the effects of compressibility.

Ristorcelli, J. R.↗
NASA NTRSDate not supplied

EMU Ventilation Loop Simulation and Assessment of Contamination of the EMU Sublimator Hydrophilic Coating

Controlling moisture is a critical function performed in the space suit during extra vehicular activity (EVA) missions. Currently, this function is carried out by a sublimator, which is coated with a hydrophilic material. Tests have shown that siloxane compounds can off gas from the helmet absorption pad extension (HAP-E) and helmet absorption band (HAB) when they are at elevated temperatures and initial testing indicated that siloxanes could damage the hydrophilic coating. Therefore, more representative testing was needed. Important test considerations were source and destination kinetics, system geometries, scrubbing efficacy by the activated charcoal from the Metox canister, and other benefits that derive from active condensation occurring on the hydrophilic coating during exposure. A ventilation loop that was originally constructed to test hardware for control of CO2 and trace contaminants under EVA conditions was modified to conduct the siloxane off gassing tests. A vacuum compatible acrylic chamber was installed to simulate the suit volume. A steel beaker located in the chamber contained the HAP-E and HAB materials and was heated with two band heaters. The flow exiting the chamber passed through a LiOH cannister, which contained carbon beds to simulate their function in the Metox. The flow then was directed into a smaller volume that contained a coupon, which was coated with the hydrophilic material and was maintained at 10°C (50°F). Moisture was injected upstream of the suit volume simulator to maintain the dew point between 12.7°C (55°F) and 18.3°C (65°F). The pressure was controlled to between 3.8 and 4.4 psia, and the flow was held at 6 acfm. The system was run for an equivalent of 25 ten-hour EVAs, where the HAP-E/HAB materials were changed out after each EVA. The effect of the potential off gassing on the coupon was measured after each EVA period.

PLSS↗
DOE OSTI2023

Origin of Rapid Delithiation In Secondary Particles Of LiNi 0.8 Co 0.15 Al 0.05 O 2 and LiNi y Mn z Co 1− y − z O 2 Cathodes

Abstract Most research on the electrochemical dynamics in materials for high‐energy Li‐ion batteries has focused on the global behavior of the electrode. This approach is susceptible to misleading analyses resulting from idiosyncratic kinetic conditions, such as surface impurities inducing an apparent two‐phase transformation within LiNi 0.8 Co 0.15 Al 0.05 O 2 . Here, nano‐focused X‐ray probes are used to measure delithiation operando at the scale of secondary particle agglomerates in layered cathode materials during charge. After an initial latent phase, individual secondary particles undergo rapid, stochastic, and largely uniform delithiation, which is in contrast with the gradual increase in cell potential. This behavior reproduces across several layered oxides. Operando X‐ray microdiffraction (‐XRD) leverages the relationship between Li content and lattice parameter to further reveal that rate acceleration occurs between Li‐site fraction ( x Li ) ≈0.9 and ≈0.5 for LiNi 0.8 Co 0.15 Al 0.05 O 2 . Physics‐based modeling shows that, to reproduce the experimental results, the exchange current density ( i 0 ) must depend on x Li , and that i 0 should increase rapidly over three orders of magnitude at the transition point. The specifics and implications of this jump in i 0 are crucial to understanding the charge‐storage reaction of Li‐ion battery cathodes.

25 ENERGY STORAGE↗
DOE OSTI2020

A New General Paradigm for Understanding and Preventing Li Metal Penetration through Solid Electrolytes

The use of lithium (Li) or sodium (Na) metal anodes together with highly ion-conductive solid electrolytes (SEs) could provide batteries with a step improvement in volumetric and gravimetric energy densities. Unfortunately, these SEs face significant technical challenges, in large part because Li and Na dendrites can penetrate through SEs, leading to short circuits. The ability of such a soft material (Li or Na metal) to penetrate through ceramic is surprising from the point of view of models widely used in the Li-battery field. We introduce a concept, new to the battery field, for preventing penetration of lithium dendrites through SEs by putting the SE surfaces into a state of residual compressive stress. For a sufficiently high compressive stress, cracks have difficulty forming, and cracks that do form are forced to close, inhibiting dendrite penetration. This approach is widely used to solve commercially important stress corrosion cracking problems in metals and static fatigue problems in ceramics and glasses (e.g., Gorilla Glass). However, the technique will not be useful for SEs if the Li-ion transport rate through a SE is substantially reduced when the SE is under compression. Our molecular dynamics calculations for Li-ion transport through a common SE demonstrate that the introduction of even very high residual compressive stresses (~10 GPa) has only a modest effect on Li-ion transport kinetics, suggesting that this approach is viable and capable of providing a new paradigm for developing high-performance and mechanically stable SEs.

25 ENERGY STORAGE↗
NASA NTRS2016

Micro Acoustic Resonant Chambers for Heating/Agitating/Mixing (MARCHAM)

A variety of applications require the mixing and/or heating of a slurry made from a powder/fluid mixture. One of these applications, Sub Critical Water Extraction (SCWE), is a process where water and an environmental powder sample (sieved soil, drill cuttings, etc.) are heated in a sealed chamber to temperatures greater than 200 degrees Celsius by allowing the pressure to increase, but without reaching the critical point of water. At these temperatures, the ability of water to extract organics from solid particulate increases drastically. This paper describes the modeling and experimentation on the use of an acoustic resonant chamber which is part of an amino acid detection instrument called Astrobionibbler [Noell et al. 2014, 2015]. In this instrument we use acoustics to excite a fluid- solid fines mixture in different frequency/Amplitude regimes to accomplish a variety of sample processing tasks. Driving the acoustic resonant chamber at lower frequencies can create circulation patterns in the fluid and mixes the liquid and fines, while driving the chamber at higher frequencies one can agitate the fluid and powder and create a suspension. If one then drives the chamber at high amplitude at resonance heating of the slurry occurs. In the mixing and agitating cell the particle levitation force depends on the relative densities and compressibility’s of the particulate and fluid and on the kinetic and potential energy densities associated with the velocity and pressure fields [Glynne-Jones, Boltryk and Hill 2012] in the cell. When heating, the piezoelectric transducer and chamber is driven at high power in resonance where the solid/fines region is modelled as acoustic transmission line with a large loss component. In this regime, heat is pumped into the solution/fines mixture and rapidly heats the sample. We have modeled the piezoelectric transducer/chamber/ sample using Mason’s equivalent circuit. In order to assess the validity of the model we have built and tested a variety of chambers. This paper describes the experimental results which are in general agreement with theory within the limitations of the modeling.

fluid particle interaction↗
DOE OSTI2026

Impact of atmospheric turbulence on performance and loads of wind turbines: knowledge gaps and research challenges

Wind energy harvesting from the atmosphere takes place in the atmospheric boundary layer. The boundary layer shear and buoyancy create three-dimensional turbulent eddies spanning a range of scales that form a continuous forward cascade of kinetic energy to the smallest scales of motion where energy is dissipated. Large-scale atmospheric circulations modulate the boundary layer turbulence, characterized by coherence and intermittency. As wind turbines grow in size and the integrated control of both turbines and wind farms spans greater distances, the relationship between the scales of atmospheric turbulence and the design and operation of wind energy facilities has entered new territory. The boundary layer turbulence impacts both wind turbine power production and turbine loads. Optimizing wind turbine and wind farm performance requires an understanding of how turbulence affects both wind turbine efficiency and reliability. While the characteristics of atmospheric boundary layer turbulence have been observed and studied in detail over the last few decades, there are still significant gaps in our understanding of the impact of turbulence on wind power resources and wind farm operations. This paper outlines the current state of turbulence research relevant to wind energy applications and points to gaps in our knowledge that need to be addressed to effectively utilize wind resources.

Kosović, Branko [Johns Hopkins Univ., Baltimore, M↗
DOE OSTI2021

Energy barriers for Dy and H penetrating graphene on 6 H -SiC(0001) and freestanding bilayer graphene from first-principles calculations

Currently, intercalation of foreign guest atoms into two-dimensional (2D) layered van der Waals materials is an active research area motivated in part by the development of next-generation energy-storage technologies and optoelectronic devices. One such extensively studied 2D material is the graphene-on-SiC system. To realize and control the desired intercalated structures, it is fundamentally important to understand the kinetic process of intercalation. For the intercalation of a guest atom into graphene layers on SiC substrate, a critical kinetic parameter is the energy barrier of a guest atom penetrating the perfect graphene top layer into the gallery under it. However, accurate theoretical calculations for such penetration barriers are unavailable in literature. From our first-principles density functional theory calculations, we obtain the global energy barriers of 3.47 and 1.80 eV for single Dy and H atoms penetrating the graphene top layer on a graphene buffer layer supported by a Si-terminated 6H-SiC(0001) substrate, respectively. Furthermore, for comparison as well as for examining the lateral strain effects, we also obtain the global barriers of 5.05 and 1.50 eV for single Dy and H atoms penetrating freestanding bilayer graphene with a tensile strain of about 8.8% to match our model for supported graphene, as well as the global barriers of 7.21 and 4.18 eV for penetrating unstrained freestanding bilayer graphene, respectively. From corresponding minimum energy paths with multiple energy minima and saddle points, we can also obtain various local energy barriers and the global backward barrier from the graphene gallery back to the top surface.

2D materials↗