Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Analysis”

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

Numerical Analysis of Detonability Assessment in a Natural Gas-Air Fueled Rotating Detonation Engine

Rotating Detonation Engines (RDE) offer an alternative combustion strategy to replace conventional constant pressure combustion with a process that could produce a pressure gain without the use of a mechanical compressor. Recent numerical and experimental publications that consider air as the oxidizer have primarily focused on the ability of these annular combustors to sustain a stable continuous detonation wave when fueled by hydrogen. However, for this to be a viable consideration for the land-based power generation it is necessary to explore the ability to detonate natural gas and air within the confines of the annular geometry of an RDE. Previous studies on confined detonations have expressed the importance of permitting detonation cells to fully form within the combustor in order to achieve stability. This poses a challenge for natural gas–air fueled processes as their detonation cell size can be quite large even at moderate pressures. Despite the practical importance, only a few studies are available on natural gas detonations for air-breathing RDE applications. Moreover, the extreme thermodynamic condition (high temperature inside the combustor) allows limited accessibility inside the combustor for detailed experimental instrumentations, providing mostly single-point data. Recent experimental studies at the National Energy Technology Laboratory (NETL) have reported detonation failure at higher methane concentration in an air-breathing RDE fueled by natural gas-hydrogen fuel blends. This encourages to perform a detailed numerical investigation on the wave characteristics of detonation in a natural gas-air fueled RDE to understand the various aspects of instability associated with the natural gas-air detonation.This study is a numerical consideration of a methane-air fueled RDE with varying operating conditions to ascertain the ability to achieve a stable, continuous detonation wave. The simulations have been performed in a 2D unwrapped RDE geometry using the open-source CFD library “OpenFOAM” employing an unsteady pressure-based compressible reactive flow solver with a <em>k</em>–<em>ε</em> turbulence model in a structured rectangular grid system. Both reduced and detailed chemical kinetic models have been used to assess the effect of the chemistry on the detonation wave characteristics and the underlying flow features. A systematic grid sensitivity study has been conducted with various grid sizes to quantify the weakly stable overdriven detonation on a coarse mesh and oscillating features at fine mesh resolutions.The main focus of the current study is to investigate the effects of operating injection pressure on detonation wave characteristics of an air-breathing Rotating Detonation Engine (RDE) fueled with natural gas-hydrogen fuel blends. Wave speeds, peak pressures and temperatures, and dominant frequencies have been computed from the time histories. The flow structures were then visualized using 2D contours of temperature and species concentration.

Saha, Pankaj↗

Numerical Analysis of Liquid Metal MHD Flow and Heat Transfer for Open-Surface Li Divertor in FNSF

Within the ongoing U.S.-based program on the development of liquid metal plasma-facing components, numerical simulations and analyses are performed to address the feasibility of the open-surface Li divertor. In the previous scoping studies (Smolentsev, 2021), heat-removal capabilities of the divertor were assessed using a simplified flow model for a slug-type velocity profile and constant flow thickness. Here, new analyses take into account forces acting on the flowing Li layer. Three reduced-order mathematical models are applied under the conditions of the U.S. Fusion Nuclear Science Facility (FNSF) to access magnetohydrodynamic (MHD) flow development effects, velocity distribution, and surface waves: 1) fully developed MHD flow; 2) quasi-2-D developing MHD flow; and 3) multiphase MHD flow. The obtained results for MHD flows and the surface heat flux computed with the plasma code scrape-off layer plasma simulation for ITER (SOLPS-ITER) are then used as input data to compute the temperature distribution in the divertor by solving the convection–diffusion energy equation.

Smolentsev, Sergey↗

Numerical Analysis on Hydraulic Power Take-Off for Wave Energy Converter and Power Smoothing Methods: Preprint

