Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulator Study”

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 829 records · Page 46

Alpha image reconstruction (AIR): A new iterative CT image reconstruction approach using voxel-wise alpha blending

Purpose: Iterative image reconstruction gains more and more interest in clinical routine, as it promises to reduce image noise (and thereby patient dose), to reduce artifacts, or to improve spatial resolution. Among vendors and researchers, however, there is no consensus of how to best achieve these aims. The general approach is to incorporatea priori knowledge into iterative image reconstruction, for example, by adding additional constraints to the cost function, which penalize variations between neighboring voxels. However, this approach to regularization in general poses a resolution noise trade-off because the stronger the regularization, and thus the noise reduction, the stronger the loss of spatial resolution and thus loss of anatomical detail. The authors propose a method which tries to improve this trade-off. The proposed reconstruction algorithm is called alpha image reconstruction (AIR). One starts with generating basis images, which emphasize certain desired image properties, like high resolution or low noise. The AIR algorithm reconstructs voxel-specific weighting coefficients that are applied to combine the basis images. By combining the desired properties of each basis image, one can generate an image with lower noise and maintained high contrast resolution thus improving the resolution noise trade-off. Methods: All simulations and reconstructions are performed in native fan-beam geometry. A water phantom with resolution bar patterns and low contrast disks is simulated. A filtered backprojection (FBP) reconstruction with a Ram-Lak kernel is used as a reference reconstruction. The results of AIR are compared against the FBP results and against a penalized weighted least squares reconstruction which uses total variation as regularization. The simulations are based on the geometry of the Siemens Somatom Definition Flash scanner. To quantitatively assess image quality, the authors analyze line profiles through resolution patterns to define a contrast factor for contrast-resolution plots. Furthermore, the authors calculate the contrast-to-noise ratio with the low contrast disks and the authors compare the agreement of the reconstructions with the ground truth by calculating the normalized cross-correlation and the root-mean-square deviation. To evaluate the clinical performance of the proposed method, the authors reconstruct patient data acquired with a Somatom Definition Flash dual source CT scanner (Siemens Healthcare, Forchheim, Germany). Results: The results of the simulation study show that among the compared algorithms AIR achieves the highest resolution and the highest agreement with the ground truth. Compared to the reference FBP reconstruction AIR is able to reduce the relative pixel noise by up to 50% and at the same time achieve a higher resolution by maintaining the edge information from the basis images. These results can be confirmed with the patient data. Conclusions: To evaluate the AIR algorithm simulated and measured patient data of a state-of-the-art clinical CT system were processed. It is shown, that generating CT images through the reconstruction of weighting coefficients has the potential to improve the resolution noise trade-off and thus to improve the dose usage in clinical CT.

60 APPLIED LIFE SCIENCES↗

Study on Conversion from an Existing Alternating Current Medium-Voltage Line to a Direct Current Line: A Case of Load or Source Expansion

Simultaneous increase of power generation through renewable systems and load due owing to the rapid growth in electrified road transportation systems, leads to new challenges in the design, operation, and protection of distribution-level electrical power systems. This paper addresses challenges related to the integration of new sources or loads to a distribution grid through long alternating current lines. The International Council on Large Electric Systems medium-voltage benchmark model developed for North America was used for simulation studies. Use cases were developed to illuminate the issues related to long alternating current lines and the usefulness of medium-voltage direct current lines in a distribution network.

Jaldanki, Sreenivasa↗

Patient motion effects on the quantification of regional myocardial blood flow with dynamic PET imaging

