Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphysics”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A Novel Partitioned Approach for Reduced Order Model—Finite Element Model (ROM-FEM) and ROM-ROM Coupling

Partitioned methods allow one to build a simulation capability for coupled problems by reusing existing single-component codes. In so doing, partitioned methods can shorten code development and validation times for multiphysics and multiscale applications. In this work, we consider a scenario in which one or more of the “codes” being coupled are projection-based reduced order models (ROMs), introduced to lower the computational cost associated with a particular component. We simulate this scenario by considering a model interface problem that is discretized independently on two non-overlapping subdomains. Here we then formulate a partitioned scheme for this problem that allows the coupling between a ROM “code” for one of the subdomains with a finite element model (FEM) or ROM “code” for the other subdomain. The ROM “codes” are constructed by performing proper orthogonal decomposition (POD) on a snapshot ensemble to obtain a low-dimensional reduced order basis, followed by a Galerkin projection onto this basis. The ROM and/or FEM “codes” on each subdomain are then coupled using a Lagrange multiplier representing the interface flux. To partition the resulting monolithic problem, we first eliminate the flux through a dual Schur complement. Application of an explicit time integration scheme to the transformed monolithic problem decouples the subdomain equations, allowing their independent solution for the next time step. We show numerical results that demonstrate the proposed method’s efficacy in achieving both ROM-FEM and ROM-ROM coupling.

de Castro, Amy↗
DOE OSTI2014

Application of the Ta liner technique to produce Ca beams at INFN-Legnaro National Laboratories (INFN-LNL)

The ECR ion sources are able to produce a wide variety of highly charged metallic ion beams thanks to the development of different techniques (ovens, sputtering, direct insertion, metal ions from volatile compounds (MIVOC)). In the case of the ovens, the sticking of the hot vapors on the surface of the plasma chamber leads to high material consumption rates. For elements like Ca, a tantalum liner inserted inside the chamber can be used to limit this phenomenon. The modeling of temperature distribution inside the chamber with and without the liner was carried out with COMSOL-multiphysics code. Results of simulation and the comparison with experiments performed at INFN-Legnaro National Laboratories with Ca beams are discussed.

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Full-wave feasibility study of magnetic diagnostic based on O-X mode conversion and oblique reflectometry imaging

An innovative millimeter wave diagnostic is proposed to measure the local magnetic field and the edge current as a function of the minor radius in the pedestal region. The idea behind such diagnostic is to localize and characterize a direction of reduced reflectivity at the O-mode cutoff layer. We modeled the wave scattering and mode-conversion processes by means of the finite-element COMSOL Multiphysics code in two dimensions (2D). Sensitivity studies were performed for parameters mocking up DIII-D plasmas. Simulations confirmed the presence of a minimum in reflectivity of an externally injected O-mode beam, and confirmed that this minimum depends on the magnetic field at the cutoff, as expected from the OX mode conversion physics. This study gives confidence in the feasibility of the diagnostic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Full wave simulation of waves in ECRIS plasmas based on the finite element method

This paper describes the modeling and the full wave numerical simulation of electromagnetic waves propagation and absorption in an anisotropic magnetized plasma filling the resonant cavity of an electron cyclotron resonance ion source (ECRIS). The model assumes inhomogeneous, dispersive and tensorial constitutive relations. Maxwell's equations are solved by the finite element method (FEM), using the COMSOL Multiphysics{sup ®} suite. All the relevant details have been considered in the model, including the non uniform external magnetostatic field used for plasma confinement, the local electron density profile resulting in the full-3D non uniform magnetized plasma complex dielectric tensor. The more accurate plasma simulations clearly show the importance of cavity effect on wave propagation and the effects of a resonant surface. These studies are the pillars for an improved ECRIS plasma modeling, that is mandatory to optimize the ion source output (beam intensity distribution and charge state, especially). Any new project concerning the advanced ECRIS design will take benefit by an adequate modeling of self-consistent wave absorption simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Measurement strategy for rectangular electrical capacitance tomography sensor

To investigate the influence of the measurement strategy for the rectangular electrical capacitance tomography (ECT) sensor, a Finite Element Method (FEM) is utilized to create the model for simulation. The simulation was carried out using COMSOL Multiphysics(trade mark, serif) and Matlab(trade mark, serif). The length-width ratio of the rectangular sensing area is 5. Twelve electrodes are evenly arranged surrounding the pipe. The covering ratio of the electrodes is 90%. The capacitances between different electrode pairs are calculated for a bar distribution. The air of the relative permittivity 1.0 and the material of the permittivity 3.0 are used for the calibration. The relative permittivity of the second phase is 3.0. The noise free and noise data are used for the image reconstruction using the Linear Back Projection (LBP). The measurement strategies with 1-, 2- and 4- electrode excitation are compared using the correlation coefficient. Preliminary results show that the measurement strategy with 2-electrode excitation outperforms other measurement strategies with 1- or 4-electrode excitation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Optothermal transport behavior in whispering gallery mode optical cavities