One of the primary challenges for wave energy converter (WEC) systems is the fluctuating nature of wave resources, which require the WEC components to be designed to handle loads (i.e., torques, forces, and powers) that are many times greater than the average load. This approach requires a much greater power take-off (PTO) capacity than the average power output and indicates a higher cost for the PTO. Moreover, additional design requirements, such as battery storage, are needed, particularly for practical electrical grid connection, and can be a problem for sensitive equipment (e.g., radar, computing devices, and sensors). Therefore, it is essential to investigate potential methodologies to reduce the overall power fluctuation while trying to optimize the power output from WECs. In this study, a detailed hydraulic PTO model was developed and coupled with a time-domain hydrodynamics model (WEC-Sim) to evaluate the PTO efficiency for WECs and the trade-off between power output and fluctuation using different power smoothing methods, including energy storage, pressure relief mechanism, and a powerbased setpoint control method. The study also revealed that the maximum power fluctuation for WECs can be significantly reduced by one order of magnitude when these power smoothing methods are applied.

16 TIDAL AND WAVE POWER↗

Numerical analysis of the unintegrated double gluon distribution

We present detailed numerical analysis of the unintegrated double gluon distribution which includes the dependence on the transverse momenta of partons. The unintegrated double gluon distribution was obtained following the Kimber-Martin-Ryskin method as a convolution of the perturbative gluon splitting function with the collinear integrated double gluon distribution and the Sudakov form factors. We analyze the dependence on the transverse momenta, longitudinal momentum fractions and hard scales. We find that the unintegrated gluon distribution factorizes into a product of two single unintegrated gluon distributions in the region of small values of x, provided the splitting contribution is included and the momentum sum rule is satisfied.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A numerical analysis of a magnetocaloric refrigerator with a 16-layer regenerator

A numerical analysis was conducted to study a room temperature magnetocaloric refrigerator with a 16-layer parallel plates active magnetic regenerator (AMR). Sixteen layers of LaFeMnSiH having different Curie temperatures were employed as magnetocaloric material (MCM) in the regenerator. Measured properties data was used. A transient one dimensional (1D) model was employed, in which a unique numerical method was developed to significantly accelerate the simulation speed of the multi-layer AMR system. As a result, the computation speed of a multi-layer AMR case was very close to the single-layer configuration. The performance of the 16-layer AMR system in different frequencies and utilizations has been investigated using this model. To optimize the layer length distribution of the 16-layer MCMs in the regenerator, a set of 137 simulations with different MCM distributions based on the Design of Experiments (DoE) method was conducted and the results were analyzed. The results show that the 16-layer AMR system can operate up to 84% of Carnot cycle COP at a temperature span of 41 K, which cannot be obtained using an AMR with fewer layers. Here, the DoE results indicate that for a 16-layer AMR system, the uniform distribution is very close to the optimized design.

36 MATERIALS SCIENCE↗

Centrifuge modeling and numerical analysis on lateral performance of mono-bucket foundation for offshore wind turbines

The lateral loading performance of the mono-bucket foundation for offshore wind turbines is investigated in this study through a series of centrifuge tests and numerical analysis. Six mono-bucket models are considered to assess the influences of aspect ratio under the ultimate condition. The FE model is validated against the centrifuge test results, and the failure mechanism of the foundation is revealed. The parametric study and sensitivity analysis of the soil parameters are performed in the numerical simulations. A modified theoretical calculation method is proposed to estimate the ultimate lateral capacity of the suction bucket foundation by solving the limit equilibrium equations. A correction term is added in the calculation using the “m” method. The position of the rotation center is calculated, demonstrating a better accuracy compared to the previously reported geometric method. The modified calculation method is verified against the field test, the centrifuge test, and the laboratory test. The method is proved to be applicable in estimating the lateral performance of the mono-bucket foundation in a general form. Finally, this study aims to provide well-documented data and design references for practical engineers.

17 WIND ENERGY↗

Numerical Analysis of Novel Plate Type Heat Exchanger with Oval-Twisted Channels

