Search NASA⌕ Search

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.

At least 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modelling the development of capillary pressure in freshly 3D-printed concrete elements

3D concrete printing is a promising technology recently developed to automate construction. Since no formwork is used in this technology to support and protect fresh concrete, there are two aspects which considerably accelerate the development of capillary pressure in 3D-printed concrete in comparison to conventionally placed concrete: i) high stiffness of 3D-printed needed to provide sufficient buildability, and ii) very early and fast evaporation of pore water. Accelerated development of capillary pressure may lead to severe plastic shrinkage cracking in 3D-printed elements and, hence, need to be mitigated. This investigation aims at providing a poromechanical model for capillary pressure development in 3D-printed elements. To simulate the development of capillary pressure and plastic shrinkage, environmental factors, material properties, and element geometry need to be considered as a whole. The model inputs – coefficient of permeability, static bulk modulus, air entry pressure and chemical shrinkage rate – were determined experimentally. The model was validated for two fine-grained concretes. Both 3D-printed materials yielded faster capillary pressure increase in comparison to cast concrete, while partial substitution of cement with silica fume further accelerated the capillary pressure development. Furthermore, due to the lower permeability of the mixture containing silica fume, the gradient of capillary pressure between 3D-printed layers increased, as did the gradient of plastic shrinkage.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Massively parallel modeling and inversion of electrical resistivity tomography data using PFLOTRAN

Abstract. Electrical resistivity tomography (ERT) is a broadly accepted geophysical method for subsurface investigations. Interpretation of field ERT data usually requires the application of computationally intensive forward modeling and inversion algorithms. For large-scale ERT data, the efficiency of these algorithms depends on the robustness, accuracy, and scalability on high-performance computing resources. In this regard, we present a robust and highly scalable implementation of forward modeling and inversion algorithms for ERT data. The implementation is publicly available and developed within the framework of PFLOTRAN, an open-source, state-of-the-art massively parallel subsurface flow and transport simulation code. The forward modeling is based on a finite-volume discretization of the governing differential equations, and the inversion uses a Gauss–Newton optimization scheme. To evaluate the accuracy of the forward modeling, two examples are first presented by considering layered (1D) and 3D earth conductivity models. The computed numerical results show good agreement with the analytical solutions for the layered earth model and results from a well-established code for the 3D model. Inversion of ERT data, simulated for a 3D model, is then performed to demonstrate the inversion capability by recovering the conductivity of the model. To demonstrate the parallel performance of PFLOTRAN's ERT process model and inversion capabilities, large-scale scalability tests are performed by using up to 131 072 processes on a leadership class supercomputer. These tests are performed for the two most computationally intensive steps of the ERT inversion: forward modeling and Jacobian computation. For the forward modeling, we consider models with up to 122 ×106 degrees of freedom (DOFs) in the resulting system of linear equations and demonstrate that the code exhibits almost linear scalability on up to 10 000 DOFs per process. On the other hand, the code shows superlinear scalability for the Jacobian computation, mainly because all computations are fairly evenly distributed over each process with no parallel communication.

58 GEOSCIENCES↗
DOE OSTI2014

Height control of self-assembled quantum dots by strain engineering during capping

Strain engineering during the capping of III-V quantum dots has been explored as a means to control the height of strained self-assembled quantum dots. Results of Kinetic Monte Carlo simulations are confronted with cross-sectional Scanning Tunnel Microscopy (STM) measurements performed on InAs quantum dots grown by molecular beam epitaxy. We studied InAs quantum dots that are capped by In{sub x}Ga{sub (1−x)}As layers of different indium compositions. Both from our realistic 3D kinetic Monte Carlo simulations and the X-STM measurements on real samples, a trend in the height of the capped quantum dot is found as a function of the lattice mismatch between the quantum dot material and the capping layer. Results obtained on additional material combinations show a generic role of the elastic energy in the control of the quantum dot morphology by strain engineering during capping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Theory Aided Design of Active and Durable Nanoscale Cathode Catalysts (Final Report)

