SEARCH · Search NASA

Results for “finite element Method”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66

All sources NASA NTRS DOE OSTI

DOE OSTI2010

Collocation Method for Numerical Solution of Coupled Nonlinear Schroedinger Equation

The coupled nonlinear Schroedinger equation models several interesting physical phenomena presents a model equation for optical fiber with linear birefringence. In this paper we use collocation method to solve this equation, we test this method for stability and accuracy. Numerical tests using single soliton and interaction of three solitons are used to test the resulting scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Interface nano-confined acoustic waves in polymeric surface phononic crystals

The impulsive acoustic dynamics of soft polymeric surface phononic crystals is investigated here in the hypersonic frequency range by near-IR time-resolved optical diffraction. The acoustic response is analysed by means of wavelet spectral methods and finite element modeling. An unprecedented class of acoustic modes propagating within the polymer surface phononic crystal and confined within 100 nm of the nano-patterned interface is revealed. The present finding opens the path to an alternative paradigm for characterizing the mechanical properties of soft polymers at interfaces and for sensing schemes exploiting polymers as embedding materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Thermal management and prototype testing of Compton scattering X-ray beam position monitor for the Advanced Photon Source Upgrade

Accurate and stable x-ray beam position monitors (XBPMs) are key elements in obtaining the desired user beam stability in the Advanced Photon Source Upgrade. In the next-generation XBPMs for the canted-undulator front ends, where two undulator beams are separated by 1.0 mrad, the lower beam power (<10 kW) per undulator allows us to explore lower-cost solutions based on Compton scattering from a diamond placed edge-on to the x-ray beam. Because of the high peak power density of the x-ray beams, this diamond experiences high temperatures and has to be clamped to a water-cooled heat spreader using thermal interface materials (TIMs), which play a key role in reducing the temperature of the diamond. To evaluate temperature changes through the interface via thermal simulations, the thermal contact resistance (TCR) of TIMs at an interface between two solid materials under even contact pressure must be known. This paper addresses the TCR measurements of several TIMs, including gold, silver, pyrolytic graphite sheet, and 3D graphene foam. In addition, a prototype of a Compton-scattering XBPM with diamond blades was installed at APS Beamline 24-ID-A in May 2015 and has been tested. This study presents the design of the Compton-scattering XBPM, and compares thermal simulation results obtained for the diamond blade of this XBPM by the finite element method with in situ empirical measurements obtained by using reliable infrared technology.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Gunn threshold voltage characterization in GaAs devices with wedge-shaped tapering

Here, we fabricate gallium arsenide-based devices with a wedge-shaped tapering region connected to a rectangular-shaped region and measure the threshold voltage required to trigger the Gunn effect. The threshold voltage reduction is attributed to the focusing of the electric field toward the narrower end of the device and is effective when the device has a steep enough tapering. We also model the electric field profile for the tapered devices using an intuitive graphical approach and the finite element method and provide estimates for the threshold voltages of tapered devices. Finally, we compare the estimates to the measured values and provide possible reasons for the discrepancies. We believe the capability of threshold voltage reduction with the wedge-shaped tapering design could be useful in device applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Droplet formation simulation using mixed finite elements

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Solid-state synthesis of UV-plasmonic Cr2N nanoparticles

Materials that exhibit plasmonic response in the UV region can be advantageous for many applications, such as biological photodegradation, photocatalysis, disinfection, and bioimaging. Transition metal nitrides have recently emerged as chemically and thermally stable alternatives to metal-based plasmonic materials. However, most free-standing nitride nanostructures explored so far have plasmonic responses in the visible and near-IR regions. Herein, we report the synthesis of UV-plasmonic Cr2N nanoparticles using a solid-state nitridation reaction. The nanoparticles had an average diameter of 9 ± 5 nm and a positively charged surface that yields stable colloidal suspension. The particles were composed of a crystalline nitride core and an amorphous oxide/oxynitride shell whose thickness varied between 1 and 7 nm. Calculations performed using the finite element method predicted the localized surface plasmon resonance (LSPR) for these nanoparticles to be in the UV-C region (100–280 nm). While a distinctive LSPR peak could not be observed using absorbance measurements, low-loss electron energy loss spectroscopy showed the presence of surface plasmons between 80 and 250 nm (or ∼5 to 15 eV) and bulk plasmons centered around 50–62 nm (or ∼20 to 25 eV). Plasmonic coupling was also observed between the nanoparticles, resulting in resonances between 250 and 400 nm (or ∼2.5 to 5 eV).

