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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Modeling Planar Electrodes and Zero–Gap Membrane Electrode Assemblies for CO 2 Electrolysis

Multiphysics modeling enables probing of conditions inside a CO 2 electrolyzer that are difficult to measure, such as local concentrations and pH, as well as rapid testing of possible design changes. A one-dimensional model for a zero-gap membrane electrode assembly (MEA) CO 2 electrolyzer was developed with the assumption that catalyst layers interact with the membrane ionomer such that the ionomer affects the underlying kinetics. The kinetics for bicarbonate reacting to form hydrogen are fit using a planar electrode model for silver with an ionomer coating. The MEA model results are validated against experimental studies for current density and product selectivity. Flooding of the cathode is modeled using saturation curves, and results show that blocked pores in the microporous layer play a significant role in limiting the mass transport at high potentials (>2.8 V). Sensitivity studies showed that CO Faradaic efficiency can be increased by decreasing catalyst layer thickness and porosity, and decreasing KHCO 3 concentration.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2016

Microscale modeling of effective mechanical and electrical properties of textiles

Here, a computational framework for assisting in the development of novel textiles is presented. Electronic textiles are key in the rapidly growing field of wearable electronics for both consumer and military uses. There are two main challenges to the modeling of electronic textiles: the discretization of the textile microstructure and the interaction between electromagnetic and mechanical fields. A director-based beam formulation with an assumed electrical current is used to discretize the fabric at the level of individual fibrils. The open-source package FEniCS was used to implement the finite element model. Contact integrals were added into the FEniCS framework so that multiphysics contact laws can be incorporated in the same framework, leveraging the code generation and automated differentiation capabilities of FEniCS to produce the tangents needed by the implicit solution method. The computational model is used to construct and determine the mechanical, thermal, and electrical properties of a representative volume elements of a plain woven textile. Dynamic relaxation to solve the mechanical fields and the electrical and thermal fields is solved statically for a given mechanical state. The simulated electrical responses are fit to a simplified Kirchhoff network model to determine effective resistances of the textile.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The strip method for shape derivatives

Abstract A major challenge in shape optimization is the coupling of finite element method (FEM) codes in a way that facilitates efficient computation of shape derivatives. This is particularly difficult with multiphysics problems involving legacy codes, where the costs of implementing and maintaining shape derivative capabilities are prohibitive. The volume and boundary methods are two approaches to computing shape derivatives. Each has a major drawback: the boundary method is less accurate, while the volume method is more invasive to the FEM code. We introduce the strip method , which computes shape derivatives on a strip adjacent to the boundary. The strip method makes code coupling simple. Like the boundary method, it queries the state and adjoint solutions at quadrature nodes, but requires no knowledge of the FEM code implementations. At the same time, it exhibits the higher accuracy of the volume method. As an added benefit, its computational complexity is comparable to that of the boundary method, that is, it is faster than the volume method. We illustrate the benefits of the strip method with numerical examples.

Hardesty, Sean↗
DOE OSTI2025

A Fundamentally New Coupled Approach to Contact Mechanics via the Dirichlet‐Neumann Schwarz Alternating Method

Contact phenomena are crucial for understanding the behavior of mechanical systems. However, existing computational approaches for simulating mechanical contact often face numerical challenges, such as inaccurate physical predictions, energy conservation errors, and unwanted oscillations. Here, we introduce an alternative technique for simulating dynamic contact based on the non‐overlapping Schwarz alternating method, originally developed for domain decomposition. In multibody contact scenarios, this method treats each body as a separate, non‐overlapping domain and prevents interpenetration using an alternating Dirichlet–Neumann iterative process. This approach has a strong theoretical foundation, eliminates the need for contact constraints, and offers flexibility, making it ideal for multiscale and multiphysics applications. We conducted a numerical comparison between the Schwarz method and traditional methods, such as the Lagrange multiplier and penalty methods, focusing on a benchmark impact problem. Our results indicate that the Schwarz alternating method outperforms traditional methods in several key areas: it provides more accurate predictions for various measurable quantities and demonstrates exceptional energy conservation capabilities. To address unwanted oscillations in contact velocities and forces, we explored various algorithms and stabilization techniques, ultimately opting for the naïve‐stabilized Newmark scheme for its simplicity and effectiveness. Additionally, we validated the efficiency of the Schwarz method in a three‐dimensional impact problem, highlighting its inherent capacity to accommodate different mesh topologies, time‐integration schemes, and time steps for each interacting body.

Schwarz alternating method↗
DOE OSTI2016

Modeling quantum cascade lasers: Coupled electron and phonon transport far from equilibrium and across disparate spatial scales