One of the greatest technological hurdles to deployment of fuel cells relates to the sluggish activity, low durability and the high cost of the catalysts that are currently employed. For automotive PEM fuel cells to become commercially viable, the Pt-specific power density would need to be reduced to less than 0.2gPt/kW (at cell voltages >0.65 V). This would require the Pt loadings to be less than 0.15 mg Pt /cm 2 MEA within the membrane electrode assembly. This could be achieved by enhancing the catalytic activity at the cathode, thus lowering its overpotential. Various different Pt-alloys have shown 2-4 times enhanced activities over Pt alone but still suffer some of the same durability issues as that of the pure Pt. There is a general loss of active Pt due to dissolution and sintering. While there have been a number of elegant fundamental experimental and theoretical studies on ideal single crystal Pt and Pt alloy surfaces which have helped to elucidate the factors that control the activity, there have been very few fundamental studies focused on the stability, reactivity and durability of well-defined Pt nanoparticles. We carried out ab initio density functional theory together with a novel double reference method that we developed to simulate constant potential electrochemical systems in order to model the electrocatalytic reduction of oxygen over model Pt alloy surfaces and nanoparticles. These simulations were used to probe the factors that control the electrocatalytic activity, guide the potential selection of new materials and test their stability under reaction conditions. Ab initio calculations were used to determine the reaction energies and activation barriers for a comprehensive array of different elementary adsorption, desorption, surface reaction and diffusion steps over Pt and Pt alloys involved in the electroreduction of oxygen as a function of surface coverage, the electrochemical potential and temperature. The calculations were used to simulate the potential dependent adsorption and surface reaction energies along with activation barriers in order to determine the kinetics for different surface structures and structural features (step edge and corner sites) to provide the necessary input for kinetic Monte Carlo simulations to follow the rates of reaction. More coarse-grained simulated annealing methods were used to help establish the lowest energy structures and morphologies for different Pt and Pt-alloy nanoparticles that form. The results from the DFT calculations were used to establish an ab initio-derived kinetic database that was used in both 2D and 3D kinetic Monte Carlo simulations that were used to follow the electrocatalytic performance over different particle sizes, shapes and compositions. The ab initio calculations together with the kinetic Monte Carlo simulations were used to complete the following objectives: 1) Determine the controlling elementary reaction pathways and intrinsic kinetics involved in ORR and their potential dependent behavior; 2) Establish the influence of the extrinsic reaction environment including surface structure, alloy composition and spatial arrangement and the humidity on ORR kinetics; 3) Elucidate the effects of particle size and morphology as well as the atomic structure and composition of nanoparticles of Pt, Pt 3 Co, Pt 3 Ni, Pt 3 Fe and other Pt-alloys on the electrocatalytic activity and reaction selectivity; 5) Construct surface Pourbaix phase diagrams for different Pt and Pt-alloys in order to map out the most stable surface phases for different material compositions and reaction conditions; 6) Elucidate the mechanisms that control dissolution and re-deposition of Pt.

08 HYDROGEN↗
DOE OSTI2021

Quantifying the unknown impact of segmentation uncertainty on image-based simulations

Abstract Image-based simulation, the use of 3D images to calculate physical quantities, relies on image segmentation for geometry creation. However, this process introduces image segmentation uncertainty because different segmentation tools (both manual and machine-learning-based) will each produce a unique and valid segmentation. First, we demonstrate that these variations propagate into the physics simulations, compromising the resulting physics quantities. Second, we propose a general framework for rapidly quantifying segmentation uncertainty. Through the creation and sampling of segmentation uncertainty probability maps, we systematically and objectively create uncertainty distributions of the physics quantities. We show that physics quantity uncertainty distributions can follow a Normal distribution, but, in more complicated physics simulations, the resulting uncertainty distribution can be surprisingly nontrivial. We establish that bounding segmentation uncertainty can fail in these nontrivial situations. While our work does not eliminate segmentation uncertainty, it improves simulation credibility by making visible the previously unrecognized segmentation uncertainty plaguing image-based simulation.

42 ENGINEERING↗
DOE OSTI2016

Quantitative three-dimensional characterization of block copolymer directed self-assembly on combined chemical and topographical prepatterned templates

