SEARCH · Engineering Papers

Results for “spherical fitting”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

516 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Statistical characterization of experimental magnetized liner inertial fusion stagnation images using deep-learning-based fuel–background segmentation

Significant variety is observed in spherical crystal x-ray imager (SCXI) data for the stagnated fuel–liner system created in Magnetized Liner Inertial Fusion (MagLIF) experiments conducted at the Sandia National Laboratories Z-facility. As a result, image analysis tasks involving, e.g., region-of-interest selection (i.e. segmentation), background subtraction and image registration have generally required tedious manual treatment leading to increased risk of irreproducibility, lack of uncertainty quantification and smaller-scale studies using only a fraction of available data. We present a convolutional neural network (CNN)-based pipeline to automate much of the image processing workflow. This tool enabled batch preprocessing of an ensemble of N scans = 139 SCXI images across N exp = 67 different experiments for subsequent study. The pipeline begins by segmenting images into the stagnated fuel and background using a CNN trained on synthetic images generated from a geometric model of a physical three-dimensional plasma. The resulting segmentation allows for a rules-based registration. Our approach flexibly handles rarely occurring artifacts through minimal user input and avoids the need for extensive hand labelling and augmentation of our experimental dataset that would be needed to train an end-to-end pipeline. Here we also fit background pixels using low-degree polynomials, and perform a statistical assessment of the background and noise properties over the entire image database. Our results provide a guide for choices made in statistical inference models using stagnation image data and can be applied in the generation of synthetic datasets with realistic choices of noise statistics and background models used for machine learning tasks in MagLIF data analysis. We anticipate that the method may be readily extended to automate other MagLIF stagnation imaging applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1984

Impact experiments in viscous fluid media

Available phase and group velocity data are inverted by a procedure which includes the effects of transverse anisotropy, anelastic dispersion, sphericity, and gravity. The resulting models, for average Earth, average ocean, and oceanic regions divided according to the age of the ocean floor, are quite different from previous results which ignore the above effects. The models show a low-veocity zone with age dependent anisotropy and velocities higher than derived in previous surface wave studied. The correspondence between the anisotropy variation with age and a physical model based on flow aligned olivine is suggestive. For most of the Earth SHSV in the vicinity of the low-velocity zone. Near the East Pacific Rise, however, SVSH at depth, consistent with ascending flow. Anisotropy is as important as temperature in causing radial and lateral variations in velocity. The models have a high velocity nearly isotropic layer at the top of the mantle that thickens with age. This layer defines the LID, or seismic lithosphere. In the Pacific, the LID thickens with age to a maximum thickness of about 50 km. This thickness is comparable to the thickness of the elastic lithosphere. The LID thickness is thinner than derived using isotropic or pseudo-isotropic procedures A new model for Average Earth is obtained which includes a thin LID. This model extends the fit of a P.R.E.M. type model to shorter period surface waves.

Greeley, R.↗

NASA NTRS2009

A Model of Thermal Conductivity for Planetary Soils: Theory for Cemented Soils - 2

A numerical model of heat conduction through particulate media made of spherical grains cemented by various bonding agents is presented. The pore-filling gas conductivity, volume fraction, and thermal conductivity of the cementing phase are tunable parameters. Cement fractions <0.001-0.01% in volume have small effects on the soil bulk thermal conductivity. A significant conductivity increase (factor 3-8) is observed for bond fractions of 0.01 to 1% in volume. In the 1 to 15% bond fraction domain, the conductivity increases continuously but less intensely (25-100% conductivity increase compared to a 1% bond system). Beyond 15% of cements, the conductivity increases vigorously and the bulk conductivity rapidly approaches that of bedrock. The composition of the cements (i.e. conductivity) has little influence on the bulk thermal inertia of the soil, especially if the volume of bond <10%. These results indicate that temperature measurements are sufficient to detect cemented soils and quantify the amount of cementing phase, but the mineralogical nature of the bonds and the typical grain size are unlikely to be determined from orbit. On Mars, a widespread surface unit characterized by a medium albedo (0.19-0.26) and medium/high thermal inertia (200-600 J s(0.5)/sq m/K) has long been hypothesized to be associated with a duricrust. The fraction of cement required to fit the thermal data is less than approx.1-5% by volume. This small amount of material is consistent with orbital observations, confirming that soil cementation is an important factor controlling the thermal inertia of the Martian surface

Piqueux, S.↗

DOE OSTI2014

Pressure solution creep of random packs of spheres