Quantum cascade lasers (QCLs) are high-power coherent light sources in the midinfrared and terahertz parts of the electromagnetic spectrum. They are devices in which the electronic and lattice systems are far from equilibrium, strongly coupled to one another, and the problem bridges disparate spatial scales. Here, we present our ongoing work on the multiphysics and multiscale simulation of far-from-equilibrium transport of charge and heat in midinfrared QCLs.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Ultrasound Manipulation and Extrusion of Active Nanorods

Abstract Synthetic self‐propelled nano and microparticles have a growing appeal for targeted drug delivery, collective functionality, and manipulation at the nanoscale. However, it is challenging to control their positions and orientations under confinement, e.g., in microchannels, nozzles, and microcapillaries. This study reports on the synergistic effect of acoustic and flow‐induced focusing in microfluidic nozzles. In a microchannel with a nozzle, the balance between the acoustophoretic forces and the fluid drag due to streaming flows generated by the acoustic field controls the microparticle's dynamics. This study manipulates the positions and orientations of dispersed particles and dense clusters inside the channel at a fixed frequency by tuning the acoustic intensity. The main findings are: first, this study successfully manipulates the positions and orientations of individual particles and dense clusters inside the channel at a fixed frequency by tuning the acoustic intensity. Second, when an external flow is applied, the acoustic field separates and selectively extrudes shape‐anisotropic passive particles and self‐propelled active nanorods. Finally, the observed phenomena are explained by multiphysics finite‐element modeling. The results shed light on the control and extrusion of active particles in confined geometries and enable applications for acoustic cargo (e.g., drug) delivery, particle injection, and additive manufacturing via printed self‐propelled active particles.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Thermodynamic Modeling of Nuclear Fuel Materials

Thermodynamics is essential, fundamental information for understanding the behavior of any system. For nuclear materials, the situation is even more critical. The fuel in particular becomes highly multicomponent with burnup and is constantly evolving from beginning to end of life. Calculation of Phase Diagram (CALPHAD) modeling gives a continuous thermodynamic description of a system and allows coupling of the information to other phenomena facilitating multiphysics simulations.

Mcmurray, Jake W.↗
DOE OSTI2017

Coupled Physics Simulation of Expansive Reactions in Concrete with the Grizzly Code

The Grizzly code is being developed under the US Department of Energy’s Light Water Sustainability program as a tool to model aging mechanisms and their effects on the integrity of critical nuclear power plant components. An important application for Grizzly is the modeling of aging in concrete structures, which can be due to a number of mechanisms. Initial focus in this area has been on modeling expansive reactions due to alkali-silica reactions or radiation-induced volumetric expansion. Grizzly is an inherently multiphysics modeling platform that naturally permits including the effects of multiple coupled physics in a simulation. Models have been developed for transport of heat and moisture in concrete, and these have been coupled and used as inputs to models for expansive reactions. This paper summarizes this capability as well as work currently underway to couple it with mechanical deformation models that account for nonlinear constitutive behavior.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Generalized thermo-mechanical framework for heterogeneous materials through asymptotic homogenization

Abstract A fundamental understanding of the interaction between microstructure and underlying physical mechanisms is essential, especially for developing more accurate multi-physics models for heterogeneous materials. Effects of microstructure on the material response at the macroscale are modeled by using the generalized thermomechanics. In this study, strain gradient theory is employed as a higher-order theory on the macroscale with thermodynamics modeled as a first-order theory on the microscale. Hence, energy depends only on the temperature such that we circumvent an extension of Fourier’s law and analyze the “simplest” thermo-mechanical model in strain gradient elasticity. Developing multiphysics models for heterogeneous materials is indeed a challenge and even this “simplest” model in generalized thermomechanics creates dozens of parameters to be determined. We develop a thermo-mechanical framework, in which microstructure is modeled as a periodic structure and through asymptotic homogenization approach, higher-order parameters at macroscopic scale are calculated. To illustrate the importance of higher-order parameters in overall thermo-mechanical response of a heterogeneous materials, finite element method (FEM) is employed with the aid of open-source codes (FEniCS). Verification example of a bulk system and several case studies of porous structures demonstrate how such numerical framework can be beneficial in the design of materials with tailored microstructures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Fast and accurate reduced-order modeling of a MOOSE-based additive manufacturing model with operator learning

