Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMSOL simulations”

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 19 records

3D COMSOL Simulations for Thermal Deflection of HFIR Fuel Plate in the "Cheverton-Kelley" Experiments

Three dimensional simulation capabilities are currently being developed at Oak Ridge National Laboratory using COMSOL Multiphysics, a finite element modeling software, to investigate thermal expansion of High Flux Isotope Reactor (HFIR) s low enriched uranium fuel plates. To validate simulations, 3D models have also been developed for the experimental setup used by Cheverton and Kelley in 1968 to investigate the buckling and thermal deflections of HFIR s highly enriched uranium fuel plates. Results for several simulations are presented in this report, and comparisons with the experimental data are provided when data are available. A close agreement between the simulation results and experimental findings demonstrates that the COMSOL simulations are able to capture the thermal expansion physics accurately and that COMSOL could be deployed as a predictive tool for more advanced computations at realistic HFIR conditions to study temperature-induced fuel plate deflection behavior.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

COMSOL Simulations for Steady State Thermal Hydraulics Analyses of ORNL s High Flux Isotope Reactor

Simulation models for steady state thermal hydraulics analyses of Oak Ridge National Laboratory s High Flux Isotope Reactor (HFIR) have been developed using the COMSOL Multiphysics simulation software. A single fuel plate and coolant channel of each type of HFIR fuel element was modeled in three dimensions; coupling to adjacent plates and channels was accounted for by using periodic boundary conditions. The standard k- turbulence model was used in simulating turbulent flow with conjugate heat transfer. The COMSOL models were developed to be fully parameterized to allow assessing impacts of fuel fabrication tolerances and uncertainties related to low enriched uranium (LEU) fuel design and reactor operating parameters. Heat source input for the simulations was obtained from separate Monte Carlo N Particle calculations for the axially non-contoured LEU fuel designs at the beginning of the reactor cycle. Mesh refinement studies have been performed to calibrate the models against the pressure drop measured across the HFIR core.

CFD↗

COMSOL Simulation of Optical Whispering Gallery Mode (WGM) Resonator

We simulate whispering-gallery mode coupling in a Si$_3$N$_4$--LiNbO$_3$ waveguide--resonator system using 2D and 3D COMSOL models. Results show that the TM$_{00}$ mode couples far more efficiently than TE$_{00}$ due to weaker confinement and stronger evanescent overlap, with critical coupling reached near a 0.20~µm gap. The study highlights polarization dependence, geometric sensitivity, and birefringence as key factors shaping resonant behavior.

42 ENGINEERING↗

COMSOL Simulation of Optical Whispering Gallery Mode (WGM) Resonator

We simulate whispering-gallery mode coupling in a Si$_3$N$_4$--LiNbO$_3$ waveguide--resonator system using 2D and 3D COMSOL models. Results show that the TM$_{00}$ mode couples far more efficiently than TE$_{00}$ due to weaker confinement and stronger evanescent overlap, with critical coupling reached near a 0.20~µm gap. The study highlights polarization dependence, geometric sensitivity, and birefringence as key factors shaping resonant behavior.

Xiaoyu, Chen [Fermilab]↗

Polarity effects and apparent ion recombination in microionization chambers

