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

Engineering on-chip nanoporous gold material libraries via precision photothermal treatment [Precision Photothermal Annealing of Nanoporous Gold Thin Films for the Microfabrication of a Single-ship Material Libraries]

Single-chip material libraries of thin films of nanostructured materials are a promising approach for high throughput studies of structure-property relationship in the fields of physics and biology. Nanoporous gold (np-Au), produced by an alloy corrosion process, is a nanostructured material of specific interest in both these fields. One attractive property of np-Au is its self-similar coarsening behavior by thermally induced surface diffusion. However, traditional heat application techniques for the modification of np-Au are bulk processes that cannot be used to generate a library of different pore sizes on a single chip. Laser micromachining offers an attractive solution to this problem by providing a means to apply energy with high spatial and temporal resolution. In our present study we use finite element multiphysics simulations to predict the effects of laser mode (continuous-wave vs. pulsed) and supporting substrate thermal conductivity on the local np-Au film temperatures during photothermal annealing and subsequently investigate the mechanisms by which the np-Au network is coarsening. Our simulations predict that continuous-wave mode laser irradiation on a silicon supporting substrate supports the widest range of morphologies that can be created through the photothermal annealing of thin film np-Au. Using this result we successfully fabricate a single-chip material library consisting of 81 np-Au samples of 9 different morphologies for use in increased throughput material interaction studies.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A closer look at the physical and optical properties of gold nanostars: an experimental and computational study

A combined experimental and computational study was carried out to design a semi-empirical method to determine the volume, surface area, and extinction coefficients of gold nanostars. The values obtained were confirmed by reconstructing the nanostar 3D topography through high-tilt TEM tomography and introducing the finite elements in COMSOL Multiphysics through which we have also calculated the morphology-dependent extinction coefficient. We have, for the first time, modeled the heat losses of a real, experimentally synthesized nanostar, and found the plasmon resonances to be in excellent agreement with those obtained experimentally. Furthermore, we believe that our approach could substantially improve the applicability of this remarkable nanomaterial.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Modeling gas-diffusion electrodes for CO 2 reduction

CO 2 reduction conducted in electrochemical cells with planar electrodes immersed in an aqueous electrolyte is severely limited by mass transport across the hydrodynamic boundary layer. This limitation can be minimized by use of vapor-fed, gas-diffusion electrodes (GDEs), enabling current densities that are almost two orders of magnitude greater at the same applied cathode overpotential than what is achievable with planar electrodes in an aqueous electrolyte. The addition of porous cathode layers, however, introduces a number of parameters that need to be tuned in order to optimize the performance of the GDE cell. In this study, we develop a multiphysics model for gas diffusion electrodes for CO 2 reduction and used it to investigate the interplay between species transport and electrochemical reaction kinetics. The model demonstrates how the local environment near the catalyst layer, which is a function of the operating conditions, affects cell performance. Finally, we also examine the effects of catalyst layer hydrophobicity, loading, porosity, and electrolyte flowrate to help guide experimental design of vapor-fed CO 2 reduction cells.

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

Towards membrane-electrode assembly systems for CO 2 reduction: a modeling study

Membrane-electrode assemblies (MEAs) are an attractive cell design for the electrochemical reduction of CO 2 because they exhibit low ohmic loss and high energy efficiency. We describe here the development and application of a multiphysics model to investigate the fundamental limitations of two MEA designs: one with gaseous feeds at both the anode and cathode (full-MEA), and the other with an aqueous anode feed (KHCO 3 or KOH exchange solution) and a gaseous cathode feed (exchange-MEA). The total current density for the three cases follows the order: KOH-MEA > KHCO 3 -MEA > full-MEA. This trend is established by examining the distribution of the applied voltage. We show that the main charge-carrying species are carbonate anions for an MEA that uses an anion-exchange membrane (AEM). The amount of CO 2 consumed but not converted to CO decreases with increasing current densities above 100 mA cm –2 for a full-MEA, but converges to 50% for exchange-MEAs. The full-MEA becomes limited by ohmic resistance as the membrane dehydrates with increasing cell temperature, and eventually becomes limited due to water mass transport. The exchange-MEAs can maintain membrane hydration and the local ion concentration at the anode, but are limited by salt precipitation at the cathode, as well as a higher tendency to flood. Lastly, we explore the effects of temperature and discuss the possibility of increasing water supply to the full-MEA to improve its performance at elevated temperatures. The MEA model and the understanding of MEA performance for the electrochemical reduction of CO 2 presented in this study should help guide the design of next-generation CO 2 reduction cells.

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

Direct observation of the local microenvironment in inhomogeneous CO 2 reduction gas diffusion electrodes via versatile pOH imaging

In this study, we report how the micrometer-scale morphology of a carbon dioxide reduction (CO 2 R) gas diffusion electrode (GDE) affects the mass transport properties and with it, the local CO 2 R performance. We developed a technique to probe the microenvironment in a CO 2 R GDE via local pOH imaging with time- and three-dimensional spatial, micrometer-scale resolution. The local activity of hydroxide anions (OH - ), represented by the pOH value, around a GDE in contact with an aqueous electrolyte is a crucial parameter that governs the catalytic activity and CO 2 R selectivity. Here, we use fluorescence confocal laser scanning microscopy (CLSM) to create maps of the local pOH around a copper GDE by combining two ratiometric fluorescent dyes, one of which is demonstrated as a pOH sensor for the first time in this work. We observe that the local pOH decreases when current is applied due to the creation of OH - as a byproduct of CO 2 R. Interestingly, the pOH is lower inside microtrenches compared to the electrode surface and decreases further as trenches become more narrow due to enhanced trapping of OH - . We support our experimental results with multiphysics simulations that correlate exceptionally well with measurements. These simulations additionally suggest that the decreased pOH inside microcavities in the surface of a CO 2 R GDE leads to locally enhanced selectivity towards multicarbon (C 2+ ) products. This study suggests that narrow microstructures on the length scale of 5 μm in a GDE surface serve as local CO 2 R hotspots, and thus highlights the importance of a GDE's micromorphology on the CO 2 R performance.

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

