Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiphysics simulation”

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 415 records · Page 23

High-order, linearly stable, partitioned solvers for general multiphysics problems based on implicit–explicit Runge–Kutta schemes

This work introduces a general framework for constructing high-order, linearly stable, partitioned solvers for multiphysics problems from a monolithic implicit–explicit Runge–Kutta (IMEX-RK) discretization of the semi-discrete equations. The generic multiphysics problem is modeled as a system of n systems of partial differential equations where the i th subsystem is coupled to the other subsystems through a coupling term that can depend on the state of all the other subsystems. This coupled system of partial differential equations reduces to a coupled system of ordinary differential equations via the method of lines where an appropriate spatial discretization is applied to each subsystem. The coupled system of ordinary differential equations is taken as a monolithic system and discretized using an IMEX-RK discretization with a specific implicit–explicit decomposition that introduces the concept of a predictor for the coupling term. We propose four coupling predictors that enable the monolithic system to be solved in a partitioned manner, i.e., subsystem-by-subsystem, and preserve the IMEX-RK structure and therefore the design order of accuracy of the monolithic scheme. The four partitioned solvers that result from these predictors are high-order accurate, allow for maximum re-use of existing single-physics software, and two of the four solvers allow the subsystems to be solved in parallel at a given stage and time step. We also analyze the stability of a coupled, linear model problem with a specific coupling structure and show that one of the partitioned solvers achieves unconditional linear stability for this problem, while the others are unconditionally stable only for certain values of the coupling strength. We demonstrate the performance of the proposed partitioned solvers on several classes of multiphysics problems including a simple linear system of ODEs, advection–diffusion–reaction systems, fluid–structure interaction problems, and particle-laden flows, where we verify the design order of the IMEX schemes and study various stability properties.

97 MATHEMATICS AND COMPUTING↗

Outdoor performance of a tandem InGaP/Si photovoltaic luminescent solar concentrator

Here, we report the design, fabrication and outdoor characterization of a tandem luminescent solar concentrator/Si multi-junction photovoltaic module. Our tandem LSC/Si device consists of an InGaP LSC functioning as a top cell and a passivated contact Si bottom cell. The LSC comprises of an InGaP microcell array coupled to a polymer waveguide, loaded with CdSe/CdS core-shell quantum dot luminophores. The light trapping efficiency of the LSC waveguide is enhanced by encapsulation with photoluminescence trapping mirrors consisting of dielectric multilayer thin films. We demonstrate the performance of the LSC/Si device through a series of outdoor tests under various irradiance conditions conducted at the National Renewable Energy Laboratory. We report the first outdoor testing data of an LSC/Si tandem module, displaying maintained performance across varied diffusivity conditions for the LSC component. Finally, we model the tandem module performance using a ray optic simulation-based multiphysics model and forecast a pathway for high efficiency tandem LSC/Si module performance.

14 SOLAR ENERGY↗

Ripple Effects: Bed Form Morphodynamics Cascading Into Hyporheic Zone Biogeochemistry

The water quality and ecosystem health of river corridors depend on the biogeochemical processes occurring in the hyporheic zones (HZs) of the beds and banks of rivers. HZs in riverbeds often form because of bed forms. Despite widespread and persistent variation in river flow, how the discharge- and grain size-dependent geometry of bed forms and how bed form migration collectively and systematically affects hyporheic exchange flux, solute transport, and biogeochemical reaction rates are unknown. As such, we investigated these linked processes through morphodynamically consistent multiphysics numerical simulation experiments. Several realistic ripple geometries based on bed form stability criteria using mean river flow velocity and median sediment grain size were designed. Ripple migration rates were estimated based primarily on the river velocity. The ripple geometries and migration rates were used to drive hyporheic flow and reactive transport models which quantified HZ nitrogen transformation. Results from fixed bed form simulations were compared with matching migrating bed form scenarios. We found that the turnover exchange due to ripple migration has a large impact on reactant supply and reaction rates. The nitrate removal efficiency increased asymptotically with Damköhler number for both mobile and immobile ripples, but the immobile ripple always had a higher nitrate removal efficiency. Since moving ripples remove less nitrogen, and may even be net nitrifying at times, consideration for bed form morphodynamics may therefore lead to reduction of model-based estimates of denitrification. The connection between nitrate removal efficiency and Damköhler number can be integrated into frameworks for quantifying transient, network-scale, HZ nitrate dynamics.