Over the past century, whispering gallery mode optical cavities have enabled numerous advances in science and engineering, such as discoveries in quantum mechanics and non-linear optics, as well as the development of optical gyroscopes and add drop filters. One reason for their widespread appeal is their ability to confine light for long periods of time, resulting in high circulating intensities. However, when sufficiently large amounts of optical power are coupled into these cavities, they begin to experience optothermal or photothermal behavior, in which the optical energy is converted into heat. Above the optothermal threshold, the resonance behavior is no longer solely defined by electromagnetics. Previous work has primarily focused on the role of the optothermal coefficient of the material in this instability. However, the physics of this optothermal behavior is significantly more complex. In the present work, we develop a predictive theory based on a generalizable analytical expression in combination with a geometry-specific COMSOL Multiphysics finite element method model. The simulation couples the optical and thermal physics components, accounting for geometry variations as well as the temporal and spatial profile of the optical field. To experimentally verify our theoretical model, the optothermal thresholds of a series of silica toroidal resonant cavities are characterized at different wavelengths (visible through near-infrared) and using different device geometries. The silica toroid offers a particularly rigorous case study for the developed optothermal model because of its complex geometrical structure which provides multiple thermal transport paths.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Determination of surface electric charge profile in pyroelectric crystals

Pyroelectric crystals are used to produce high energy self-focused electron beams. Here, an experimental analysis in combination with simulation studies will be reported to investigate possible sources of this effect. In the experiments, the surface of crystal was divided into six separated parts and the rate of surface electric charge production was measured accordingly. A non-steady and spatially non-uniform distribution of the surface charge generation was observed, in which it tends to a uniform distribution in the course of experiment. The obtained surface electric charges from the experiments were used to simulate the electric field and potential around the crystal by COMSOL Multiphysics. It was observed that emitted electrons from the crystal surface were focused, and the non-uniformity in spatial charge is responsible for this phenomenon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Characteristics of temperature rise in variable inductor employing magnetorheological fluid driven by a high-frequency pulsed voltage source

A variable inductor with magnetorheological (MR) fluid has been successfully applied to power electronics applications; however, its thermal characteristics have not been investigated. To evaluate the performance of the variable inductor with respect to temperature, we measured the characteristics of temperature rise and developed a numerical analysis technique. The characteristics of temperature rise were determined experimentally and verified numerically by adopting a multiphysics analysis technique. In order to accurately estimate the temperature distribution in a variable inductor with an MR fluid-gap, the thermal solver should import the heat source from the electromagnetic solver to solve the eddy current problem. To improve accuracy, the B–H curves of the MR fluid under operating temperature were obtained using the magnetic property measurement system. In addition, the Steinmetz equation was applied to evaluate the core loss in a ferrite core. The predicted temperature rise for a variable inductor showed good agreement with the experimental data and the developed numerical technique can be employed to design a variable inductor with a high-frequency pulsed voltage source.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Progress towards developing neutron tolerant magnetostrictive and piezoelectric transducers

Current generation light water reactors (LWRs), sodium cooled fast reactors (SFRs), small modular reactors (SMRs), and next generation nuclear plants (NGNPs) produce harsh environments in and near the reactor core that can severely tax material performance and limit component operational life. To address this issue, several Department of Energy Office of Nuclear Energy (DOE-NE) research programs are evaluating the long duration irradiation performance of fuel and structural materials used in existing and new reactors. In order to maximize the amount of information obtained from Material Testing Reactor (MTR) irradiations, DOE is also funding development of enhanced instrumentation that will be able to obtain in-situ, real-time data on key material characteristics and properties, with unprecedented accuracy and resolution. Such data are required to validate new multi-scale, multi-physics modeling tools under development as part of a science-based, engineering driven approach to reactor development. It is not feasible to obtain high resolution/microscale data with the current state of instrumentation technology. However, ultrasound-based sensors offer the ability to obtain such data if it is demonstrated that these sensors and their associated transducers are resistant to high neutron flux, high gamma radiation, and high temperature. To address this need, the Advanced Test Reactor National Scientific User Facility (ATR-NSUF) is funding an irradiation, led by PSU, at the Massachusetts Institute of Technology Research Reactor to test the survivability of ultrasound transducers. As part of this effort, PSU and collaborators have designed, fabricated, and provided piezoelectric and magnetostrictive transducers that are optimized to perform in harsh, high flux, environments. Four piezoelectric transducers were fabricated with either aluminum nitride, zinc oxide, or bismuth titanate as the active element that were coupled to either Kovar or aluminum waveguides and two magnetostrictive transducers were fabricated with Remendur or Galfenol as the active elements. Pulse-echo ultrasonic measurements of these transducers are made in-situ. This paper will present an overview of the test design including selection criteria for candidate materials and optimization of test assembly parameters, data obtained from both out-of-pile and in-pile testing at elevated temperatures, and an assessment based on initial data of the expected performance of ultrasonic devices in irradiation conditions.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2015

