Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMSOL 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 73 records · Page 4

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↗

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↗

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↗

Experimental and Mathematical Simulation Study on the Granulation of a Modified Basic Oxygen Furnace Steel Slag

Basic oxygen furnace (BOF) steel slag is a major waste of the steelmaking industry. Utilization of BOF slag contributes to the sustainability of the steel industry by alleviating its environmental impact. Vitrification during cooling is an effective method to promote the cementitious activity of slags with the aim to apply slags in high value-added applications. In the present study, an Al{sub 2}O{sub 3}- and SiO{sub 2}-modified BOF slag was water granulated at a pilot scale. The amorphous and mineral fractions were measured quantitatively. The critical cooling rate to vitrify the modified slag was calculated through the Time-Temperature-Transformation diagrams constructed using in situ confocal scanning laser microscopy. To provide an insight into the crystallization behavior during the granulation process, a mathematical model was developed. The model was validated by comparing the amorphous fraction obtained from experiments with that from simulation. Temperature profiles of the slag particles with varied sizes were calculated with the aid of COMSOL Multiphysics software. The effect of particle size on the vitrified fraction was discussed in detail and the temperature gradient from surface to center of the particle was identified. The results provide novel fundamental understanding of the vitrification process regarding to the slag valorization, which will help the industrial implementation of granulation for pollution remediation and other engineered environmental systems such as chemical industry and pharmaceutical industry.

36 MATERIALS SCIENCE↗

Research on the Corrosion Behavior of X70 Pipeline Steel Under Coupling Effect of AC + DC and Stress

The coupling effect of stray AC (alternating current) + DC (direct current) and tensile stress-induced X70 pipeline corrosion is complicated with its mechanism uncovered. In this work, the corrosion behavior of X70 pipeline steel under three different corrosion (including an individual AC corrosion, mixed AC + 150 MPa tensile stress corrosion, coupling effect of AC + a 30 A/m{sup 2} DC interference and 150 MPa tensile stress corrosion) conditions in Dagang soil simulated solution was investigated by COMSOL Multiphysics software simulation, electrochemical tests, scanning electron microscope and other techniques. The results indicated that with the increase in AC density, the open-circuit potential of the X70 steel samples was negatively shifted and corrosion rate increased in three different corrosion conditions. The corrosion rates of X70 steel samples under individual AC interference were almost unchanged when 150.9 MPa tensile stress was applied. When applying 30 A/m{sup 2} DC interference, the corrosion rate of X70 steel increases significantly. Therefore, the factors affecting corrosion performance of X70 pipeline steel were DC > AC > stress. Under individual AC interference, uniform corrosion occurred on the X70 steel at an AC current density of less than 30 A/m{sup 2}, while the pitting corrosion occurred at an AC current density of greater than 100 A/m{sup 2}. The application of tensile stress and DC had a certain inhibitory effect on pitting corrosion.

36 MATERIALS SCIENCE↗

Solving macroscopic and microscopic pin-fin heat sink problems by adapted spectral method

In this paper, simulation of heat transfer in a heat sink with macroscopic and microscopic scales when one pin-fin is added to the system is investigated by the proposed spectral method. In the microscopic problems, heat transfer model uses dual-phase lag formulas in contrast with macroscopic problems when Fourier law is used to formulate the governing equation. In macroscopic problem, the results are compared with COMSOL multiphysics software results and a good agreement between the results are shown. In microscopic problems, 3D Gaussian heat source is used and boundary conditions obey the Newton law. Comparisons show the efficiency of the current method, while the results are compared with existed literature. It is shown that: by adding one pin fin to the heat sink problem, temperature of the system will decrease. Numerical results are given and graphs of the temperature for various kinds of pin-fin heat sink problems are depicted.

97 MATHEMATICS AND COMPUTING↗

Model-based optimization strategies for direct hydrogenation of carbon dioxide to dimethyl ether

Here, this study discusses model-based optimization strategies for CO 2 hydrogenation to dimethyl ether (DME) over aCuZnZr (CZZ) and ferrite (FER) mixed catalyst system in a packed-bed reactor configuration. A two-dimensional axisymmetric, nonisothermal packed-bed reactor model was developed using COMSOL Multiphysics 6.2 software. The model solves two-dimensional (radial and axial) heat and mass transport equations in the packed-bed and integrates intraparticle diffusion and heat transfer in a 1D approach. This powerful feature differs from a traditional porous media approach and takes into account any heat and mass transfer limitations that may exist. Analysis shows that the heat transfer limitations are negligible, but strong internal mass transfer limitations were observed at 10 ≤ WHSV ≤ 90 h -1 on the FER catalyst and at 240 °C. The optimum catalyst composition (i.e., mixing ratio) varies depending on the operating regime. The FER catalyst weight in the mixture can be as low as 5 wt.%, but the ideal composition de-pends on the internal mass transfer limitation and its relationship with the operating regime (i.e., weight hourly space velocity, temperature). A catalyst composition of 80 wt.% CZZ and 20 wt.% FER was suggested; this composition can provide high CO 2 conversion and DME production rates at a wide range of temperatures and flow rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimization of direct air capture processes using reactive transport models of adsorption-desorption cycles

