Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “V CODES”

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,261 records · Page 70

Measuring magnetic fields in laser-driven coils with dual-axis proton deflectometry

Abstract By driving hot electrons between two metal plates connected by a wire loop, high power lasers can generate multi-tesla, quasi-static magnetic fields in miniature coil targets. Many experiments involving laser-coil targets rely on proton deflectometry directed perpendicular to the coil axis to extract a measurement of the magnetic field. In this paper, we show that quantitative measurements using perpendicular probing are complicated by the presence of GV m −1 electric fields in the target that develop on sub-ns timescales. Probing parallel to the coil axis with fiducial grids is shown to reliably separate the electric and magnetic field measurements, giving current estimates of I ≈ 5 kA in 1 mm- and 2 mm-diameter wire loops. An analytic model of proton deflection in electric and magnetic fields is used to benchmark results from the particle-in-cell code and help deconvolve the magnetic and electric field deflections. Results are used to motivate a new experimental scheme that combines a single-plate target with axial proton probing and direct current measurements. This scheme has several important advantages over the traditional target and diagnostic set-up, enabling the robust measurement of coil magnetic fields and plasma properties, as well as making it easier to validate different theoretical models at a range of laser intensities.

Bradford, P. (ORCID:0000000296186342)↗

Computation of hypersonic laminar viscous flow past spinning sharp and blunt cones at high angle of attack

Computational results, obtained with a parabolic Navier-Stokes marching code, are presented for hypersonic viscous flow past spinning sharp and blunt cones at angle of attack. The code takes into account the asymmetries in the flow field resulting from spinning motion and computes the asymmetric shock shape, crossflow and streamwise shear, heat transfer, crossflow separation, and vortex structure. The Magnus force and moments are also computed. Comparisons are made with other theoretical analyses based on boundary-layer and boundary-region equations, and an anomaly is discovered in the displacement thickness contribution to the Magnus force when compared with boundary-layer results. In addition, a new criterion for defining crossflow separation behind spinning bodies is introduced which generalizes the Moore-Rott-Sears criterion for two-dimensional unsteady separation. A condition which characterizes the onset of separation in the flow field is defined.

Agarwal, R.↗

Water Vapor on Titan: The Stratospheric Vertical Profile from Cassini/CIRS Infrared Spectra

Water vapor in Titan's middle atmosphere has previously been detected only by disk-average observations from the Infrared Space Observatory (Coustenis et al., 1998). We report here the successful detection of stratospheric water vapor using the Cassini Composite Infrared Spectrometer (CIRS, Flasar et al., 2004) following an earlier null result (de Kok et al., 2007a). CIRS senses water emissions in the far-infrared spectral region near 50 microns, which we have modeled using two independent radiative transfer and inversion codes (NEMESIS, Irwin et al 2008 and ART, Coustenis et al., 2010). From the analysis of nadir spectra we have derived a mixing ratio of (0.14 plus or minus 0.05) ppb at 100 km, corresponding to a column abundance of approximately (3.7 plus or minus 1.3) x 10(exp 14) moles per square centimeter. Using limb observations, we obtained mixing ratios of (0.13 plus or minus 0.04) ppb at 125 km and (0.45 plus or minus 0.15) ppb at 225 km of altitude, confirming that the water abundance has a positive vertical gradient as predicted by photochemical models. In the latitude range (80 deg. S - 30 deg. N) we see no evidence for latitudinal variations in these abundances within the error bars.

Cottini, V.↗

Wing design for a civil tiltrotor transport aircraft

