Search NASA⌕ Search

SEARCH · Search NASA

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 235 records · Page 13

Multicolor eclipse studies of UU Aquarii. 1: Observations and system parameters

A study of the eclipses in UU Aqr from multicolor high-speed photometry is presented. A revised ephemeris for the times of minimum and an upper limit for orbital period variations are obtained. We use measurements of contact phases in the eclipse light curve to derive the binary geometry and to estimate masses and relevant dimensions. We find a mass ratio of q = 0.30 +/- 0.07 and an inclination of i = 78 deg +/- 2 deg. The masses of the component stars are M(sub 1) = 0.67 +/- 0.14 solar mass and M(sub 2) = 0.20 +/- 0.07 solar mass. Our photometric model predicts K(sub 1) = 84 +/- 26 km/s, which is approximately 30% smaller than the velocity amplitude obtained from the emission lines. From the white dwarf fluxes we estimate T(sub wd) approximately = 34,000 K and a distance of d = 270 +/- 50 pc if the inner disk is opaque. UU Aqr has long term brightness variations of approximately = 0.3 m on timescales of approximately 4 yr. The system was in a 'high' state in 1989 and 1990 and in a 'low' state in 1988 and 1992. The high state results from an increase in the brightness of the outer and cooler parts of the disk, mainly due to the appearance of a bright spot at disk rim. Based on the smooth and gradual eclipse shape and on the absence of a prominent hump in the light curve we suggest that UU Aqr is a high mass-transfer nova-like system with a relatively bright and optically thick accretion disk. We find no perceptible eclipse in the H-alpha emission line. The fluxes at mid-eclipse can be fitted by a compostion of a late-type spectrum plus an optically thin hydrogen emission-line spectrum. These evidences suggest that the emission lines are formed in an extended region only partially occulted during eclipse.

Baptista, R.↗

A study of overflow simulations using MPAS-Ocean: Vertical grids, resolution, and viscosity

MPAS-Ocean is used to simulate an idealized, density-driven overflow using the dynamics of overflow mixing and entrainment (DOME) setup. Additionally, numerical simulations are carried out using three of the vertical coordinate types available in MPAS-Ocean, including z-star with partial bottom cells, z-star with full cells, and sigma coordinates. The results are first benchmarked against other models, including the MITgcm’s z-coordinate model and HIM’s isopycnal coordinate model, which are used to set the base case used for this work. A full parameter study is presented that looks at how sensitive overflow simulations are to vertical grid type, resolution, and viscosity. Horizontal resolutions with 50 km grid cells are under-resolved and produce poor results, regardless of other parameter settings. Vertical grids ranging in thickness from 15 m to 120 m were tested. A horizontal resolution of 10 km and a vertical resolution of 60 m are sufficient to resolve the mesoscale dynamics of the DOME configuration, which mimics real-world overflow parameters. Mixing and final buoyancy are least sensitive to horizontal viscosity, but strongly sensitive to vertical viscosity. This suggests that vertical viscosity could be adjusted in overflow water formation regions to influence mixing and product water characteristics. In conclusion, the study shows that sigma coordinates produce much less mixing than z-type coordinates, resulting in heavier plumes that go further down slope. Sigma coordinates are less sensitive to changes in resolution but as sensitive to vertical viscosity compared to z-coordinates.

58 GEOSCIENCES↗

Parameter optimization for transitions between memory states in small arrays of Josephson junctions

Coupled arrays of Josephson junctions possess multiple stable zero voltage states. Such states can store information and consequently can be utilized for cryogenic memory applications. Basic memory operations can be implemented by sending a pulse to one of the junctions and studying transitions between the states. In order to be suitable for memory operations, such transitions between the states have to be fast and energy efficient. Here in this article we employed simulated annealing, a stochastic optimization algorithm, to study parameter optimization of array parameters which minimizes times and energies of transitions between specifically chosen states that can be utilized for memory operations (Read, Write, and Reset). Simulation results show that such transitions occur with access times on the order of 10–100 ps and access energies on the order of 10 -19 –5×10 -18 J. Numerical simulations are validated with approximate analytical results.

97 MATHEMATICS AND COMPUTING↗

Wave-Particle Interactions Associated with Nongyrotropic Distribution Functions: A Hybrid Simulation Study

Nongyrotropic plasma distribution functions can be formed in regions of space where guiding center motion breaks down as a result of strongly curved and weak ambient magnetic fields. Such are the conditions near the current sheet in the Earth's middle and distant magnetotail, where observations of nongyrotropic ion distributions have been made. Here a systematic parameter study of nongyrotropic proton distributions using electromagnetic hybrid simulations is made. We model the observed nongyrotropic distributions by removing a number of arc length segments from a cold ring distribution and find significant differences with the results of simulations that initially have a gyrotropic ring distribution. Model nongyrotropic distributions with initially small perpendicular thermalization produce growing fluctuations that diffuse the ions into a stable Maxwellian-like distribution within a few proton gyro periods. The growing waves produced by nongyrotropic distributions are similar to the electromagnetic proton cyclotron waves produced by a gyrotropic proton ring distribution in that they propagate parallel to the background magnetic field and occur at frequencies on the order of the proton gyrofrequency, The maximum energy of the fluctuating magnetic field increases as the initial proton distribution is made more nongyrotropic, that is, more highly bunched in perpendicular velocity space. This increase can be as much as twice the energy produced in the gyrotropic case.

Convery, P. D.↗

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↗