Search NASA⌕ Search

SEARCH · Search NASA

Results for “Embedded boundary”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Investigation of the Rocket Induced Flow Field in a Rectangular Duct

Several tests were performed on a one-sixth scale Rocket Based Combined Cycle (RBCC) engine model at the University of Alabama in Huntsville. The UAH RBCC facility consists of a rectangular duct with a vertical strut mounted in the center. The scaled strut consists of two supersonic rocket nozzles with an embedded vertical turbine between the rocket nozzles. The tests included mass flow, flow visualization and horizontal pressure traverses. The mass flow test indicated a c:hoked condition when the rocket chamber pressure is between 200 psi and 300 psi. The flow visualization tests narrowed the rocket chamber pressure range from, 250 psi to 300 psi. Also, from this t.est, an assumption of a minimum <area choke was made. Horizontal pressure traverses were performed at locations both upstream and downstream of the nozzle exit plane. These traverses aided in examining the induced airflow due to the ejector effect and the growth of the: boundary layer in the strut sidewall gap.

Landrum, D. Brian↗
NASA NTRS1972

Theoretical study of stress concentrations at circular holes and inclusions in strain hardening materials.

Nonlinear boundary value problems of an infinite elastic-plastic plate with a circular hole subjected to pure tension and pure shear at infinity are solved by a method involving Fourier series and finite difference. On the basis of these solutions, the validity of Neuber's relationship between the stress and strain concentration factors for the plane stress problems is examined and a generalized Stowell formula for the stress concentration factor is proposed for problems in which the applied loading may be pure shear as well as pure tension and, furthermore, other stress states. By the same method of solution, the stress distributions around a rigid circular cylindrical inclusion embedded in an infinite rigid-plastic matrix subjected to uniform transverse pure shear and tension are obtained.

Huang, W. C.↗
NASA NTRS2014

Using a Commercial Ethernet PHY Device in a Radiation Environment