The goal of this research is the proper tailoring of the civil tiltrotor's composite wing-box structure leading to a minimum-weight wing design. With focus on the structural design, the wing's aerodynamic shape and the rotor-pylon system are held fixed. The initial design requirement on drag reduction set the airfoil maximum thickness-to-chord ratio to 18 percent. The airfoil section is the scaled down version of the 23 percent-thick airfoil used in V-22's wing. With the project goal in mind, the research activities began with an investigation of the structural dynamic and aeroelastic characteristics of the tiltrotor configuration, and the identification of proper procedures to analyze and account for these characteristics in the wing design. This investigation led to a collection of more than thirty technical papers on the subject, some of which have been referenced here. The review of literature on the tiltrotor revealed the complexity of the system in terms of wing-rotor-pylon interactions. The aeroelastic instability or whirl flutter stemming from wing-rotor-pylon interactions is found to be the most critical mode of instability demanding careful consideration in the preliminary wing design. The placement of wing fundamental natural frequencies in bending and torsion relative to each other and relative to the rotor 1/rev frequencies is found to have a strong influence on the whirl flutter. The frequency placement guide based on a Bell Helicopter Textron study is used in the formulation of frequency constraints. The analysis and design studies are based on two different finite-element computer codes: (1) MSC/NASATRAN and (2) WIDOWAC. These programs are used in parallel with the motivation to eventually, upon necessary modifications and validation, use the simpler WIDOWAC code in the structural tailoring of the tiltrotor wing. Several test cases were studied for the preliminary comparison of the two codes. The results obtained so far indicate a good overall agreement between the two codes.

Rais-Rohani, Masoud↗

Deformation of the $0^+_{1,2}$ states in 110 Cd from low-energy Coulomb excitation

Electromagnetic properties of 110 Cd were studied via low-energy Coulomb excitation with 32 S and 14 N beams. Magnitudes and relative signs of eight E2 matrix elements, including quadrupole moments of the $2^+_1$ and $2^+_2$ states, were determined using the least-squares code GOSIA. From those, quadrupole deformation parameters of the $0^+_{1,2,3}$ states were inferred, providing for the first time conclusive evidence for the non-axial character of the ground state in 110 Cd. The experimental results were compared with new calculations using the general quadrupole collective Bohr Hamiltonian model with SLy4 and UNEDF0 interactions. The non-axiality of the ground state is reproduced by the present calculations, independently of the interaction used.

Collective models↗

Towards Z-Next: The Integration of Theory, Experiments, and Computational Simulation in a Bayesian Data Assimilation Framework

Making reliable predictions in the presence of uncertainty is critical to high-consequence modeling and simulation activities, such as those encountered at Sandia National Laboratories. Surrogate or reduced-order models are often used to mitigate the expense of performing quality uncertainty analyses with high-fidelity, physics-based codes. However, phenomenological surrogate models do not always adhere to important physics and system properties. This project develops surrogate models that integrate physical theory with experimental data through a maximally-informative framework that accounts for the many uncertainties present in computational modeling problems. Correlations between relevant outputs are preserved through the use of multi-output or co-predictive surrogate models; known physical properties (specifically monotoncity) are also preserved; and unknown physics and phenomena are detected using a causal analysis. By endowing surrogate models with key properties of the physical system being studied, their predictive power is arguably enhanced, allowing for reliable simulations and analyses at a reduced computational cost.

97 MATHEMATICS AND COMPUTING↗

Computation of nonequilibrium three-dimensional inviscid flow over blunt-nosed bodies flying at supersonic speeds

A computer code based on the method of characteristics is described and applied to the study of two- and three-dimensional chemical nonequilibrium flow over sharp and blunt-nosed bodies. Nonequilibrium flow over a wedge is used to show the approach to equilibrium flow, and to demonstrate the nature of the reaction zone behind the bow shock wave. The structure and development of a blunt-body entropy layer in nonequilibrium flow is examined for a blunt cone at zero incidence. Three-dimensional computations for the space shuttle body at 30 deg angle of attack are presented. A nondimensional scaling parameter, the Damkoehler number, which is the ratio of flow time to chemical reaction time, is calculated and its significance discussed.

Rakich, J. V.↗

Bar Code Events in the Juno-UVS Data: Signature ~10 MeV Electron Microbursts at Jupiter

One of the most intriguing discoveries of Juno is the quasi-systematic detection of upgoing electrons above the auroral regions. Here we discuss a by-product of the most energetic component of this population: a contamination resembling bar codes in the Juno-UVS images. This pattern is likely caused by bursts of ∼10 MeV electrons penetrating the instrument. These events are mostly detected when Juno’s magnetic footprint is located poleward of the main emission relative to the magnetic pole. The signal is not periodic, but the bursts are typically 0.1–1 s apart. They are essentially detected when Juno-UVS is oriented toward Jupiter, indicating that the signal is due to upgoing electrons. The event detections occur between 1 and 7 Jovian radii above the 1-bar level, suggesting that the electron acceleration takes place close to Jupiter and is thus both strong and brief.

