SEARCH · Engineering Papers
Results for “Numerical 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.
Assessment of uncertainty in the numerical simulation of solar irradiance over inclined PV panels: New algorithms using measurements and modeling tools
Development of accurate transposition models to simulate plane-of-array (POA) irradiance from horizontal measurements or simulations is a complex process mainly because of the anisotropic distribution of diffuse solar radiation in the atmosphere. The limited availability of reliable POA measurements at large temporal and spatial scales leads to difficulties in the comprehensive evaluation of transposition models. This paper proposes new algorithms to assess the uncertainty of transposition models using both surface-based observations and modeling tools. We reviewed the analytical derivation of POA irradiance and the approximation of isotropic diffuse radiation that simplifies the computation. Two transposition models are evaluated against the computation by the rigorous analytical solution. We proposed a new algorithm to evaluate transposition models using the clear-sky measurements at the National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory (SRRL) and a radiative transfer model that integrates diffuse radiances of various sky-viewing angles. We found that the radiative transfer model and a transposition model based on empirical regressions are superior to the isotropic models when compared to measurements. We further compared the radiative transfer model to the transposition models under an extensive range of idealized conditions. Our results suggest that the empirical transposition model has slightly higher cloudy-sky POA irradiance than the radiative transfer model, but performs better than the isotropic models under clear-sky conditions. Significantly smaller POA irradiances computed by the transposition models are observed when the photovoltaics (PV) panel deviates from the azimuthal direction of the sun. The new algorithms developed in the current study have opened the door to a more comprehensive evaluation of transposition models for various atmospheric conditions and solar and PV orientations.
Side-jet effects in high-pressure turbulent flows: Direct Numerical Simulation of nitrogen injected into carbon dioxide
Not Available
Numerical simulation of the coalescence of two bubbles in an ultrasound field
Not Available
Phase contrast imaging measurements and numerical simulations of turbulent density fluctuations in gas-fuelled ECRH discharges in Wendelstein 7-X
The fundamental nature of turbulent density fluctuations in standard Wendelstein 7-X (W7-X) stellarator discharges is investigated experimentally via phase contrast imaging (PCI) in combination with gyrokinetic simulations with the code GENE. We find that density fluctuations are ion-temperature-gradient-driven and radially localised in the outer half of the plasma. It is shown that the line-integrated PCI measurements cover the right range of wavenumbers and a favourable toroidal and poloidal location to capture some of the strongest density fluctuations in W7-X. Due to the radial localisation of fluctuations, measured wavenumber–frequency spectra exhibit a dominant phase velocity, which can be related to the E × B rotation velocity at the radial position of a well in the neoclassical radial electric field. The match is robust against variations of heating power and line-integrated density, which is partly due to the localisation of fluctuations and partly due to effects of the radial gradient in the E × B velocity profile on the wavenumber–frequency spectrum. The latter effect is studied with a newly built synthetic PCI diagnostic and global gyrokinetic simulations with GENE-3D.
Theoretical Model and Numerical Simulation of Adsorption and Deformation in Flexible Metal–Organic Frameworks
We report that metal–organic frameworks (MOFs) have recently attracted considerable attention as new nanoporous materials with applications to separation of fluid mixtures, catalysis, gas capture and storage, and drug delivery. A subclass of MOFs, the MIL-53 (with Al or Cr) family, exhibits a complex structural phase transition, often referred to as the breathing transition, which occurs when gas molecules are adsorbed in its pore space. In this phenomenon, the material morphology oscillates between two distinct phases, usually referred to as the narrow-pore (np) and large-pore (lp) phases. We describe a statistical mechanical model based on the energetics of the system that couples the adsorbates with the host and the ensuing structural deformation of the materials. The Hamiltonian of the system consists of the elastic energy of the MOF, the fluid phase energy, and the interaction energy between the gas and MOFs. Minimizing the total energy with respect to both the gas density and the displacement field in the material yields the governing equations for both. The model is used to study the adsorption of CO 2 and CH 4 in MIL-53(Al). Here, the results demonstrate that all reported experimental features of the phenomenon are produced by the model. In particular, consistent with experimental data, the model predicts that the material experiences contraction in the np phase, whereas it undergoes swelling in the lp phase, and that the two phases coexist in the transition zone. The model is, however, completely general and is applicable to any type of guest-responsive material that undergoes deformation as a result of its interaction with the guest molecules.
The Mechanism of Methane Gas Migration Through the Gas Hydrate Stability Zone: Insights From Numerical Simulations
Free gas migration through the gas hydrate stability zone (GHSZ) and subsequent gas seepage at the seabed are characteristic features in marine gas hydrate provinces worldwide. The biogenic or thermogenic gas is typically transported along faults from deeper sediment strata to the GHSZ. Several mechanisms have been proposed to explain free gas transport through the GHSZ. While inhibition of hydrate formation by elevated salinities and temperatures have been addressed previously in studies simulating unfocused, area-wide upward advection of gas, which is not adequately supported by field observations, the role of focused gas flow through chimney-like structures has been underappreciated in this context. Our simulations suggest that gas migration through the GHSZ is, fundamentally, a result of methane gas supply in excess of its consumption by hydrate formation. The required high gas flux is driven by local overpressure, built up from gas accumulating below the base of the GHSZ that fractures the overburden when exceeding a critical pressure, thereby creating the chimney-like migration pathway. Initially rapid hydrate formation raises the temperature in the chimney structure, thereby facilitating further gas transport through the GHSZ. As a consequence, high hydrate saturations form preferentially close to the seafloor, where temperatures drop to bottom water values, producing a prominent subsurface salinity peak. Over time, hydrates form at a lower rate throughout the chimney structure, while initial temperature elevation and salinity peak dissipate. Thus, our simulations suggest that the near-surface salinity peak and elevated temperatures are a result of transient high-flux gas migration through the GHSZ.
Parameterizations of Entrainment-Mixing Mechanisms and their Effects on Cloud Droplet Spectral Width Based on Numerical Simulations
Entrainment-mixing mechanisms significantly affect cloud droplet number concentration, radius, and spectral shape. Quantitative examination of entrainment-mixing effects on cloud droplet spectral width is lacking. Here, we examine the effects of entrainment-mixing processes on cloud microphysics by 12,218 different setups, each simulated 10 times using the Explicit Mixing Parcel Model (EMPM) driven by the observational data from the Third Tibetan Plateau Atmospheric Scientific Experiment (TIPEX-III) campaign. Parameterizations of entrainment-mixing mechanisms are developed by relating homogeneous mixing degree to transition scale number that depends on the dissipation rate and droplet evaporation time scale. The correlation between relative dispersion of cloud droplet size distribution and homogeneous mixing degree changes from negative to positive with the decreasing homogeneous mixing degree. The different relationships are closely related to the competition between complete and partial droplet evaporation and the number concentration of small droplets, which are quantitatively described by two newly introduced dimensionless numbers. The competition is significantly affected by relative humidity and mixing fraction of entrained air as well as turbulence dissipation rate, but not much by cloud droplet number concentration. Especially, when relative humidity and dissipation rate are high, there is only a negative correlation. This study sheds new light on generalizing the homogeneous/inhomogeneous concept by considering relative dispersion, and also provides parameterizations of entrainment-mixing processes and relative dispersion for atmospheric models.
Numerical Simulations of the Geospace Response to the Arrival of an Idealized Perfect Interplanetary Coronal Mass Ejection
Abstract Previously, Tsurutani and Lakhina (2014, https://doi.org/10.1002/2013GL058825 ) created estimates for a “perfect” interplanetary coronal mass ejection and performed simple calculations for the response of geospace, including . In this study, these estimates are used to drive a coupled magnetohydrodynamic‐ring current‐ionosphere model of geospace to obtain more physically accurate estimates of the geospace response to such an event. The sudden impulse phase is examined and compared to the estimations of Tsurutani and Lakhina (2014, https://doi.org/10.1002/2013GL058825 ). The physics‐based simulation yields similar estimates for Dst rise, magnetopause compression, and equatorial values as the previous study. However, results diverge away from the equator. values in excess of 30 nT/s are found as low as magnetic latitude. Under southward interplanetary magnetic field conditions, magnetopause erosion combines with strong region one Birkeland currents to intensify the response. Values obtained here surpass those found in historically recorded events and set the upper threshold of extreme geomagnetically induced current activity at Earth.
Spontaneous imbibition in a square tube with corner films: theoretical model and numerical simulation
Spontaneous imbibition in an angular tube with corner films is a fundamental problem in many scientific and engineering processes. In this paper, a modified interacting capillary bundle model is developed to describe the liquid imbibition dynamics in a square tube with corner films. The square tube is decomposed into several interacting sub-capillaries and the local capillary pressure in each sub-capillary is derived based on the specific shape of its meniscus. The conductance of each sub-capillary is calculated using single-phase lattice Boltzmann simulation. The modified interacting capillary bundle model and color-gradient lattice Boltzmann method are used to simulate the liquid imbibition dynamics in the square tube with different fluid properties. The predictions by the modified interacting capillary bundle model match well with the lattice Boltzmann simulation results for different conditions, demonstrating the accuracy and robustness of the interacting capillary bundle model to describe the imbibition dynamics with corner films. In addition, the interacting capillary bundle model is helpful to investigate the mechanisms during spontaneous imbibition and the influences of fluid viscosity, surface tension, wetting phase contact angle and gravity on imbibition dynamics. Finally, a universal scaling law of imbibition dynamics for the main meniscus is developed and the scaling law for arc meniscus is also analyzed.
Numerical simulation of ion charge breeding in electron beam ion source
The Electron Beam Ion Source particle-in-cell code (EBIS-PIC) tracks ions in an EBIS electron beam while updating electric potential self-consistently and atomic processes by the Monte Carlo method. Recent improvements to the code are reported in this paper. The ionization module has been improved by using experimental ionization energies and shell effects. The acceptance of injected ions and the emittance of extracted ion beam are calculated by extending EBIS-PIC to the beam line transport region. An EBIS-PIC simulation is performed for a Cs charge-breeding experiment at BNL. The charge state distribution agrees well with experiments, and additional simulation results of radial profiles and velocity space distributions of the trapped ions are presented.
Numerical simulation of temperature field, microstructure evolution and mechanical properties of HSS during hot stamping
The hot stamping of boron steels is widely used to produce ultra high strength automobile components without any spring back. The ultra high strength of final products is attributed to the fully martensitic microstructure that is obtained through the simultaneous forming and quenching of the hot blanks after austenization. In the present study, a mathematical model incorporating both heat transfer and the transformation of austenite is presented. A FORTRAN program based on finite element technique has been developed which permits the temperature distribution and microstructure evolution of high strength steel during hot stamping process. Two empirical diffusion-dependent transformation models under isothermal conditions were employed respectively, and the prediction capability on mechanical properties of the models were compared with the hot stamping experiment of an automobile B-pillar part.
A numerical simulation study of gallium-phosphide/silicon heterojunction passivated emitter and rear solar cells
The performance of passivated emitter and rear (PERC) solar cells made of p-type Si wafers is often limited by recombination in the phosphorus-doped emitter. To overcome this limitation, a realistic PERC solar cell is simulated, whereby the conventional phosphorus-doped emitter is replaced by a thin, crystalline gallium phosphide (GaP) layer. The resulting GaP/Si PERC cell is compared to Si PERC cells, which have (i) a standard POCl{sub 3} diffused emitter, (ii) a solid-state diffused emitter, or (iii) a high efficiency ion-implanted emitter. The maximum efficiencies for these realistic PERC cells are between 20.5% and 21.2% for the phosphorus-doped emitters (i)–(iii), and up to 21.6% for the GaP emitter. The major advantage of this GaP hetero-emitter is a significantly reduced recombination loss, resulting in a higher V{sub oc}. This is so because the high valence band offset between GaP and Si acts as a nearly ideal minority carrier blocker. This effect is comparable to amorphous Si. However, the GaP layer can be contacted with metal fingers like crystalline Si, so no conductive oxide is necessary. Compared to the conventional PERC structure, the GaP/Si PERC cell requires a lower Si base doping density, which reduces the impact of the boron-oxygen complexes. Despite the lower base doping, fewer rear local contacts are necessary. This is so because the GaP emitter shows reduced recombination, leading to a higher minority electron density in the base and, in turn, to a higher base conductivity.
Numerical simulation of the sustaining discharge in radio frequency hollow cathode discharge in argon
In this paper, a two-dimensional fluid model was developed to study the radio frequency (RF) hollow cathode discharge (HCD) in argon at 1 Torr. The evolutions of the particle density distribution and the ionization rate distribution in RF HCD at 13.56 MHz indicate that the discharge mainly occurs inside the hollow cathode. The spatio-temporal distributions of the ionization rate and the power deposition within the hollow cathode imply that sheath oscillation heating is the primary mechanism to sustain the RF HCD, whereas secondary electron emission plays a negligible role. However, as driving frequency decreases, secondary electron heating becomes a dominant mechanism to sustain the discharge in RF hollow cathode.