Abstract We performed numerical calculations of compaction in aggregates of spherical grains, using Lehner and Leroy's (2004, hereinafter LL) constitutive model of pressure solution at grain contacts. That model is founded on a local definition of the thermodynamic driving force and leads to a fully coupled formulation of elastic deformation, dissolution, and diffusive transport along the grain boundaries. The initial geometry of the aggregate was generated by random packing of spheres with a small standard deviation of the diameters. During the simulations, isostatic loading was applied. The elastic displacements at the contacts were calculated according to Digby's (1981) nonlinear contact force model, and deformation by dissolution was evaluated using the LL formulation. The aggregate strain and porosity were tracked as a function of time for fixed temperature, applied effective pressure, and grain size. We also monitored values of the average and standard deviation of total load at each contact, the coordination number for packing, and the statistics of the contact dimensions. Because the simulations explicitly exclude processes such as fracturing, plastic flow, and transport owing to surface curvature, they can be used to test the influence of relative changes in the kinetics of dissolution and diffusion processes caused by contact growth and packing rearrangements. We found that the simulated strain data could be empirically fitted by two successive power laws of the form, ε x ∝ t ξ , where ξ was equal to 1 at very early times, but dropped to as low as 0.3 at longer times. The apparent sensitivity of strain rate to stress found in the simulations was much lower than predicted from constitutive laws that assume a single dominant process driven by average macroscopic loads. Likewise, the apparent activation enthalpy obtained from the simulated data was intermediate between that assumed for dissolution and diffusion, and, further, tended to decrease with time. These results are similar to the experimental observations of Visser et al.'s (2012), who used an aggregate geometry and physical conditions closely resembling the present numerical simulations.

Bernabé, Y.↗

NASA NTRS1990

Velocity-Resolved Far-Infrared Spectra of [Fe 2]: Evidence for Mixing and Clumping in SN 1987A

We present approx. 400 km/s resolution profiles of the 17.94 and 25.99 micron [Fe II] transitions from SN 1987A at t approx. 400 days after core collapse. These observations used the facility cooled grating spectrometer aboard NASA's Kuiper Airborne Observatory. The two profiles are similar and have FWHM line widths of approx. 2700 km/s. The higher signal-to-noise 18 micron profile is somewhat asymmetric, falling off more steeply on the redshifted side than on the blue. Gaussian fits to the profiles yield an average centroid velocity of 280 +/- 140 km/s relative to the Large Magellanic Cloud. The wings of the profiles extend to velocities is approx. greater than 3000 km/s. This shows that a significant fraction of the iron has been mixed outward into the hydrogen-rich envelope, which has a minimum expansion velocity of 2100-2400 km/s. Both profiles also contain an unresolved 3-5 sigma emission feature on the redshifted wing at nu(LSR) approx. + 3900 km/s. We interpret this feature as emission from a high-velocity clump of material containing approx. 3% of the total iron mass. The total line flux of the 26 micron ground-state transition yields an optically thin, singly ionized iron mass of 0.026 solar mass, relatively independent of the assumed temperature. This is significantly less than the 0.06 Me of Fe+ determined from the decline of the optical light curve and the ionization of measured nickel lines, implying that the iron transitions still have appreciable optical depth. However, because of the small change in the 26 micron line flux from our measurement at 250 days, and the similarity of our profiles to the 1.26 micron [Fe II] profile, most of the emission is believed to originate from optically thin material with a temperature of 4406 +/- 400 K. A comparison of the data with spherically symmetric models indicates a power-law density exponent of -3.2 +/- 1.1 and a minimum expansion velocity of 650 +/- 650 km/s for this optically thin component. The [Fe II] line fluxes and profiles also imply that the remainder of the material has high optical depth and is distributed in clumps throughout the ejecta, rather than being concentrated at low velocities in the center of a smooth density distribution.

Haas, Michael R.↗

DOE OSTI2014

Fusion of {sup28}Si + {sup28,30}Si: Different Trends at Sub-Barrier Energies