Microbursts↗

Overview of T and D–T results in JET with ITER-like wall

In 2021 JET exploited its unique capabilities to operate with T and D–T fuel with an ITER-like Be/W wall (JET-ILW). This second major JET D–T campaign (DTE2), after DTE1 in 1997, represented the culmination of a series of JET enhancements—new fusion diagnostics, new T injection capabilities, refurbishment of the T plant, increased auxiliary heating, in-vessel calibration of 14 MeV neutron yield monitors—as well as significant advances in plasma theory and modelling in the fusion community. DTE2 was complemented by a sequence of isotope physics campaigns encompassing operation in pure tritium at high T-NBI power. Carefully conducted for safe operation with tritium, the new T and D–T experiments used 1 kg of T (vs 100 g in DTE1), yielding the most fusion reactor relevant D–T plasmas to date and expanding our understanding of isotopes and D–T mixture physics. Furthermore, since the JET T and DTE2 campaigns occurred almost 25 years after the last major D–T tokamak experiment, it was also a strategic goal of the European fusion programme to refresh operational experience of a nuclear tokamak to prepare staff for ITER operation. The key physics results of the JET T and DTE2 experiments, carried out within the EUROfusion JET1 work package, are reported in this paper. Progress in the technological exploitation of JET D–T operations, development and validation of nuclear codes, neutronic tools and techniques for ITER operations carried out by EUROfusion (started within the Horizon 2020 Framework Programme and continuing under the Horizon Europe FP) are reported in (Litaudon et al Nucl. Fusion accepted), while JET experience on T and D–T operations is presented in (King et al Nucl. Fusion submitted).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A rotor-mounted digital instrumentation system for helicopter blade flight research measurements

A rotor mounted flight instrumentation system developed for helicopter rotor blade research is described. The system utilizes high speed digital techniques to acquire research data from miniature pressure transducers on advanced rotor airfoils which are flight tested on an AH-1G helicopter. The system employs microelectronic pulse code modulation (PCM) multiplexer digitizer stations located remotely on the blade and in a hub mounted metal canister. As many as 25 sensors can be remotely digitized by a 2.5 mm thick electronics package mounted on the blade near the tip to reduce blade wiring. The electronics contained in the canister digitizes up to 16 sensors, formats these data with serial PCM data from the remote stations, and transmits the data from the canister which is above the plane of the rotor. Data are transmitted over an RF link to the ground for real time monitoring and to the helicopter fuselage for tape recording. The complete system is powered by batteries located in the canister and requires no slip rings on the rotor shaft.

Knight, V. H., Jr.↗

Relativistic tidal interaction of a white dwarf with a massive black hole

We compute encounters of a realistic white dwarf model with a massive black hole in the regime where relativistic effects are important, using a three-dimensional, finite-difference, Eulerian, piecewise parabolic method (PPM) hydrodynamical code. Both disruptive and nondisruptive encounters are considered. We identify and discuss relativistic effects important for the problem: relativistic shift of the pericenter distance, time delay, relativistic precession, and the tensorial structure of the tidal forces. In the nondisruptive case, stripping of matter takes place. In the surface layers of the surviving core, complicated hydrodynamical phenomena are revealed. In both disruptive and nondispruptive encounters, material flows out in the form of two thin, S-shaped, supersonic jets. Our results provide realistic initial conditions for the subsequent investigation of the dynamics of the debris in the field of the black hole. We evaluate the critical conditions for complete disruption of the white dwarf, and compare our results with the corresponding results for nonrelativistic encounters.

Frolov, V. P.↗

Simulation and Experimental Measurements of Inductively Coupled CF4 and CF4/Ar Plasmas

The recently developed code SEMS (semiconductor equipment modeling software)is applied to the simulation of CF4 and CF4/Ar inductively coupled plasmas (ICP). This work builds upon the earlier nitrogen, transformer coupled plasma (TCP) SEMS research by demonstrating its accuracy for more complex reactive mixtures, moving closer to the realization of a virtual plasma reactor. Attention is given to the etching of and/or formation of carbonaceous films on the quartz dielectric window and diagnostic aperatures. The simulations are validated through comparisons with experimental measurements using FTIR (Fourier Transform Infrared) and UV absorption spectroscopy for CFx and SiFx neutral radicals, QMS (quadrupole mass spectrometry) for the ions, and Langmuir probe measurements of electron number density and temperature in an ICP GEC reference cell.

