Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter studies”

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 253 records · Page 14

Target geometrical effects on the stagnation layer formed by colliding a pair of laser produced copper plasmas

We present the results of a time and space resolved optical-spectroscopic study of colliding plasmas formed at the front surfaces of flat and inclined Cu slab targets as a function of both the distance and the wedge angle between them for angles ranging from 100° to 180° (laterally colliding plasmas). The key parameters studied are stagnation layer density, temperature, duration, and kinetics of atomic/ionic spatial distributions and all have been found to vary significantly with wedge angle. It is found that the density and temperature of the stagnation layer decrease with increasing wedge angle. It is also found that the larger the wedge angle, the tighter and more well defined the stagnation layer formed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-F-T-55: Reproducibility of Interstitial HDR Brachytherapy Plans

Purpose: Treating gynecological cancers with interstitial high-dose-rate (HDR) brachytherapy requires precise reconstruction of catheter positions to obtain accurate dosimetric plans. In this study, we investigated the degree of reproducibility of dosimetric plans for Syed HDR brachytherapy. Methods: We randomly selected five patients having cervix-vaginal cancer who were recently treated in our clinic with interstitial HDR brachytherapy with a prescription dose of 25–30 Gy in five fractions. Interstitial needles/catheters were placed under fluoroscopic guidance and intra-operative 3T MRI scan was performed to confirm the desired catheter placement for adequate target volume coverage. A CT scan was performed and fused with the MRI for delineating high-risk CTV (HR-CTV), intermediate-risk CTV (IR-CTV) and OARs. HDR treatment plans were generated using Oncentra planning software. A single plan was used for all five fractions of treatment for each patient. For this study, we took the original clinical plan and removed all the reconstructed catheters from the plan keeping the original contours unchanged. Then, we manually reconstructed all the catheters and entered the same dwell time from the first original clinical plan. The dosimetric parameters studied were: D90 for HR-CTV and IR-CV, and D2cc for bladder, rectum, sigmoid and bowel. Results: The mean of absolute differences in dosimetric coverage (D90) were (range): 1.3% (1.0–2.0%) and 2.0% (0.9–3.6%) for HR-CTV and IR-CTV, respectively. In case of OARs, the mean of absolute variations in D2cc were (range): 4.7% (0.7–8.9%) for bladder, 1.60% (0.3–3.2%) for rectum, 1.6% (0–3.9%) for sigmoid, and 1.8% (0–5.1%) for bowel. Conclusion: Overall, the reproducibility of interstitial HDR plans was within clinically acceptable limit. Observed maximum variation in D2cc for bladder. If number of catchers and dwell points were relatively low or any one catheter was heavily loaded, then reproducibility of the plan was more sensitive to the accuracy of catheter reconstruction.

60 APPLIED LIFE SCIENCES↗

Sample variance in weak lensing: How many simulations are required?

Constraining cosmology using weak gravitational lensing consists of comparing a measured feature vector of dimension N b with its simulated counterpart. An accurate estimate of the N b × N b feature covariance matrix C is essential to obtain accurate parameter confidence intervals. When C is measured from a set of simulations, an important question is how large this set should be. To answer this question, we construct different ensembles of N r realizations of the shear field, using a common randomization procedure that recycles the outputs from a smaller number N s ≤ N r of independent ray-tracing N-body simulations. We study parameter confidence intervals as a function of (N s , N r ) in the range 1 ≤ N s ≤ 200 and 1 ≤ N r ≲ 105. Previous work [S. Dodelson and M. D. Schneider, Phys. Rev. D 88, 063537 (2013)] has shown that Gaussian noise in the feature vectors (from which the covariance is estimated) lead, at quadratic order, to an O(1/N r ) degradation of the parameter confidence intervals. Using a variety of lensing features measured in our simulations, including shear-shear power spectra and peak counts, we show that cubic and quartic covariance fluctuations lead to additional O(1/N 2 r ) error degradation that is not negligible when N r is only a factor of few larger than N b . We study the large N r limit, and find that a single, 240 Mpc/h sized 512 3 -particle N-body simulation (N s = 1) can be repeatedly recycled to produce as many as N r = few × 10 4 shear maps whose power spectra and high-significance peak counts can be treated as statistically independent. Lastly, a small number of simulations (N s = 1 or 2) is sufficient to forecast parameter confidence intervals at percent accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