Characterization of the three-dimensional (3D) structure in directed self-assembly (DSA) of block copolymers is crucial for understanding the complex relationships between the guiding template and the resulting polymer structure so DSA could be successfully implemented for advanced lithography applications. Here, we combined scanning transmission electron microscopy (STEM) tomography and coarse-grain simulations to probe the 3D structure of P2VP- b -PS- b -P2VP assembled on prepatterned templates using solvent vapor annealing. The templates consisted of nonpreferential background and raised guiding stripes that had PS-preferential top surfaces and P2VP-preferential sidewalls. The full 3D characterization allowed us to quantify the shape of the polymer domains and the interface between domains as a function of depth in the film and template geometry and offered important insights that were not accessible with 2D metrology. Sidewall guiding was advantageous in promoting the alignment and lowering the roughness of the P2VP domains over the sidewalls, but incommensurate confinement from the increased topography could cause roughness and intermittent dislocations in domains over the background region at the bottom of the film. The 3D characterization of bridge structures between domains over the background and breaks within domains on guiding lines sheds light on possible origins of common DSA defects. The positional fluctuations of the PS/P2VP interface between domains showed a depth-dependent behavior, with high levels of fluctuations near both the free surface of the film and the substrate and lower fluctuation levels in the middle of the film. As a result, this research demonstrates how 3D characterization offers a better understanding of DSA processes, leading to better design and fabrication of directing templates.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The GBS code for tokamak scrape-off layer simulations

We describe a new version of GBS, a 3D global, flux-driven plasma turbulence code to simulate the turbulent dynamics in the tokamak scrape-off layer (SOL), superseding the code presented by Ricci et al. (2012) [14]. The present work is driven by the objective of studying SOL turbulent dynamics in medium size tokamaks and beyond with a high-fidelity physics model. We emphasize an intertwining framework of improved physics models and the computational improvements that allow them. The model extensions include neutral atom physics, finite ion temperature, the addition of a closed field line region, and a non-Boussinesq treatment of the polarization drift. GBS has been completely refactored with the introduction of a 3-D Cartesian communicator and a scalable parallel multigrid solver. We report dramatically enhanced parallel scalability, with the possibility of treating electromagnetic fluctuations very efficiently. The method of manufactured solutions as a verification process has been carried out for this new code version, demonstrating the correct implementation of the physical model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

Virtual tool mark generation for efficient striation analysis in forensic science

In 2009, a National Academy of Sciences report called for investigation into the scienti c basis behind tool mark comparisons (National Academy of Sciences, 2009). Answering this call, Chumbley et al. (2010) attempted to prove or disprove the hypothesis that tool marks are unique to a single tool. They developed a statistical algorithm that could, in most cases, discern matching and non-matching tool marks made at di erent angles by sequentially numbered screwdriver tips. Moreover, in the cases where the algorithm misinterpreted a pair of marks, an experienced forensics examiner could discern the correct outcome. While this research served to con rm the basic assumptions behind tool mark analysis, it also suggested that statistical analysis software could help to reduce the examiner's workload. This led to a new tool mark analysis approach, introduced in this thesis, that relies on 3D scans of screwdriver tip and marked plate surfaces at the micrometer scale from an optical microscope. These scans are carefully cleaned to remove noise from the data acquisition process and assigned a coordinate system that mathematically de nes angles and twists in a natural way. The marking process is then simulated by using a 3D graphics software package to impart rotations to the tip and take the projection of the tip's geometry in the direction of tool travel. The edge of this projection, retrieved from the 3D graphics software, becomes a virtual tool mark. Using this method, virtual marks are made at increments of 5 and compared to a scan of the evidence mark. The previously developed statistical package from Chumbley et al. (2010) performs the comparison, comparing the similarity of the geometry of both marks to the similarity that would occur due to random chance. The resulting statistical measure of the likelihood of the match informs the examiner of the angle of the best matching virtual mark, allowing the examiner to focus his/her mark analysis on a smaller range of angles. Preliminary results are quite promising. In a study with both sides of 6 screwdriver tips and 34 corresponding marks, the method distinguished known matches from known non-matches with zero false positive matches and only two matches mistaken for non-matches. For matches, it could predict the correct marking angle within 5-10 . Moreover, on a standard desktop computer, the virtual marking software is capable of cleaning 3D tip and plate scans in minutes and producing a virtual mark and comparing it to a real mark in seconds. These results support several of the professional conclusions of the tool mark analysis com- munity, including the idea that marks produced by the same tool only match if they are made at similar angles. The method also displays the potential to automate part of the comparison process, freeing the examiner to focus on other tasks, which is important in busy, backlogged crime labs. Finally, the method o ers the unique chance to directly link an evidence mark to the tool that produced it while reducing potential damage to the evidence.