Designing optimal core–shell MOFs for direct air capture

Metal–organic frameworks (MOFs), along with other novel adsorbents, are frequently proposed as candidate materials to selectively adsorb CO 2 for carbon capture processes. However, adsorbents designed to strongly bind CO 2 nearly always bind H 2 O strongly (sometimes even more so). Given that water is present in significant quantities in the inlet streams of most carbon capture processes, a method that avoids H 2 O competition for the CO 2 binding sites would be technologically valuable. Here, in this paper, we consider a novel core–shell MOF design strategy, where a high-CO 2 -capacity MOF “core” is protected from competitive H 2 O-binding via a MOF “shell” that has very slow water diffusion. We consider a high-frequency adsorption/desorption cycle that regenerates the adsorbents before water can pass through the shell and enter the core. To identify optimal core–shell MOF pairs, we use a combination of experimental measurements, computational modeling, and multiphysics modeling. Our library of MOFs is created from two starting MOFs-UiO-66 and UiO-67-augmented with 30 possible functional group variations, yielding 1740 possible core–shell MOF pairs. After defining a performance score to rank these pairs, we identified 10 core–shell MOF candidates that significantly outperform any of the MOFs functioning alone.

42 ENGINEERING↗
DOE OSTI2026

In situ catalyst activation and regeneration enable energy-efficient high-current CO 2 reduction to ethanol-rich C 2+ mixtures

Electrochemical conversion of dissolved CO 2 in bicarbonate electrolytes, i.e., bicarbonate electrolysis, offers distinct advantages over gas diffusion electrode systems by enabling direct utilization of the CO 2 capture electrolyte while bypassing the energy-intensive CO 2 release step. However, bicarbonate electrolysis faces challenges such as CO 2 mass-transfer limitation, local pH-driven CO 2 depletion, and high cathodic potentials. The higher potential often causes catalyst surface reorganization, leading to a gradual loss of active sites and variations in selectivity during CO 2 reduction. Here, we report a directed, in situ activation and regeneration method that allows precatalysts to equilibrate under dynamic (pulsed) electrolysis conditions. We demonstrate in situ activation of a scalable Cu 2 O/Cu mesh that, under short-width (t = 4 s) pulsed electrolysis, provides stable mixed oxidation states of Cu, favoring the formation of an ethanol-rich crude mixture. The pulsed electrolysis waveform, consisting of six distinct segments, is tuned to form Cu + oxides, which are then reduced to generate local alkaline conditions favoring C–C coupling. This synergistic effect results in FEs of 73% for C2+ products and 39% for ethanol at an applied current density of −150 mA cm −2 and a cathodic potential of −1.45 V (vs. RHE). The overall half-cell energy efficiency is ∼30% for C 2+ products. The in situ Raman experiments confirm the role of pCO 2 R in dynamically regenerating Cu+-containing surface species during pulsed operation, thereby steering selectivity towards C 2+ products. A comprehensive multiscale, multiphysics model is developed to investigate the dynamic behavior of copper surface species (Cu, Cu + , and Cu 2+ ) and local microenvironmental conditions during the pCO 2 R. The results reveal that the coexistence of different copper oxidation states, especially the Cu+ intermediate, is critical in steering selectivity towards multicarbon (C 2+ ) products. The dynamic modulation of surface redox states via tailored pulsing strategies favors C–C coupling pathways by inducing localized alkaline conditions and stabilizing reactive intermediates. This work establishes a predictive modeling platform that links pulse waveform design with mechanistic insights into catalyst state evolution and product selectivity. Overall, this study provides valuable insights into the synergistic effect of in situ activation of pre-catalysts and pulsed electrolysis for higher selectivity towards C 2+ products.

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

An analysis of parameters affecting ampacity in aircraft bipolar MVDC power cables via coupled electrical, thermal, and computational fluid dynamic modelling

Abstract The next generation of aircraft, including more electric aircraft and all‐electric aircraft (AEA), requires electric power systems with high power density and low system mass specifications. Increasing the voltage of the system to the range of a few kV, medium voltage (MV), is a reasonable approach to achieving high‐power‐density and low‐system‐mass EPSs for aircraft applications. Higher voltages, however, pose many challenges for aviation MV power cables such as arcs and arc tracking, partial discharges (PDs), and thermal management. In this regard, thermal management is more challenging since heat transfer by convection is greatly reduced at wide‐body aircraft's cruising altitudes due to the reduced air pressure. In this paper, a finite element method (FEM) model is developed in COMSOL Multiphysics for an aircraft bipolar MVDC (±5 kV) power cable. Using the model, the maximum permissible cable current at a low pressure of 18.8 kPa (at an altitude of 12.2 km from sea level, the usual cruising altitude for wide‐body aircraft) is calculated. Also, an analytical model is developed based on analytical and proven empirical correlations governing conductive, radiative, and convective heat transfers at the steady state to estimate the ampacity of the bipolar cable system at reduced pressure. It was shown that the proposed analytical model can be used for atmospheric pressure and systems with a larger number of poles, expanding its range of applications. The results of the FEM and analytical models correlate at wide ranges of parameters such as ambient temperature, duct size, distance between the positive and negative pole cables, and the overall diameter of the cables. The influence of horizontal and vertical arrangement of poles is included in the analytical model. The results of this study can be used to design bipolar MVDC power cable systems for the envisaged wide‐body AEA.

Azizi, Arian↗
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 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

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↗