The novel heat exchanger (HX) designs implementing geometrical and surface modification combinations are expected to perform better than traditional HX technologies. Applied enhancement techniques seek to achieve (1) higher overall heat transfer performance, (2) increased compactness, and (3) simplified or comparable manufacturability. One innovative enhancement technique is to generate swirling flow vortices induced by channel or tube twisting. The turbulence generated by the twisted cross-section greatly enhances the heat transfer rate with minimal increases in pressure drop. Plate-type HXs and Printed Circuit Heat Exchangers (PCHXs) are compact designs that achieve high heat transfer rates per unit volume by utilizing several small channels, which maximizes the heat transfer surface area between the hot and cold fluids. Currently, advanced manufacturing technologies enable the design and fabrication of compact-type units with complicated channel geometries to achieve the highest performance and meet the compactness criteria of innovative HX technology. The proposed HX design concept combines the compactness of plate-type HXs and twisted channels, which provide additional turbulence and flow swirl enhancement. The plate-type oval-twisted HX (PTOTHX) is a crossflow configuration, with 16 short channels on one side (for hot fluid) and 8 long channels on the perpendicular side (for cold fluid). The inlet plenums have flow guide vanes to redirect flow and produce uniformity across the various flow paths. The compact size and purportedly improved heat transfer performance of the PTOTHX investigated herein prove its viability in various applications. Some notable potential nuclear applications of the PTOTHX include reactor core, spent fuel cooling, and residual heat dissipation. To establish a reference case, circular channels (PTCHX) are also considered in the present study for comparison with (PTOTHX). This paper aims to outline the numerical analysis procedures for determining the viability of the PTOTHX by comparing its heat transfer performance with the PTCHX units. The computational study used STAR-CCM+, a commercial computational fluid dynamics (CFD) code. Sensitivity analysis and model selection studies are conducted to determine the appropriate mesh density and turbulence model to provide the reported results. Numerical analyses comparing the Nusselt number (Nu) of the PTOTHX design with a comparable HX unit, including a circular cross-section and no twisting (PTCHX), show an overall heat transfer performance increase of 29-55% for balanced flow and 29-59% for imbalanced flow. Oval-cross-sectional twisted channels induce swirling flow vortices, enhancing the working fluid's convective heat transfer capabilities.

25 ENERGY STORAGE↗

Numerical analysis of gas exhaust in Wendelstein 7-X using the direct simulation Monte Carlo method

The present work is focused on a 3D numerical assessment of the Wendelstein 7-X (W7-X) particle exhaust. For all the numerical simulations the direct simulation Monte Carlo solver of the DIVGAS workflow, has been employed. The complex 3D geometry of the sub-divertor region includes the pumping gap panel, supporting structures, cooling pipes as well as the cryo-vacuum pump. All the considered flow simulations correspond to the Standard magnetic configuration of W7-X. The main conclusions, which can be extracted from the present numerical analysis could be summarized as follows; The coupling between EMC3-EIRENE and DIVGAS, which considers the fact that the incoming neutral particle flux at the sub-divertor is based on realistic plasma background, has been demonstrated. Three plasma scenarios have been considered, for which is clearly seen that by increasing the heating power, the neutral pressure as well as the resulting pumping efficiency is increased. The obtained numerical results of the neutral pressure in the sub-divertor lie within a more general scan matrix, which assumes a wider range of incoming particle flux, namely 10 19 –10 24 (s −1 ). It has been observed that, the sub-divertor neutral pressure is proportional to the incoming neutral particle flux, with the effective pumping speed to be a constant of proportionality. The influence of switching off the cryo-vacuum pump on the sub-divertor pressure is rather modest and a weak increase of the neutral pressure in the sub-divertor is expected. Correlations of the sub-divertor pressure with the total incoming particle flux as well as the individual pumped flux at each of the AEH and AEP sections have been deduced. Moreover, it has been demonstrated that the influence of the incoming neutral particle flux on the albedo coefficient at the AEH and AEP pumping gaps is rather weak. All the above numerical findings will actively support the optimization of the W7-X particle exhaust, in view of future experimental campaigns.

DSMC method↗

Numerical analysis of ELM stability with rotation and ion diamagnetic drift effects in JET