Purpose: Patient motion is a common problem during dynamic positron emission tomography (PET) scans for quantification of myocardial blood flow (MBF). The purpose of this study was to quantify the prevalence of body motion in a clinical setting and evaluate with realistic phantoms the effects of motion on blood flow quantification, including CT attenuation correction (CTAC) artifacts that result from PET–CT misalignment. Methods: A cohort of 236 sequential patients was analyzed for patient motion under resting and peak stress conditions by two independent observers. The presence of motion, affected time-frames, and direction of motion was recorded; discrepancy between observers was resolved by consensus review. Based on these results, patient body motion effects on MBF quantification were characterized using the digital NURBS-based cardiac-torso phantom, with characteristic time activity curves (TACs) assigned to the heart wall (myocardium) and blood regions. Simulated projection data were corrected for attenuation and reconstructed using filtered back-projection. All simulations were performed without noise added, and a single CT image was used for attenuation correction and aligned to the early- or late-frame PET images. Results: In the patient cohort, mild motion of 0.5 ± 0.1 cm occurred in 24% and moderate motion of 1.0 ± 0.3 cm occurred in 38% of patients. Motion in the superior/inferior direction accounted for 45% of all detected motion, with 30% in the superior direction. Anterior/posterior motion was predominant (29%) in the posterior direction. Left/right motion occurred in 24% of cases, with similar proportions in the left and right directions. Computer simulation studies indicated that errors in MBF can approach 500% for scans with severe patient motion (up to 2 cm). The largest errors occurred when the heart wall was shifted left toward the adjacent lung region, resulting in a severe undercorrection for attenuation of the heart wall. Simulations also indicated that the magnitude of MBF errors resulting from motion in the superior/inferior and anterior/posterior directions was similar (up to 250%). Body motion effects were more detrimental for higher resolution PET imaging (2 vs 10 mm full-width at half-maximum), and for motion occurring during the mid-to-late time-frames. Motion correction of the reconstructed dynamic image series resulted in significant reduction in MBF errors, but did not account for the residual PET–CTAC misalignment artifacts. MBF bias was reduced further using global partial-volume correction, and using dynamic alignment of the PET projection data to the CT scan for accurate attenuation correction during image reconstruction. Conclusions: Patient body motion can produce MBF estimation errors up to 500%. To reduce these errors, new motion correction algorithms must be effective in identifying motion in the left/right direction, and in the mid-to-late time-frames, since these conditions produce the largest errors in MBF, particularly for high resolution PET imaging. Ideally, motion correction should be done before or during image reconstruction to eliminate PET-CTAC misalignment artifacts.

60 APPLIED LIFE SCIENCES↗

A Design Study of the Parallel-Slit Ring Collimator for Fast Neutron Emission Tomography of Spent Fuel

For the past 18 months, Oak Ridge National Laboratory has been developing a passive fast-neutron emission tomography capability. The goal of this development is a new passive capability to detect diversion of single fuel pins (rods) from spent or used nuclear fuel assemblies by performing neutron emission tomography. At present, a laboratory prototype imager is under construction. The purpose of the prototype imager is to demonstrate imaging capability sufficient to resolve individual fuel pins in a spent nuclear fuel assembly. This report documents the results of simulation studies to optimize the design of the prototype imager.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Fast Responding Voltage Regulator and Dynamic VAR Compensator