In this study, we develop and implement a reactive transport model in COMSOL Multiphysics® to address the challenges of direct air carbon capture. The model is validated against experimental data and used to simulate the cyclic steady state of the adsorption-desorption process. The optimization of this model is achieved through advanced trust-region methods integrated with Gaussian Processes. Key decision variables, including adsorption and desorption times, desorption temperature and pressure, input velocity, bed porosity, column length, and radius were optimized to minimize the capture cost. After optimization, a sensitivity analysis revealed the complex interplay between the decision variables and their effect on the specific energy and cost of removing the CO 2 . We optimized the capture cost while taking into account the trade-off between energy consumption and productivity. The resulting minimum capture cost was determined to be 265.2 $/t-CO 2 , which aligns with expected values reported in the literature. Numerical results suggest the effectiveness of the optimization strategies applied, and underscore the importance of simultaneous decision variable selection in improving the performance in direct air capture processes. We also extend the modeling approach to a 2D axisymmetric model to better visualize CO₂ uptake and temperature profiles, revealing significant radial gradients during the regeneration step. As a main drawback, this enhanced model comes with a computational cost approximately 40 times higher than that of the 1D model.

Adsorption-desorption process↗

Effects of fiber dosage, loading orientation and stress on frequency response of enhanced Carbon Nano-Fiber Aggregates

Nanomaterial additives have been broadly used in cement and concrete-based sensors to measure the stress and strain in the structure. Most of the studies focused on measuring the resistivity of sensors with Direct Current (DC) or relatively low-frequency. Here in this work, the impedance responses of enhanced Carbon Nano-Fiber Aggregates (CNFAs) at different frequencies are rigorously studied and utilized as a tool for real-time structural health monitoring (SHM). The CNFA impedance measurements at different measurement frequencies provide an extra dimension (frequency dimension) to allow the characterization of different phenomena and multimodal measurements simultaneously. First, a comprehensive physical model is established to understand the frequency responses of the CNFA impedance. Second, COMSOL Multiphysics is used to simulate the frequency responses of the CNFA impedance, which provides more insights into the electric field and current density distribution. The impedance of CNFAs decreased with the increased dosage of carbon nanofibers (CNFs). Third, experimental studies are reported in detail. The electrical impedance variation (EZV) of the CNFA in parallel orientation is 151% higher than that of perpendicular orientation at the frequency of 2612 Hz and uniaxial compression stress of 4.65 MPa. The strain-EZV curve obtained from the response spectrum is linear for frequencies ranging from 98 Hz to 463.9 kHz. The strain-EZV curve for 5625 Hz has a linear fit with the gauge factor of 147.78. At 1000 Hz, the CNFA exhibited a repeatable behavior up to 9.35 MPa and detection limit up to the stress of 18.62 MPa.

36 MATERIALS SCIENCE↗

Coupling thermal energy storage with a thermally anisotropic building envelope for building demand-side management across various US climate conditions

Here, the thermally anisotropic building envelope (TABE) is a novel active building envelope that enhances energy efficiency and thermal comfort in buildings by transferring heat and cold between building envelopes and hydronic loops. When coupled with thermal energy storage (TES) units, the TABE + TES enables the storage of both heat and cold energy captured by the TABE roof or exterior walls. This stored energy can be later released by the TABE floor for indoor heating and cooling, benefiting both the grid and the end user. This paper evaluates the merits of TABE + TES for building demand-side management across various US climate conditions, focusing on peak load shaving, annual energy savings, and cost savings under time-of-use (TOU) electric rate schedules. Simulations were conducted by integrating time-of-day–informed, rule-based control strategies in MATLAB, TABE components and TES units in COMSOL Multiphysics, and whole-building energy analysis in EnergyPlus. A case study using the US Department of Energy’s prototype single-family detached house model in Birmingham, Alabama; Los Angeles, California; Oak Ridge, Tennessee; and Denver, Colorado, showed that the TABE + TES system achieved (1) 70 % peak load shaving in Los Angeles and Denver and 20 % in Birmingham and Oak Ridge; (2) significant peak electricity savings of 351–497 kWh, reducing peak energy consumption by 38 %–78 %; and (3) annual heating cost savings of 0.79 $\$$/m2–1.17 $\$$/m 2 and cooling cost savings of 0.60 $\$$/m 2 –1.17 $\$$/m 2 using a normal utility rate or low-TOU rate. The benefits of employing the TABE + TES system are even more significant under high TOU rates.