Hash, D. B.↗

Experimental spectrum averaged cross sections (SACS) in 252 Cf(sf) neutron field and its impact on the evaluation of neutron standards

A new evaluation of spectrum averaged cross sections (SACS) of 235 U, 238 U, and 239 Pu measured in the 252 Cf(sf) reference neutron field is presented and found to be consistent with the original Mannhart SACS evaluation in IRDF-2002. The comprehensive vetted experimental database that includes SACS ratio and absolute SACS measurements of major actinides is being used to update the SACS database employed as input of the new GMApy code to derive the Neutron Standards. An update of the current neutron standards based on Time Projection Chamber (TPC) shape data, a new comprehensive uncertainty quantification, and revised SACS experimental database is proposed which result in a 0.7% increase of the evaluated 239 Pu(n,f)/ 235 U(n,f) crosssection ratio in the 1–5 MeV energy region. The increase is due to a 0.3% reduction of the standard 235 U(n,f) cross section and a 0.4% increase of the 239 Pu(n,f) reference cross section in the 1–5 MeV energy region. Those changes are well within estimated USU fission cross-section uncertainties of 1.2%, but are relevant for the evaluated mean values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Design and execution of a verification, validation, and uncertainty quantification plan for a numerical model of left ventricular flow after LVAD implantation

