Engineering Papers⌕ Search

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.

At least 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

A novel Lagrangian approach for the stable numerical simulation of fault and fracture mechanics

The simulation of the mechanics of geological faults and fractures is of paramount importance in several applications, such as ensuring the safety of the underground storage of wastes and hydrocarbons or predicting the possible seismicity triggered by the production and injection of subsurface fluids. However, the stable numerical modeling of ground ruptures is still an open issue. The present work introduces a novel formulation based on the use of the Lagrange multipliers to prescribe the constraints on the contact surfaces. The variational formulation is modified in order to take into account the frictional work along the activated fault portion according to the principle of maximum plastic dissipation. The numerical model, developed in the framework of the Finite Element method, provides stable solutions with a fast convergence of the non-linear problem. The stabilizing properties of the proposed model are emphasized with the aid of a realistic numerical example dealing with the generation of ground fractures due to groundwater withdrawal in arid regions. - Highlights: • A numerical model is developed for the simulation of fault and fracture mechanics. • The model is implemented in the framework of the Finite Element method and with the aid of Lagrange multipliers. • The proposed formulation introduces a new contribution due to the frictional work on the portion of activated fault. • The resulting algorithm is highly non-linear as the portion of activated fault is itself unknown. • The numerical solution is validated against analytical results and proves to be stable also in realistic applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Numerical Simulation with Experimental Validation of the Draping Behavior of Woven Fabrics

Woven fabric composites are extensively used in molding complex geometrical shapes due to their high conformability compared to other fabrics. Preforming is an important step in the overall process. In this step, the two-dimensional fabric is draped to become the three-dimensional shape of the part prior to resin injection. During preforming, the orientation of the tows may change significantly compared to the initial orientations. Accurate prediction of the tow orientations after molding is important for evaluating the structural performance of the final part. This paper investigates the fiber angle changes for carbon fiber woven fabrics during draping over a truncated pyramid tool designed and fabricated at the General Motors Research Labs. This aspect of study is a subset of the broad study conducted under the purview of a Department of Energy project funded to GM in developing state of the art computational tools for integrated manufacturing and structural performance prediction of carbon fiber composites. Fabric bending, picture frame testing, and bias-extension evaluations were carried out to determine the material parameters for these fabrics. The PAM-FORM computer program was used to model the draping behavior of these fabrics. Following deformation, fiber angle changes at different locations on the truncated pyramid were measured experimentally. The predicted angles matched the experimental results well as measured along the centerline and at several different locations on the deformed fabric. Details of the test methods used as well as the numerical results with various simulation parameters will be provided.

42 ENGINEERING↗
DOE OSTI2016

Numerical simulations of a horizontal axis water turbine designed for underwater mooring platforms

In order to extend the operational life of Underwater Moored Platforms (UMPs), a horizontal axis water turbine is designed to supply energy for the UMPs. The turbine, equipped with controllable blades, can be opened to generate power and charge the UMPs in moored state. Three-dimensional Computational Fluid Dynamics (CFD) simulations are performed to study the characteristics of power, thrust and the wake of the turbine. Particularly, the effect of the installation position of the turbine is considered. Simulations are based on the Reynolds Averaged Navier-Stokes (RANS) equations and the shear stress transport k-ω turbulent model is utilized. The numerical method is validated using existing experimental data. The simulation results show that this turbine has a maximum power coefficient of 0.327 when the turbine is installed near the tail of the UMP. The flow structure near the blade and in the wake are also discussed.

13 HYDRO ENERGY↗
DOE OSTI2014

Numerical simulation of shock initiation of Ni/Al multilayered composites

The initiation of chemical reaction in cold-rolled Ni/Al multilayered composites by shock compression is investigated numerically. A simplified approach is adopted that exploits the disparity between the reaction and shock loading timescales. The impact of shock compression is modeled using CTH simulations that yield pressure, strain, and temperature distributions within the composites due to the shock propagation. The resulting temperature distribution is then used as initial condition to simulate the evolution of the subsequent shock-induced mixing and chemical reaction. To this end, a reduced reaction model is used that expresses the local atomic mixing and heat release rates in terms of an evolution equation for a dimensionless time scale reflecting the age of the mixed layer. The computations are used to assess the effect of bilayer thickness on the reaction, as well as the impact of shock velocity and orientation with respect to the layering. Computed results indicate that initiation and evolution of the reaction are substantially affected by both the shock velocity and the bilayer thickness. In particular, at low impact velocity, Ni/Al multilayered composites with thick bilayers react completely in 100 ms while at high impact velocity and thin bilayers, reaction time was less than 100 μs. Quantitative trends for the dependence of the reaction time on the shock velocity are also determined, for different bilayer thickness and shock orientation.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Numerical simulation of the dynamics of a liquid crystal in the case of plane strain using GPUs