ZND Verification Tests for Reactive Burn Models in FLAG

The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Fundamental Aeroelastic Analysis of an Urban Air Mobility Rotor

A computational investigation was conducted of the aeroelastic characteristics of a representative rotor typical of aircraft for Urban Air Mobility. A rotorcraft comprehensive analysis was used to examine the aeroelastic characteristics of the model, both in hover and in forward flight. Throughout the study, parameters examined included rotor rotational speed, forward flight speed, blade torsion stiffness, and chordwise locations of the blade elastic axis and center of gravity. Three sets of airfoil properties were analyzed: NACA 0012, VR8, and VR12. For hover, the effect of dynamic inflow was studied via a comparison of uniform inflow and dynamic inflow models. Thrust coefficient and blade tip twist were analyzed to ascertain performance characteristics, and eigen analysis was employed to determine stability. For the forward flight cases, rpm and forward flight speed were swept separately, and the effect of each on performance and stability was assessed. Thrust coefficient and mean drag coefficient were used to evaluate performance of the rotor, and eigen analysis was invoked to study rotor stability. Forward flight stability results generated via a constant coefficient analysis were compared with results generated via Floquet Theory.

Fundamental Aeroelastic Analysis↗

A study of eigenvalue sensitivity for hydrodynamic stability operators

The eigenvalue sensitivity for hydrodynamic stability operators is investigated. Classical matrix perturbation techniques as well as the concept of epsilon-pseudospectra are applied to show that parts of the spectrum are highly sensitive to small perturbations. Applications are drawn from incompressible plane Couette flow, trailing line vortex flow, and compressible Blasius boundary-layer flow. Parameter studies indicate a monotonically increasing effect of the Reynolds number on the sensitivity. The phenomenon of eigenvalue sensitivity is due to the nonnormality of the operators and their discrete matrix analogs and may be associated with large transient growth of the corresponding initial value problem.

Schmid, Peter J.↗

Polarized neutron scattering study of the multiple order parameter system NdB 4 [Polarized neutron scattering study on multiple order parameter system NdB 4 ]

Neutron polarization analysis has been carried out in order to clarify the magnetic structures of multiple order parameter f-electron system NdB 4 . We confirmed the noncollinear “all-in all-out” structure (Γ 4 ) of the in-plane moment, which is in good agreement with our previous neutron powder diffraction study. We found that the magnetic moment along the c-axis m c showed diagonally antiferromagnetic structure (Γ 10 ), inconsistent with previously reported “vortex” structure (Γ 2 ). The microscopic mixture of these two structures with → q 0 =(0,0,0) appears in phase II and remains stable in phases III and IV, where an incommensurate modulation coexists. The unusual magnetic ordering is phenomenologically understood via Landau theory with the primary order parameter Γ 4 coupled with higher-order secondary order parameter Γ 10 . The magnetic moments were estimated to be 1.8 ± 0.2 and 0.2 ± 0.05μ B at T = 7.5K for Γ 4 and Γ 10 , respectively. We also found a long-period incommensurate modulation of the → q 1 =(0,0,1/2) antiferromagnetic structure of mc with the propagation → q s1 =(0.14,0.14,0.1) and → q s2 =(0.2,0,0.1) in phase III and IV, respectively. The amplitude of sinusoidal modulation was about m c =1.0 ± 0.2μ B at T=1.5 K. The local (0,0,1/2) structure consists of in-plane ferromagnetic and out-of-plane antiferromagnetic coupling of m c , opposite to the coexisting Γ 10 . The mc of Γ 10 is significantly enhanced up to 0.6μ B at T=1.5 K, which is accompanied by the incommensurate modulations. As a result, the Landau phenomenological approach indicates that the higher-order magnetic and/or multipole interactions based on the pseudoquartet f-electron state play important roles.

