Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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

Striations in moderate pressure dc driven nitrogen glow discharge

Plasma stratification has been studied for more than a century. Despite the many experimental studies reported on this topic, theoretical analyses and numerical modeling of this phenomenon have been mostly limited to rare gases. In this work, a one-dimensional fluid model with detailed kinetics of electrons and vibrationally excited molecules is employed to simulate moderate-pressure (i.e. a few Torrs) dc discharge in nitrogen in a 15.5 cm long tube of radius 0.55 cm. The model also considers ambipolar diffusion to account for the radial loss of ions and electrons to the wall. The proposed model predicts self-excited standing striations in nitrogen for a range of discharge currents. The impact of electron transport parameters and reaction rates obtained from a solution of local two-term and a multi-term Boltzmann equation on the predictions are assessed. In-depth kinetic analysis indicates that the striations result from the undulations in electron temperature caused due to the interaction between ionization and vibrational reactions. Furthermore, the vibrationally excited molecules associated with the lower energy levels are found to influence nitrogen plasma stratification and the striation pattern strongly. A balance between ionization processes and electron energy transport allows the formation of the observed standing striations. Simulations were conducted for a range of discharge current densities from ~0.018 to 0.080 mA cm –2 , for an operating pressure of 0.7 Torr. Parametric studies show that the striation length decreases with increasing discharge current. Furthermore, the predictions from the model are compared against experimental measurements and are found to agree favorably.

42 ENGINEERING↗
DOE OSTI2024

Simulation of metal nanoparticles growth in methane atmosphere of arc discharge: comparison to experiment

A direct current arc discharge in a methane atmosphere is a scalable and sustainable method to produce metal-carbon core–shell nanoparticles and single-walled carbon nanotubes, where a metal catalyst can be continuously supplied through evaporation of an anode made from the catalyst material. The size of catalyst particles is of critical importance as it can affect the synthesis yield and properties of nanotubes and core–shell nanoparticles. This study presents a numerical model describing the formation and growth of metal particles for the conditions representative of the arc discharge with an evaporating iron anode at near-atmospheric pressure of a methane-rich atmosphere. The model incorporates carbon adsorption to the metal surface and explains the limiting effect of carbon coverage on the size of metal nanoparticles. The predicted particle sizes are compared with experimental observations. The model also predicts higher concentrations of metal particles with the increasing partial pressure of methane.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Development of a predictive capability of short-pulse laser-driven broadband x-ray radiography

High intensity, short-pulse laser interaction with a solid metal target produces broadband hard x-rays potentially for various applications of x-ray radiography. In this work, experimental benchmarking of numerical modelling for short-pulse laser-driven broadband x-ray radiography is presented. Angular dependent x-ray spectra are first calculated with a hybrid particle-in-cell code, Large Scale Plasma (LSP), using fast electron parameters inferred from an analysis of measured bremsstrahlung signals. Subsequently, a calculated x-ray spectrum in the direction of radiography is used in photon transport calculations using a Monte Carlo code, Particle and Heavy Ion Transport code System (PHITS), to simulate a radiographic image including a modelled 3D test object, an x-ray attenuation filter and an image plate detector. Simulated radiographic images are compared with measurements obtained in an experiment using a 50-TW Leopard short-pulse laser at the University of Nevada Reno. Results show that simulations reproduce the experimental images well for three different attenuation filters (plastic, aluminium, and brass), while 1D transmission profiles for the plastic and aluminium filters are quantitatively in good agreement. The modelling approach established in this work could be used as a predictive tool to simulate radiographic images of complex 3D solid objects at any arbitrary angular position or to optimize experimental components such as the source spectrum, x-ray attenuation filters and a detector type depending on a radiographic object without carrying out radiographic experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Toroidal plasma acceleration due to NBI fast ion losses in LTX- β