Measurement and simulation of thermoelectric efficiency for single leg

Thermoelectric efficiency measurements were carried out on n-type bismuth telluride legs with the hot-side temperature at 100 and 150°C. The electric power and heat flow were measured individually. Water coolant was utilized to maintain the cold-side temperature and to measure heat flow out of the cold side. Leg length and vacuum pressure were studied in terms of temperature difference across the leg, open-circuit voltage, internal resistance, and heat flow. Finite-element simulation on thermoelectric generation was performed in COMSOL Multiphysics, by inputting two-side temperatures and thermoelectric material properties. The open-circuit voltage and resistance were in good agreement between the measurement and simulation. Much larger heat flows were found in measurements, since they were comprised of conductive, convective, and radiative contributions. Parasitic heat flow was measured in the absence of bismuth telluride leg, and the conductive heat flow was then available. Finally, the maximum thermoelectric efficiency was derived in accordance with the electric power and the conductive heat flow.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Oscillatory instability development in extraction system of a negative ion source

Conditions of oscillatory instability development in the extraction system of a negative hydrogen ion source based on a volume-produced plasma are studied. Such an ion source is characterized by the presence of the parent gas in the extraction system due to the leakage from the gas-discharge chamber. The secondary electrons in the area of the ion-optical system become the reason of oscillation appearance and possible beam current modulation. Analytically the range of the stable beam propagation is found. The instability increment is shown to be rather small. Maximum increment of the oscillations corresponds to the beam velocity equal to the thermal velocity of plasma electrons. The group velocity of the oscillations is close to the beam velocity so the oscillations are convective. Simulation of the low energy beam propagation is performed in COMSOL Multiphysics, the beam current modulation being observed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Mixed Pierce-two-stream instability development in an extraction system of a negative ion source

Mixed Pierce-two-stream instability may occur in an extraction system of a negative ion source based on a volume-produced plasma. The reasons for instability development are discussed. Analytically the conditions of unstable beam propagation are determined. The instability threshold is shown to be increased compared with the pure Pierce instability. The influence of inclined perturbations on the instability behavior is investigated. The numerical calculations are performed in COMSOL Multiphysics. The simulation results confirm the existence of such a mixed instability appearance that develops due to both the electrons of the external circuit and the background positive ions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Interaction between the lower hybrid wave and density fluctuations in the scrape-off layer

In the present paper, the perturbation of the launched power spectrum of the Lower Hybrid wave at the separatrix by electron density fluctuations in the scrape-off layer is investigated. Considering a slab geometry with magnetic field lines parallel to the toroidal direction, the full wave equation is solved using Comsol Multiphysics® for a fully active multi-junction like LH antenna made of two modules. When electron density fluctuations are incorporated in the dielectric tensor over a thin perturbed layer in front of the grill, it is shown that the power spectrum may be strongly modified from the antenna mouth to the plasma separatrix as the wave propagates. The diffraction effect leads to the appearance of multiple satellite lobes with randomly varying positions, a feature consistent with the recently developed model that has been applied successfully to high density discharges on the Tokamak Tore Supra corresponding to the large spectral gap regime [Decker J. et al. Phys. Plasma 21 (2014) 092504]. The perturbation is found to be maximum for the Fourier components of the fluctuating spectrum in the vicinity of the launched LH wavelength.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Geometric consideration of nanostructures for energy storage systems

Battery performance and its fade are determined by various aspects such as the transport of ions and electrons through heterogeneous internal structures; kinetic reactions at the interfaces; and the corresponding interplay between mechanical, chemical, and thermal responses. The fundamental factor determining this complex multiscale and multiphysical nature of a battery is the geometry of active materials. In this work, we systematically consider the tradeoffs among a selection of limiting geometries of media designed to store ions or other species via a diffusion process. Specifically, we begin the investigation by considering diffusion in spheres, rods, and plates at the particle level, in order to assess the effects of geometry, diffusivity, and rate on capacity. Then, the study is extended to considering of the volume fraction and particle network, as well as kinetics at the interface with electrolyte. Our study suggests that, in terms of overall bulk level material performance, thin film batteries may generate the highest energy density with high power capability when they are implemented at nanoscales or with highly diffusion materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Retention modeling for ultra-thin density of Cu-based conductive bridge random access memory (CBRAM)