36 MATERIALS SCIENCE↗

Momentum injection by supernovae in the interstellar medium

Supernova (SN) explosions deposit prodigious energy and momentum in their environments, with the former regulating multiphase thermal structure and the latter regulating turbulence and star formation rates in the interstellar medium (ISM). However, systematic studies quantifying the impact of SNe in realistic inhomogeneous ISM conditions have been lacking. Using three-dimensional hydrodynamic simulations, we investigate the dependence of radial momentum injection on both physical conditions (considering a range of mean density n{sub 0} = 0.1–100 cm{sup −3}) and numerical parameters. Our inhomogeneous simulations adopt two-phase background states that result from thermal instability in atomic gas. Although the supernova remnant (SNR) morphology becomes highly complex for inhomogeneous backgrounds, the radial momentum injection is remarkably insensitive to environmental details. For our two-phase simulations, the final momentum produced by a single SN is given by 2.8×10{sup 5} M{sub ⊙} km s{sup −1}n{sub 0}{sup −0.17}. This is only 5% less than the momentum injection for a homogeneous environment with the same mean density, and only 30% greater than the momentum at the time of shell formation. The maximum mass in hot gas is also quite insensitive to environmental inhomogeneity. Strong magnetic fields alter the hot gas mass at very late times, but the momentum injection remains the same. Initial experiments with multiple spatially correlated SNe show a momentum per event nearly as large as single-SN cases. We also present a full numerical parameter study to assess convergence requirements. For convergence in the momentum and other quantities, we find that the numerical resolution Δ and the initial size of the SNR r{sub init} must satisfy Δ,r{sub init}<r{sub sf}/3, where the shell formation radius is given by r{sub sf}=30 pc n{sub 0}{sup −0.46} for two-phase models (or 30% smaller for a homogeneous medium).

79 ASTRONOMY AND ASTROPHYSICS↗

Numerical Studies of Optimization and Aberration Correction Methods for the Preliminary Demonstration of the Parametric Ionization Cooling (PIC) Principle in the Twin Helix Muon Cooling Channel

Muon colliders have been proposed for the next generation of particle accelerators that study high-energy physics at the energy and intensity frontiers. In this paper we study a possible implementation of muon ionization cooling, Parametric-resonance Ionization Cooling (PIC), in the twin helix channel. The resonant cooling method of PIC offers the potential to reduce emittance beyond that achievable with ionization cooling with ordinary magnetic focusing. We examine optimization of a variety of parameters, study the nonlinear dynamics in the twin helix channel and consider possible methods of aberration correction.

43 PARTICLE ACCELERATORS↗

Simulation study of interactions of Space Shuttle-generated electron beams with ambient plasmas

This report summarizes results obtained through the support of NASA Grant NAGW-1936. The objective of this report is to conduct large scale simulations of electron beams injected into space. The topics covered include the following: (1) simulation of radial expansion of an injected electron beam; (2) simulations of the active injections of electron beams; (3) parameter study of electron beam injection into an ionospheric plasma; and (4) magnetosheath-ionospheric plasma interactions in the cusp.

Lin, Chin S.↗

The chemistry of interstellar argonium and other probes of the molecular fraction in diffuse clouds

We present a general parameter study in which the abundance of interstellar argonium (ArH{sup +}) is predicted using a model for the physics and chemistry of diffuse interstellar gas clouds. Results have been obtained as a function of UV radiation field, cosmic-ray ionization rate, and cloud extinction. No single set of cloud parameters provides an acceptable fit to the typical ArH{sup +}, OH{sup +}, and H{sub 2}O{sup +} abundances observed in diffuse clouds within the Galactic disk. Instead, the observed abundances suggest that ArH{sup +} resides primarily in a separate population of small clouds of total visual extinction of at most 0.02 mag per cloud, within which the column-averaged molecular fraction is in the range 10{sup −5}--10{sup −2}, while OH{sup +} and H{sub 2}O{sup +} reside primarily in somewhat larger clouds with a column-averaged molecular fraction ∼0.2. This analysis confirms our previous suggestion that the argonium molecular ion is a unique tracer of almost purely atomic gas.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimization of duct acoustic liners of finite length