The recent Lithium Tokamak Experiment-Beta (LTX- β ) upgrade includes the addition of neutral beam injection (NBI) in the same direction as the plasma current (co- I P ) and a new toroidal Mirnov array for MHD characterization. In initial NBI experiments, a spontaneously rotating n = 1 MHD mode is seen to accelerate during NBI in the counter-beam direction, accompanied by a rise in electron density consistent with the beam-injected inventory but without a clear increase in plasma pressure. Together with analytic and numerical modeling of beam optics and fast ion confinement, these observations indicate the prompt loss of all or nearly all beam ions. However, the same modeling also suggests that planned upgrades to the Ohmic heating system should provide the fast ion confinement necessary for beam heating and core fueling. Here, a simple analytic model relates the momentum confinement time τ Φ to the observed evolution of mode rotation due to the combination of NBI momentum coupling, fast ion loss $\vec{J}\times\vec{B}$, and anomalous viscous torques, yielding τ Φ values consistent with past measurements of electron energy confinement time τ E,e .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Investigation of the role of hydrogen molecules in 1D simulation of divertor detachment

The role of neutral and charged hydrogenic molecules in detached regimes of tokamak plasmas is investigated using a simplified 1D numerical model. Using MAST Upgrade like conditions, simulations are implemented to study the rollover of target flux Γ in upstream density scan. It is found that if H 2 and H 2 + are considered in simulations a lower target temperature and a larger upstream density will be required to trigger divertor detachment under the same input power and particle flux, and the critical detachment threshold (the critical ratio of upstream static pressure to the power entering the recycling region) is found to be $^{ρ_{up}}_{ρ_{recl}}$ ~ 8.1 NMW -1 at rollover. Molecule–plasma interactions are found to be as crucial as atom–plasma interactions during divertor detachment, both of which account for the majority of plasma momentum loss in the cases studied here. Further analysis of the momentum loss decomposition shows molecule-plasma elastic collisions dominate molecule-plasma interactions, while molecular charge exchange cannot effectively reduce plasma momentum. In terms of H alpha emission, a strong rise of H alpha signal is found to be due to molecular excitation channels when the upstream density further increases after rollover.

SD1D↗
DOE OSTI2022

Impact of CH 4 addition on the electron properties and electric field dynamics in a Ar nanosecond-pulsed dielectric barrier discharge

Non-equilibrium plasmas derive their low temperature reactivity from producing and driving energetic electrons and active species under large electric fields. Therefore, the impact of reactants on the plasma properties including electron number density, electric field, and electron temperature is critical for applications such as plasma CH 4 reforming. Due to experimental complexity, electron properties and the electric field are rarely measured together in the same discharge. Here, we combine time-resolved Thomson scattering and electric field induced second harmonic generation (EFISH) to probe electron temperature, electron density, and electric field strength in a 60 Torr CH 4 /Ar nanosecond-pulsed dielectric barrier discharge (ns-DBD) while varying the CH 4 mole fraction from 0 to 8%. These measurements are compared to a 1-D numerical model to benchmark its predictions and identify areas of uncertainty. Nonlinear coupling between CH 4 addition, electron temperature, electron density, and the electric field was directly observed. Contrary to previous measurements in He, the electron temperature increased with CH 4 mole fraction. This rise in electron temperature is identified as electron heating by residual electric fields that increased with larger CH 4 mole fraction. Moreover, the electron number density has been found to decrease rapidly with the increase of methane mole fraction. Comparison of these measurements with the model yielded better agreement at higher CH 4 mole fractions and with the usage of ab initio calculated Ar electron-impact cross-sections from the B-spline R-matrix (BSR) database. Furthermore, the calculated plasma properties are shown to be sensitive to the residual surface charge implanted on the quartz dielectric surfaces. Without considering surface charge in the simulations, the calculated electric field profiles agreed well with the measurements, but the electron properties were underpredicted by more than a factor of three. Therefore, measurements of either the electric field or electron properties measurements alone are insufficient to fully validate modelling predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Measurement and reduction of Ar metastable densities by nitrogen admixing in electron beam-generated plasmas

Abstract Electron beam (e-beam) generated plasmas are useful for material processing applications such as deposition and etching because the plasmas deliver a large fluence of very low energy of ions to surfaces. Metastable species produced in the beam-region can also transport significant energy to the plasma periphery and surfaces. In this work, we have investigated the spatially resolved density of metastable Ar 1s 5 species produced in an Ar and Ar/N 2 e-beam generated plasma at pressures of 60–67 mTorr using laser-induced fluorescence (LIF). The experiments provide the first direct measure of absolute density and reduction of Ar 1s 5 in an e-beam generated plasma when argon is diluted with nitrogen. These results are consistent with previous predictions of numerical modeling and measurements using optical emission spectroscopy. The present spatially resolved LIF measurements directly quantify the reduction of Ar 1s 5 in the e-beam generated plasma by nitrogen admixing. This reduction was observed in the region of the electron beam and in the plasma periphery, where substrates are usually placed. For example, up to a threefold reduction of the density of Ar 1s 5 was measured when the argon background was diluted with 15.5% nitrogen at pressure of 60 mTorr. Ar 1s 5 reduction is attributed to excitation exchange with nitrogen molecules as well as the cooling of plasma electrons via inelastic collisions with nitrogen molecules.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Single-atom heat engine as a sensitive thermal probe