Purpose: Microchambers demonstrate anomalous voltage-dependent polarity effects. Existing polarity and ion recombination correction factors do not account for these effects. As a result, many commercial microchamber models do not meet the specification of a reference-class ionization chamber as defined by the American Association of Physicists in Medicine. The purpose of this investigation is to determine the cause of these voltage-dependent polarity effects. Methods: A series of microchamber prototypes were produced to isolate the source of the voltage-dependent polarity effects. Parameters including ionization-chamber collecting-volume size, stem and cable irradiation, chamber assembly, contaminants, high-Z materials, and individual chamber components were investigated. Measurements were performed with electrodes coated with graphite to isolate electrode conductivity. Chamber response was measured as the potential bias of the guard electrode was altered with respect to the collecting electrode, through the integration of additional power supplies. Ionization chamber models were also simulated using COMSOL Multiphysics software to investigate the effect of a potential difference between electrodes on electric field lines and collecting volume definition. Results: Investigations with microchamber prototypes demonstrated that the significant source of the voltage-dependent polarity effects was a potential difference between the guard and collecting electrodes of the chambers. The voltage-dependent polarity effects for each prototype were primarily isolated to either the guard or collecting electrode. Polarity effects were reduced by coating the isolated electrode with a conductive layer of graphite. Polarity effects were increased by introducing a potential difference between the electrodes. COMSOL simulations further demonstrated that for a given potential difference between electrodes, the collecting volume of the chamber changed as the applied voltage was altered, producing voltage-dependent polarity effects in the chamber response. Ionization chamber measurements and COMSOL simulations demonstrated an inverse relationship between the chamber collecting volume size and the severity of voltage-dependent polarity effects on chamber response. The effect of a given potential difference on chamber polarity effects was roughly ten times greater for microchambers as compared to Farmer-type chambers. Stem and cable irradiations, chamber assembly, contaminants, and high-Z materials were not found to be a significant source of the voltage-dependent polarity effects. Conclusions: A potential difference between the guard and collecting electrodes was found to be the primary source of the voltage-dependent polarity effects demonstrated by microchambers. For a given potential difference between electrodes, the relative change in the collecting volume is smaller for larger-volume chambers, illustrating why these polarity effects are not seen in larger-volume chambers with similar guard and collecting electrode designs. Thus, for small-volume chambers, it is necessary to reduce the potential difference between the guard and collecting electrodes in order to reduce polarity effects for reference dosimetry measurements.

60 APPLIED LIFE SCIENCES↗

Comsol multiphysics simulation of the 1958 Y-12 criticality accident - 14654

This paper presents the simulation of the 1958 Y-12 criticality excursion using the commercial multiphysics code COMSOL Multiphysics. ICNC 2011 reported the initial effort to simulate the Y-12 accident power history using a fully coupled neutronic and computational fluid dynamic (CFD) approach with the Imperial College London FETCH code. That simulation included the initial spike and a subsequent power history 30 seconds afterwards. A time-dependent lumped parameter excursion model was later developed at Imperial in November 2014 and taken to 20 min. The COMSOL simulation is viewed as a hybrid treatment of these efforts. The use of COMSOL as a platform for criticality excursion simulation had been previously demonstrated with SILENE step reactivity insertion benchmarks. The Y-12 accident represented a unique time dependent reactivity insertion mechanism by dilution of the fissile concentration over time that facilitated the excursion termination. The COMSOL model includes a moving mesh vertical boundary, compartment radiolytic gas model, coupled turbulent fluid flow and conjugate heat transfer. The results are compared with both FETCH and lumped parameter simulations. The COMSOL power history results are in agreement with the reported and reconstructed initial spike fissions, power and total fissions over time. The simulated power trace shows explicit subsequent spikes approximately 10-20 seconds apart from the initial spike, over the first three minutes of the excursion. This result qualitatively confirms the radiation dose rate trace from a far field monitor located well beyond the affected facility. Additionally, the COMSOL model shows the onset of boiling after 3 minutes, as postulated in the literature. The combination of boiling and fissile solution dilution contributed to the shutdown mechanism. (authors)

42 ENGINEERING↗

Acoustic Modeling of Novel Over-The-Rotor Acoustic Liner Concepts Using COMSOL Multiphysics Simulation Software

This paper presents simulation predictions of novel over-the-rotor acoustic liner configurations compared with experimental data. This is a follow on to the investigation on the effect of grazing flow on grooved over-the-rotor casing treatments. The purpose of this modeling investigation stems from the determination that traditional liner modeling techniques are inadequate for nontraditional internal liner cavity geometries. A series of these expansion chamber type liner design concepts were previously developed by Glenn Research Center and tested using the Langley Research Center Normal Incidence Tube. The measured normalized impedance and absorption coefficient is compared with that obtained from the COMSOL Multiphysics® simulation. Additionally, using a parametric case study, another series of simulations were conducted to determine what additional benefits, if any, this genre of liner has to offer the acoustics community. A total of 53 permutations of this concept was investigated. The results indicate that with longer internal fins, the resonance is sifted to lower frequencies.

Acoustics↗

Energy performance analysis of smart wall system with switchable insulation and thermal storage capacity