Lining design procedures generally ignore the effects of discontinuous impedance changes on reflection and modal redistribution of incident energy. The present study involves identification of the mechanisms of optimal suppression by segmented liners, including impedance discontinuity effects. The mathematical model involves rectangular ducts in the absence of mean flow. Segmented acoustic fields are represented by duct eigenfunction expansions, and mode matching is used to ensure continuity of the total field. Extensive parameter studies were used to identify attenuation mechanisms and to define preliminary liner configurations. These are used as starting values in conjugate-gradient optimization studies. Optimal segmented liners are presented and it is shown that, if the sound source is sufficiently well-defined, conventional infinite duct optimum attenuation rates can be improved upon by taking advantage of wall impedance discontinuities. Finally, follow-on work which will involve experimental verification and extension to the case with mean flow in the duct is described.

Beckemeyer, R. J.↗

Simulating the Cranfield geological carbon sequestration project with high-resolution static models and an accurate equation of state

In this study, a field-scale carbon dioxide (CO 2 ) injection pilot project was conducted as part of the Southeast Regional Sequestration Partnership (SECARB) at Cranfield, Mississippi. We present higher-order finite element simulations of the compositional two-phase CO 2 -brine flow and transport during the experiment. High- resolution static models of the formation geology in the Detailed Area Study (DAS) located below the oil- water contact (brine saturated) are used to capture the impact of connected flow paths on breakthrough times in two observation wells. Phase behavior is described by the cubic-plus-association (CPA) equation of state, which takes into account the polar nature of water molecules. Parameter studies are performed to investigate the importance of Fickian diffusion, permeability heterogeneity, relative permeabilities, and capillarity. Simulation results for the pressure response in the injection well and the CO 2 breakthrough times at the observation wells show good agreement with the field data. For the high injection rates and short duration of the experiment, diffusion is relatively unimportant (high P clet numbers), while relative permeabilities have a profound impact on the pressure response. High-permeability pathways, created by fluvial deposits, strongly affect the CO 2 transport and highlight the importance of properly characterizing the formation heterogeneity in future carbon sequestration projects.

54 ENVIRONMENTAL SCIENCES↗

Continued Development of the Brine Evaporation Bag

The Brine Evaporation Bag (BEB) is a membrane-based bag system for the dewatering of brine. Previous studies showed the ability of the BEB to dewater brine at low temperatures with a 96 percent mass reduction. Additionally, a microgravity flight showed the BEB is microgravity compatible. Current work focuses on the effects of temperature, vacuum, purge gas flow rate, membrane area, and membrane permeability on the rate of dewatering within a vacuum oven configured to mimic the Heat Melt Compactor. Within this study, it was found that changing the temperature or level of vacuum would change the rate of dewatering. The purge gas, membrane area, and membrane permeability did not affect the dewater rate. The reason for this behavior may be that the dewatering is heat transfer limited, and out of all the parameters studied, only the temperature and vacuum have an effect on the heat transfer rate. The ISS (International Space Station) produces brine at a rate of 1.2 liters per day. This initial study showed that it is possible to remove water from a BEB at a rate of 1.6 liters per day in this breadboard configuration; even at moderate temperatures. Development of a dedicated BEB Evaporator will be discussed. In addition, it is further postulated that a specifically designed BEB Evaporator would result in an increased dewatering rate allowing for even lower operating temperatures or faster dewatering rates.

Recycling↗

Low-density photodissociation regions