54 ENVIRONMENTAL SCIENCES↗

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 ↗

Vibration energy harvesting in a small channel fluid flow using piezoelectric transducer

This work is aimed at developing a way to harvest energy from a fluid stream with the application of piezoelectric transducers in a small channel. In this COMSOL Multiphysics based simulation study, it is attempted to harvest energy from the abundant renewable source of energy available in the form of kinetic energy of naturally occurring flow of fluids. The strategy involves harnessing energy from a fluid-actuator through generation of couples, eddies and vortices, resulting from the stagnation and separation of flow around a semi-circular bluff-body attached to a cantilever beam containing a piezoceramic layer. Fluctuation of fluidic pressure impulse on the beam due to vortex shedding and varying lift forces causes the flexible cantilever beam to oscillate in the direction normal to the fluid flow in a periodic manner. The periodic application and release of a mechanical strain upon the beam effected a generation of electric potential within the piezoelectric layer, thus enabling extraction of electrical energy from the kinetic energy of the fluid. The piezoelectric material properties and transducer design are kept unchanged throughout the study, whereas the configuration is tested with different fluids and varying flow characteristics. The size and geometry of the obstructing entity are systematically varied to closely inspect the output from different iterations and for finding the optimum design parameters. The intermittent changes in the generated forces and subsequent variation in the strain on the beam are also monitored to find definitive relationship with the electrical energy output.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development and Verification of an Improved Wake-Added Turbulence Model in FAST.Farm

We introduce a generalized wake-added turbulence (WAT) model in the multiphysics, multiturbine simulation tool FAST.Farm. The WAT model introduces additional small-scale turbulence that represents the breakdown of vortical structures and shear layers in the wake. The article describes the development, implementation, calibration, and verification of the model. The novelties of the model include support for wake asymmetry, buildup of WAT across the wind farm, and secondary effects of wake-induced turbulence (e.g., wake meandering) driven by smaller-scale turbulence structures that arise from wake breakdown. Large-eddy simulations were run to support the calibration of the WAT parameters and verification of the model. Previous studies hypothesized that the lack of WAT modeling was the source of underprediction of fatigue loads, in particular for cases at low turbulence intensities and/or stable atmospheric boundary layers. This study confirms that the newly implemented WAT model enhances the loads predictions in these cases.

17 WIND ENERGY↗

Capacitor Design for Self-Resonant Coils for Long-Distance Wireless Power Transfer System

In this paper, an integrated capacitor design is proposed for higher-order resonant tank topologies for self-resonant coils, such as series, parallel, LCC, LLC, etc. The capacitor is one of the large, lossy, and thermally vulnerable components of a high-frequency resonant tank, and designing a high-voltage, thermally stable resonant capacitor can be highly challenging. Designing the extremely high-voltage capacitor as an integral part of the coil reduces the size and complexity of the coil assembly. This paper proposes a low-loss PCB-based high-voltage capacitor design to achieve that target, which can be implemented as an integral part of the coil. The proposed capacitor designs are simulated using Multiphysics FEA and tested experimentally. A 23 kV, 133 pF capacitor prototype was built and tested as part of a 1 kW long-distance wireless charging system. The test results verify the capacitor’s voltage, current, and thermal resiliency performance.

Mohammad, Mostak [ORNL] (ORCID:0000000256388783)↗

Community Organizations: Changing the Culture in Which Research Software Is Developed and Sustained