Stability to the type-I edge localized mode (ELM) in JET plasmas was investigated numerically by analyzing the stability to a peeling-ballooning mode with the effects of plasma rotation and ion diamagnetic drift. The numerical analysis was performed by solving the extended Frieman–Rotenberg equation with the MINERVA-DI code. To take into account these effects in the stability analysis self-consistently, the procedure of JET equilibrium reconstruction was updated to include the profiles of ion temperature and toroidal rotation, which are determined based on the measurement data in experiments. With the new procedure and MINERVA-DI, it was identified that the stability analysis including the rotation effect can explain the ELM trigger condition in JET with ITER like wall (JET-ILW), though the stability in JET with carbon wall (JET-C) is hardly affected by rotation. The key difference is that the rotation shear in JET-ILW plasmas analyzed in this study is larger than that in JET-C ones, the shear which enhances the dynamic pressure destabilizing a peeling-ballooning mode. In addition, the increase of the toroidal mode number of the unstable MHD mode determining the ELM trigger condition is also important when the plasma density is high in JET-ILW. Though such modes with high toroidal mode number are strongly stabilized by the ion diamagnetic drift effect, it was found that plasma rotation can sometimes overcome this stabilizing effect and destabilizes the peeling-ballooning modes in JET-ILW.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical analysis of air dehumidification through electrospray-enhanced vortical flow cyclone separator

The dehumidification of air is a highly energy-intensive process, typically requiring vapor compression systems to cool the air significantly below its dew point temperature. Here, this study numerically investigates air dehumidification achieved by electrostatic and dielectrophoresis forces in a converging duct and centrifugal principles in cyclone separators. An electrospray was employed to energize water droplets, which then captured water vapor via dielectrophoresis. A novelty of the present study is the inclusion of a cyclone separator positioned downstream of the electrospray. This combination increased the dehumidification efficiency and controlled the flow of the condensate. By integrating electrospray technology and applying an electric field, this study examines the impact of particle separation efficiency, pressure drops, and cyclone performance on air dehumidification. The control parameters for air dehumidification were investigated. Three cyclones were modelled for airflow rates of 5 cubic feet per minute [cfm] (0.0024 m 3 /s), 38 cfm (0.0179 m 3 /s), and 200 cfm (0.094 m 3 /s). The results indicated air dehumidification was up to 5.7 % under environmental conditions, similar to those used to validate the CFD model. The degree of air dehumidification was primarily influenced by the electric charge acquired by the water droplets, their flight trajectory, and the vorticity intensity within the cyclone. Air dehumidification could increase by up to 13 % by increasing the number of droplets injected into the air stream, the flight time, and the electrical charge. The breakdown electric field strength threshold near the electrospray was a limiting factor for further enhancing dehumidification performance. This study examines the optimization of an electro-assisted dehumidification system that includes cyclone-based droplet separation through numerical analysis. The results provide insights into the design of energy-efficient air treatment technologies.

42 ENGINEERING↗

Numerical analysis and wave tank validation on the optimal design of a two-body wave energy converter

To improve the performance of a ‘point absorber’ type wave energy converter (WEC), an additional submerged body can be deployed. The submerged body can be used to increase the equivalent excitation force on the WEC, as well as provide resonance tuning. This paper presents numerical analysis and experimental validation of a two-body point absorber type WEC using a mechanical motion rectifier (MMR) based power takeoff. The two-body point absorber consists of a floating buoy connected to a neutrally buoyant submerged body via the power-takeoff. The mechanical motion rectifier and a ball screw translate the relative heave motion of the two bodies into unidirectional rotation, which in turn spins a generator. Frequency domain analysis suggests there is an optimal submerged body mass for maximum WEC power absorption. Regular wave simulations in the time domain are compared to the results obtained in the frequency domain. While the time domain and frequency domain results predict the same optimal mass ratio, time domain analysis provides a more complex and accurate power result. To validate the time domain model, experimental wave tank testing is conducted using a 1:30 scale model WEC. The experiment shows the two-body WEC can produce twice the amount of power as the single-body WEC with same floating buoy and can be further increased by PTO design and power electronics optimization. Furthermore, wave tank testing shows the two-body WEC has a capture width ratio up to 58% at 59 kW/m and 51% at 36 kW/m.

16 TIDAL AND WAVE POWER↗

Analytical and numerical analysis of the slope of von Mises planar trusses