A theoretical parameter study of the temperature and chemical structure of photodissociation regions (PDRs) and their resultant spectrum is presented. The models, which solve for the chemical and thermal steady state of a one-dimensional slab exposed to an incident FUV flux, relate the observed line and continuum emission from these regions to physical parameters of interest. Intensities of numerous fine-structure lines, the molecular rotational transitions (C-12)O J = 1-0 and H2 0-0 S(0), and the dust continuum intensities at 60 and 100 microns are predicted, and the column density of FUV-pumped vibrationally excited H2 is estimated.

Hollenbach, David J.↗

On the decay of outbursts in dwarf novae nad X-ray novae

We perform computations using a time-dependent model for the accretion disk limit-cycle mechanism to examine the decay of the optical light following the peak of a dwarf nova outburst. We present the results of a parameter study of the physical input variables which affect the decay rate. In the model, the decay is brought about by a cooling transition front which begins at large radii in the disk and moves inward. The nature of the decay is strongly influenced by the radial dependence of the accretion disk viscosity parameter alpha. To obtain exponential decays for typical dwarf nova parameters, we require alpha proportional to r(exp epsilon(sub 0)), where epsilon(sub 0) approximately = 0.3-0.4. The exact value of epsilon(sub 0) which produces exponential decays depends on factors such as the mass of the accreting star and the inner radius of the accretion disk. Therefore, the observed ubiquity of exponential decays in two different types of systems (dwarf novae and X-ray novae) leads us to believe that alpha is an unnatural scaling for the viscosity. The physics of the cooling transition front must be self-regulating in that the timescale (-parital derivative of lnSigma(r)/partial derivative +)(exp -1) (where Sigma is the surface density) for mass extraction across the front remains constant. This may be consistent with a scaling alpha proportional to (h/r)(exp n), where h is the local disk semi-thickness and n approximately 1-2. As regards the speed of the cooling front, we find v(sub F)(r) proportional to r(exp p), where p approximately 3 at large radii, with an abrupt transition to p approximately 0 at some smaller radius. The r(exp 3) dependence is much steeper than has been found by previous workers and appears to result from the strong variation of specific heat within the cooling front when the front resides at a large radius in the disk. The outflow of disk material across the cooling front causes a significant departure of dln T(sub dff0/dln r from the standard value of -0.75 (expected from steady state accretion) within about 0.2 dex in radius of the break associated with the cooling front -- T(sub eff) aproximately 10(exp 3.9) K (r/10(exp 10 cm)) (exp -0.1). These effects should be observable with eclipse mapping. Finally, it appears that the relatively slow decay rate for the optical flux in the 1975 outburst of A0620-00 can be accounted for if the primary is a approximately 10 Solar mass black hole.

Cannizzo, John K.↗

Characterizing the effects of drive asymmetries, component offsets, and joint gaps in double shell capsule implosions

This work provides a numerical study of how double shell capsule deformations caused by drive asymmetries and fabrication imperfections affect implosion symmetry and neutron yield. Hydrodynamics simulations are performed in two dimensions and focus on low-mode deformations that are caused by corresponding asymmetries in the Hohlraum drive, component offsets, and ablator joint gaps. By providing a parameter study of these features, our goal is to understand the dominant sources for inner shell deformation and yield degradation. The discussed capsules are composed of an aluminum ablator with a chromium inner shell. The latter encloses a carbon-deuterium foam ball that serves as fuel. We find that for clean capsules, even-numbered low-mode asymmetries in the drive are imprinted on the ablator and smoothly transferred to the inner shell during shell collision. The resulting deformation of the inner shell is more pronounced with larger fuel radius, while the yield is inversely proportional to the amplitude of the drive asymmetry and varies by factors ≤4 in comparison with clean simulations. Capsule component offsets in the vertical direction and ablator thickness nonuniformity result in p1-type deformations of the imploding inner shell. Finally, joint gaps have the largest effect in deforming the ablator and inner shell and degrading yield. While small gap widths (1 μm) result in prolate inner shells, larger gap widths (4 μm) cause an oblate deformation. More importantly, capsules with a small outer gap (1 μm) experience a dramatic drop in yield, typically <3% of a clean simulation.

Goodarzi, S. (ORCID:000000018556489X)↗