We investigate the effect of Cu concentration On-state resistance retention characteristics of W/Cu/Ti/HfO{sub 2}/Pt memory cell. The development of RRAM device for application depends on the understanding of the failure mechanism and the key parameters for device optimization. In this study, we develop analytical expression for cations (Cu{sup +}) diffusion model using Gaussian distribution for detailed analysis of data retention time at high temperature. It is found that the improvement of data retention time depends not only on the conductive filament (CF) size but also on Cu atoms concentration density in the CF. Based on the simulation result, better data retention time is observed for electron wave function associated with Cu{sup +} overlap and an extended state formation. This can be verified by analytical calculation of Cu atom defects inside the filament, based on Cu{sup +} diffusion model. The importance of Cu diffusion for the device reliability and the corresponding local temperature of the filament were analyzed by COMSOL Multiphysics simulation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Analysis of hydrogen plasma in a microwave plasma chemical vapor deposition reactor

The aim of this work is to build a numerical model of hydrogen plasma inside a microwave plasma chemical vapor deposition system. This model will help in understanding and optimizing the conditions for the growth of carbon nanostructures. A 2D axisymmetric model of the system is implemented using the finite element high frequency Maxwell solver and the heat transfer solver in COMSOL Multiphysics. The system is modeled to study variation in parameters with reactor geometry, microwave power, and gas pressure. The results are compared with experimental measurements from the Q-branch of the H{sub 2} Fulcher band of hydrogen using an optical emission spectroscopy technique. The parameter γ in Füner's model is calibrated to match experimental observations at a power of 500 W and 30 Torr. Good agreement is found between the modeling and experimental results for a wide range of powers and pressures. The gas temperature exhibits a weak dependence on power and a strong dependence on gas pressure. The inclusion of a vertical dielectric pillar that concentrates the plasma increases the maximum electron temperature by 70%, the maximum gas temperature by 50%, and the maximum electron number density by 70% when compared to conditions without the pillar at 500 W and 30 Torr. Experimental observations also indicate intensified plasma with the inclusion of a pillar.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Design and fabrication of non silicon substrate based MEMS energy harvester for arbitrary surface applications

Internet of Things (IoT) uses MEMS sensor nodes and actuators to sense and control objects through Internet. IOT deploys millions of chemical battery driven sensors at different locations which are not reliable many times because of frequent requirement of charging & battery replacement in case of underground laying, placement at harsh environmental conditions, huge count and difference between demand (24 % per year) and availability (energy density growing rate 8% per year). Energy harvester fabricated on silicon wafers have been widely used in manufacturing MEMS structures. These devices require complex fabrication processes, costly chemicals & clean room. In addition to this silicon wafer based devices are not suitable for curved surfaces like pipes, human bodies, organisms, or other arbitrary surface like clothes, structure surfaces which does not have flat and smooth surface always. Therefore, devices based on rigid silicon wafers are not suitable for these applications. Flexible structures are the key solution for this problems. Energy transduction mechanism generates power from free surrounding vibrations or impact. Sensor nodes application has been purposefully selected due to discrete power requirement at low duty cycle. Such nodes require an average power budget in the range of about 0.1 microwatt to 1 mW over a period of 3-5 seconds. Energy harvester is the best alternate source in contrast with battery for sensor node application. Novel design of Energy Harvester based on cheapest flexible non silicon substrate i.e. cellulose acetate substrate have been modeled, simulated and analyzed on COMSOL multiphysics and fabricated using sol-gel spin coating setup. Single cantilever based harvester generates 60-75 mV peak electric potential at 22Hz frequency and approximately 22 µW power at 1K-Ohm load. Cantilever array can be employed for generating higher voltage by replicating this structure. This work covers design, optimization, fabrication of harvester and schottky diodes based voltage multiplier.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Investigation of static and dynamic behavior of functionally graded piezoelectric actuated Poly-Si micro cantilever probe

In the present paper a novel functionally graded piezoelectric (FGP) actuated Poly-Si micro cantilever probe is proposed for atomic force microscope. The shear piezoelectric coefficient d{sub 15} has much higher value than coupling coefficients d{sub 31} and d{sub 33}, hence in the present work the micro cantilever beam actuated by d{sub 15} effect is utilized. The material properties are graded in the thickness direction of actuator by a simple power law. A three dimensional finite element analysis has been performed using COMSOL Multiphysics® (version 4.2) software. Tip deflection and free vibration analysis for the micro cantilever probe has been done. The results presented in the paper shall be useful in the design of micro cantilever probe and their subsequent utilization in atomic force microscopes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