The objectives of this project were to develop a dynamic VAR compensator (DVC) for voltage regulation through VAR support to demonstrate the ability to achieve greater levels of voltage control on electricity distribution networks, and faster response compared to existing grid technology. The goal of the project was to develop a prototype Fast Dynamic VAR Compensator (Fast DVC) hardware device, and this was achieved. In addition to developing the dynamic VAR compensator device, Varentec in partnership with researchers at North Carolina State University (NCSU) successfully met the objectives to model the potential positive impact of such DVCs on representative power networks. This modeling activity validated the ability of distributed dynamic VAR compensators to provide fast voltage regulation and reactive power control required to respond to grid disturbances under high penetration of fluctuating and intermittent distributed energy resources (DERs) through extensive simulation studies. Specifically the following tasks were set to be accomplished: 1) Development of dynamic VAR compensator to support dynamic voltage variations on the grid through VAR control 2) Extensive testing of the DVC in the lab environment 3) Present the operational DVC device to the DOE at Varentec’s lab 4) Formulation of a detailed specification sheet, unit assembly document, test setup document, unit bring-up plan, and test plan 5) Extensive simulations of the DVC in a system with high PV penetration. Understanding the operation with many DVC on a single distribution system 6) Creation and submittal of quarterly and final reports conveying the design documents, unit performance data, modeling simulation charts and diagrams, and summary explanations of the satisfaction of program goals. This report details the various efforts that led to the development of the Fast DVC as well as the modeling & simulation results. The report begins with the introduction in Section II which outlines the problems associated with increasing penetration of renewable resources on the grid and what are the challenges in solving these issues. Section III provides a brief summary of the initial concepts that were explored before the final Fast DVC prototype was developed. Section IV and V present the hardware, controls, communication, and mechanical architecture and features of the Fast DVC. Section VI presents the procedure for assembling the Fast DVC. The test plan and lab setup are expounded in Section VII and VIII. Section IX presents the results from testing the Fast DVC in the lab, several results are presented along with some pictures of the developed unit. Section X details the pilot demonstration activities. Finally, section XI showcases the efforts associated with the simulation and modeling of the Fast DVC in various scenarios with DERs and its impact in all these scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Impact of bowtie filter and object position on the two-dimensional noise power spectrum of a clinical MDCT system

Purpose: Noise characteristics of clinical multidetector CT (MDCT) systems can be quantified by the noise power spectrum (NPS). Although the NPS of CT has been extensively studied in the past few decades, the joint impact of the bowtie filter and object position on the NPS has not been systematically investigated. This work studies the interplay of these two factors on the two dimensional (2D) local NPS of a clinical CT system that uses the filtered backprojection algorithm for image reconstruction. Methods: A generalized NPS model was developed to account for the impact of the bowtie filter and image object location in the scan field-of-view (SFOV). For a given bowtie filter, image object, and its location in the SFOV, the shape and rotational symmetries of the 2D local NPS were directly computed from the NPS model without going through the image reconstruction process. The obtained NPS was then compared with the measured NPSs from the reconstructed noise-only CT images in both numerical phantom simulation studies and experimental phantom studies using a clinical MDCT scanner. The shape and the associated symmetry of the 2D NPS were classified by borrowing the well-known atomic spectral symbols s, p, and d, which correspond to circular, dumbbell, and cloverleaf symmetries, respectively, of the wave function of electrons in an atom. Finally, simulated bar patterns were embedded into experimentally acquired noise backgrounds to demonstrate the impact of different NPS symmetries on the visual perception of the object. Results: (1) For a central region in a centered cylindrical object, an s-wave symmetry was always present in the NPS, no matter whether the bowtie filter was present or not. In contrast, for a peripheral region in a centered object, the symmetry of its NPS was highly dependent on the bowtie filter, and both p-wave symmetry and d-wave symmetry were observed in the NPS. (2) For a centered region-ofinterest (ROI) in an off-centered object, the symmetry of its NPS was found to be different from that of a peripheral ROI in the centered object, even when the physical positions of the two ROIs relative to the isocenter were the same. (3) The potential clinical impact of the highly anisotropic NPS, caused by the interplay of the bowtie filter and position of the image object, was highlighted in images of specific bar patterns oriented at different angles. The visual perception of the bar patterns was found to be strongly dependent on their orientation. Conclusions: The NPS of CT depends strongly on the bowtie filter and object position. Even if the location of the ROI with respect to the isocenter is fixed, there can be different symmetries in the NPS, which depend on the object position and the size of the bowtie filter. For an isolated off-centered object, the NPS of its CT images cannot be represented by the NPS measured from a centered object.

60 APPLIED LIFE SCIENCES↗

Chain and Segmental Dynamics in Polymer–Nanoparticle Composites with High Nanoparticle Loading