Chemistry↗

DOE OSTI2023

Capacitance–voltage modeling of mid-wavelength infrared nBn detectors

Capacitance–voltage measurements are a powerful technique to determine doping profiles of semiconductor pn junctions and Schottky barrier diodes. The measurements were recently extended to III-V-based mid-wavelength nBn infrared detectors, and absorber doping densities have been extracted using the widely accepted Schottky approximation, where the potential drop across the device is assumed to be across the depleting absorber layer. However, this approach is limited to when the absorber region of the nBn is under high reverse bias and thus is only able to extract the absorber region doping profile. Here, we introduce a semi-analytical model that is capable of extracting barrier dopant polarity, doping concentration, and thickness, as well as contact and absorber layer doping concentrations, all from a capacitance–voltage measurement. Rather than solely considering the potential drop across the depleting layers, it considers the potential drop across the accumulating layer as well. This negative charge accumulation occurs for the contact and absorber layers in the case of reverse and forward biases, respectively. This allows for a single model to be applied to a capacitance–voltage curve at forward and reverse biases and it can provide regions of bias where the absorber transitions from depletion to accumulation. We compare the agreement of the semianalytical model with modeling results from commercially available finite element method software and experimental capacitance–voltage data. Finally, we show that the method is consistent with the Schottky approximation of extracting absorber doping densities at high reverse bias and discuss the model's limitations.

Physics↗

DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗

DOE OSTI2015

DNS of turbulent flow with hemispherical wall roughness

The present study of the effect of roughness density on the mean flow turbulence parameters is motivated by the need for new generation of boundary conditions for multiphase computational multiphase fluid dynamics (CMFD) models applied to boiling flows. Effect of roughness element density on the turbulent flow in a channel is quantified through direct numerical simulations (DNSs). The Navier--Stokes equations are solved using finite element method and bubbles are approximated as rigid near-hemispherical obstacles at the wall. Six different cases were analysed including channel flow with smooth wall and channel flow with rough wall for five different bubble nucleation site densities. Friction factor and the law of the wall was calculated and compared with the previously published results. Existing correlations for nucleating bubble site density dependency on a wall heat flux were used to obtain a relation between the heat flux and the friction factor, leading to the law of the wall dependency on the heat flux. This separate effect study provides new guidelines on how the heat flux in subcooled boiling regime affects the turbulence behaviour near the wall and guides the computational fluid dynamics model development for boiling two-phase flows.

Mechanics↗

DOE OSTI2014

Inter-strand current sharing and ac loss measurements in superconducting YBCO Roebel cables

A Roebel cable, one twist pitch long, was modified from its as-received state by soldering copper strips between the strands to provide inter-strand connections enabling current sharing. Various DC transport currents (representing different percentages of its critical current) were applied to a single strand of such a modified cable at 77 K in a liquid nitrogen bath. Simultaneous monitoring of I–V curves in different parts of the strand as well as in its interconnections with other strands was made using a number of sensitive Keithley nanovoltmeters in combination with a multichannel high-speed data acquisition card, all controlled via LabView software. Current sharing onset was observed at about 1.02 of strand I c . At a strand current of 1.3I c about 5% of the current was shared through the copper strip interconnections. A finite element method modeling was performed to estimate the inter-strand resistivities required to enable different levels of current sharing. The relative contributions of coupling and hysteretic magnetization (and loss) were compared, and for our cable and tape geometry, and at dB/dt=1 T s -1 , and our inter-strand resistance of 0.77 mΩ, (enabling a current sharing of 5% at 1.3I c ) the coupling component was 0.32% of the hysteretic component. However, inter-strand contact resistance values of 100–1000 times smaller (close to those of NbTi and Nb 3 Sn based accelerator cables) would make the coupling components comparable in size to the hysteretic components.

36 MATERIALS SCIENCE↗

DOE OSTI2015

Inter-strand current sharing and ac loss measurements in superconducting YBCO Roebel cables