Software is the key crosscutting technology that enables advances in mathematics, computer science, and domain-specific science and engineering to achieve robust simulations and analysis for science, engineering, and other research fields. However, software itself has not traditionally received focused attention from research communities; rather, software has evolved organically and inconsistently, with its development largely as by-products of other initiatives. Moreover, challenges in scientific software are expanding due to disruptive changes in computer hardware, increasing scale and complexity of data, and demands for more complex simulations involving multiphysics, multiscale modeling and outer-loop analysis. In recent years, community members have established a range of grass-roots organizations and projects to address these growing technical and social challenges in software productivity, quality, reproducibility, and sustainability. This article provides an overview of such groups and discusses opportunities to leverage their synergistic activities while nurturing work toward emerging software ecosystems.

97 MATHEMATICS AND COMPUTING↗

Direct Measurements of Ionic Transport Behavior of Dual-Layer Porous Electrodes

To improve power and cycling performance of lithium-ion batteries, dual-layer or porosity-gradient electrodes have been proposed. By engineering a higher porosity close to the separator, the intention is to improve ion transport where it is most impactful. In this research, MacMullin numbers of two dual-layer anode samples are tested using an impedance measurement technique developed previously. To characterize the microstructure of each layer independently, we developed an improved transmission-line model that accounts for each layer's properties.Virtual experiments using COMSOL Multiphysics to simulate impedance measurements are used to examine and improve the accuracy of the numerical inversion procedure. The results for the two dual-layer anodes studied show that MacMullin numbers follow expected trends, though the anodes are quite different from each other.

Dual-Layer Electrodes↗

Kelvin

Kelvin is a multiphysics sintering simulation software that uses the material point method to simulate the thermomechanical behavior of metal parts.

Billings, JayJ↗

The Alamo multiphysics solver for phase field simulations with strong-form mechanics and block structured adaptive mesh refinement

Alamo is a high-performance scientific code that uses block-structured adaptive mesh refinement to solve such problems as: the ignition and burn of solid rocket propellant, plasticity, damage and fracture in materials undergoing loading, and the interaction of compressible flow with eroding solid materials. Alamo is powered by AMReX, and provides a set of unique methods, models, and algorithms that enable it to solve solid-mechanics problems (coupled to other physical behavior such as fluid flow or thermal diffusion) using the power of block-structured adaptive mesh refinement.

36 MATERIALS SCIENCE↗

Towards a supported common NEAMS software stack

The NEAMS IPSC's are developing multidimensional, multiphysics, multiscale simulation codes based on first principles that will be capable of predicting all aspects of current and future nuclear reactor systems. These new breeds of simulation codes will include rigorous verification, validation and uncertainty quantification checks to quantify the accuracy and quality of the simulation results. The resulting NEAMS IPSC simulation codes will be an invaluable tool in designing the next generation of Nuclear Reactors and also contribute to a more speedy process in the acquisition of licenses from the NRC for new Reactor designs. Due to the high resolution of the models, the complexity of the physics and the added computational resources to quantify the accuracy/quality of the results, the NEAMS IPSC codes will require large HPC resources to carry out the production simulation runs.

97 MATHEMATICS AND COMPUTING↗

Enabling Predictive Simulation and UQ of Complex Multiphysics PDE Systems by the Development of Goal-Oriented Variational Sensitivity Analysis and A Posteriori Error Estimation Methods