25 ENERGY STORAGE↗

A distillation column for hydrogen isotope removal from liquid lithium

Recovery of tritium from plasma-facing components in fusion devices will be vital to future full-scale operation. Liquid, low-Z materials have demonstrated many inherent advantages over solid first wall materials. To this end, a thermal treatment method in the form of a distillation column for extraction of hydrogen isotopes from liquid lithium has been designed, developed, and constructed at the Center for Plasma-Material Interactions at the University of Illinois at Urbana-Champaign. Use of induction heating and lithium condensation stages are the two qualities that set this design apart from other thermal treatment systems. Induction heating capabilities were modeled using the COMSOL Multiphysics software, which were validated when commissioning the physical heater module. Proof-of-concept tests were performed in the prototype column, which were undertaken as batch processes to investigate the efficacy with which the column could remove hydrogen gas from lithium-rich and lithium hydride-rich samples. All of the tests reported used lithium hydride as a surrogate for lithium deuteride and lithium tritide. The design process and results from the initial tests will be discussed, along with the envisioned placement of this treatment scheme in a fully-functional lithium loop.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Construction of 3D MHD pressure drop correlation and flow characterization in the contraction region of a fusion blanket manifold

Inlet and outlet manifolds are typical components of liquid metal (LM) blanket designs of a fusion power reactor to be used to distribute the LM flow into breeding channels and collect it at the exit of the blanket. High pressure loss in the magnetohydrodynamic (MHD) flows featuring abrupt geometrical changes is one of the main feasibility issues of such designs. Recently, optimization studies were conducted to construct 3D MHD pressure drop correlations for a LM flow in an electrically insulating manifold with gradual expansion. Here, the 3D computational approach developed in that study is applied to the outlet manifold featuring gradual contraction. A systematic analysis was performed with a total number of 135 flow cases computed with COMSOL Multiphysics for Hartmann numbers 1000 < Ha < 10,000, Reynolds numbers 100 < Re < 12,000, and contraction angles 45° < θ < 75° for a fixed contraction ratio of 4. The effects of Ha, Re and θ on the flow recirculation, development length and the total pressure drop were carefully examined. A linear regression analysis was used to determine the power rule of pressure drop coefficient k related to Ha and Re, demonstrating a good match with the Ludford layer theory. Eventually, a correlation for the 3D MHD pressure drop coefficient was constructed as a function of Ha, Re and θ. Further, the results were compared against the inlet manifold. It was found that the flow in the inlet manifold exhibits larger recirculation zones. In the investigated range of Ha, Re and θ, the pressure drop coefficient k of the LM MHD flow in the gradual contraction is only slightly lower (< 8 %) than that in the gradual expansion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational examination of two-phase microchannel heat transfer correlations with conjugate heat spreading

Microchannel flow boiling is an attractive thermal management strategy for ever-growing volumetric heat dissipation demands associated with electronic systems. Due to difficulties related to measurement at the microscales, the majority of researchers have chosen to study relatively simple situations with uniform heat flux. However, many applications involve local hotspots which give rise to highly non-uniform heat flux and temperature gradients due to heat spreading. This necessitates the consideration of conjugate heat transfer for accurate analysis. The current work is aimed at investigating conjugate heat transfer in a two-phase microchannel array. Experimental data was collected on R134a flow boiling heat transfer for very small hydraulic diameter (<100 µm) silicon microchannels with very high heat fluxes (>1 kW cm -2 ) applied via platinum strip heaters with a footprint size of 1 cm × 1 mm. The collected experimental data was then combined with detailed computational modeling utilizing finite element modeling with COMSOL Multiphysics and MATLAB to examine the applicability of five published heat transfer correlations for use in determining local heat transfer coefficients. A two-phase correlation from Agostini and Bontemps, developed for markedly different test parameters, provided the best computational agreement with an RMS temperature difference from experiment of 3.3 °C and a predicted peak perimeter-averaged heat transfer coefficient of 116 MW m -2 K -1 . Modeling confirms the presence of highly non-uniform local heat flux and correspondingly non-uniform local heat transfer coefficient. Finally, the results of this study make clear the need for better micro-scale two-phase correlations developed to predict local heat transfer coefficients at these small scales and high local heat fluxes.

42 ENGINEERING↗

Cooling performance of an active-passive hybrid composite phase change material (HcPCM) finned heat sink: Constant operating mode