42 ENGINEERING↗
NASA NTRS1993

Experimental and numerical study of swept ramp injection into a supersonic flowfield

Time-averaged measurements of pressure, temperature, velocity, and injectant mole fraction are presented using the planar laser-induced iodine fluorescence (PLIIF) techniques in the complex 3D compressible flowfield around a swept ramp fuel injector. Within the range of thermodynamic conditions present in the test case studied, the technique's accuracy is estimated to be 4 percent for pressure, 4 percent for temperature, 4 percent for velocity, and 3 percent for injectant mole fraction. Comparisons to numerical simulations using the SPARK 3D Navier-Stokes computer code with an algebraic turbulence model are made at the centerplane of the flowfield as well as on three crossflow planes downstream of the injector, and good agreement is found. A weak asymmetry in the incoming flowfield appears to be amplified by boundary layer separation occuring when the the ramp generated shock reflects off the tunnel walls. In the near field of the injector, laminar calculations show better agreement to the measurements than turbulent calculations.

Donohue, James M.↗
DOE OSTI2016

WE-H-BRA-04: Biological Geometries for the Monte Carlo Simulation Toolkit TOPASNBio

Purpose: New advances in radiation therapy are most likely to come from the complex interface of physics, chemistry and biology. Computational simulations offer a powerful tool for quantitatively investigating radiation interactions with biological tissue and can thus help bridge the gap between physics and biology. The aim of TOPAS-nBio is to provide a comprehensive tool to generate advanced radiobiology simulations. Methods: TOPAS wraps and extends the Geant4 Monte Carlo (MC) simulation toolkit. TOPAS-nBio is an extension to TOPAS which utilizes the physics processes in Geant4-DNA to model biological damage from very low energy secondary electrons. Specialized cell, organelle and molecular geometries were designed for the toolkit. Results: TOPAS-nBio gives the user the capability of simulating biological geometries, ranging from the micron-scale (e.g. cells and organelles) to complex nano-scale geometries (e.g. DNA and proteins). The user interacts with TOPAS-nBio through easy-to-use input parameter files. For example, in a simple cell simulation the user can specify the cell type and size as well as the type, number and size of included organelles. For more detailed nuclear simulations, the user can specify chromosome territories containing chromatin fiber loops, the later comprised of nucleosomes on a double helix. The chromatin fibers can be arranged in simple rigid geometries or within factual globules, mimicking realistic chromosome territories. TOPAS-nBio also provides users with the capability of reading protein data bank 3D structural files to simulate radiation damage to proteins or nucleic acids e.g. histones or RNA. TOPAS-nBio has been validated by comparing results to other track structure simulation software and published experimental measurements. Conclusion: TOPAS-nBio provides users with a comprehensive MC simulation tool for radiobiological simulations, giving users without advanced programming skills the ability to design and run complex simulations.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2006

Functional brain imaging of a complex navigation task following one night of total sleep deprivation