We propose employing a quantum heat engine as a sensitive probe for thermal baths. In particular, we study a single-atom Otto engine operating in an open thermodynamic cycle. Owing to its cyclic nature, the engine is capable of translating small temperature differences between two baths into a macroscopic oscillation in a flywheel. We present analytical and numerical modeling of the quantum dynamics of the engine and estimate it to be capable of detecting temperature differences as small as 2 μK. This sensitivity can be further improved by utilizing quantum resources such as squeezing of the ion motion. The proposed scheme does not require quantum state initialization and is able to detect small temperature differences in a wide range of base temperatures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Symmetrically tuned large-volume conic shell-cavities for axion searches

In an earlier paper [1], a new class of thin-shell cavities were proposed to evade the steep frequency scaling of conventional axion haloscopes. In this follow-up work, we see that a generalized conic geometry enables robust frequency-tuning for these large-volume cm-wave cavities. The frequency-defining dimension of a conic shell-cavity changes symmetrically and uniformly during tuning, maintaining a high axion coupling efficiency (the "form factor") to an external solenoid field. It is further shown that such tunable geometry is not restricted to circular cones. A general prescription for arbitrary volume-filling conic shell-cavities is developed and direct solutions are obtained for the created numerical models. Additionally, the largest of the realized designs is a meandering "brain" cavity that is tunable over a frequency range of 20%. The scan rate of this cavity is three orders of magnitude larger than that of a scaled cylindrical cavity used in the current generation experiments. The prospects for such a large improvement in the scan rate should motivate R& D efforts in fabrication and other implementation techniques. If these engineering challenges can be met, cavity-based axion haloscopes can stay competitive at frequencies higher than a few GHz. We propose an experimental configuration at 20 GHz (~80 μeV) using an array of brain cavities and compare it with other proposals for similar frequencies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Application of simulations and theory of sub-cyclotron frequency modes to DIII-D

Theoretical and numerical study of Alfvén eigenmodes (AEs) in the sub-cyclotron frequency range explains the observed frequency scaling with beam parameters in experiments across different devices. Global Alfvén (GAEs) and compressional Alfvén eigenmodes (CAEs) are frequently excited during neutral beam injection (NBI) in the National Spherical Torus Experiment (NSTX/NSTX-U), as well as other beam-heated devices such as MAST and DIII-D. These modes are driven unstable through the Doppler shifted cyclotron resonance with the NBI ions and can be excited in ITER due to super-Alfvénic velocities and anisotropy of the fast ion distribution. Numerical models and theory for sub-cyclotron frequency modes, previously developed and used to study the excitation of GAEs/CAEs in the NSTX(-U), have been successfully applied to explain the DIII-D observations. Simulations for DIII-D demonstrate that the modes with ω/ω ci ≈ 0.6, previously mis-identified as compressional Alfvén eigenmodes, have shear polarization and should be identified as the GAEs. Simulation results match the observed frequencies and estimated toroidal mode numbers in DIII-D experiments. Good agreement is found between scaling of the experimentally observed GAE frequencies with NBI parameters for NSTX, NSTX-U and DIII-D and scaling predicted by the theory and simulations. Here, the DIII-D simulation results also demonstrate that GAEs can be excited for larger aspect-ratio devices and smaller (sub- Alfvénic) beam injection velocities, therefore they can be unstable in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Resonance in radio frequency sheath admittance and enhanced impurity emission near the ion cyclotron frequency