It was demonstrated that a posteriori analyses in general and in particular one that uses adjoint methods can accurately and efficiently compute numerical error estimates and sensitivity for critical Quantities of Interest (QoIs) that depend on a large number of parameters. Activities include: analysis and implementation of several time integration techniques for solving system of ODEs as typically obtained from spatial discretization of PDE systems; multirate integration methods for ordinary differential equations; formulation and analysis of an iterative multi-discretization Galerkin finite element method for multi-scale reaction-diffusion equations; investigation of an inexpensive postprocessing technique to estimate the error of finite element solution of the second-order quasi-linear elliptic problems measured in some global metrics; investigation of an application of the residual-based a posteriori error estimates to symmetric interior penalty discontinuous Galerkin method for solving a class of second order quasi-linear elliptic problems; a posteriori analysis of explicit time integrations for system of linear ordinary differential equations; derivation of accurate a posteriori goal oriented error estimates for a user-defined quantity of interest for two classes of first and second order IMEX schemes for advection-diffusion-reaction problems; Postprocessing finite element solution; and A Bayesian Framework for Uncertain Quantification of Porous Media Flows.

97 MATHEMATICS AND COMPUTING↗

Performance Improvements for the Griffin Transport Solvers

Griffin is a Multiphysics Object-Oriented Simulation Environment based reactor multiphysics analysis application jointly developed by Idaho National Laboratory and Argonne National Laboratory. Griffin includes a variety of deterministic radiation transport solvers for fixed source, k-eigenvalue, adjoint, and subcritical multiplication, as well as transient solvers for point-kinetics, improved quasi-static, and spatial dynamics. A code assessment performed in FY-20 identified two significant issues with the transport solvers in Griffin: first, the primary heterogeneous SN (discrete ordinates) transport solver based on continuous finite element methods required significant mesh refinement and higher memory usage compared to solvers based on the method of characteristic for equivalent accuracy. Second, the homogeneous PN (spherical harmonics expansion) transport solver did not adequately support polynomial refinement, which is a feature usually required for problems with spatial homogenization and pronounced streaming, typical in fast or gas-cooled reactor systems. To address the first issue, the development effort focused on the more promising discontinuous finite element method (DFEM)-based SN transport solver in Griffin. The addition of an asynchronous parallel transport sweeper and coarse mesh finite difference (CMFD) acceleration have rendered a superior heterogeneous SN transport capability for multiphysics problems that requires far less computing resources in terms of both CPU time and memory usage. This is demonstrated with typical thermal- and fast-spectrum reactor benchmark problems, including 2D Transient Reactor Test, 3D Advanced Burner Test Reactor (ABTR), and 2D and 3D Empire microreactor. For the second issue, the development effort focused on a new transport solver based on the hybrid finite element PN method (HFEM-PN), equivalent to the variational nodal method, as well as a new diffusion solver based on HFEM-Diffusion. This solver is intended for homogenized domains with multiphysics coupling (i.e., supports mesh displacement, seamless temperature feedback, etc.). Initial calculations with the HFEM-Diffusion implementation show very good parallel efficiency for the residual evaluations with the 2D ABTR benchmark. A future development effort will be centered on further improvements to the CMFD, HFEM-PN, and DFEM diffusion solvers to ensure Griffin meets performance and software quality assurance requirements for advanced reactor design and analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Updated SAM Model for the Molten Salt Reactor Experiment (MSRE)

The development of reference standard problems based on prototypical reactor designs is of particular importance to verify the adequacy of computer codes and evaluation models for specific reactor types. To support the multiphysics coupled simulation of molten-salt-fueled reactor (MSR) using SAM and Griffin computer codes for safety and licensing analysis, much efforts have been put into enhancing code capabilities and developing reference models for the MSR primary loop in SAM. In this work, a previously developed Molten Salt Reactor Experiment (MSRE) primary loop model was updated to include a two-dimensional (2-D) core region and external core components in one-dimension (1-D) or zero-dimension (0-D). To ensure accurate feedback calculation in multi-physics simulations, the delayed neutron precursor tracking model and solid graphite model were added in the SAM model. In addition, the 2-D and 1-D domains are tightly coupled using the recently developed single-solve approach in SAM. The updated model has been tested under both steady-state and transient scenarios to demonstrate its potential for the multi-physics simulation of MSRE with coupled SAM and Griffin.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