This work involved placing a commercial Ethernet PHY on its own power boundary, with limited current supply, and providing detection methods to determine when the device is not operating and when it needs either a reset or power-cycle. The device must be radiation-tested and free of destructive latchup errors. The commercial Ethernet PHY's own power boundary must be supplied by a current-limited power regulator that must have an enable (for power cycling), and its maximum power output must not exceed the PHY's input requirements, thus preventing damage to the device. A regulator with configurable output limits and short-circuit protection (such as the RHFL4913, rad hard positive voltage regulator family) is ideal. This will prevent a catastrophic failure due to radiation (such as a short between the commercial device's power and ground) from taking down the board's main power. Logic provided on the board will detect errors in the PHY. An FPGA (field-programmable gate array) with embedded Ethernet MAC (Media Access Control) will work well. The error detection includes monitoring the PHY's interrupt line, and the status of the Ethernet's switched power. When the PHY is determined to be non-functional, the logic device resets the PHY, which will often clear radiation induced errors. If this doesn't work, the logic device power-cycles the FPGA by toggling the regulator's enable input. This should clear almost all radiation induced errors provided the device is not latched up.

Parks, Jeremy↗
NASA NTRS2006

Plasma Circulation in the Magnetosphere

We investigate the global structure and dynamics of plasma circulation produced by prototypical solar wind disturbances of the interplanetary magnetic field and dynamic pressure. We track the global circulation and energization of solar wind, polar wind, and auroral wind plasmas throughout the magnetosphere, until they precipitate or escape into the downstream solar wind. We use the full equations of motion of the plasma ions within fields produced by a global MHD simulation of the dynamic solar wind interaction. We use the dynamic hot plasma density and Poynting energy flux specified at the inner boundary of the MHD simulation as drivers of conjugate ion outflow fluxes using local empirical relations obtained from the FAST and Polar missions. Birkeland currents computed by the MHD code are used to derive a field-parallel potential drop from a Knight-like relation [as modified by Lyons and Evans, 1980]. This potential drop is applied to each ion as an initial bulk energy, added to a thermal heating driven by the locally incident Poynting flux. The solar wind pressure increase case (B(sub Y) = 5; B(sub z) = 0 nT) produces an immediate substorm owing to compression of pre-existing plasmas. The SB(sub z), interval (embedded in NB(sub z)) produces a substorm after about one hour of development. Both disturbances enhance the auroral wind flux and heavy ion pressure of the magnetosphere substantially, with complex dynamic structuring by auroral acceleration vortexes and dynamic reconnection. Comparisons are made with observations during disturbed periods including the Halloween 2003 super-storm and other periods.

Moore, T. E.↗
NASA NTRS2011

Crater Flux Transfer Events: Highroad to the X Line?

We examine Cluster observations of a so-called magnetosphere crater FTE, employing data from five instruments (FGM, CIS, EDI, EFW, and WHISPER), some at the highest resolution. The aim of doing this is to deepen our understanding of the reconnection nature of these events by applying recent advances in the theory of collisionless reconnection and in detailed observational work. Our data support the hypothesis of a stratified structure with regions which we show to be spatial structures. We support the bulge-like topology of the core region (R3) made up of plasma jetting transverse to reconnected field lines. We document encounters with a magnetic separatrix as a thin layer embedded in the region (R2) just outside the bulge, where the speed of the protons flowing approximately parallel to the field maximizes: (1) short (fraction of a sec) bursts of enhanced electric field strengths (up to approximately 30 mV/m) and (2) electrons flowing against the field toward the X line at approximately the same time as the bursts of intense electric fields. R2 also contains a density decrease concomitant with an enhanced magnetic field strength. At its interface with the core region, R3, electric field activity ceases abruptly. The accelerated plasma flow profile has a catenary shape consisting of beams parallel to the field in R2 close to the R2/R3 boundary and slower jets moving across the magnetic field within the bulge region. We detail commonalities our observations of crater FTEs have with reconnection structures in other scenarios. We suggest that in view of these properties and their frequency of occurrence, crater FTEs are ideal places to study processes at the separatrices, key regions in magnetic reconnection. This is a good preparation for the MMS mission.

Farrugia, C. J.↗
NASA NTRS2016

Aerosol Meteorology of Maritime Continent for the 2012 7SEAS Southwest Monsoon Intensive Study - Part 2: Philippine Receptor Observations of Fine-Scale Aerosol Behavior

The largest 7 Southeast Asian Studies (7SEAS) operations period within the Maritime Continent (MC) occurred in the August-September 2012 biomass burning season. Data included were observations aboard the MY Vasco, dispatched to the Palawan Archipelago and Sulu Sea of the Philippines for September 2012. At these locations, the Vasco observed MC smoke and pollution entering the southwest monsoon (SWM) monsoonal trough. Here we describe the research cruise findings and the finer-scale aerosol meteorology of this convectively active region. This 2012 cruise complemented a 2-week cruise in 2011 and was generally consistent with previous findings in terms of how smoke emission and transport related to monsoonal flows, tropical cyclones (TC), and the covariance between smoke transport events and the atmosphere's thermodynamic structure. Biomass burning plumes were usually mixed with significant amounts of anthropogenic pollution. Also key to aerosol behavior were squall lines and cold pools propagating across the South China Sea (SCS) and scavenging aerosol particles in their path. However, the 2012 cruise showed much higher modulation in aerosol frequency than its 2011 counterpart. Whereas in 2011 large synoptic-scale aerosol events transported high concentrations of smoke into the Philippines over days, in 2012 measured aerosol events exhibited a much shorter-term variation, sometimes only 312h. Strong monsoonal flow reversals were also experienced in 2012. Nucleation events in cleaner and polluted conditions, as well as in urban plumes, were observed. Perhaps most interestingly, several cases of squall lines preceding major aerosol events were observed, as opposed to 2011 observations where these lines largely scavenged aerosol particles from the marine boundary layer. Combined, these observations indicate pockets of high and low particle counts that are not uncommon in the region. These perturbations are difficult to observe by satellite and very difficult to model. Indeed, the Navy Aerosol Analysis and Prediction System (NAAPS) simulations captured longer period aerosol events quite well but largely failed to capture the timing of high-frequency phenomena. Ultimately, the research findings of these cruises demonstrate the real world challenges of satellite-based missions, significant aerosol life cycle questions such as those the future Aerosol/Clouds/Ecosystems (ACE) will investigate, and the importance of small-scale phenomena such as sea breezes, squall lines, and nucleation events embedded within SWM patterns in dominating aerosol life cycle and potential relationships to clouds.

Reid, Jeffrey S.↗
NASA NTRS1993

The Sm-Nd systematics of silicate inclusions in iron meteorites: Results from Caddo (IAB)

The timing of events leading to the formation of silicate-rich and metal-rich regions in planetesimals remains an important problem in the study of planetary formation and differentiation in the early solar system. The IAB irons are especially important as they are considered to represent a magmatic differentiation series. Iron meteorites present a particular challenge for chronological studies, due to the relative paucity of phases serving as hosts for radioactive parent-daughter nuclides. Recent work using the Re-Os system, following on the pioneering work by Herr et al. and Luck and Allegre, appears promising, but investigators up to now have concentrated on whole rock isochrons. Silicate clasts enclosed within iron meteorites can provide information about the chronology and thermal history of irons. Extensive work on Rb-Sr, K-Ar, and I-Xe has been reported on silicate inclusions in iron meteorites. We report the initial results from our Sm-Nd study of an inclusion with the Caddo IAB iron, the first Sm-Nd isotopic study of a silicate clast embedded within an iron meteorite. Our results include measurements of the standard long-lived Sm-147/Nd-143 (tau = 152 AE) system, as well as the shorter-lived SM-146/Nd-142 (tau = 0.149 AE) system, which has been shown to be very useful in deciphering the history of the early solar system. The Caddo silicate clast was described by Palme et al., who kindly provided us with a major part of the inclusion. The inclusion is coarse-grained consisting predominantly of olivine, clinopyroxene, and plagioclase, with lesser amounts of orthopyroxene, Fe-Ni metal, sulfide, and phosphate. The relatively large grain size (up to 3 mm) and 120 degree grain boundaries suggest extensive metamorphism at high temperatures. Based on study of a thin section, there is evidence for metal invading along grain boundaries in some regions of the inclusion, suggesting that the Fe-Ni metal was molten when the silicate clast was incorporated. Metamorphic recrystallization may have occurred during this event.

Stewart, Brian W.↗
NASA NTRS1991

Application of the boundary element method to transient heat conduction

An advanced boundary element method (BEM) is presented for the transient heat conduction analysis of engineering components. The numerical implementation necessarily includes higher-order conforming elements, self-adaptive integration and a multiregion capability. Planar, three-dimensional and axisymmetric analyses are all addressed with a consistent time-domain convolution approach, which completely eliminates the need for volume discretization for most practical analyses. The resulting general purpose algorithm establishes BEM as an attractive alternative to the more familiar finite difference and finite element methods for this class of problems. Several detailed numerical examples are included to emphasize the accuracy, stability and generality of the present BEM. Furthermore, a new efficient treatment is introduced for bodies with embedded holes. This development provides a powerful analytical tool for transient solutions of components, such as casting moulds and turbine blades, which are cumbersome to model when employing the conventional domain-based methods.

Dargush, G. F.↗
DOE OSTI2024

Nonsteady Load Responses of Wind Turbines to Atmospheric and Mountain-Generated Turbulence Eddies, With Impacts on the Main Bearing: A Validation Study

Previous computational and field experiments identify three characteristic time scales in the aerodynamic responses of utility-scale wind turbine loads to atmospheric boundary layer (ABL) turbulence: a 30-90 second time scale for the passage of high/low speed "streaks" through the rotor plane, the blade and rotor rotation time scales (approximately 1 to 5 seconds), and a sub-second time scale created by blade rotation through gradients within eddy coherent structure. In the current study we compare aerodynamic load responses from daytime ABL turbulence quantified with large-eddy simulation and a actuator line model of the NREL 5 MW wind turbine with analysis of field data from the NREL/GE 1.5 MW wind turbine 5 kilometers east of the Rocky Mountain Front Range in Colorado. In addition, we contrast the responses to the passage of the mountain-generated eddies embedded within the westerly winds with the ABL eddies embedded within northerly/southerly winds. These analyses are in context with the nonsteady forcing of the main bearing by the aerodynamic generation of nontorque bending moments on the main shaft. Potentially relevant to main bearing failure mechanisms, both computational and field data show that the magnitudes of turbulence-generated nontorque bending moments, that we show generate nonsteady force on the main bearing, are of order, and often larger than, torque (which underlies power). However, the temporal variations in these two responses are uncorrelated, implying that the aerodynamic mechanisms that drive power and main bearing response are fundamentally different. We find this to be the case in the field with both mountain-generated eddies (westerly winds) and ABL-generated eddies (northerly/southerly winds). Whereas the time and length scales are comparable, the mountain eddies were somewhat more energetic than the northerly/southerly ABL eddies. Interestingly, however, the fluctuations in nontorque bending moment that force the main bearing were found to be stronger when forced by the ABL eddies than the mountain eddies. The field studies validate the key results from the computational study and show even stronger response in the nontorque bending moment than in the computer simulations. In all cases, the torque and nontorque bending moments are temporally uncorrelated, torque and power are driven by time variations in rotor-averaged horizontal wind velocity and nontorque bending moments are driven by time changes in the degree of nonuniformity in the distribution of velocity over the rotor plane. Thus the results generalize the mechanisms underlying nonsteady aerodynamic forcing to classes of turbulence eddy types with strength of order or stronger than ABL eddies with transverse scale of order the wind turbine rotor. These include atmospheric turbulence eddies, topography-generated turbulence eddies and, by extension, impacts of turbine-wake-scale turbulence eddies on downstream wind turbine rotors.

17 WIND ENERGY↗
NASA NTRS1991

Lunar regolith densification

Core tube samples of the lunar regolith obtained during the Apollo missions showed a rapid increase in the density of the regolith with depth. Various hypotheses have been proposed for the possible cause of this phenomenon, including the densification of the loose regolith material by repeated shaking from the seismic tremors which have been found to occur at regular monthly intervals when the moon and earth are closest to one another. A test bed was designed to study regolith densification. This test bed uses Minnesota Lunar Simulant (MLS) to conduct shaking experiments in the geotechnical centrifuge with an inflight shake table system. By reproducing realistic in-situ regolith properties, the experiment also serves to test penetrator concepts. The shake table system was designed and used for simulation experiments to study effects of earthquakes on terrestrial soil structures. It is mounted on a 15 g-ton geotechnical centrifuge in which the self-weight induced stresses are replicated by testing an n-th scale model in a gravity field which is n times larger than Earth's gravity. A similar concept applies when dealing with lunar prototypes, where the gravity ratio required for proper simulation of lunar gravity effects is that between the centrifugal acceleration and the lunar gravity. Records of lunar seismic tremors, or moonquakes, were obtained. While these records are being prepared for use as the input data to drive the shake table system, records from the El Centro earthquake of 1940 are being used to perform preliminary tests, using a soil container which was previously used for earthquake studies. This container has a laminar construction, with the layers free to slide on each other, so that the soil motion during the simulated earthquake will not be constrained by the otherwise rigid boundaries. The soil model is prepared by pluviating the MLS from a hopper into the laminar container to a depth of 6 in. The container is mounted on the shake table and the centrifuge is operated to generate an acceleration of 10 times Earth's gravity or 60 times the lunar gravity, thus simulating a lunar regolith thickness of 30 ft. The shake table is then operated using the scaled 'moonquake' as the input motion. One or more model moonquakes are used in each experiment, after which the soil is analyzed for its density profile with depth. This is accomplished by removing from the soil bed a column of soil contained within a thin rubber sleeve which has been previously embedded vertically in the soil during pluviation. This column of soil is transferred to a gamma ray device, in which the gamma ray transmission transversely through the soil is measured and compared with standard calibration samples. In this manner, the density profile can be determined. Preliminary results to date are encouraging, and the Center plans to study the effects of duration of shaking, intensity of the shaking motion, and the frequency of the motion.

Ko, Hon-Yim↗
NASA NTRS1999

Flow Trajectories in Friction Stir Welding

In the Friction Stir Welding (FSW) process, a rotating, shouldered tool with a threaded pin is inserted under pressure along the seam of two pieces of metal which are tightly clamped together, and secured against a rigid anvil underneath. The rotating pin travels along the seam and through a combination of pressure and friction heating produces a zone of plastic deformation around the pin within the workpiece on either side of the seam. As the pin is moved in the direction of welding, the plasticised material moves around the tool and bonds together behind it. The elements of the material flow behavior are a combination of three elements. There is a rotational transport of material being carried around the tool, extrusion of material being forced around the pin on both sides into the cavity created behind it, and a lifting and dropping of material as it is stirred and mixed by the rotating action of the pin. It was assumed that rotational motion of the plastic zone is the primary mechanism for transport of material around the welding tool. A kinematic mathematical model was used to compute trajectories of material movement for various distributions of rotational slip within the plastic zone. These trajectories were then compared with the results of an experiment that produced radiographs of markers embedded in a workpiece that was welded with the FSW process. It was assumed that the copper wire markers retained their original length as the aluminum material flowed around them. The kinematic model included a constraint so that the displacements were such that the total length of the wire markers did not increase after deformation. There was good agreement between the calculated trajectories for the case of localized slip at the outer surface of the plastic deformation zone and the radiographs of the copper wire markers. The trajectories differed markedly from the radiograph traces when a distributed slip zone was assumed. It was concluded that the flow field could be described generally by a rotating plug model with slip occurring on the surface at the outer boundary of the plastic zone.

Bernstein, Edward L.↗
NASA NTRS1989

Results of the first statistical study of Pioneer Venus Orbiter plasma observations in the distant Venus tail - Evidence for a hemispheric asymmetry in the pickup of ionospheric ions

Pioneer Venus Orbiter plasma and magnetometer observations from the first nine tail seasons of crossings of the Venus wake are used to study ion pickup in the far wake of an unmagnetized object embedded in the solar wind. This statistical study treats all of the plasma spectra containing pickup ions in the vicinity of the Venus tail. A hemispheric asymmetry is found in the pickup of ionospheric ions, with approximately four times more O(+) events observed in the 'northern' magnetic hemisphere (where the induced electric field points outward) than in the 'southern' magnetic hemisphere. Out of a total of 167 large O(+) events, 125, or 75 percent, occurred in the northern hemisphere when position is calculated in terms of Venus radii, and 129 or 77 percent occurred in the 'northern' hemisphere when position is expressed in gyroradii. This hemispheric asymmetry in ion pickup is consistent with the prediction of the Cloutier et al. (1974) mass loading model for Venusian ions above the ionopause boundary.

Intriligator, Devrie S.↗
NASA NTRS2020

Evaluating Precipitation Features and Rainfall Characteristics in a Multi‐Scale Modeling Framework

Cloud and precipitation systems are simulated with a multi‐scale modeling framework (MMF) and compared over the Tropics and Subtropics against the Tropical Rainfall Measuring Mission (TRMM) Radar‐defined Precipitation Features (RPFs) product. A methodology, in close analogy to the TRMM RPFs, is developed to produce simulated precipitation features (PFs) from the output of the embedded two‐dimensional (2D) cloud‐resolving models (CRMs) within an MMF. Despite the limitations of 2D CRMs, the simulated population distribution, horizontal and vertical structure of PFs, and the geographical location and local rainfall contribution of mesoscale convective systems (MCSs) are in good agreement with the TRMM observations. However, some model discrepancies are found and can be identified and quantified within the PF distributions. Using model biases in relative population and rainfall contributions, PFs can be characterized into four size categories: small, medium to large, very large, and extremely large. Four different major mechanisms might account for the model biases in each different category: (1) the two‐dimensionality of the CRMs, (2) a positive convection‐wind‐evaporation feedback loop, (3) an artificial dynamic constraint in a bounded CRM domain with cyclic boundaries, and (4) the limited CRM domain size. The second and fourth mechanisms tend to contribute to the excessive tropical precipitation biases commonly found in most MMFs, whereas the other mechanisms reduce rainfall contributions from small and very large PFs. MMF sensitivity experiments with various CRM domain sizes and grid spacings showed that larger domains (higher resolutions) tend to shift PF populations toward larger (smaller) sizes.

Jiun-Dar Chern↗
DOE OSTI2024

Elucidating the role of Cr migration in Ni-Cr exposed to molten FLiNaK via multiscale characterization

Structural materials used in nuclear reactor environments are subjected to coupled extremes such as high temperature, irradiation, and corrosion which act in concert to degrade their functional performance. Connecting alloy microstructure such as grain boundaries, and accumulating point defects with corrosion attack and pore morphology is critical to understanding underlying mechanisms. We uncover the compositional variations and morphology at multiple length scales in corrosion-damaged Ni-Cr alloy after exposure to oxidants in molten fluoride salts. A complex network of dense corrosion pores is detected by surface-level SEM observations. The corrosion pores take on a 1-dimensional morphology and are enriched with Ni and depleted of Cr 1–5 µm from the pore surface. STEM-EDS and 4D-STEM strain maps acquired simultaneously highlight the local variations in composition and structure of a ≤ 200 nm Cr-rich layer identified from a cross-section taken at the bottom of an isolated corrosion pore between the Ni-Cr alloy matrix and the embedded salt. However, the absence of an observed interface between the Ni-Cr alloy matrix and the FCC-structured Cr-rich layer suggests that Cr plating from the salt did not transpire. These findings support a proposed Cr lattice diffusion mechanism rather than Cr-precipitation from the salt to accommodate temperature transient conditions during sample cooling. Through the development of a 1D phase field model, these results are rationalized by formation energies for the Ni- and Cr-oxidation into the molten salt. This study reveals the locally altered microstructure caused by high temperature corrosion in non-steady-state molten salt nuclear reactor environments.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Convection and electrodynamic signatures in the vicinity of a Sun-aligned arc: Results from the Polar Acceleration Regions and Convection Study (Polar ARCS)

An experimental campaign designed to study high-latitude auroral arcs was conducted in Sondre Stromfjord, Greenland, on February 26, 1987. The Polar Acceleration Regions and Convection Study (Polar ARCS) consisted of a coordinated set of ground-based, airborne, and sounding rocket measurements of a weak, sun-aligned arc system within the duskside polar cap. A rocket-borne barium release experiment, two DMSP satellite overflights, all-sky photography, and incoherent scatter radar measurements provided information on the large-scale plasma convection over the polar cap region while a second rocket instrumented with a DC magnetometer, Langmuir and electric field probes, and an electron spectrometer provided measurements of small-scale electrodynamics. The large-scale data indicate that small, sun-aligned precipitation events formed within a region of antisunward convection between the duskside auroral oval and a large sun-aligned arc further poleward. This convection signature, used to assess the relationship of the sun-aligned arc to the large-scale magnetospheric configuration, is found to be consistent with either a model in which the arc formed on open field lines on the dusk side of a bifurcated polar cap or on closed field lines threading an expanded low-latitude boundary layer, but not a model in which the polar cap arc field lines map to an expanded plasma sheet. The antisunward convection signature may also be explained by a model in which the polar cap arc formed on long field lines recently reconnected through a highly skewed plasma sheet. The small-scale measurements indicate the rocket passed through three narrow (less than 20 km) regions of low-energy (less than 100 eV) electron precipitation in which the electric and magnetic field perturbations were well correlated. These precipitation events are shown to be associated with regions of downward Poynting flux and small-scale upward and downward field-aligned currents of 1-2 micro-A/sq m. The paired field-aligned currents are associated with velocity shears (higher and lower speed streams) embedded in the region of antisunward flow.

Weiss, L. A.↗
NASA NTRS1995

Formation Of the Giant Planets By Concurrent Accretion Of Solids And Gas

New numerical simulations of the formation of the giant planets are presented, in which for the first time both the gas and planetesimal accretion rates are calculated in a self-consistent, interactive fashion. The simulations combine three elements: 1) three-body accretion cross-sections of solids onto an isolated planetary embryo, 2) a stellar evolution code for the planet's gaseous envelope, and 3) a planetesimal dissolution code within the envelope, used to evaluate the planet's effective capture radius and the energy deposition profile of accreted material. Major assumptions include: The planet is embedded in a disk of gas and small planetesimals with locally uniform initial surface mass density, and planetesimals are not allowed to migrate into or out of the planet's feeding zone. All simulations are characterized by three major phases. During the first phase, the planet's mass consists primarily of solid material. The planetesimal accretion rate, which dominates that of gas, rapidly increases owing to runaway accretion, then decreases as the planet's feeding zone is depleted. During the second phase, both solid and gas accretion rates are small and nearly independent of time. The third phase, marked by runaway gas accretion, starts when the solid and gas masses are about equal. It is engendered by a strong positive feedback on the gas accretion rates, driven by the rapid contraction of the gaseous envelope and the rapid expansion of the outer boundary, which depends on the planet's total mass. The overall evolutionary time scale is generally determined by the length of the second phase. The actual rates at which the giant planets accreted small planetesimals is probably intermediate between the constant rates assumed in most previous studies and the highly variable rates that we have used. Within the context, of the adopted model of planetesimal accretion, the joint constraints of the time scale for dissipation of the solar nebula and the current high-Z masses of the giant planets lead to estimates of the initial surface density (sigma(sub init)) of planetesimals in the outer region of the solar nebula. The results show sigma(sub init) approx. = 10 g/sq cm near Jupiter's orbit and that sigma(sub init) proportional to alpha(sup -2), where alpha is the distance from the Sun. These values are a factor of 3 - 4 times as high as that of the "minimum mass" solar nebula at Jupiter's distance and a factor of 2 - 3 times as high it Saturn's distance. Our estimates for the formation time of Jupiter and Saturn are 1 - 10 million years while those for Uranus fall in the range of 2 - 16 million years. These estimates follow from the properties of our Solar System and do not necessarily apply to giant planets in other planetary systems.