Background: The fusion excitation function of the system Si-28 + Si-28 at energies near and below the Coulomb barrier is known only down to similar or equal to 15 mb. This precludes any information on both coupling effects on sub-barrier cross sections and the possible appearance of hindrance. For Si-28 + Si-30 even if the fusion cross section is measured down to similar or equal to 50 mu b, the evidence of hindrance is marginal. Both systems have positive fusion Q values. While Si-28 has a deformed oblate shape, Si-30 is spherical. Purpose: We investigate 1. the possible influence of the different structure of the two Si isotopes on the fusion excitation functions in the deep sub-barrier region and 2. whether hindrance exists in the Si + Si systems and whether it is strong enough to generate an S-factor maximum, thus allowing a comparison with lighter heavy-ion systems of astrophysical interest. Methods: Si-28 beams from the XTU Tandem accelerator of the INFN Laboratori Nazionali di Legnaro were used. The setup was based on an electrostatic beam separator, and fusion evaporation residues (ER) were detected at very forward angles. Angular distributions of ER were measured. Results: Fusion cross sections of Si-28 + Si-28 have been obtained down to similar or equal to 600 nb. The slope of the excitation function has a clear irregularity below the barrier, but no indication of a S-factor maximum is found. For Si-28 + Si-30 the previous data have been confirmed and two smaller cross sections have been measured down to similar or equal to 4 mu b. The trend of the S-factor reinforces the previous weak evidence of hindrance. Conclusions: The sub-barrier cross sections for Si-28 + Si-28 are overestimated by coupled-channels calculations based on a standard Woods-Saxon potential, except for the lowest energies. Calculations using the M3Y+repulsion potential are adjusted to fit the Si-28 + Si-28 and the existing Si-30 + Si-30 data. An additional weak imaginary potential (probably simulating the effect of the oblate Si-28 deformation) is required to fit the low-energy trend of Si-28 + Si-28. The parameters of these calculations are applied to predict the ion-ion potential for Si-28 + Si-30. Its cross sections are well reproduced by also including one-and successive two-neutron transfer channels, besides the low-lying surface excitations.

Montagnoli, G.↗

DOE OSTI2015

The spherically symmetric gravitational collapse of a clump of solids in a gas

In the subject of planetesimal formation, several mechanisms have been identified that create dense particle clumps in the solar nebula. The present work is concerned with the gravitational collapse of such clumps, idealized as being spherically symmetric. Fully nonlinear simulations using the two-fluid model are carried out (almost) up to the time when a central density singularity forms. We refer to this as the collapse time. The end result of the study is a parametrization of the collapse time, in order that it may be compared with timescales for various disruptive effects to which clumps may be subject in a particular situation. An important effect that determines the collapse time is that as the clump compresses, it also compresses the gas due to drag. This increases gas pressure, which retards particle collapse and can lead to oscillation in the size and density of the clump. In the limit of particles perfectly coupled to the gas, the characteristic ratio of gravitational force to gas pressure becomes relevant and defines a two-phase Jeans parameter, J{sub t}, which is the classical Jeans parameter with the speed of sound replaced by an effective wave speed in the coupled two-fluid medium. The parameter J{sub t} remains useful even away from the perfect coupling limit because it makes the simulation results insensitive to the initial density ratio of particles to gas (Φ{sub 0}) as a separate parameter. A simple ordinary differential equation model is developed. It takes the form of two coupled non-linear oscillators and reproduces key features of the simulations. Finally, a parametric study of the time to collapse is performed and a formula (fit to the simulations) is developed. In the incompressible limit J{sub t}→0, collapse time equals the self-sedimentation time, which is inversely proportional to the Stokes number. As J{sub t} increases, the collapse time decreases with J{sub t} and eventually becomes approximately equal to the dynamical time. Values of collapse time versus clump size are given for a minimum-mass solar nebula. Finally, the timescale of clump erosion due to turbulent strain is estimated.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Isospin-invariant Skyrme energy-density-functional approach with axial symmetry

Density functional theory (DFT) is the microscopic tool of choice to describe properties of nuclei over the entire nuclear landscape, with a focus on medium-mass and heavy complex systems. Modern energy density functionals (EDFs) often offer a level of accuracy typical of phenomenological approaches based on parameters locally fitted to the data. It is clear, however, that in order to achieve high quality of predictions to guide spectroscopic studies, current functionals must be improved, especially in the isospin channel. In this respect, experimental studies of short-lived nuclei far from stability offer a unique test of isospin aspects of the many-body theory. In this study, we develop the isospin-invariant Skyrme-EDF method by considering local densities in all possible isospin channels and proton-neutron ($p-n$) mixing terms as mandated by the isospin symmetry. The EDF employed has the most general form that depends quadratically on the isoscalar and isovector densities. We test and benchmark the resulting $p-n$ EDF approach, and study the general properties of the new scheme by means of the cranking in the isospin space. We extend the existing axial DFT solver HFBTHO to the case of isospin-invariant EDF approach with all possible $p-n$ mixing terms. Explicit expressions have been derived for all the densities and potentials that appear in the isospin representation. In practical tests, we consider the Skyrme EDF SkM* and, as a first application, concentrate on Hartree-Fock aspects of the problem, i.e., pairing has been disregarded. Calculations have been performed for the ($A = 78, T ≃ 11), (A = 40, T ≃ 8$), and ($A = 48, T ≃ 4$) isobaric analog chains. Isospin structure of self-consistent $p-n$ mixed solutions has been investigated with and without the Coulomb interaction, which is the sole source of isospin symmetry breaking in our approach. The extended axial HFBTHO solver has been benchmarked against the symmetry-unrestricted HFODD code for deformed and spherical states. We developed and tested a general isospin-invariant Skyrme-EDF framework. The new approach permits spin-isospin densities that may give rise to hitherto unexplored modes in the excitation spectrum. The new formalism has been tested in the Hartree-Fock limit. A systematic comparison between HFODD and HFBTHO results show a maximum deviation of about 10 keV on the total binding energy for deformed nuclei when the Coulomb term is included. Without this term, the results of both solvers agree down to a ~10 eV level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Skyrme pseudopotentials at next-to-next-to-leading order: Construction of local densities and first symmetry-breaking calculations