Left ventricular assist devices (LVADs) are implantable pumps that act as a life support therapy for patients with severe heart failure. Despite improving the survival rate, LVAD therapy can carry major complications. Particularly, the flow distortion introduced by the LVAD in the left ventricle (LV) may induce thrombus formation. While previous works have used numerical models to study the impact of multiple variables in the intra-LV stagnation regions, a comprehensive validation analysis has never been executed. The main goal of this work is to present a model of the LV-LVAD system and to design and follow a verification, validation and uncertainty quantification (VVUQ) plan based on the ASME V&V40 and V&V20 standards to ensure credible predictions. The experiment used to validate the simulation is the SDSU cardiac simulator, a bench mock-up of the cardiovascular system that allows mimicking multiple operation conditions for the heart-LVAD system. The numerical model is based on Alya, the BSC’s in-house platform for numerical modelling. Alya solves the Navier-Stokes equation with an Arbitrary Lagrangian-Eulerian (ALE) formulation in a deformable ventricle and includes pressure-driven valves, a 0D Windkessel model for the arterial output and a LVAD boundary condition modeled through a dynamic pressure-flow performance curve. The designed VVUQ plan involves: (a) a risk analysis and the associated credibility goals; (b) a verification stage to ensure correctness in the numerical solution procedure; (c) a sensitivity analysis to quantify the impact of the inputs on the four quantities of interest (QoIs) (average aortic root flow $Q^{avg}_{Ao}$, maximum aortic root flow $Q^{max}_{Ao}$, average LVAD flow $Q^{avg}_{VAD}$, and maximum LVAD flow $Q^{max}_{VAD}$; (d) an uncertainty quantification using six validation experiments that include extreme operating conditions. Numerical code verification tests ensured correctness of the solution procedure and numerical calculation verification showed a grid convergence index (GCI)95% <3.3%. The total Sobol indices obtained during the sensitivity analysis demonstrated that the ejection fraction, the heart rate, and the pump performance curve coefficients are the most impactful inputs for the analysed QoIs. The Minkowski norm is used as validation metric for the uncertainty quantification. It shows that the midpoint cases have more accurate results when compared to the extreme cases. The total computational cost of the simulations was above 100 [core-years] executed in around three weeks time span in Marenostrum IV supercomputer. This work details a novel numerical model for the LV-LVAD system, that is supported by the design and execution of a VVUQ plan created following recognised international standards. We present a methodology demonstrating that stringent VVUQ according to ASME standards is feasible but computationally expensive.

59 BASIC BIOLOGICAL SCIENCES↗

Calculations of the flow past bluff bodies, including tilt-rotor wing sections at alpha = 90 deg

An attempt was made to model in two dimensions the effects of rotor downwash on the wing of the tilt-rotor aircraft and to compute the drag force on airfoils at - 90 deg angle of attack, using a well-established Navier-Stokes code. However, neither laminar nor turbulent calculations agreed well with drag and base-pressure measurements at high Reynolds numbers. Therefore, further efforts were concentrated on bluff-body flows past various shapes at low Reynolds numbers, where a strong vortex shedding is observed. Good results were obtained for a circular cylinder, but the calculated drag of a slender ellipse at right angles to the freestream was significantly higher than experimental values reported in the literature for flat plates. Similar anomalous results were obtained on the tilt-rotor airfoils, although the qualitative effects of flap deflection agreed with the wind tunnel data. The ensemble of results suggest that there may be fundamental differences in the vortical wakes of circular cylinders and noncircular bluff bodies.

Raghavan, V.↗

Turbulent Diffusivities for Momentum, Heat, Salt and Passive Scalars

A program began ten years ago to build a turbulence model to describe high Reynolds numbers flows. Specifically, the aim was to devise a model that would satisfy two basic conditions: reproduce available turbulence data (laboratory, numerical simulations, etc.) concerning shear driven flows, buoyancy driven flows, 2D turbulence, freely decaying turbulence etc., and yet be manageable so as to be used, for example, in GCM'S. The model was presented in a series of papers that have appeared in Physics of Fluids since 1996. A total of about 80 turbulence statistics were reproduced. The model has no adjustable parameters. The next step was to apply the model to construct the vertical diffusities K for moment K(sub m), temperature K(sub h), salt K(sub s), and passive scalars K(sub c). First, we took K(sub c) = K(sub h) and tested the model using the GFDL ocean model. The results for the profiles of T and S vs. depth are indistinguishable from those derided using the latest model, the KPP model by the NCAR group. Presently, we are running the same GFDL code relaxing the assumptive K(sub c) = K(sub h). Indeed, the turbulence model yelds a salt diffusivity that depends on R(sub i) and R(sub rho) = ((beta)(delta)S/(alpha)(delta)T/(delta)z) in such a way that K(sub s) may be quite different from K(sub h). Salt fingers and double diffusivity laboratory data are reproduced. Results from the ocean model will be available shortly. Finally, we are trying to derive the horizontal diffusivities with the goal of providing a physically acceptable representation of mesoscale eddies. The recently suggested GMW parameterization has improved several O-GCM results and the goal here is to try to derive/justify it from a turbulence model and/or propose improvements/modifications. Theoretical work is in progress.

Canuto, V. M.↗

Integration of scanning probe microscope with high-performance computing: Fixed-policy and reward-driven workflows implementation

The rapid development of computation power and machine learning algorithms has paved the way for automating scientific discovery with a scanning probe microscope (SPM). The key elements toward operationalization of the automated SPM are the interface to enable SPM control from Python codes, availability of high computing power, and development of workflows for scientific discovery. Here, we build a Python interface library that enables controlling an SPM from either a local computer or a remote high-performance computer, which satisfies the high computation power need of machine learning algorithms in autonomous workflows. We further introduce a general platform to abstract the operations of SPM in scientific discovery into fixed-policy or reward-driven workflows. Furthermore, our work provides a full infrastructure to build automated SPM workflows for both routine operations and autonomous scientific discovery with machine learning.

47 OTHER INSTRUMENTATION↗

Simulating energetic particle losses in JET plasmas with a reverse integration biasing scheme

We report that an integrated energetic particle transport model has been constructed in JET plasmas constrained by experimental fast ion loss measurements. The model incorporates a synthetic fast ion loss detector identical to JET’s thin-foil Faraday cup fast ion loss detector array. The loss model combines analyses from the TRANSP and ORBIT-kick codes with enhanced features for producing the synthetic diagnostic. Extensions to the ORBIT code framework allow a full-orbit representation within the vacuum region that can map particles directly to an installed detector geometry. Since synthetic fast ion loss detectors are plagued by weak loss statistics, a novel reverse integration biasing scheme has been implemented to boost computational efficiency. The model is validated against experimental loss measurements induced by long-lived kink modes and is found to be in good agreement. This confirms the development of a fully integrated transport/loss model which can be quantitatively verified against experiment allowing for future validation and predictive studies. The model is particularly useful for more complicated plasma scenarios that involve multiple fast ion species such as JET’s 2021 DT-campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