Parallel computational algorithm is worked out for numerical implementation of two-dimensional dynamic model of a liquid crystal, which takes into account the mutual influence of three physical effects – transfer of acoustic energy due to translational motion, viscoelastic rotation of the particles under the action of tangential stresses, anisotropic thermal expansion and compression. The algorithm is based on the Godunov gap decay method, the Ivanov scheme with controlled dissipation of energy and the splitting method with respect to spatial variables. The CUDA technology for computer systems with graphic accelerators is used. Results of computations demonstrating the efficiency of proposed method and algorithm are represented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Numerical Simulations of Cathodoluminescence Measurements in Thin-Film Solar Cells

Cathodoluminescence (CL) measurements can be applied to assess grain-boundary (GB) and grain-interior (GI) recombination in thin-film solar cell materials and made quantitative if we can develop CL models that account for material and measurement complexities. Recently, we developed a three-dimensional numerical CL model, based in MATLAB, that simulates the GI CL intensity as a function of four parameters: grain size, GI lifetime, and GB and surface recombination velocities. The model assumes that GB electrostatic potentials are screened by the high excesscarrier densities used in the CL measurement such that transport is governed by ambipolar diffusion. Here, we develop models to address directly GB potentials and their effects on these measurements. First, we transfer the MATLAB-based model to COMSOL software, and then introduce shallow donors to the GBs to produce potentials. We also develop a two-dimensional model in COMSOL to simulate CL GB contrast with GB potentials. Simulations indicate that GB potentials can increase or decrease CL intensities relative to the zero-potential case. However, the high electron-beam currents typically applied in CL measurements minimize the impact of GB potentials.

14 SOLAR ENERGY↗
DOE OSTI2020

Numerical simulation of high-efficiency, scalable, all-back-contact Cd(Se,Te) solar cells

In this work, all-back-contact thin film photovoltaic devices in lattice back contact (LBC) configuration are studied by device simulation to set the requirements for achieving efficiency > 20% with lateral contact spacing > 50 µm. The numerical device model is first validated against published data for an LBC perovskite device and extended to predict the performance of Cd(Se,Te)-based polycrystalline devices. Recent advances in the electronic properties of Cd(Se,Te) films have made high efficiency LBC devices feasible with industrially-scalable contact dimensions. Constraints on bulk, interface, and contact properties are quantified and discussed

14 SOLAR ENERGY↗
DOE OSTI2013

Numerical simulation study of positron production by intense laser-accelerated electrons

Positron production by ultra-intense laser-accelerated electrons has been studied with two-dimensional particle-in-cell and Monte Carlo simulations. The dependence of the positron yield on plasma density, plasma length, and converter thickness was investigated in detail with fixed parameters of a typical 100 TW laser system. The results show that with the optimal plasma and converter parameters a positron beam containing up to 1.9 × 10{sup 10} positrons can be generated, which has a small divergence angle (10°), a high temperature (67.2 MeV), and a short pulse duration (1.7 ps)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analyzing Relationships between Tropical Cyclone Intensity and Rain Rate over the Ocean Using Numerical Simulations