Study Objectives: To assess the cerebral effects associated with sleep deprivation in a simulation of a complex, real-world, high-risk task. Design and Interventions: A two-week, repeated measures, cross-over experimental protocol, with counterbalanced orders of normal sleep (NS) and total sleep deprivation (TSD). Setting: Each subject underwent functional magnetic resonance imaging (fMRI) while performing a dual-joystick, 3D sensorimotor navigation task (simulated orbital docking). Scanning was performed twice per subject, once following a night of normal sleep (NS), and once following a single night of total sleep deprivation (TSD). Five runs (eight 24s docking trials each) were performed during each scanning session. Participants: Six healthy, young, right-handed volunteers (2 women; mean age 20) participated. Measurements and Results: Behavioral performance on multiple measures was comparable in the two sleep conditions. Neuroimaging results within sleep conditions revealed similar locations of peak activity for NS and TSD, including left sensorimotor cortex, left precuneus (BA 7), and right visual areas (BA 18/19). However, cerebral activation following TSD was substantially larger and exhibited higher amplitude modulations from baseline. When directly comparing NS and TSD, most regions exhibited TSD>NS activity, including multiple prefrontal cortical areas (BA 8/9,44/45,47), lateral parieto-occipital areas (BA 19/39, 40), superior temporal cortex (BA 22), and bilateral thalamus and amygdala. Only left parietal cortex (BA 7) demonstrated NS>TSD activity. Conclusions: The large network of cerebral differences between the two conditions, even with comparable behavioral performance, suggests the possibility of detecting TSD-induced stress via functional brain imaging techniques on complex tasks before stress-induced failures.

Strangman, Gary↗
DOE OSTI2018

Prediction of residual stresses in a multipass pipe weld by a novel 3D finite element approach

Due to enormous computation cost, current residual stress simulation of multipass girth welds are mostly performed using two-dimensional (2D) axisymmetric models. The 2D model can only provide limited estimation on the residual stresses by assuming its axisymmetric distribution. In this study, a highly efficient thermal-mechanical finite element code for three-dimensional (3D) model has been developed based on high performance Graphics Processing Unit (GPU) computers. Our code is further accelerated by considering the unique physics associated with welding processes that are characterized by steep temperature gradient and a moving arc heat source. It is capable of modeling large-scale welding problems that cannot be easily handled by the existing commercial simulation tools. To demonstrate the accuracy and efficiency, our code was compared with a commercial software by simulating a 3D multi-pass girth weld model with over 1 million elements. Our code achieved comparable solution accuracy with respect to the commercial one but with over 100 times saving on computational cost. Moreover, the three-dimensional analysis demonstrated more realistic stress distribution that is not axisymmetric in hoop direction.

Huang, Hui↗
NASA NTRS1991

Investigation of vortex development on a pitching slender body of revolution

A computational study of the unsteady flow about a pitching 3.5 caliber tangent ogive forebody is presented. The flow is simulated using the full 3D unsteady Navier-Stokes equations and a time-accurate implicit algorithm. Comparison to available experimental data for a steady 20 deg case is presented as validation. Effects of grid resolution and a comparison of solutions using full Navier-Stokes and the thin-layer approximation are included. The forebody is simulated in a 'pitch-up to 20 deg and hold' maneuver, and two different pitch axis locations are used in the study. Examination of the unsteady vorticity field for the pitch-up cases reveals the formation of strong shear layers as the body decelerates, and their roll-up into vortical structures in a process similar to that observed in airfoil dynamic stall. Pronounced vortex/surface interactions are seen which produce multiple secondary separation regions, ejection of vorticity from the surface, and embedded regions of high suction.

Stanek, M. J.↗
DOE OSTI2019

Ensemble Shear Strength, Stability, and Permeability of Mixed Mineralogy Fault Gouge Recovered From 3D Granular Models

We conduct numerical shear experiments on mixtures of quartz and talc gouge using a three-dimensional (3D) distinct element model. A modified slip-weakening constitutive law is applied at contacts. We perform velocity-stepping experiments on both uniform and layered mixtures of quartz and talc analogs. We separately vary the proportion of talc in the uniform mixtures and talc layer thickness in the layered mixtures. Shear displacements are cycled through velocities of 1 and 10 μm/s. We follow the resulting evolution of ensemble shear strength, slip stability, and permeability of the gouge mixture and explore the mesoscopic mechanisms. Simulation results show that talc has a strong weakening effect on shear strength—a thin shear-parallel layer of talc (three particles wide) can induce significant weakening. However, the model offsets laboratory-derived strong weakening effects of talc observed in uniform mixtures, implying the governing mechanisms may be the shear localization effect of talc, which is enhanced by its natural platy shape or preimposed layered structure. Ensemble stability (a – b ) can be enhanced by increasing talc content in uniform talc-quartz mixtures. Reactivation-induced permeability increase is amplified with increased quartz content before the maturation of shear localization. Postmaturation permeability enhances on velocity upsteps and diminishes on velocity downsteps. Talc enhances compaction at velocity downsteps, potentially reducing fault permeability. Evolution trends of stability relating to the composition and structure of the fault gouge are straightforwardly obtained from the 3D simulation. Local friction evolution indicates that talc preferentially organizes and localizes in the shear zone, dominating the shear strength and frictional stability of faults