There is an ongoing quest to improve on the spectroscopic quality of nuclear energy density functionals (EDFs) of the Skyrme type through extensions of its traditional form. One direction for such activities is the inclusion of terms of higher order in gradients in the EDF. We report on exploratory symmetry-breaking calculations performed for an extension of the Skyrme EDF that includes central terms with four gradients at next-to-next-to-leading order (N2LO) and for which the high-quality parametrization SN2LO1 has been constructed recently. Up to now, the investigation of such functionals with higher-order terms was limited to infinite matter and spherically symmetric configurations of singly and doubly magic nuclei. We address here nuclei and phenomena that require us to consider axial and nonaxial deformation, both for reflection-symmetric and also reflection-asymmetric shapes, as well as the breaking of time-reversal invariance. Achieving these calculations demanded a number of formal developments. These all resulted from the formulation of the N2LO EDF requiring the introduction of new local densities with additional gradients that are not present in the EDF at NLO. Their choice is not unique, but can differ in the way the gradients are coupled. While designing a numerical implementation of N2LO EDFs in Cartesian three-dimensional coordinate-space representation, we have developed a novel definition and a new unifying notation for normal and pair densities that contain gradients at arbitrary order. Besides having mnemonic advantages, the new notation allows for the easy identification of redundancies and reducibilities in a given set of local densities, and the new definition makes it straightforward to construct densities that automatically adopt the symmetries of the many-body state they are constructed from. The resulting scheme resolves several issues with some of the choices that have been made for local densities in the past, in particular when breaking time-reversal symmetry. Guided by general practical considerations, we propose an alternative form of the N2LO contribution to the Skyrme EDF that is built from a different set of densities. It has exactly the same physics content, but is much more efficient to handle in formal discussions and, compared to the original formulation, leads to a substantial reduction of computational cost and memory requirements in deformed codes. As representative examples for the performance of SN2LO1, we have chosen the ground states of even-even Kr and Nd isotopes, the fission barrier of 240 Pu as well as the superdeformed rotational band of 194 Hg. Overall, for the nuclei and phenomena studied here, the SN2LO1 parametrization does not yet present a systematic improvement over standard NLO parametrizations. This finding calls for improved fit protocols that better discriminate between NLO and N2LO terms and better exploit the unique features of the additional degrees of freedom offered by the latter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

CLASH-VLT: constraints on f (R) gravity models with galaxy clusters using lensing and kinematic analyses

We perform a maximum likelihood kinematic analysis of the two dynamically relaxed galaxy clusters MACS J1206.2-0847 at z =0.44 and RXC J2248.7-4431 at z =0.35 to determine the total mass profile in modified gravity models, using a modified version of the MAMPOSSt code of Mamon, Biviano and Bou and apos;e. Our work is based on the kinematic and lensing mass profiles derived using the data from the Cluster Lensing And Supernova survey with Hubble (hereafter CLASH) and the spectroscopic follow-up with the Very Large Telescope (hereafter CLASH-VLT). We assume a spherical Navarro-Frenk-White (NFW hereafter) profile in order to obtain a constraint on the fifth force interaction range λ for models in which the dependence of this parameter on the environment is negligible at the scale considered (i.e. λ= const ) and fixing the fifth force strength to the value predicted in f (R) gravity. We then use information from lensing analysis to put a prior on the other NFW free parameters. In the case of MACSJ 1206 the joint kinematic+lensing analysis leads to an upper limit on the effective interaction range λ≤1.61 mpc at Δχ{sup 2}=2.71 on the marginalized distribution. For RXJ 2248 instead a possible tension with the ΛCDM model appears when adding lensing information, with a lower limit λ≥0.14 mpc at Δχ{sup 2}=2.71. This is consequence of the slight difference between the lensing and kinematic data, appearing in GR for this cluster, that could in principle be explained in terms of modifications of gravity. We discuss the impact of systematics and the limits of our analysis as well as future improvements of the results obtained. This work has interesting implications in view of upcoming and future large imaging and spectroscopic surveys, that will deliver lensing and kinematic mass reconstruction for a large number of galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Pn-FEM Method for Solving the First-Order Neutron Transport Equation