In this study, the relationship between tropical cyclone (TC) intensity and rain rate over the ocean is investigated using a full-physics numerical model (WRF) and a physics-based TC rainfall model (TCR). TC intensity is found to be nearly linearly correlated with the average rain rate in the inner core [~0.97 (mm h –1 m –2 )/(m s –1 )], while the correlation is weak at outer radii. This difference is induced because TC intensity is significantly correlated with both the vertical velocity and specific humidity in the inner core but is not significantly correlated with the vertical velocity in the outer radii. Further investigation shows that the intensity–rain-rate relationship at the outer radii is influenced by the TC evolution stage. The rain rate for the outer radii is positively correlated with TC intensity for nondecaying TCs, while this correlation is reduced for decaying TCs due to systematic downdrafts in the outer radii. In the context of climate change, the sensitivity of the TC rain rate to sea surface temperature (SST) is found to be +9% per 1 K increase of SST, roughly the product of the sensitivity of TC intensity to SST (+3%) and the Clausius-Clapeyron scaling (+7%). Coupled with synthetic storms, evolution of the TC rain rate over the twenty-first century under the SSP5-8.5 scenario is projected by the TCR (calibrated with the WRF simulations). The annual increase rates of averaged TC rain rate are 0.17% and 0.20% for the inner core and outer radii, respectively, larger than the annual increase rate of TC intensity (0.046%) but comparable to that of cube of intensity (0.18%).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Light hydrocarbon and noble gas migration as an analogue for potential CO 2 leakage: numerical simulations and field data from three hydrocarbon systems

Abstract Sites where anthropogenic CO 2 captured from industrial sources is stored in deep geological formations for climate change mitigation are required to show secure retention of the injected CO 2 . Monitoring, reporting, and verification (MRV) plans are needed to indicate that no CO 2 release has occurred. We explored the degree to which direct comparison between a surface anomaly and reservoir geochemistry using various geochemical parameters can be used for attribution. We used data collected on light hydrocarbons and noble gases throughout the sedimentary column at three CO 2 enhanced oil recovery (EOR) sites to understand the processes that may cause fluid evolution. Light hydrocarbon and noble gases were sampled from reservoirs, gas‐bearing intervals above reservoirs, and groundwater. Vertical profiles indicated that lighter components are relatively enriched during migration (i.e. chromatographic separation). Static and numerical models were designed to simulate episodic gas migration and geochemical alteration of these geochemical parameters from solubility and sorption. The effects of hydrocarbon solubility were minimal (Bernard ratio changes within 5.2%) although field data were within the range of expected alterations from sorption. Forward models of CO 2 migration and noble gas interactions showed that CO 2 stripping causes an enrichment of crustal noble gases. In areas where natural fluxes of CO 2 from depth are non‐existent, the occurrence of crustal noble gas signature may distinguish fugitive CO 2 from the reservoir from natural near‐surface sources, and could be considered to explain apparent fluid anomalies. © 2019 Society of Chemical Industry and John Wiley & Sons, Ltd.

Anderson, Jacob↗
DOE OSTI2019

Numerical Simulation of the Three-Phase Flow of a Bubble Interacting with the Steel–Slag Interface During the Secondary Refining Process

Gas stirring is an effective method for the secondary refining process of steel. The microphenomena among the phases of an argon bubble floating in steel and motion across the steel–slag interface were simulated and discussed with a sharp interface tracking method. A water-model experiment was set to validate the numerical model, and the experimental results showed agreement with the numerical model. A critical size for bubbles to pass the interface was found using the model.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Quasioptical modeling of wave beams with and without mode conversion: III. Numerical simulations of mode-converting beams

This work continues a series of papers where we propose an algorithm for quasioptical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. In the second paper, we presented a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams. Here, we report quasioptical simulations of mode-converting beams for the first time. We also demonstrate that PARADE can model splitting of two-mode beams. The numerical results produced by PARADE show good agreement with those of one-dimensional full-wave simulations and also with conventional ray tracing (to the extent that one-dimensional and ray-tracing simulations are applicable).

ECRH↗
DOE OSTI2019

Quasioptical modeling of wave beams with and without mode conversion. III. Numerical simulations of mode-converting beams

This work continues a series of papers where we propose an algorithm for quasi-optical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. In the second paper, we presented a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams. Here, we report quasi-optical simulations of mode-converting beams for the first time. Here, we also demonstrate that PARADE can model splitting of two-mode beams. The numerical results produced by PARADE show good agreement with those of one-dimensional full-wave simulations and also with conventional ray tracing (to the extent that one-dimensional and ray-tracing simulations are applicable).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

2D numerical simulation of unsteady flows for large scale floods prediction in real time