58 GEOSCIENCES↗
DOE OSTI2022

Disc tearing leads to low and high frequency quasi-periodic oscillations in a GRMHD simulation of a thin accretion disc

Black hole X-ray binaries (BHXRBs) display a wide range of variability phenomena, from long duration spectral state changes to short-term broad-band variability and quasi-periodic oscillations (QPOs). A particularly puzzling aspect is the production of QPOs, which – if properly understood – could be used as a powerful diagnostic tool of black hole accretion and evolution. In this work, we analyse a high-resolution 3D general relativistic magnetohydrodynamic simulation of a geometrically thin accretion disc, which is tilted by 65° with respect to the black hole spin axis. We find that the Lense–Thirring torque from the rapidly spinning 10 M ⊙ black hole causes several sub-discs to tear off within ~10–20 gravitational radii. Tearing occurs in cycles on time-scales of seconds. During each tearing cycle, the inner sub-disc precesses for 1–5 periods before it falls into the black hole. We find a precession frequency of ~3 Hz, consistent with observed low-frequency QPOs. In addition, we find a high frequency QPO (HFQPO) with centroid frequency of ~55 Hz in the power spectra of the mass-weighted radius of the inner disc and the radial mass flux. This signal is caused by radial epicyclic oscillations of a dense ring of gas at the tearing radius, which suggests a corresponding modulation of the X-ray light curve and may thus explain some of the observed HFQPOs.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRSDate not supplied

Progressive Damage Response of 3D Woven Composites via the Multiscale Recursive Micromechanics Solution with Tailored Fidelity

Progressive failure simulations have been performed for orthogonal 3D woven composites consisting of RTM6 resin matrix and AS4 carbon fibers. The Multiscale Recursive Micromechanics approach has been used, which, while being computationally efficient, captures the primary effects of the microstructure at each considered length scale. This approach also enables use of any micromechanics theory at any length scale, and herein, the fidelity of the chosen theories across the scales has been tailored to strike a balance with computational efficiency. The Mori-Tanaka method is employed at the lowest length scale, the Generalized Method of Cells is used at intermediate scales, and the High-Fidelity Generalized Method of Cells is used at the highest woven composite repeating unit cell scale. Furthermore, two different damage models, also with different levels of fidelity and efficiency, have been used for the resin material at the lowest length scale. Results for the mechanical behavior in response to loading in various directions are compared for the two damage models and with available test data.

3D woven↗
NASA NTRS2014

Numerical Modeling of Ocular Dysfunction in Space

Upon introduction to microgravity, the near-loss of hydrostatic pressure causes a marked cephalic (headward) shift of fluid in an astronaut's body. The fluid shift, along with other factors of spaceflight, induces a cascade of interdependent physiological responses which occur at varying time scales. Long-duration missions carry an increased risk for the development of the Visual Impairment and Intracranial Pressure (VIIP) syndrome, a spectrum of ophthalmic changes including posterior globe flattening, choroidal folds, distension of the optic nerve sheath, kinking of the optic nerve and potentially permanent degradation of visual function. In the cases of VIIP found to date, the initial onset of symptoms occurred after several weeks to several months of spaceflight, by which time the gross bodily fluid distribution is well established. We are developing a suite of numerical models to simulate the effects of fluid shift on the cardiovascular, central nervous and ocular systems. These models calculate the modified mean volumes, flow rates and pressures that are characteristic of the altered quasi-homeostatic state in microgravity, including intracranial and intraocular pressures. The results of the lumped models provide initial and boundary data to a 3D finite element biomechanics simulation of the globe, optic nerve head and retrobulbar subarachnoid space. The integrated set of models will be used to investigate the evolution of the biomechanical stress state in the ocular tissues due to long-term exposure to microgravity.

intraocular pressure↗
DOE OSTI2024

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