A deterministic method was investigated and developed for the photon-neutron transport calculation for fusion reactor blanket. Firstly, considering the fact that the second order transport equation cannot describe void region accurately due to the appearance of zero total cross section in the denominator, the first order equation was selected as the starting point of the solution. The angular discretization uses the spherical harmonics (Pn) method. Although discrete ordinates (Sn) method has been most commonly applied due to its efficiency, it suffers from the ray effect. To describe unstructured geometry, the finite element method was used to discrete spatial variable. It is well known that Bubnov-Galerkin discretization of the first-order equation is inappropriate which will create oscillations. To obtain robust solutions, Petrov- Galerkin methods have been developed such as the Streamline-Upwind Petrov-Galerkin (SUPG) method and the Generalized Least-Squares (GLS) method. These methods require a stabilization term added to its test functions, which will incur additional computational cost. Another finite element formulation that usually used for first order equation is the Discontinuous-Galerkin (DG) method. In this method, each element has a unique and separate set of basis function. Then the adjacent elements can be connected by using the in-coming flux. Its main drawback is the number of unknowns, much larger than Continuous-Galerkin (CG) method. The Sub-Grid Scale (SGS) finite element method combines the advantages of CG and DG methods together. It inserts the DG method in each element to correct the CG solution, and provides numerically stable solution of first-order equation without increasing the number of unknowns. The boundary conditions and interface conditions for Pn method is a hard nut to crack. We use Riemann method to deal with this problem. It decomposes the neutron angular flux into in-coming and out-going information through an interface, and applies boundary condition to the in-coming one. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2016

Effect of reconstruction methods and x-ray tube current–time product on nodule detection in an anthropomorphic thorax phantom: A crossed-modality JAFROC observer study

Purpose: To evaluate nodule detection in an anthropomorphic chest phantom in computed tomography (CT) images reconstructed with adaptive iterative dose reduction 3D (AIDR{sup 3D}) and filtered back projection (FBP) over a range of tube current–time product (mAs). Methods: Two phantoms were used in this study: (i) an anthropomorphic chest phantom was loaded with spherical simulated nodules of 5, 8, 10, and 12 mm in diameter and +100, −630, and −800 Hounsfield units electron density; this would generate CT images for the observer study; (ii) a whole-body dosimetry verification phantom was used to ultimately estimate effective dose and risk according to the model of the BEIR VII committee. Both phantoms were scanned over a mAs range (10, 20, 30, and 40), while all other acquisition parameters remained constant. Images were reconstructed with both AIDR{sup 3D} and FBP. For the observer study, 34 normal cases (no nodules) and 34 abnormal cases (containing 1–3 nodules, mean 1.35 ± 0.54) were chosen. Eleven observers evaluated images from all mAs and reconstruction methods under the free-response paradigm. A crossed-modality jackknife alternative free-response operating characteristic (JAFROC) analysis method was developed for data analysis, averaging data over the two factors influencing nodule detection in this study: mAs and image reconstruction (AIDR{sup 3D} or FBP). A Bonferroni correction was applied and the threshold for declaring significance was set at 0.025 to maintain the overall probability of Type I error at α = 0.05. Contrast-to-noise (CNR) was also measured for all nodules and evaluated by a linear least squares analysis. Results: For random-reader fixed-case crossed-modality JAFROC analysis, there was no significant difference in nodule detection between AIDR{sup 3D} and FBP when data were averaged over mAs [F(1, 10) = 0.08, p = 0.789]. However, when data were averaged over reconstruction methods, a significant difference was seen between multiple pairs of mAs settings [F(3, 30) = 15.96, p < 0.001]. Measurements of effective dose and effective risk showed the expected linear dependence on mAs. Nodule CNR was statistically higher for simulated nodules on images reconstructed with AIDR{sup 3D} (p < 0.001). Conclusions: No significant difference in nodule detection performance was demonstrated between images reconstructed with FBP and AIDR{sup 3D}. mAs was found to influence nodule detection, though further work is required for dose optimization.

60 APPLIED LIFE SCIENCES↗