Reprint of "CFD simulations of circulating fluidized bed risers, part II, evaluation of differences between 2D and 3D simulations
Explore the source record for details and available documents.
Circulating fluidized bed↗
SEARCH · Search NASA
Results for “3D 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.
Correction of Cloud Optical Thickness Retrievals from Nadir Reflectances in the Presence of 3D Radiative Effects. Part I: Concept and Tests on 3D RT Simulations
3D effects cause substantial underestimation of cloud optical thickness (COT) in airborne and satellite retrievals based on 1D radiative transfer computations (such as in the case of widely used bispectral technique). For a single-layer isolated cloud we propose a simple linear correction of the retrieved COT with the renormalization factor dependent on the cloud’s aspect ratio (the ratio between vertical and horizontal dimensions of the cloud). This is an empirical assumption which we successfully test using synthetic 3D RT data. We introduce a heuristic “block model” of 3D radiative effects and show that the functional form of the renormalization factor is consistent with the process of radiation escape from cloud sides in an essentially 3D geometry. We also extend the block model to the case of single-layer broken cloud field with radiative interaction between the neighboring clouds. In this case the renormalization factor depends also on the distance between clouds.
cloud optical thickness↗
3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses
Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.
42 ENGINEERING↗
Towards a self-consistent model of the convective core boundary in upper main sequence stars: I. 2.5D and 3D simulations
There is strong observational evidence that the convective cores of intermediate-mass and massive main sequence stars are substantially larger than those predicted by standard stellar-evolution models. However, it is unclear what physical processes cause this phenomenon or how to predict the extent and stratification of stellar convective boundary layers. Convective penetration is a thermal-timescale process that is likely to be particularly relevant during the slow evolution on the main sequence. We use our low-Mach-number S EVEN -L EAGUE H YDRO code to study this process in 2.5D and 3D geometries. Starting with a chemically homogeneous model of a 15 M ⊙ zero-age main sequence star, we construct a series of simulations with the luminosity increased and opacity decreased by the same factor, ranging from 10 3 to 10 6 . After reaching thermal equilibrium, all of our models show a clear penetration layer; its thickness becomes statistically constant in time and it is shown to converge upon grid refinement. The penetration layer becomes nearly adiabatic with a steep transition to a radiative stratification in simulations at the lower end of our luminosity range. This structure corresponds to the adiabatic ‘step overshoot’ model often employed in stellar-evolution calculations. The simulations with the highest and lowest luminosity differ by less than a factor of two in the penetration distance. The high computational cost of 3D simulations makes our current 3D data set rather sparse. Depending on how we extrapolate the 3D data to the actual luminosity of the initial stellar model, we obtain penetration distances ranging from 0.09 to 0.44 pressure scale heights, which is broadly compatible with observations.
79 ASTRONOMY AND ASTROPHYSICS↗
The dark side of photovoltaic — 3D simulation of glare assessing risk and discomfort
Photovoltaic (PV) systems form an important force in the implementation of renewable energies, but as we all know, the force has always its dark side. Besides efficiency considerations and discussions about architectures of power distribution networks, the increasing numbers of installations of PV systems for implementing renewable energies have secondary effects. PV systems can generate glare due to optical reflections and hence might be a serious concern. On the one hand, glare could affect safety, e.g. regarding traffic. On the other hand, glare is a constant source of discomfort in vicinities of PV systems. Hence, assessment of glare is decisive for the success of renewable energies near municipalities and traffic zones for the success of solar power. Several courts decided on the change of PV systems and even on their de-installation because of glare effects. Thus, location-based assessments are required to limit potential reflections and to avoid risks for public infrastructure or discomfort of residents. The question arises on how to calculate reflections accurately according to the environment's topography. Our approach is founded in a 3D-based simulation methodology to calculate and visualize reflections based on the geometry of the environment of PV systems. This computational model is implemented by an interactive tool for simulation and visualization. Hence, project planners receive flexible assistance for adjusting the parameters of solar panels amid the planning process and in particular before the installation of a PV system. - Highlights: • Solar panels cause glare that impacts neighborhoods and traffic infrastructures. • Glare might cause disability and discomfort. • 3D environment for the calculation of glare • Interactive tool to simulate and visualize reflections • Impact assessment of solar power plant farms.
54 ENVIRONMENTAL SCIENCES↗
3D simulations of vertical displacement events in tokamaks: A benchmark of M3D-C 1 , NIMROD, and JOREK
In recent years, the nonlinear 3D magnetohydrodynamic codes JOREK, M3D-C1, and NIMROD developed the capability of modeling realistic 3D vertical displacement events (VDEs) including resistive walls. In this paper, a comprehensive 3D VDE benchmark is presented between these state-of-the-art codes. The simulated case is based on an experimental NSTX plasma but with a simplified rectangular wall. There are differences between the physics models and numerical methods, and the VDE evolution leads to sensitivities on the initial conditions that cannot be avoided as can be done in edge localized modes (ELM) and sawtooth simulations (due to the non-cyclical nature of VDEs). Nonetheless, the comparison serves to quantify the level of agreement in the relevant quantities used to characterize disruptions, such as the 3D wall forces and energy decay. The results bring confidence regarding the use of the mentioned codes for disruption studies, and they distinguish aspects that are specific to the models used (e.g., reduced vs full MHD models). Here, the simulations show important 3D fea- tures for a NSTX plasma, such as the self-consistent evolution of the halo current and the origin of the wall forces. In contrast to other reduced MHD models based on an ordering in the aspect ratio, the ansatz-based JOREK reduced MHD model allows capturing many aspects of the 3D dynamics even in the spherical tokamak limit considered here.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Formation and evolution of protostellar accretion discs – II. From 3D simulation to a simple semi-analytic model of Class 0/I discs
Here we use a 3D radiative non-ideal magnetohydrodynamic simulation to investigate the formation and evolution of a young protostellar disc from a magnetized pre-stellar core. The simulation covers the first ${\sim }10\, {\rm kyr}$ after protostar formation and shows a massive, weakly magnetized disc with radius that initially grows and then saturates at ${\sim }30\, {\rm au}$. The disc is gravitationally unstable with prominent large-amplitude spiral arms. We use our simulation results and a series of physical arguments to construct a predictive and quantitative physical picture of Class 0/I protostellar disc evolution from several aspects, including (i) the angular-momentum redistribution in the disc, self-regulated by gravitational instability to make most of the disc marginally unstable; (ii) the thermal profile of the disc, well-approximated by a balance between radiative cooling and accretion heating; and (iii) the magnetic-field strength and magnetic-braking rate inside the disc, regulated by non-ideal magnetic diffusion. Using these physical insights, we build a simple 1D semi-analytic model of disc evolution. We show that this 1D model, when coupled to a computationally inexpensive simulation for the evolution of the surrounding pseudo-disc, can be used reliably to predict disc evolution in the Class 0/I phase. The predicted long-term evolution of disc size, which saturates at ${\sim }30\, {\rm au}$ and eventually shrinks, is consistent with a recent observational survey of Class 0/I discs. Such hierarchical modelling of disc evolution circumvents the computational difficulty of tracing disc evolution through Class 0/I phase with direct, numerically converged simulations.
79 ASTRONOMY AND ASTROPHYSICS↗
3D simulations of inertial confinement fusion implosions part 2: systematic flow anomalies and impact of low modes on performances in OMEGA experiments
We present the first 3D radiation-hydrodynamic simulations of directly driven inertial confinement fusion implosions with an inline package for polarized crossed-beam energy transfer, which were used to assess the impact of the current distributed polarization rotators (DPRs) on OMEGA as well as other known sources of asymmetry. Applied to OMEGA implosions, the simulations predict bang times with no need for ad hoc multipliers, as well as yields—if you separately account for the impacts of imprint and fuel age. The magnitude of the flow is well reproduced when the low mode sources are large, whereas the modeling of the stalk is thought to be required to match the flow magnitude in the remaining cases. For the cases explored in more detail, polarized cross-beam energy transfer (CBET)—the only known systematic drive asymmetry, brought the results closest to the measured flow vectors. The remaining discrepancies are shown to be stemming from the limited knowledge of the laser pointing modes. For typical current levels of beam mispointing, power imbalance, target offset, and asymmetry caused by polarized CBET, low modes degrade the yield by more than 40%. The current strategy of attempting to compensate the mode-1 asymmetry with a preimposed target offset recovers only about one-third of the losses caused by the low modes due to the dynamic nature of the multiple asymmetries and the presence of low modes other than l = 1. Furthermore, addressing the root causes of the drive asymmetries is apt to be more beneficial. To that end, one possible solution to the specific issue of polarized CBET (10 µm DPRs) is shown to work well.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗
Preliminary 3D Simulations of Anomalous Thermocouple Coupling to Surface Electric Fields
Anomalous measurements from subsurface thermocouples have been observed during re-entry of the Orion capsule (EFT-1 plugs). These measurements are considered anomalous because the amplitude and frequency of the inferred temperature modulations are too large to be due to plausible heating processes. Consequently, it is believed the anomalies are caused by electrical interference generated by the re-entry plasma. Quantifying these effects requires several steps: modeling the external plasma flows, modeling the plasma-surface interaction, and finally modeling the surface-thermocouple coupling. This paper addresses this last step using 3D finite element simulations of a wire embedded in a weakly conductive medium.
thermocouple anomaly↗
3D simulations of inertial confinement fusion implosions part 1: inline modeling of polarized cross beam energy transfer and subsequent drive anomalies on OMEGA and NIF
Inertial confinement fusion experiments are sensitive to cross-beam energy transfer (CBET), a laser-plasma instability that redistributes laser energy in the coronal plasma through self-generated ion acoustic wave (IAW) gratings. The detailed CBET coupling depends on the polarization state of the crossing wavefields. CBET itself can also scramble the beam polarizations by inducing ellipticity through the IAW grating, and rotating the seed polarization toward that of the pump. We develop a ray-based model that describes the polarized CBET coupling and that is compatible with the framework of 3D inline radiative hydrodynamics simulations. The model is implemented in the ASTER/IFRIIT code and verified against an academic test case and an offline polarized CBET post-processor. It is then applied to the detailed configuration of the distributed polarization rotator system on OMEGA, where results highlight how polarized CBET induces significant low modes in the collisional absorption source term. Finally, the modeling is applied to a simple indirect-drive configuration, comparing CBET calculations with 96 unpolarized or polarized beams with 24 unpolarized quads. It is shown that these cases produce similar power amplification per cone of beams grouped with similar polar angles. However, the 96 beam geometry itself is found to reduce azimuthal variations in quad power after the interaction and favors beams with larger polar angles within the cones, an effect that is amplified by the polarized CBET. Application of the model to inline calculations of OMEGA implosions are presented in a companion paper.
3D inline modeling↗
3D Simulations of the Early Martian Hydrological Cycle Mediated by a H2-CO2 Greenhouse
For decades, the scientific community has been trying to reconcile abundant evidence for fluvialactivity on Noachian and early Hesperian Mars with the faint young Sun and reasonableconstraints on ancient atmospheric pressure and composition. Recently, the investigation of H2-CO2 collision-induced absorption has opened up a new avenue to warm Noachian Mars. We usethe ROCKE-3D global climate model to simulate plausible states of the ancient Martian climatewith this absorptive warming and reasonable constraints on surface paleopressure. We find that1.5-2 bar CO2-dominated atmospheres with ≥3% H2 can produce global mean surfacetemperatures above freezing, while also providing sufficient warming to avoid surfaceatmospheric CO2 condensation at 0°-45° obliquity. Simulations conducted with both moderntopography and a paleotopography, before Tharsis formed, highlight the importance of Tharsis asa cold trap for water on the planet. Additionally, we find that low obliquity (modern and 0°) ismore conducive to rainfall over valley network locations than high (45°) obliquity. PLAIN LANGUAGE SUMMARY Much evidence tells us that ancient Mars had liquid water on its surface. But reconciling thatwith the fainter young Sun and reasonable constraints on Mars' early atmosphere is challenging.We use a 3D global climate model with an atmosphere containing some amount of hydrogen (apossible greenhouse gas) to study Mars' early hydrological cycle and compare it to the geologicalevidence of surface liquid water billions of years ago. We find that hydrogen and carbon dioxidetogether in an atmosphere twice as thick as modern Earth can warm early Mars above thefreezing point of water and that a low axial tilt produces rainfall patterns that best match thegeologic evidence
Scott D Guzewich↗
R Coronae Borealis Star Evolution: Simulating 3D Merger Events to 1D Stellar Evolution Including Large-scale Nucleosynthesis
R Coronae Borealis (RCB) stars are rare hydrogen-deficient carbon-rich variable supergiants thought to be the result of dynamically unstable white dwarf mergers. In this article, we attempt to model RCB stars through all the relevant timescales by simulating a merger event in Octo-tiger, a 3D adaptive mesh refinement (AMR) hydrodynamics code, and mapping the post-merger object into MESA, a 1D stellar evolution code. We then post-process the nucleosynthesis on a much larger nuclear reaction network to study the enhancement of s-process elements. We present models that match observations or previous studies in most surface abundances, isotopic ratios, early evolution, and lifetimes. We also observe similar mixing behavior to previous modeling attempts that result in the partial He-burning products visible on the surface in observations. However, we do note that our subsolar models lack any enhancement in s-process elements, which we attribute to a lack of hydrogen in the envelope. We also find that the 16 O/ 18 O isotopic ratio is very sensitive to initial hydrogen abundance and increases outside of the acceptable range with a hydrogen mass fraction greater than 10 -4 .
79 ASTRONOMY AND ASTROPHYSICS↗
3D Simulations of gas injection on callovo-oxfordian claystone assuming spatial heterogeneity and anisotropy
A series of gas injection tests on Callovo-Oxfordian (COx) claystone from the Bure underground research laboratory (URL) in France were carried out at the British Geological Survey (BGS). The tests were performed using a triaxial apparatus specifically designed to capture small volumetric strains induced by the injected gas flow and consequent material dilatancy. The long-duration experiments were monitored throughout. Measurements also included pressure, stresses (axial and radial stresses prescribed for each test stage), rate of gas inflow, gas outflow volume as well as pore-pressures observed at various points of the sample. A coupled hydro-gas-mechanical 3D numerical model has been developed to simulate the tests. Initial permeability is assumed heterogeneous throughout the specimen and embedded fractures are incorporated in the formulation. Gas pressure-induced deformations during the test lead to variations of permeability due to changes in matrix porosity and, especially, fracture aperture as well as fracture orientation due to material anisotropy. A programme of sensitivity analyses involving the variation of different aspects and parameters of the model contributes to a better understanding of the phenomena and highlights its complexity. The model is able to reproduce the observed behaviour of the tests.
58 GEOSCIENCES↗
3D SIMULATIONS OF REALISTIC POWER HALOS IN MAGNETOHYDROSTATIC SUNSPOT ATMOSPHERES: LINKING THEORY AND OBSERVATION
The well-observed acoustic halo is an enhancement in time-averaged Doppler velocity and intensity power with respect to quiet-Sun values that is prominent for the weak and highly inclined field around the penumbra of sunspots and active regions. We perform 3D linear wave modeling with realistic distributed acoustic sources in a magnetohydrostatic sunspot atmosphere and compare the resultant simulation enhancements with multiheight SDO observations of the phenomenon. We find that simulated halos are in good qualitative agreement with observations. We also provide further proof that the underlying process responsible for the halo is the refraction and return of fast magnetic waves that have undergone mode conversion at the critical a = c atmospheric layer. In addition, we also find strong evidence that fast Alfvén mode conversion plays a significant role in the structure of the halo, taking energy away from photospheric and chromospheric heights in the form of field-aligned Alfvén waves. This conversion process may explain the observed “dual-ring” halo structure at higher (>8 mHz) frequencies.
79 ASTRONOMY AND ASTROPHYSICS↗
Relevence of 3D Simulations and Sandwich Core Topology for the Modeling of Honeycomb Core Sandwich Panels Undergoing Mode I Delamination
A recently developed cohesive zone traction-separation law, which includes the effects of fiber bridging in a novel way, is extended from 2D to 3D. The proposed cohesive model is applied to low fidelity (i.e. homogenized core representation) and high fidelity (i.e. directly accounting for the core topology) finite element models of a composite panel comprised of carbon fiber reinforced plastic facesheets and a honeycomb sandwich core. This enables the investigation of 2D to 3D parameter transferability, width-dependent effects such as thumbnail-shaped crack growth, and the verification of plane strain / plane stress assumptions. A pronounced curvature of the initial interface-related crack front is observed, while the bridging-related crack front is straight. Furthermore, it is found that the cohesive parameters can easily be transferred from 2D to 3D under plane stress assumptions, but not under plane strain assumptions. The numerical predictions are compared to experimental load-displacement and R-curves.
Composites↗
RELAP5-3D Simulation of PG-27 Test at the HTTF Facility
The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a quarter scaled integral-effect test facility designed to examine the transient phenomena occurring in a Modular Hight Temperature Gas-cooled Reactor (MHTGR). It is specifically designed to handle the pressurized and depressurized conduction cooldown (PCC and DCC) events. The PCC event is an accident scenario in which there is a loss of forced convection of the coolant through the system. However, the pressure boundary remains intact. This study focuses on the Reactor Excursions and Leak Analysis Program (RELAP)5-3D simulation of the PG-27 test (PCC phenomena) and the sensitivity of those simulations from the material properties and primary helium mass flow rate standpoint. The aim is to capture a reasonable one-dimensional picture of a three dimensional phenomena. The simulation will help in the code assessment and in understanding the various unknown parameters involved. The presented results will demonstrate the code is able to capture the detailed one-dimensional physics accurately. However, input parameters to the code (i.e., the primary mass flow rate, refractory material thermal conductivity, and heat capacity at elevated temperatures) impact the fuel temperature the most. This work will discuss the sensitivity of simulations to the identified input parameters.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
3D Simulations of the Richtmyer-Meshkov Instability with Re-Shock
We present results of inviscid simulations, in three dimensions, of Richtmyer-Meshkov instability for high incident shock Mach number. The growth rate of a single harmonic perturbation is quantified and compared with the results of a 2D calculation. Upon re-shock, the perturbation amplitude undergoes a phase reversal while the mean velocity of the interface is zero. Before re-shock the normalized growth rate of a 2D and 3D interface are nearly the same, but the growth rate after re-shock is significantly larger for the 3D than the 2D case. We also examine the evolution of multiple harmonic perturbations. Computational and parallelization issues of the simulation code will also be briefly discussed. The computations were done on the T3E at Pittsburgh Supercomputing Center.
Meiron, Daniel I.↗