Ion cyclotron resonance heating (ICRH) is of considerable interest among all auxiliary heating techniques, because it transfers power directly to ions, targets the high-density core, and involves the cheapest radio-frequency (RF) components. During ICRH operation, RF sheaths form on the ICRH antenna itself, nearby hardware, and far-field surfaces. These sheaths are associated with large hot-spot formation and impurity emissions. This work presents high-resolution numerical modelling of RF sheaths in nuclear fusion scenarios using hPIC2, a Debye-scale particle-in-cell code. The modelling reveals a new RF sheath phenomenon which occurs when the RF is resonant with a species ion cyclotron frequency or harmonic in an oblique magnetic field. This resonance of the RF sheath causes modifications of the energy-angle distributions of the ions impacting on the walls, with consequent increase in wall impurity emission. Due to the 1/ R scaling of the magnetic field in a tokamak, such cases are possible on divertor or vacuum chamber surfaces, depending on the geometry of the tokamak and the RF heating scenario. A simple physics interpretation using a driven damped harmonic oscillator is proposed, where the ion sheath transit time plays the role of damping the RF sheath admittance. Critically, the resonance leads to increased ion heat flux at the wall as well as increased physical sputtering, despite a lack of increase in RF rectified sheath potential.

ICRH antenna↗
DOE OSTI2022

Results from the FOCAL Experiment Campaign 1: Turbine Control Co-Design

This paper summarizes the scaling approach and results from the first experimental campaign of the Floating Offshore-wind and Controls Advanced Laboratory (FOCAL) Experimental Program, which is focused on developing public data sets to support the validation of numerical models for floating wind control co-design. A 1:70 Froude-scale performance-matched model of the IEA Wind 15-MW reference turbine was designed, built, and characterized on a fixed base. The National Renewable Energy Laboratory's Reference Open-Source COntroller (ROSCO) was integrated into the scale testing to examine the impact advanced turbine and generator control strategies have on turbine performance and loads, enabling validation of control strategies for floating wind design optimization.

17 WIND ENERGY↗
DOE OSTI2024

Wind Turbine Aeroelastic Stability in OpenFAST

Wind turbines are growing in size and increasingly suffer from aeroelastic instabilities. Unfortunately, numerical models often show inconsistent results during verification studies. We address this gap by first introducing novel linearization capabilities within the open-source aero-hydro-servo-elastic framework OpenFAST. Next, a code-to-code benchmark study is presented that compares modal parameters between OpenFAST and HAWCStab2 for a land-based version of the International Energy Agency 15-MW reference wind turbine modeled with quasi-steady aerodynamics. The two solvers are in strong agreement except for discrepancies in the second rotor flapwise modes. The differences are attributed to the torsional flexibility of the tower, which is assumed torsionally stiff in the OpenFAST model. Work is ongoing to close this modeling gap. The aeroelastic stability of a low-specific-power land-based wind turbine is also investigated. The impact of design choices is discussed, high-lighting how narrow the margins are between a stable design and an unstable design.

17 WIND ENERGY↗
DOE OSTI2021

Increasing freshwater and dissolved organic carbon flows to Northwest Alaska’s Elson lagoon

Manifestations of climate change in the Arctic are numerous and include hydrological cycle intensification and permafrost thaw, both expected as a result of atmospheric and surface warming. Across the terrestrial Arctic dissolved organic carbon (DOC) entrained in arctic rivers may be providing a carbon subsidy to coastal food webs. Yet, data from field sampling is too often of limited duration to confidently ascertain impacts of climate change on freshwater and DOC flows to coastal waters. This study applies numerical modeling to investigate trends in freshwater and DOC exports from land to Elson Lagoon in Northwest Alaska over the period 1981-2020. While the modeling approach has limitations, the results point to significant increases in freshwater and DOC exports to the lagoon over the past four decades. The model simulation reveals significant increases in surface, subsurface (suprapermafrost), and total freshwater exports. Significant increases are also noted in surface and subsurface DOC production and export, influenced by warming soils and associated active-layer thickening. The largest changes in subsurface components are noted in September, which has experienced a ~50% increase in DOC export emanating from suprapermafrost flow. Direct coastal suprapermafrost freshwater and DOC exports in late summer more than doubled between the first and last five years of the simulation period, with a large anomaly in September 2019 representing a more than fourfold increase over September direct coastal export during the early 1980s. These trends highlight the need for dedicated measurement programs that will enable improved understanding of climate change impacts on coastal zone processes in this data sparse region of Northwest Alaska.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A novel construct for scaling groundwater–river interactions based on machine-guided hydromorphic classification