Here, the present study explores a hybrid thermal management technology based on air cooling and hybrid composite phase change material (HcPCM) filled finned heat sink for cooling performance of lower to medium heat flux dissipating electronic devices. Two-dimensional numerical simulations are conducted to study the conjugate heat transfer effects of three types of finned heat sink: air-cooled finned heat sink, HcPCM-cooled finned heat sink, and hybrid (air-HcPCM) cooled finned heat sink. A heat sink with a constant volume faction of plate-fins is designed in all cases and simultaneous effects of hybrid nanoparticles and air are investigated to keep the heat sink temperature at safe operating conditions between 40–60°C. The effect of air is incorporated into the heat sink by applying the convective heat transfer coefficient of h c = 10–100W /m 2 .K which tends to create the natural convection and forced convection heat transfer characteristics. The heat flux is varied from 25–40kW /m 2 in the current study. The hybrid nanoparticles of carbon additives (GO and MWCNTs) are dispersed into the RT-35HC, used as a PCM, with a volume fraction of 0% to 6%. Transient simulations are carried out using COMSOL Multiphysics to solve the governing equations for PCM based conjugate heat transfer model. The results showed that forced convection heat transfer improved the cooling performance of the hybrid heat sink compared to natural convection heat transfer. The addition of nanoparticles further enhanced thermal enhancement and uniform melting distribution of PCM inside the finned heat sink. The h c between 30 to 50W/m 2 .K shows optimized values for forced convection heat transfer operating conditions. The volume fraction of 2% of GO+MWCNTs nanoparticles in recommended or optimum concentration for uniform melting of PCM inside the finned heat sink.

36 MATERIALS SCIENCE↗

Modeling and experimental validation of a W f /W-fabrication by chemical vapor deposition and infiltration

Tungsten (W) has a unique combination of excellent thermal properties, low sputter yield, low hydrogen retention, and acceptable activation. Therefore, W is presently the main candidate for the first wall material in future fusion devices. However, its intrinsic brittleness and its further embrittlement during operation bears the risk of a sudden and catastrophic component failure. As a countermeasure, tungsten fiber-reinforced tungsten (W f /W) with extrinsic toughening is being developed. A possible synthesis route is chemical vapor deposition (CVD) using heated W fabrics as substrate. The challenge is that the growing CVD-W can isolate domains from precursor access leading to strength-reducing pores. To deepen the process understanding and to optimize the CVD parameters, models were developed with COMSOL Multiphysics and validated experimentally. W deposition rate equations as function of the temperature and the partial pressures of the precursors H 2 and WF 6 were experimentally validated in previous work. In the present article, these equations are applied to obtain partial pressures within the CVD reactor. The results are taken as input for transient simulations in the microscale, in which W coatings, growing onto multiple adjacent W fibers, were simulated via mesh deformation and remeshing. The surface-to-surface contact of the W coatings and the corresponding potential pore formation were simulated by implementing sophisticated deposition rate stop conditions. Within the measuring uncertainties of ≃ ±1%, the models are validated successfully by experimental comparison regarding the deposition rate, pore structure, and relative densities ranging from 0.6 to 0.9.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Heat pulse testing at monitoring wells to estimate subsurface fluid velocities in geological CO 2 storage

Monitoring the injected CO 2 during geological CO 2 storage (GCS) is essential to assure containment and identify CO 2 leakage. Here in this work, a new approach is introduced to estimate the evolution of the downhole fluid velocity at a monitoring well and identify CO 2 arrival time using in-well heat pulse/tracer test. The proposed technique involves using a downhole heater to generate a series of heat pulses and measuring their corresponding temperature response. The surface temperature of the downhole heater is controlled by the supplied electrical power and the heat loss by convection to the surroundings. Convective heat transfer is well described using Newton's law of cooling in which the temperature difference between the heater and the surrounding fluids drives the heat transfer, for which the convection heat transfer coefficient (h) controls the magnitude of heat loss. Among various factors that control h, it depends on the type of the flowing fluid and its velocity. Through analyzing the measured temperature at different heat pulses, the changes in h - due to mobilization of the in-situ brine or CO 2 arrival - can be estimated. Consequently, the velocity of the flowing fluid across the heater can be obtained. Since heat transfer by convection is sensitive to the type of the surrounding fluid, intrusion of CO 2 can be detected from the relatively higher surface temperature obtained at CO 2 arrival. Churchill and Bernstein (1977)'s correlation is adopted to estimate the change of fluid velocity in terms of the change in h. To demonstrate the validity of the proposed technique, the results are applied and validated against those of COMSOL Multiphysics simulation tool for single-phase brine (before CO 2 arrival) and single-phase CO 2 (after CO 2 arrival). The observed temperature heating is sensitive to the flowing fluid velocity and fluid type. The temperature signal observed at CO 2 arrival is large and easily detectable using temperature monitoring tool which provides reliable indication for tracking CO 2 arrival at monitoring wells compared with passive temperature monitoring. The results obtained using the proposed technique agree very well with the numerical results obtained from the simulation tool with a maximum estimation error of 7 percent.

02 PETROLEUM↗