The ever-increasing global energy demand and the issues caused by population growth and unsustainable energy resource usages have several environmental and economic impacts. The on-demand capability of dynamic wall systems with switchable insulation systems can contribute toward energy efficiency and reduced electric cost using "building-as-a-battery." In this paper, the performance of an exterior envelope system that employs a switchable insulation system is investigated. A COMSOL simulation was used to study the envelope performance under three switchable insulation locations in the wall system. The on and off switching cycle included insulating or conducting the exterior, interior, or both sides of the mass wall system. To validate the simulation work, an experimental test was conducted in climate chamber on a 4 × 8 in. wall system with identical wall components used in the COMSOL simulation. Results show a good agreement between the experiment and simulation results. The wall system with switchable insulation placed on the exterior side provided the highest interior inward heat flux when compared with the interior or split insulation. Also, an exterior switchable insulation system in addition to a thermal mass maintained a positive heat flow for more hours when the outdoor temperature was lower than the interior temperature.

Iffa, Emishaw↗

An improved back projection algorithm of ultrasound tomography

Binary logic back projection algorithm is improved in this work for the development of fast ultrasound tomography system with a better effect of image reconstruction. The new algorithm is characterized by an extra logical value ‘2’ and dual-threshold processing of collected raw data. To compare with the original algorithm, a numerical simulation was conducted by the verification of COMSOL simulations formerly, and then a set of ultrasonic tomography system is established to perform the experiments of one, two and three cylindrical objects. The object images are reconstructed through the inversion of signals matrix acquired by the transducer array after a preconditioning, while the corresponding spatial imaging errors can obviously indicate that the improved back projection method can achieve modified inversion effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A 1-D Model of the 4 Bed Molecular Sieve of the Carbon Dioxide Removal Assembly

Developments to improve system efficiency and reliability for water and carbon dioxide separation systems on crewed vehicles combine sub-scale systems testing and multi-physics simulations. This paper describes the development of COMSOL simulations in support of the Life Support Systems (LSS) project within NASA's Advanced Exploration Systems (AES) program. Specifically, we model the 4 Bed Molecular Sieve (4BMS) of the Carbon Dioxide Removal Assembly (CDRA) operating on the International Space Station (ISS).

Coker, Robert↗

RF cavity design and qualification for proton accelerator

Alvarez type Drift Tube Linac (DTL) is used for acceleration of proton beam in low energy section of beta ranging from 0.04 to 0.40. DTL is cylindrical RF cavity resonating in TM010 mode at 352.21 MHz frequency. It consists of array of drift tubes arranged ensuring that DTL centre and Drift Tube centre are concentric. The Drift Tubes also houses Permanent Magnet Quadrupole for transverse focusing of proton beam. A twelve cell prototype of DTL section is designed, developed and fabricated at Bhabha Atomic Research Centre, Trombay. Complete DTL accelerator consists of eight such DTL sections. High frequency microwave simulations are carried out in SOPRANO, vector fields and COMSOL simulation software. This prototype DTL is 1640.56 mm long cavity with 520 mm ID, 600 mm OD and consists of eleven Drift Tubes, two RF end flanges, three slug tuners, six post couplers, three RF field monitors, one RF waveguide coupler, two DN100 vacuum flanges and DTL tank platform with alignment features. Girder based Drift tube mounting arrangement utilizing uncompressing energy of disc springs for optimum combo RF-vacuum seal compression is worked out and implemented. This paper discusses design of this RF vacuum cavity operating at high accelerating field gradient in ultra-high vacuum. Detailed vacuum design and results of RF and vacuum qualifications are discussed. Results on mechanical accuracy achieved on scaled pre-prototype are also presented. Paper summarizes the engineering developments carried out for this RF cavity and brings out the future activities proposed in indigenous development of high gradient RF cavities for ion accelerators. (author)

43 PARTICLE ACCELERATORS↗

Modeling of Charge Amplification in LArTPCs

Single-phase liquid argon time projection chambers (LArTPCs) are an important part of Fermilab's neutrino experimental program. Current LArTPCs are able to detect interactions producing energy deposits above the MeV scale, and focus on neutrinos at GeV energies. The Liquid Argon Charge Amplification Device, or LArCADe, program aims to employ modified anode geometries in order to obtain electron avalanches that would help reduce detection thresholds, opening the door to new physics opportunities at the keV scale (e.g. CEvNS and dark-matter WIMP searches). To model and quantify this phenomenon, a simulation of electron transport in LAr was built, able to simulate excitations and ionizations in non-uniform fields. Using COMSOL simulations of the anode structures, the built simulation produced quantitative and realistic estimates of LAr charge amplification near different tip geometries. These results will inform the next generation of tips produced for LArCADe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