Recent experimental and simulation studies have shown that polymer–nanoparticle (NP) composites (PNCs) with ultra-high NP loading (>50%) exhibit remarkable mechanical properties and dramatic increases in polymer glass-transition temperature, viscosity, and thermal stability compared to the bulk polymer. These deviations in macroscopic properties suggest a slowdown in both segmental and chain-scale polymer dynamics due to confinement. In this work, we examine the polymer conformations and dynamics in these PNCs using molecular dynamics simulations of both unentangled and entangled coarse-grained polymers in random-close-packed NP packings with varying polymer fill fractions. We find that the changes in the polymer dynamics depend on the number of NPs in contact with a polymer segment. Using the number of polymer–NP contacts and different polymer chain conformations as criteria for categorization, we further examine the polymer dynamics at multiple length scales to show the high level of dynamic heterogeneity in PNCs with ultra-high NP loading.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A High-Fidelity Electromagnetic Transient Model of Inverter-based Resources Integrated to an IEEE-9 Bus System for Benchmarking Studies

An inverter-based resource (IBR) is a source of electricity that is asynchronously connected to the electrical grid via an electronic power converter. These power sources lack the intrinsic behaviors of the standard power plants, presenting specific challenges to system stability. These power sources lack the intrinsic behaviors of the standards power plants and their features are almost entirely defined by the control strategy used on them, presenting specific challenges to system stability as their penetration increases into the already fragile Bulk Power System (BPS). Utilizing renewable energy has a great upside but the electrical industry needs to understand the requirements like performing electromagnetic transient (EMT) simulation studies in planning and/or in interconnection studies. This paper presents the development of a high-fidelity EMT model based on commercial equipment and field parameters, which at the same time aims to study the increasing penetration of IBRs and their impact on the actual BPS. The open-source EMT model developed will be a benchmarking example that can be used to evaluate algorithms in the power grid (including evaluating grid-forming controls, power system monitoring and operation, etc.).

Martinez Montejano, Misael↗

Monte Carlo computer simulation of sedimentation of charged hard spherocylinders

In this article we present a NVT Monte Carlo computer simulation study of sedimentation of an electroneutral mixture of oppositely charged hard spherocylinders (CHSC) with aspect ratio L/σ = 5, where L and σ are the length and diameter of the cylinder and hemispherical caps, respectively, for each particle. This system is an extension of the restricted primitive model for spherical particles, where L/σ = 0, and it is assumed that the ions are immersed in an structureless solvent, i.e., a continuum with dielectric constant D. The system consisted of N = 2000 particles and the Wolf method was implemented to handle the coulombic interactions of the inhomogeneous system. Results are presented for different values of the strength ratio between the gravitational and electrostatic interactions, Γ = (mgσ)/(e{sup 2}/Dσ), where m is the mass per particle, e is the electron's charge and g is the gravitational acceleration value. A semi-infinite simulation cell was used with dimensions L{sub x} ≈ L{sub y} and L{sub z} = 5L{sub x}, where L{sub x}, L{sub y}, and L{sub z} are the box dimensions in Cartesian coordinates, and the gravitational force acts along the z-direction. Sedimentation effects were studied by looking at every layer formed by the CHSC along the gravitational field. By increasing Γ, particles tend to get more packed at each layer and to arrange in local domains with an orientational ordering along two perpendicular axis, a feature not observed in the uncharged system with the same hard-body geometry. This type of arrangement, known as tetratic phase, has been observed in two-dimensional systems of hard-rectangles and rounded hard-squares. In this way, the coupling of gravitational and electric interactions in the CHSC system induces the arrangement of particles in layers, with the formation of quasi-two dimensional tetratic phases near the surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of surface electric charge profile in pyroelectric crystals

Pyroelectric crystals are used to produce high energy self-focused electron beams. Here, an experimental analysis in combination with simulation studies will be reported to investigate possible sources of this effect. In the experiments, the surface of crystal was divided into six separated parts and the rate of surface electric charge production was measured accordingly. A non-steady and spatially non-uniform distribution of the surface charge generation was observed, in which it tends to a uniform distribution in the course of experiment. The obtained surface electric charges from the experiments were used to simulate the electric field and potential around the crystal by COMSOL Multiphysics. It was observed that emitted electrons from the crystal surface were focused, and the non-uniformity in spatial charge is responsible for this phenomenon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A study on dc fault handling of dc-dc converter based on series and parallel connected DABs in MVdc applications