36 MATERIALS SCIENCE↗

DOE OSTI2018

Characterization of the arc in crossflow using a two-temperature nonequilibrium plasma flow model

Diverse industrial applications such as circuit breakers and wire arc spraying involve the interaction between an electric arc and a stream of gas impinging perpendicular to it, a configuration commonly referred to as the arc in crossflow. The arc in crossflow is simulated using a three-dimensional time-dependent two-temperature (heavy-species and electrons) plasma flow model to better capture plasma-gas interactions and deviations from Local Thermodynamic Equilibrium (LTE). The coupled fluid-electromagnetic flow model is solved in a monolithic manner using Variational Multiscale Finite Element Method. Simulation findings are validated with experimental findings and contrasted against results obtained with a LTE model. Results from the two-temperature model corroborate experimental observations while providing quantification of the deviation between heavy-species and electron temperatures. The model is used to characterize the arc in crossflow as a function of the Reynolds and Enthalpy dimensionless numbers, which encapsulate the inter-dependence among the main parameters total current, inflow velocity, and inter-electrode spacing. Furthermore, the characterization revealed the behavior of arc shape, voltage drop, arc power, the degree of nonequilibrium, as well as the characteristic plasma front thickness, with varying controlling parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Spin wave modes in a cylindrical nanowire in crossover dipolar-exchange regime

Nanoscale magnetic systems have been studied extensively in various geometries, such as wires of different cross-sections, arrays of wires, dots, rings, etc. Such systems have interesting physical properties and promising applications in advanced magnetic devices. Uniform magnetic nanowires are the basic structures which were broadly investigated. However, some of their dynamical properties, like: (anti)crossing between the spin wave modes and impact of the magnetic field on the spin wave spectrum, still need to be exploited. We continue this research by investigation of the spin wave dynamics in solid Ni nanowire of the circular cross-section. We use two approaches: semi-analytical calculations and numerical computations based on the finite element method. We solve coupled Landau-Lifshitz and Maxwell equations and consider both dipolar and exchange interactions. We identify the dispersion brunches and their (anti)crossing by plotting the spatial profiles of spin wave amplitudes and magnetostatic potential. We also check how we can tune the spectrum of the modes by application of the external magnetic field and how it affects the modes and their dominating type of interaction.

magnonics↗

DOE OSTI2022

Three-dimensional modelling of a self-sustained atmospheric pressure glow discharge

The atmospheric pressure glow discharge (APGD) is a relatively simple and versatile plasma source used in diverse applications. Stable APGD operation at high currents, generally a challenge due to instabilities leading to glow-to-arc transition, has been demonstrated using actively-controlled cathodic cooling. This article presents the computational modelling and simulation of a self-sustained direct-current APGD in helium within a 10 mm pin-to-plate inter-electrode gap for currents ranging from 4 to 40 mA. The APGD model is comprised of the conservation equations for total mass, chemical species, momentum, thermal energy of heavy-species and of free electrons, and electric charge. The model equations are discretized using a nonlinear variational multi-scale finite element method that has demonstrated superior accuracy in other plasma flow problems, on a temporal and three-dimensional computational domain suitable to unveil the potential occurrence of instabilities. Modelling results show good agreement with experimental measurements of voltage drop and the same trend but higher values of temperature. The higher temperatures obtained by the simulations appear to be due to the absence of a near-cathode heat dissipation model. Here, the results also reveal that the distribution of electron density and of the ratio of atomic helium ions to total ions transitions from monotonically increasing away from the cathode to presenting a minimum near the centre of the gap with increasing current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Combined crystal plasticity and grain boundary modeling of creep in ferritic-martensitic steels: II. The effect of stress and temperature on engineering and microstructural properties