One predominant challenge in additive manufacturing (AM) is to achieve specific material properties by manipulating manufacturing process parameters during the runtime. Such manipulation tends to increase the computational load imposed on existing simulation tools employed in AM. The goal of the present work is to construct a fast and accurate reduced-order model (ROM) for an AM model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, ultimately reducing the time/cost of AM control and optimization processes. Our adoption of the operator learning (OL) approach enabled us to learn a family of differential equations produced by altering process variables in the laser’s Gaussian point heat source. More specifically, we used the Fourier neural operator (FNO) and deep operator network (DeepONet) to develop ROMs for time-dependent responses. Furthermore, we benchmarked the performance of these OL methods against a conventional deep neural network (DNN)-based ROM. Ultimately, we found that OL methods offer comparable performance and, in terms of accuracy and generalizability, even outperform DNN at predicting scalar model responses. The DNN-based ROM afforded the fastest training time. Furthermore, all the ROMs were faster than the original MOOSE model yet still provided accurate predictions. FNO had a smaller mean prediction error than DeepONet, with a larger variance for time-dependent responses. Unlike DNN, both FNO and DeepONet were able to simulate time series data without the need for dimensionality reduction techniques. Finally, the present work can help facilitate the AM optimization process by enabling faster execution of simulation tools while still preserving evaluation accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

An Enhanced Virtual Crack Closure Technique for Stress Intensity Factor Calculation along Arbitrary Crack Fronts and the Application in Hydraulic Fracturing Simulation

The virtual crack closure technique (VCCT) is widely used for calculating energy release rates along crack fronts and modeling the propagation of cracks in solid materials. Although the VCCT formulation for smooth crack fronts has been sufficiently addressed in the literature, the application of VCCT to a nonsmoothed crack front with sharp corners warrants further investigation. In this study, we describe an enhanced VCCT to calculate energy release rates and stress intensity factors for cracks with arbitrary shapes in 3D domains discretized on structured grids. The formulations of the enhanced VCCT were developed and implemented into a multiphysics simulation environment capable of simulating crack propagation in the framework of linear elastic fracture mechanics. Comparisons with existing analytical/numerical solutions and other VCCT approaches were performed to verify the enhanced VCCT in terms of SIF calculation along nonsmoothed crack fronts. We then applied the enhanced VCCT to a hydraulically driven penny-shaped fracture problem to further demonstrate its capability to simulate nonsmoothed fracture propagation.

02 PETROLEUM↗
DOE OSTI2024

Transition in ICF Capsule Implosions

Longstanding design and reproducibility challenges in inertial confinement fusion (ICF) capsule implosion experiments involve recognizing the need for appropriately characterized and modeled three-dimensional initial conditions and high-fidelity simulation capabilities to predict transitional flow approaching turbulence, material mixing characteristics, and late-time quantities of interest—e.g., fusion yield. We build on previous coarse graining simulations of the indirect-drive national ignition facility (NIF) cryogenic capsule N170601 experiment-a precursor of N221205 which resulted in net energy gain. We apply effectively combined initialization aspects and multiphysics coupling in conjunction with newly available hydrodynamics simulation methods, including directional unsplit algorithms and low Mach-number correction-key advances enabling high fidelity coarse grained simulations of radiation-hydrodynamics driven transition. Our presentation includes discussion of the capsule initialization and implosion dynamics, analysis of the vorticity production budget, transition signatures, quantities of interest—late-time ion temperature and fusion-neutron yield, numerical uncertainty quantification, and comparisons with NIF data.

42 ENGINEERING↗
DOE OSTI2026

Molten Halide Salt Surface Tension: Methods and Correlations

Here, this paper reviews various methods for studying surface tension and their applicability to fluoride and chloride molten salt systems, including a comparison of benefits and drawbacks. Such a comparison aids in experiment design based on desired factors such as scale, accuracy, and repeatability. A detailed review is presented for existing literature data regarding the surface tension of molten fluoride and chloride salts. These reference data were compiled and analyzed to determine cross-validated correlation equations for several alkali and alkaline earth fluoride and chloride salts as functions of temperature. These correlations are necessary for reliable multiphysics modeling approaches as well as accurate design and analysis of multiphase molten salt phenomena such as gas sparging and bubble formation/transport. This analysis supports the development of the thermophysical arm of the Molten Salt Thermal Properties Database (MSTDB-TP) managed by Oak Ridge National Laboratory.

Chloride↗
DOE OSTI2018

Electromagnetic Simulations of Newly Designed Semiconductor Bolometers for Submillimeter Observations

In the frame of the SPICA mission, new semiconductor bolometers have been designed for the POL instrument which is directly part of the HERSCHEL/PACS legacy. These detectors are cooled down to 50 mK and are sensitive to polarization for the 100, 200 and 350 μm wavelength bands. In this paper, 3D simulations have been carried out with the COMSOL Multiphysics software to design the two suspended interlaced silicon spirals which support the orthogonal absorbers (sensitive to the two polarizations) of the POL detectors. These simulations allow to optimize the design of the absorbers and to quantify the cross-polarization. In the last part of this paper, we explain how it is possible to reach the 200 and 350 μm bands from the layout of the 100-μm detector.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Air-Coupled Nondestructive Evaluation Dissected