Pollack, James B.↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
NASA NTRS1998

Method and apparatus for reducing the drag of flows over surfaces

An apparatus, and its accompanying method, for reducing the drag of flows over a surface includes arrays of small disks and sensors. The arrays are embedded in the surface and may extend above, or be depressed below, the surface, provided they remain hydraulically smooth either when operating or when inactive. The disks are arranged in arrays of various shapes, and spaced according to the cruising speed of the vehicle on which the arrays are installed. For drag reduction at speeds of the order of 30 meters/second, preferred embodiments include disks that are 0.2 millimeter in diameter and spaced 0.4 millimeter apart. For drag reduction at speeds of the order of 300 meters/second, preferred embodiments include disks that are 0.045 millimeter in diameter and spaced 0.09 millimeter apart. Smaller and larger dimensions for diameter and spacing are also possible. The disks rotate in the plane of the surface, with their rotation axis substantially perpendicular to the surface. The rotating disks produce velocity perturbations parallel to the surface in the overlying boundary layer. The sensors sense the flow at the surface and connect to control circuitry that adjusts the rotation rates and duty cycles of the disks accordingly. Suction and blowing holes can be interspersed among, or made coaxial with, the disks for creating general three-component velocity perturbations in the near-surface region. The surface can be a flat, planar surface or a nonplanar surface, such as a triangular riblet surface. The present apparatus and method have potential applications in the field of aeronautics for improving performance and efficiency of commercial and military aircraft, and in other industries where drag is an obstacle, including gas and oil delivery through long-haul pipelines.

Keefe, Laurence R.↗