Hydrologic exchange between river channels and adjacent subsurface environments is a key process that influences water quality and ecosystem function in river corridors. Predictive numerical models are needed to understand responses of river corridors to environmental change and to support sustainable watershed management. We posit that systematic hydromorphic classification provides a scaling construct that facilitates extrapolation of outputs from local-scale mechanistic models to reduced-order models applicable at reach and watershed scales. This in turn offers the potential to improve large-scale predictions of river corridor hydrobiogeochemical processes. Here we present a new machine-guided hydromorphic classification methodology that addresses the key requirements of this objective, and we demonstrate its application to a segment of the Columbia River in the northwestern United States. The resulting hydromorphic classes form spatially coherent and physically interpretable hydromorphic units that exhibit distinct behaviors in terms of distributions of subsurface residence times (a primary control on critical biogeochemical reactions). This approach forms the basis of ongoing research that is evaluating the formulation of reduced-order models and transferability of results to other river reaches and larger scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Analysis of the thermal recovery of HTS cables after an overcurrent fault

High temperature superconducting (HTS) power transmission cables are cooled to operating temperatures typically below 80 K. A valuable feature of some HTS cables is that they can limit short-circuit fault currents, which can be more than ten times the maximum operating current of a cable. When a fault occurs the HTS wire is no longer superconducting and the heating which occurs will raise the temperature above the allowable operating point. To recover from the fault, the conductor must be cooled down before returning to service. The recovery is affected by the refrigeration system's performance which varies with load and temperature. A transient numerical model of an HTS cable is applied to investigate the thermal recovery after a fault for some HTS cable system configurations and operating conditions.

Demko, Jonathan A.↗
DOE OSTI2024

Bilayer Ion Trap Design for 2D Arrays

Junctions are fundamental elements that support qubit locomotion in two-dimensional ion trap arrays and enhance connectivity in emerging trapped-ion quantum computers. In surface ion traps they have typically been implemented by shaping radio frequency (RF) electrodes in a single plane to minimize the disturbance to the pseudopotential. However, this method introduces issues related to RF lead routing that can increase power dissipation and the likelihood of voltage breakdown. Here, we propose and simulate a novel two-layer junction design incorporating two perpendicularly rotoreflected (rotated, then reflected) linear ion traps. The traps are vertically separated, and create a trapping potential between their respective planes. The orthogonal orientation of the RF electrodes of each trap relative to the other provides perpendicular axes of confinement that can be used to realize transport in two dimensions. While this design introduces manufacturing and operating challenges, as now two separate structures have to be precisely positioned relative to each other in the vertical direction and optical access from the top is obscured, it obviates the need to route RF leads below the top surface of the trap and eliminates the pseudopotential bumps that occur in typical junctions. Here in this paper the stability of idealized ion transfer in the new configuration is demonstrated, both by solving the Mathieu equation analytically to identify the stable regions and by numerically modeling ion dynamics. Our novel junction layout has the potential to enhance the flexibility of microfabricated ion trap control to enable large-scale trapped-ion quantum computing.

42 ENGINEERING↗
DOE OSTI2020

Estimating perturbed stress from 3-D borehole displacements induced by fluid injection in fractured or faulted shales

SUMMARY Hydrofracturing stress measurements in fractured and anisotropic shales are notoriously difficult, because opening of existing geological features tends to prevent the creation of a pure hydraulic fracture perpendicular to the least compressive principal stress. Here we show how adding 3-D borehole-displacement measurements while conducting the hydraulic injection test helps to better constrain the principal stress orientations and magnitudes. We developed a 3-D fully coupled hydromechanical numerical model to analyse the displacement, fluid pressure and injection flow-rate data measured during an injection pressure-step-rate test conducted to activate a faulted borehole interval in the Mont Terri Opalinus Clay (Switzerland). We find that injected fluids can only penetrate the fault when it is at or above the Coulomb failure pressure. Borehole displacement orientations are sensitive to a ∼15° variation in the stress–tensor orientation and a 1 MPa stress magnitude variation. Although some dispersion occurs while rupture is propagating along the fault plane ∼4 m away from the borehole, the maximum density of displacement orientations consistently informs about the stress orientation. Thus, an extended injection step-rate approach coupled with an accurate in situ measurement of the borehole wall displacements can be used to better constrain the local stress field perturbations in fractured shales and in heterogeneous rock in general.

Guglielmi, Yves↗