3D multiphysics simulation and analysis of a low temperature liquid metal magnetohydrodynamic power generator prototype

Magnetohydrodynamic (MHD) power generation is an interdisciplinary process, which was studied heavily in the late 1960’s but interest soon subsided. However, the current socioeconomic climate demands more diversity in renewable energy sources. MHD generators using liquid metal as the working fluid represent a viable option for this augmentation. Additionally, the steady increase in accessible computational capabilities provides a powerful new tool for interdisciplinary modeling and multiphysics simulations. These factors suggest that this is an opportune time to reevaluate MHD power generation. This paper provides mathematical background as well as a detailed multiphysics model of a low temperature liquid metal MHD power generator built in COMSOL. Simulation details are provided in addition to a comparison to experimental results. Finally, the prototype studied in this paper is the first proof of concept of a permanent magnet excited liquid metal MHD generator using Gallium as its working fluid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Finite Element Simulation of the Acoustic Pressure Inside a Beverage Container for Non-Thermal, Ultrasound-based Pasteurization

The purpose of this effort is to investigate whether large acoustic pressure waves can be transmitted inside beverage containers to enable pasteurization. Acoustic waves are known to induce large nonlinear compressive forces and shock waves in fluids, suggesting that compression waves may be capable of damaging bacteria inside beverage containers without appreciably increasingly the temperature or altering the freshness and flavor of the beverage contents. Although a combined process such as thermosonication (e.g., sonication with heating) is likely more efficient, it is instructive to compute the acoustic pressure field distribution inside the beverage container. The COMSOL simulations used two and three-dimensional models of beverage containers placed in a water bath to compute the acoustic pressure field. A limitation of these COMSOL models is that they cannot determine the bacterial lysis efficiency, rather the models provide an indirect metric of bacterial lysis based on the magnitude of the pressure field and its distribution.

42 ENGINEERING↗

Investigating the Impact of Thickness, Calendering and Channel Structures of Printed Electrodes on the Energy Density of LIBs - 3D Simulation and Validation

Current lithium ion batteries (LIBs) are expensive and bulky, limited by relatively low charging rates. To increase the rate of charging and reduce weight, thin electrodes with high energy density are required. The increase in energy density can be achieved by several techniques including boosting electrolyte transport, high loading/utilization of active material, employing high conductive electrolytes and electrodes with advanced architectures, and increasing cell temperature. In this paper, a 3D physics-based electrochemical model of LIBs is developed in COMSOL simulation software for different thickness, calendering steps as well as channel structures (conical, cylindrical) to optimize the electrode design and in turn maximize volumetric energy density. The simulation results demonstrated that calendering the electrodes with high initial porosity increases the volumetric energy density of the cell. In addition, cylindrical channel structures with relatively lower edge-to-edge distance also results in increased volumetric energy density. The simulation results of the 3D model was validated by comparing it with experimental results.

improving volumetric energy density↗

Thermal Hydraulic Characteristics of Fuel Defects in Plate Type Nuclear Research Reactors

Turbulent flow coupled with heat transfer is investigated for a High Flux Isotope Reactor (HFIR) fuel plate. The Reynolds Averaged Navier-Stokes Models are used for fluid dynamics and the transfer of heat from a thermal nuclear fuel plate using the Multi-physics code COMSOL. Simulation outcomes are compared with experimental data from the Advanced Neutron Source Reactor Thermal Hydraulic Test Loop. The computational results for the High Flux Isotope Reactor core system provide a more physically accurate simulation of this system by modeling the turbulent flow field in conjunction with the diffusion of thermal energy within the solid and fluid phases of the model domain. Recommendations are made regarding Nusselt number correlations and material properties for future thermal hydraulic modeling efforts

COMSOL LEU HFIR↗

On the Inductance of a Single-Channel Spark-Gap Switch for Pulsed Power Applications

In this paper, we derive an alternative expression for calculating the inductance of a single-channel spark-gap switch. We believe this new inductance expression provides greater clarity compared to other approaches. Interestingly, both the Rosa formula and the coaxial inductance formula, which are commonly used in such calculations, emerge as special cases of our new expression. To validate the proposed formula, we conducted COMSOL simulations that show excellent agreement with our theoretical predictions.

COMSOL↗