With the increase of direct current (dc) based power generating stations and loads, multi-terminal medium voltage dc (MT-MVdc) systems are gaining popularity. The dc sources and loads are interfaced with the MT-MVdc systems using series and parallel connected dual-active-bridge (SP-DAB). In this paper, the dc fault ride through procedure of SP-DAB meant for MT-MVdc application is analyzed. The need for dc breaker (DCB) under the requirements of MT-MVdc is established through EMT simulation studies. It is shown that an additional dc filter capacitor is essential at the output terminals of the SP-DAB to ride through the dc faults.

Jaldanki, Sreenivasa [ORNL] (ORCID:000000028122322↗

Symmetry-Breaking and Self-Sorting in Block Copolymer-Based Multicomponent Nanocomposites

Co-assembly of inorganic nanoparticles (NPs) and nanostructured polymer matrix represents an intricate interplay of enthalpic or entropic forces. Particle size largely affects the phase behavior of the nanocomposite. Theoretical studies indicate that new morphologies would emerge when the particles become comparable to the soft matrix's size, but this has rarely been supported experimentally. By designing a multicomponent blend composed of NPs, block copolymer-based supramolecules, and small molecules, a 3-D ordered lattice beyond the native BCP's morphology was recently reported when the particle is larger than the microdomain of BCP. The blend can accommodate various formulation variables. In this work, when the particle size equals the microdomain size, a symmetry-broken phase appears in a narrow range of particle sizes and compositions, which we named the "train track" structure. In this phase, the NPs aligned into a 3-D hexagonal lattice and packed asymmetrically along the c axis, making the projection of the ac and the bc plane resemble train tracks. Computational studies show that the broken symmetry reduces the polymer chain deformation and stabilizes the metastable hexagonally perforated lamellar morphology. Given the mobility of the multicomponent blend, the system shows a self-sorting behavior: segregating into two macroscopic phases with different nanostructures based on only a few nanometers NP size differences. Smaller NPs form "train track" morphology, while larger NPs form a "simple hexagon" structure, where the NPs take a symmetric hexagonal arrangement. Detailed structural evolution and simulation studies confirm the systematic-wide cooperativity across different components, indicating the strong self-regulation of the multicomponent system.

"train track" morphology↗

Review and meta-analysis of electron temperatures from high-intensity laser–solid interactions

The accelerated electron spectrum from high-intensity laser–solid interaction is often conveniently described using a Boltzmann distribution, whose temperature is known within the field as the hot-electron temperature. The importance of the electron temperature is highlighted by the sheer number of experimental and simulation studies on the subject over the past three decades. Recently, multi-kJ, multi-ps pulses have yielded electron spectra with temperatures far beyond the expected ponderomotive result. Expressions that predict the electron temperature considering laser parameters beyond intensity and wavelength have been developed, albeit using small datasets. In this review, we present what is, to the best of our knowledge, the largest dataset of electron temperatures gathered from experimental measurements and particle-in-cell simulations. This dataset allows us to compare existing analytical and empirical hot-electron temperature scaling models over a wide parameter range. We also develop new scaling models that incorporate the laser pulse duration of the laser and the plasma scale length. Three models that include pulse-duration and scale length dependence are especially successful at predicting both simulated and experimental data. The dataset will soon be made publicly available to encourage further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pressure Transient Analyses and Poroelastic Modeling of Hydraulic Fracture Dilation for Multiple Injections at the Devine Fracture Pilot Site

Our team has conducted electromagnetic (EM) surveys for the past six years to monitor hydraulic-fracture behavior at the Devine Fracture Pilot Site (DFPS). The sub-horizontal orientation of a shallow hydraulic fracture at the DFPS provides uniform access to the fracture area for interrogation and data collection. Ahmadian et al. (2023) suggested a possible correlation between spatiotemporal changes in the flow rate, bottomhole pressure (BHP), and the observed surface recorded electric field at the DFPS. In this paper, we present the development of poroelastic forward models and pressure transient analyses (PTAs) to support the development of a multiphysics inverse model for these EM surveys. First, we conducted PTAs of the shut-in periods after six injections out of 10 to determine the fracture closure pressure (FCP) or the overburden pressure used in a poroelastic fracture reopening model. Second, we developed a finite-element poroelastic model throughout five injection cycles to include the effect of the cumulative injected volumes due to the previous injections on current fracture dilation in the presence of highly permeable unpropped and propped zones adjacent to the cohesive layer that models fracture reopening. Fracture reopening in this poroelastic model is based on a calibrated traction-separation response using the bottomhole pressure collected in two injection campaigns in 2020 and 2022. We used the outcomes of a previous simulation study of the primary hydraulic-fracturing stimulation to define the dimension of an unpropped fracture zone ahead of the propped fracture area. The PTAs led to FCPs consistent with those obtained using the injection data collected at the DFPS in 2020. Further, these analyses showed that at later injections, the fracture closure occurred at a later time with respect to the shut-in time, inferring the effect of cumulative injected volumes in previous injections. The simulation results show that considering the propped and unpropped fracture zones improves our poroelastic model in predicting the injection-well BHP. The numerical simulation results demonstrate a significant excess pore pressure near the fracture because of the preceding formation loadings by the previous injections. The obtained fracture dilation area and fluid pressure distribution provide a basis to improve the development of a multiphysics inverse model. Furthermore, in an iteratively coupled scheme, this pressure distribution can be introduced into EM models to render a holistic view of the causative mechanisms for the surface signal anomalies.

02 PETROLEUM↗

Optimization of electrical treatment strategy for surface roughness reduction in conducting thin films

The surface roughness of deposited conducting thin films is responsible for various materials reliability problems in nanoelectronics and nanofabrication technologies. In this work, we report a modeling and simulation study that aims at optimizing the electrical surface treatment of deposited conducting thin films as a physical processing strategy for their surface roughness reduction. Our research is based on a continuum model of film surface morphological evolution that accounts for the residual stress in the deposited conducting thin film, the film’s wetting of the substrate layer that it is deposited on, film texture and surface diffusional anisotropy, and surface electromigration. Through systematic linear stability analysis and dynamical simulation protocols, we examine in detail the effects of film surface crystallographic orientation and applied electric field direction toward minimizing the electric field strength required for film surface smoothening. We believe that the critical electric field strength requirement for surface roughness reduction on {110}, {100}, and {111} surfaces of face-centered cubic crystalline conducting thin films exhibits a very strong dependence on the applied electric field direction, expressed as the electric field misalignment with respect to the principal residual stress directions in the film and the fast surface diffusion directions. Due to these findings, we optimize the electrical treatment strategy for surface roughness reduction of conducting thin films with respect to all relevant processing and material parameters.

42 ENGINEERING↗

Time projection chamber for GADGET II

The established Gaseous Detector with Germanium Tagging (GADGET) detection system is used to measure weak, low-energy 𝛽-delayed proton decays. It consists of the Gaseous Proton Detector equipped with a MICROMEGAS (MM) readout to detect protons and other charged particles calorimetrically, surrounded by the Segmented Germanium Array (SeGA) for high-resolution detection of prompt 𝛾 rays. To upgrade GADGET's Proton Detector to operate as a compact time projection chamber (TPC) for the detection, three-dimensional imaging and identification of low-energy 𝛽-delayed single- and multiparticle emissions mainly of interest to astrophysical studies. A new high granularity MM board with 1024 pads has been designed, fabricated, installed, and tested. A high-density data acquisition system based on generic electronics for TPCs (GET) has been installed and optimized to record and process the gas avalanche signals collected on the readout pads. The TPC's performance has been tested using a 220 Rn 𝛼-particle source and cosmic-ray muons. In addition, decay events in the TPC have been simulated by adapting the attpcroot data analysis framework. Furthermore, a novel application of two-dimensional convolutional neural networks for GADGET II event classification is introduced. The optimization of data throughput is also addressed. The GADGET II TPC is capable of detecting and identifying 𝛼 particles as well as measuring their track direction, range, and energy. The extracted energy resolution of the GADGET II TPC using P10 gas is about 5.4% at 6.288 MeV ( 220 Rn 𝛼 events), computed using charge integration. Based on a systematic simulation study, we estimated the detection efficiency of the GADGET II TPC for protons and 𝛼 particles, respectively. It has also been demonstrated that the GADGET II TPC is capable of tracking minimum-ionizing particles (i.e., cosmic-ray muons). From these measurements, the electron drift velocity was measured under typical operating conditions. In addition to being one of the first generation of micropattern gaseous detectors (MPGDs) to utilize a resistive anode applied to low-energy nuclear physics, the GADGET II TPC will also be the first TPC surrounded by a high-efficiency array of high-purity germanium 𝛾-ray detectors. As a result, the TPC of GADGET II has been designed, fabricated, and tested and is ready for operation at the Facility for Rare Isotope Beams for radioactive-beam-line experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The importance of GR for the radius of massive white dwarfs

In this work we investigate the structure of WD stars using the Tolman-Oppenheimer-Volkoff equations and compare with the Newtonian equations of gravitation in order to put in evidence the importance of the General Relativity in the study of these stars. Instead of using politropic equations of state we have solved the equations using the exact relativistic energy equation for the model of completely degenerate electron gas. We find a good fit of the TOV solution with the general EoS for the WD mass-radius diagram. We propose that our fit has to be used as relation between mass and radius for general relativistic WD instead of that Newtonian M ∼ 1/R{sup 3}, this fit is given by M = R/(a + bR + cR{sup 2} + dR{sup 3} + kR{sup 4}), where a, b, c and d are parameters and 1/k is the constant of the Newtonian mass-radius relation and it can be used in simulation study of binary systems that occurs accretion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling of plasma processes in the slowly diverging magnetic fields at the exit of an applied-field magnetoplasmadynamic thruster

The performance of plasma thrusters with applied electric and magnetic fields can be enhanced by increasing the magnetic field strength, which is applied in the thrust chamber and the exit region propulsive plume. The ejected plasma which passes through a slowly diverging magnetic field will expand but can be restricted within the magnetic nozzle fields. To examine in detail the processes that occur, a new method with Particle-in-cell calculations is applied here. A two-dimensional axisymmetric particle dynamic code is used to model an AF-MPDT (Applied-field Magnetoplasmadynamic Thruster) for which extensive experimental data are available; it used Ar propellant and had applied magnetic coils of 101.5 mm radius and 153 mm length. From the results of the simulation study, it is found that total thrust increases linearly with magnetic field strength in the range of 0–0.1 T, but it decreases with increasing applied magnetic field up to 0.6 T. Thrust efficiency is found to increase to a maximum of 8.4% when B = 0.1 T; further, the peak value of nozzle efficiency reaches 91% at a moderate magnetic field (0.3 T). In detail, it is found that distributions of plasma density (10{sup 14}–10{sup 15} m{sup −3}) that form in the magnetic nozzle demonstrate a significant pattern of concentration up to fields of B = 0.3 T where ions begin to be magnetized. However, azimuthal velocities of ions behave differently with different degrees of magnetization, i.e., weakly magnetized ions follow rotating electrons in a right-handed direction, while fully magnetized ions revolve in left-handed direction due to electromagnetic forces. Notably, a feedback effect on total magnetic field due to plasma motion identified in other studies is not found to be present in the working conditions of the AF-MPDT examined here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