In the present paper, there are presented post-critical stress states which will occur at loading by vertical shift of the top joint in the direction downwards. The formation of certain stress states depends on the size of the angle formed by a straight beam of the von Mises planar truss with horizontal plane. Numerical and analytical methods and their problems with finding the angle were described. The numerical solution applies the method of searching for a minimum of potential energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical Analysis and Flow Induced Vibration Studies in Pressurized Helium Gas Cooling Channels

Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 from the irradiation of Mo-100 targets. The NMR target consists of an Inconel window that allows the electron beam to penetrate and irradiate a stack of Mo-100 discs. The irradiation process generates large amounts of heat and is cooled using pressurized helium gas flowing at 400 g/s. LANL provides both numerical analysis support and experimental support. Computational fluid dynamics (CFD) simulations for both 300 g/s and 400 g/s are performed on the NMR target to analyze the helium flow through the NMR system. The pressure drop, the mass flow rate through each of the disc and window channels and the force applied on the first disc are some of the quantities estimated using CFD. These results will be combined with Monte Carlo N-Particle (MCNP) estimates of heat flux into the system from the electron beam to perform conjugate heat transfer analysis (CHT). The experiments presented in this report focus on the effect of the flow induced vibrations from helium flowing through the coolant gaps between the discs. There is concern that the flow induced vibrations may lead to excessive mass loss of the enriched Mo-100 discs, and the goal would be to mitigate the mass loss. This work presents benchtop tests of flow through a subset of the target stack which involves high speed imaging and displacement measurements used to estimate disc vibration, as well as mass loss.

42 ENGINEERING↗

Numerical analysis of coalescence-induced bubble departure for enhanced boiling heat transfer

Boiling heat transfer plays a crucial role in a wide range of applications, such as power generation, refrigeration, electronics cooling, and pharmaceutics. Among the various factors that influence boiling heat transfer, the dynamics of vapor bubble nucleation, growth, and departure from the heated surface stand out as particularly important. An emerging phenomenon that can promote the departure of bubbles smaller than the Fritz diameter is coalescence-induced departure. If the dynamics of this process are fully understood, then surfaces can be engineered to promote faster bubble departure and substantially increase the performance of boiling heat transfer. Further, this work expands on published results by presenting a detailed numerical analysis of bubble coalescence and departure for a range of initial bubble diameters and size ratios between coalescing bubbles. Analysis of the results is focused on explaining how the release of surface energy and bubble surface dynamics lead to bubble departure, as well as fundamentally distinguishing capillary–inertial jumping and buoyant–inertial departure mechanisms across different bubble sizes and size ratios. The results show that both the initial sizes of the coalescing bubbles and the ratio between their sizes can determine whether the merged bubble will leave the surface through capillary–inertial jumping or buoyant departure. Below a certain bubble size, the release of surface energy by the merger is not sufficient to propel the merged bubble from the surface.

42 ENGINEERING↗

Longitudinal spatial hole burning in high-power semiconductor lasers: numerical analysis

Longitudinal spatial hole burning (LSHB) in high-power semiconductor lasers is analysed by numerically solving one-dimensional (1D) rate equations. Calculations are performed for GaAs-based lasers operating at a wavelength of 1.06 μm. It is shown that the LSHB-induced decrease in output power can be accounted for by two mechanisms: build-up of spontaneous recombination and decrease in slope efficiency, equivalent to a rise in internal optical loss. We analyse the influence of different laser chip parameters on the magnitude of the LSHB effect. In particular, it is shown that to suppress LSHB it is preferable to increase the optical confinement factor Γ. We examine the relationship between LSHB and other mechanisms capable of reducing the output power. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical analysis on seismic behavior of reinforced concrete beam to concrete filled steel tubular column connections with ring-beam

This paper presents numerical study of the seismic behavior of reinforced concrete beam to concrete filled steel tube column connections with ring-beam. The material stress-strain relations, element type and boundary condition are selected, which are consistent with actual situation. Then the seismic behavior of this type of joint are researched by ABAQUS, and finite element analyses are conducted under cyclic loading. Its parameters are discussed including thickness of steel tubular column wall, sectional dimension of the ring-beam and strength of the core concrete. The results show that the ultimate capacity of the connections is improved with sectional dimension of the ring-beam increased. In the meanwhile, the influence on skeleton curve of the joints is slight of which included thickness of steel tubular column wall and strength of the core concrete.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