This work develops a two-dimensional theoretical model to simulate the behavior of a fully non-contact air-coupled nondestructive evaluation system for a thin isotropic plate. The model is divided into transmission, guided wave propagation and reception phase. The validation of the complete model was carried out by modeling the same system by means of finite element method using a Multiphysics software. In addition, the dependency of the generated Lamb waves on different transmitter’s parameters and incidence angle is thoroughly investigated. The results of the acoustic pressure excited by the transducer, the out-of-plane velocity amplitudes for the generated first antisymmetric Lamb wave mode, and the radiated pressure from the plate caused by the leaky Lamb wavefield were all compared between the two models and a reasonable degree of similarity was found.

42 ENGINEERING↗
DOE OSTI2015

Magnetic Properties of Electroformed Ni and Ni-Fe for Micromagnetic MEMS Applications

This paper reports a new technique which is based on electroforming process for high permeability Ni and Ni-Fe alloy. The magnetic properties of Ni and Ni-Fe alloy are investigated as a function of plating current density via scanning electron microscope (SEM), energy dispersive spectroscopy (EDS), and vibrating sample magnetometer (VSM). Concentration-dependent kinetics is used to formulate the deposition rate and also the same is verified by using finite element analysis tool COMSOL Multiphysics. Several experiments are performed to derive the relation between electroplating parameters and magnetic properties. The process conditions are optimized to develop high permeability (∼16000) Ni-Fe alloy (85.71–14.29 %). The higher current density shows a monotonically increasing behavior of relative permeability (μ{sub r}). In addition, the B-H curves indicate that the saturation magnetization (B{sub sat}) increases and coercivity decreases with increasing plating current density. Our experimental data confirms that the higher current density increase the softness of Ni. The observed properties of Ni and Ni-Fe based on experimental results can be used to study the magnetic behavior of micromagnetic sensors and large force actuators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Effect Laws of Different Factors on Levitation Characteristics of High-T{sub c} Superconducting Maglev System with Numerical Solutions

The interaction between the high-temperature superconducting (HTS) bulk and permanent magnet guideway (PMG) has been a topic among the existing HTS magnetic levitation systems. Up to now, the theories about the electromagnetic and force characteristics have been relatively consummated. Here, using the finite element model by COMSOL Multiphysics 5.3a, we investigated the impact of the intrinsic and external factors on the levitation characteristics of HTS maglev systems, the application laws are concluded. Factors including the critical current density and geometric dimensions of the bulk superconductor, as well as the moving velocity, are discussed. The effects of the multiple back and forth movements are also analyzed. By comparing calculation results of the distributions of the induced current and the levitation force of the bulk superconductor, we can clearly confirm the effect laws of various factors on the electromagnetic characteristics and levitation force. These observations support the physical parameters that are difficult to be tested by experiments. And they also have a significant contribution to the experimental condition selections and optimization of the high-temperature superconducting magnetic levitation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Twin-core sunflower-type photonic quasicrystal fibers incorporated gold, silver, and copper microwire: an ultrashort and broad bandwidth polarization splitter

Polarization splitters are beneficial in the field of passive optical devices which are principally used in the coherent optical communication systems. An ideal polarization splitter must possess an efficient coupling with significant optical bandwidth. Hence, we proposed a novel metal-filled polarization splitter based on twin-core sunflower-type photonic quasicrystal fiber (PQF). The characteristics of the proposed splitter such as coupling length, extinction ratio, and available optical bandwidth were investigated by using the commercially available software, COMSOL Multiphysics. These aforementioned characteristics were determined by optimizing the geometrical parameters of proposed PQF—first by using gold-filled PQF (AuPQF). With optimized geometrical parameters, the splitter achieved an ultra-short coupling length of 132 µm and an ultra-broad bandwidth of 301 nm (wavelength range 1.36–1.68 µm) with an extinction ratio below than − 20 dB which comprised E + S+C + L+U band. Secondly, a comparative study to use silver- and copper-filled PQF (AgPQF and CuPQF) was explored. Both AgPQF and CuPQF retained a short length with adequate bandwidth and a significant extinction ratio below than − 20 dB. The AgPQF has a coupling length of 137 µm with a bandwidth of 136 nm which operated at 1.44 µm wavelength, while CuPQF with 134 µm lengths with a bandwidth of 149 nm and operated at 1.50 µm wavelengths.

36 MATERIALS SCIENCE↗