This paper describes a series of physically-based crystal plasticity finite element method (CPFEM) simulations of long-term creep and creep rupture of Grade 91 steel. It is Part 2 of a two part series of papers. Part 1 describes the simulation framework; this part focuses on specific simulations and on how the predicted long-term creep properties of Grade 91 compare to the assumptions used in current high temperature design practices. This work extends the model developed in Part 1 to look at creep properties at different temperatures, principal stresses, and multiaxial stress states. The simulations show that empirically extrapolating creep rupture stresses from short-term experimental data may substantially over predict the actual long-term creep properties of Grade 91. Additionally, the CPFEM calculations predict a transition from notch strengthening creep behavior for high values of maximum principal stress and moderate notch severity to notch weakening behavior for low principal stresses and more severe notches. The latter regime better categorizes conditions in engineering components designed for long term elevated temperature use, which implies Grade 91 may be a notch weakening material in actual service. This would have a significant impact on high temperature design practices, though confirmatory test data on long-life, low stress notched specimens is difficult to obtain. Finally, one advantage of the physically-based modeling approach adopted here is that the simulation results also elucidate the microstructural mechanisms causing the macroscopic trends in engineering properties predicted by the simulations. This paper shows that the detailed micromechanical mechanisms predicted by the CPFEM simulations can be abstracted with a simple micromechanical model that can be used to both explain the detailed results and make improved predictions of engineering properties from experimental data.

36 MATERIALS SCIENCE↗

DOE OSTI2017

HIS-TORIC: extending core ICRF wave simulation to include realistic SOL plasmas

Here, this paper presents a novel approach to incorporating an arbitrarily shaped edge scrape-off-layer (SOL) plasma and an ion cyclotron range of frequency (ICRF) antenna structure into existing core ICRF wave simulation models. We partition the entire computation domain into two sub-domains: a core and an edge region. Simulations in each domain are performed separately with appropriate numerical solvers. For the core, the TORIC ICRF solver (Brambilla 1999 Plasma Phys. Control. Fusion 41 1) was modified to impose an essential (Dirichlet) boundary condition at its interface with the edge domain. In the edge, a finite element method is used to solve a cold collisional plasma model. The domains are then joined together using the continuity boundary condition for the tangential electric and magnetic fields at their interfaces (Hybrid Integration of SOL to TORIC: HIS-TORIC). The model developed here was tested using an ICRH H minority heating scenario on the Alcator C-Mod tokamak (Hutchinson et al 1994 Phys. Plasmas 1 1511). The simulated pattern of core wave propagation agrees well with a standard TORIC simulation. This approach opens the possibility of using a realistic diverted SOL plasma and a complicated 3D RF antenna together with a rigorous hot core plasma model, while requiring only minimal modification to existing RF codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

3D full wave fast wave modeling with realistic HHFW antenna geometry and SOL plasma in NSTX-U

This paper reports the significant advancement of our ability to model and to understand how RF waves interact with the SOL plasma, by developing for the full torus 3D scrape-off layer (SOL) plasma simulation together with the antenna and core plasma in the NSTX-U device in HHFW frequency regime for a cold plasma model. The present paper extends the previous modeling (Bertelli et al 2020 AIP Conf. Proc. 2254 030001) to a full 3D device geometry including realistic antenna geometry in order to capture a more realistic experimental situation including the fully 3D wave field effects and the antenna plasma interaction in the SOL plasma, and, at the same time, resulting core wave propagation. Here, the central tool of the present work is the Petra-M code, which is a newly developed state-of-the-art generic electromagnetic simulation tool for modeling RF wave propagation based on MFEM (http://mfem.org), an open source scalable C++ finite element method library.

3D RF modeling↗

DOE OSTI2023

Magnetic potential based formulation for linear and non-linear 3D RF sheath simulation

Abstract This paper reports a new numerical scheme to simulate the radio-frequency (RF) induced RF sheath, which is suitable for a large 3D simulation. In the RF sheath boundary model, the tangential component of the electric field ( E t ) is given by the gradient of a scalar electric field potential. We introduce two additional scalar potentials for the tangential components of the magnetic field, which effectively impose the normal electric displacement ( D n ) on the plasma sheath boundary condition via in-homogeneous Neumann boundary condition and constrain the tangential electric field on the surface as curl-free ( ∇ × E t = 0 ). In our approach, the non-linear sheath impedance is formulated as a natural extension of the large thickness (or asymptotic) sheath limit ( D n = 0 ), allowing for handling both asymptotic and non-linear regimes seamlessly. The new scheme is implemented using the Petra-M finite element method analysis framework and is verified with simulations in the literature. The significance of non-linearity is discussed in various plasma conditions. An application of this scheme to asymptotic RF sheath simulation on the WEST ICRF antenna side limiters is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