The challenge of finding a compromise between computational time and level of accuracy and robustness has traditionally expanded the use simplified models rather than full two-dimensional (2D) models for flood simulation. This work presents a GPU accelerated 2D shallow water model for the simulation of flood events in real time. In particular, an explicit first-order finite volume scheme is detailed to control the numerical instabilities that are likely to appear when used in complex topography. The model is first validated with the benchmark test case of the Toce River (Italy) and numerical fixes are demonstrated to be necessary. The model is next applied to reproduce real events in a reach of the Ebro River (Spain) in order to compare simulation results with field data. The second case deals with a large domain (744 km 2 ) and long flood duration (up to 20 days) allowing an analysis of the performance and speed-up achieved by different GPU devices. We report the high values of fit between observed and simulated results as well as the computational times achieved are encouraging to propose the use of the model as forecasting system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Experimental Validation of Numerical Simulations of a New-Generation NiTi Endodontic File Under Bending

The superelasticity of NiTi shape memory alloy has been used in endodontics since the 1990s. To study the mechanical behavior of endodontic instruments, a traditional approach consists in experimental investigations. However, finite element analysis constitutes another way to assess their mechanical behavior and to facilitate their design. The main aim of this study is to compare experimental and numerical bending results on different structures (NiTi wire, spreader, and instruments) to estimate the reliability of the finite element simulations. These investigations were carried out as follows. Firstly, experimental material parameters identification was performed using NiTi wires. These parameters were implemented in an appropriate NiTi model. Bending was numerically applied to the meshed structures generated by the finite element method. Experimental tests were performed on real structures under bending, with exactly the same loading, in order to compare experimental and numerical results. These results were in good agreement for each of the considered structures. This enabled the validation of the simulation results and the use of simulations to design new endodontic instruments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Estimation of water diffusion coefficient into polycarbonate at different temperatures using numerical simulation

Nowadays, many electronic systems are exposed to harsh conditions of relative humidity and temperature. Mass transport properties of electronic packaging materials are needed in order to investigate the influence of moisture and temperature on reliability of electronic devices. Polycarbonate (PC) is widely used in the electronics industry. Thus, in this work the water diffusion coefficient into PC is investigated. Furthermore, numerical methods used for estimation of the diffusion coefficient and their assumptions are discussed. 1D and 3D numerical solutions are compared and based on this, it is shown how the estimated value can be different depending on the choice of dimensionality in the model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Collisionless microinstabilities in stellarators. II. Numerical simulations

Microinstabilities exhibit a rich variety of behavior in stellarators due to the many degrees of freedom in the magnetic geometry. It has recently been found that certain stellarators (quasi-isodynamic ones with maximum-J geometry) are partly resilient to trapped-particle instabilities, because fast-bouncing particles tend to extract energy from these modes near marginal stability. In reality, stellarators are never perfectly quasi-isodynamic, and the question thus arises whether they still benefit from enhanced stability. Here, the stability properties of Wendelstein 7-X and a more quasi-isodynamic configuration, QIPC, are investigated numerically and compared with the National Compact Stellarator Experiment and the DIII-D tokamak. In gyrokinetic simulations, performed with the gyrokinetic code GENE in the electrostatic and collisionless approximation, ion-temperature-gradient modes, trapped-electron modes, and mixed-type instabilities are studied. Wendelstein 7-X and QIPC exhibit significantly reduced growth rates for all simulations that include kinetic electrons, and the latter are indeed found to be stabilizing in the energy budget. These results suggest that imperfectly optimized stellarators can retain most of the stabilizing properties predicted for perfect maximum-J configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Numerical simulation of the plasma generated by the interaction high-current electron beam with Al target

The results of one dimensional particle-in-cell simulations of the dynamics of plasma generated during the interaction of a high-energy ({<=}200 keV) and high-current ({<=}15 A/cm{sup 2}) electron beam with an aluminum target are presented. The generated target plasma is low-ionized and characterized by non-Maxwellian electron energy distribution. The density and electron temperature of the plasma, which expands toward the anode at a typical velocity of {approx}10{sup 5} cm/s, does not exceed 4 Multiplication-Sign 10{sup 14} cm{sup -3} and 1 eV, respectively, which is in satisfactory agreement with the experimental results presented in W. An et al., J. Appl. Phys. 110, 093304 (2011). The results of the simulations showed also acceleration of the ions from the target plasma toward the anode by the potential of the non-compensated space charge of the electron beam. The typical velocity of these energetic ions is {approx}10{sup 8} cm/s and depends on the electron current density and energy. These ions partially compensate the space charge of the electron beam, which leads to a decrease in the depth of the potential well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