Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical core”

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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

The spatial extent and dynamics of the Antarctic Cold Reversal

Antarctic ice cores show that a millennial-scale cooling event, the Antarctic Cold Reversal (14,700 to 13,000 years ago), interrupted the last deglaciation1,2,3. The Antarctic Cold Reversal coincides with the Bølling–Allerød warm stage in the North Atlantic, providing an example of the inter-hemispheric coupling of abrupt climate change generally referred to as the bipolar seesaw4,5,6,7,8,9. However, the ocean–atmosphere dynamics governing this coupling are debated10,11,12,13,14,15. Here we examine the extent and expression of the Antarctic Cold Reversal in the Southern Hemisphere using a synthesis of 84 palaeoclimate records. We find that the cooling is strongest in the South Atlantic and all regions south of 40° S. At the same time, the terrestrial tropics and subtropics show abrupt hydrologic variations that are significantly correlated with North Atlantic climate changes. Our transient global climate model simulations indicate that the observed extent of Antarctic Cold Reversal cooling can be explained by enhanced northward ocean heat transport from the South to North Atlantic10, amplified by the expansion and thickening of sea ice in the Southern Ocean. The hydrologic variations at lower latitudes result from an opposing enhancement of southward heat transport in the atmosphere mediated by the Hadley circulation. Our findings reconcile previous arguments about the relative dominance of ocean5,10,11 and atmospheric14,15 heat transports in inter-hemispheric coupling, demonstrating that the spatial pattern of past millennial-scale climate change reflects the superposition of both.

Geology↗
NASA NTRS1990

Lateral temperature variations at the core-mantle boundary deduced from the magnetic field

Recent studies of the secular variation of the earth's magnetic field over periods of a few centuries have suggested that the pattern of fluid motion near the surface of earth's outer core may be strongly influenced by lateral temperature variations in the lowermost mantle. This paper introduces a self-consistent method for finding the temperature variations near the core surface by assuming that the dynamical balance there is geostrophic and that lateral density variations there are thermal in origin. As expected, the lateral temperature variations are very small. Some agreement is found between this pattern and the pattern of topography of the core-mantle boundary, but this does not conclusively answer to what extent core surface motions are controlled by the mantle, rather than being determined by processes in the core.

Bloxham, Jeremy↗
DOE OSTI2014

Energy-aware Thread and Data Management in Heterogeneous Multi-core, Multi-memory Systems

By 2004, microprocessor design focused on multicore scaling—increasing the number of cores per die in each generation—as the primary strategy for improving performance. These multicore processors typically equip multiple memory subsystems to improve data throughput. In addition, these systems employ heterogeneous processors such as GPUs and heterogeneous memories like non-volatile memory to improve performance, capacity, and energy efficiency. With the increasing volume of hardware resources and system complexity caused by heterogeneity, future systems will require intelligent ways to manage hardware resources. Early research to improve performance and energy efficiency on heterogeneous, multi-core, multi-memory systems focused on tuning a single primitive or at best a few primitives in the systems. The key limitation of past efforts is their lack of a holistic approach to resource management that balances the tradeoff between performance and energy consumption. In addition, the shift from simple, homogeneous systems to these heterogeneous, multicore, multi-memory systems requires in-depth understanding of efficient resource management for scalable execution, including new models that capture the interchange between performance and energy, smarter resource management strategies, and novel low-level performance/energy tuning primitives and runtime systems. Tuning an application to control available resources efficiently has become a daunting challenge; managing resources in automation is still a dark art since the tradeoffs among programming, energy, and performance remain insufficiently understood. In this dissertation, I have developed theories, models, and resource management techniques to enable energy-efficient execution of parallel applications through thread and data management in these heterogeneous multi-core, multi-memory systems. I study the effect of dynamic concurrent throttling on the performance and energy of multi-core, non-uniform memory access (NUMA) systems. I use critical path analysis to quantify memory contention in the NUMA memory system and determine thread mappings. In addition, I implement a runtime system that combines concurrent throttling and a novel thread mapping algorithm to manage thread resources and improve energy efficient execution in multi-core, NUMA systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Properties of molecular clumps and cores in colliding magnetized flows

ABSTRACT We simulate the formation of molecular clouds in colliding flows of warm neutral medium with the adaptive mesh refinement code flash in eight simulations with varying initial magnetic field strength, between 0.01–5 μG. We include a chemical network to treat heating and cooling and to follow the formation of molecular gas. The initial magnetic field strength influences the fragmentation of the forming cloud because it prohibits motions perpendicular to the field direction and hence impacts the formation of large-scale filamentary structures. Molecular clump and core formation occurs anyhow. We identify 3D clumps and 3D cores, which are defined as connected, CO-rich regions. Additionally, 3D cores are heavily shielded. While we do not claim those 3D objects to be directly comparable to observations, this enables us to analyse their full virial state. With increasing field strength, we find more fragments with a smaller average mass; yet the dynamics of the forming clumps and cores only weakly depends on the initial magnetic field strength. The molecular clumps are mostly unbound, probably transient objects, which are weakly confined by ram pressure or thermal pressure, indicating that they are swept up by the turbulent flow. They experience significant fluctuations in the mass flux through their surface, such that the Eulerian reference frame shows a dominant time-dependent term due to their indistinct nature. We define the cores to encompass highly shielded molecular gas. Most cores are in gravitational-kinetic equipartition and are well described by the common virial parameter $\alpha _\mathrm{vir}$, while some undergo minor dispersion by kinetic surface effects.

Weis, M. (ORCID:0000000256838860)↗
DOE OSTI2022

From Protein Design to the Energy Landscape of a Cold Unfolding Protein

Understanding protein folding is crucial for protein sciences. The conformational spaces and energy landscapes of cold (unfolded) protein states, as well as the associated transitions, are hardly explored. Furthermore, it is not known how structure relates to the cooperativity of cold transitions, if cold and heat unfolded states are thermodynamically similar, and if cold states play important roles for protein function. We created the cold unfolding 4-helix bundle DCUB1 with a de novo designed bipartite hydrophilic/hydrophobic core featuring a hydrogen bond network which extends across the bundle in order to study the relative importance of hydrophobic versus hydrophilic protein–water interactions for cold unfolding. Structural and thermodynamic characterization resulted in the discovery of a complex energy landscape for cold transitions, while the heat unfolded state is a random coil. Below ~0 °C, the core of DCUB1 disintegrates in a largely cooperative manner, while a near-native helical content is retained. The resulting cold core-unfolded state is compact and features extensive internal dynamics. Below -5 °C, two additional cold transitions are seen, that is, (i) the formation of a water-mediated, compact, and highly dynamic dimer, and (ii) the onset of cold helix unfolding decoupled from cold core unfolding. Our results suggest that cold unfolding is initiated by the intrusion of water into the hydrophilic core network and that cooperativity can be tuned by varying the number of core hydrogen bond networks. Protein design has proven to be invaluable to explore the energy landscapes of cold states and to robustly test related theories.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

High pressure induced phase transition and superdiffusion in anomalous fluid confined in flexible nanopores

The behavior of a confined spherical symmetric anomalous fluid under high external pressure was studied with Molecular Dynamics simulations. The fluid is modeled by a core-softened potential with two characteristic length scales, which in bulk reproduces the dynamical, thermodynamical, and structural anomalous behavior observed for water and other anomalous fluids. Our findings show that this system has a superdiffusion regime for sufficient high pressure and low temperature. As well, our results indicate that this superdiffusive regime is strongly related with the fluid structural properties and the superdiffusion to diffusion transition is a first order phase transition. We show how the simulation time and statistics are important to obtain the correct dynamical behavior of the confined fluid. Our results are discussed on the basis of the two length scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Elastic scattering in the 12 N + 197 Au system at 𝐸 lab = 70 MeV

Halo structure is an interesting exotic configuration developed in some light weakly bound nuclei, where a valence particle orbits a nuclear core. Signatures of halo structure can be observed in the angular distributions of the elastic scattering induced by these nuclei at energies around the Coulomb barrier. There are some well-studied reactions with neutron-rich halo nuclei, such as 6 He and 11 Li . However, the information is scarce on the proton-rich side. Recent works confirm the halo structure in the 8B nuclei but still lack more experimental studies for other proton-halo candidates, such as 12 N and 17 F. In this work, we report experimental data for the elastic scattering of 12 N on 197 Au target at 𝐸 lab = 70 MeV. The 12 N is a proton-rich nucleus with proton separation energy 𝑆 𝑝 = 600 keV, which is higher than the 8 B (𝑆 𝑝 = 137 keV) and almost the same as the 17 F (𝑆 𝑝 = 601 keV). Data were obtained at the Cyclotron Institute of Texas A&M University where the 12 N radioactive beam was produced by the momentum achromatic recoil spectrometer. The calculation of the optical model was used to fit the measured angular distribution for the elastic scattering and to obtain the reaction cross section 𝜎 𝑅 . We also performed continuum discretized coupled-channel calculations to compare with the experimental data. The angular distribution of the elastic cross sections exhibits a suppression of the Fresnel peak. From the fitting of the optical model, we obtain the total reaction cross section, 𝜎 𝑅 = 1269 ± 41 mb. The agreement between the CDCC calculation and the experimental elastic scattering data is limited and the breakdown does not exhaust the measured 𝜎 𝑅 . The resulting reduced reaction cross section 𝜎 Red for 12 N + 197 Au is large and comparable to the one obtained for the 8 B + 208 Pb system. This suggests a strong decoupling of the valence proton from the core because of the low binding energy and a dynamic polarization effect. Here, more research is required to estimate the contribution of core excitation.

Elastic scattering reactions↗
DOE OSTI2022

Liquid Core Detection and Strand Condition Monitoring in a Continuous Caster Using Optical Fiber

Real-time monitoring of the liquid core position during the continuous casting of steel has been demonstrated using low-cost distributed optical-fiber-based strain sensors. These sensors were installed on the containment roll support structures in the segments of a production continuous caster to detect the position of the solid–liquid interface and monitor the strand condition during the continuous casting. Distributed Fiber Bragg Grating sensors (FBGs) were used in this work to monitor strain at six roll positions in the caster. The sensor performance was first validated by comparing optical strain measurements with conventional strain gauge measurements in the lab. Next, optical strain measurements were performed on an isolated caster segment in a segment maintenance facility using hydraulic jacks to simulate the presence of a liquid core under the roll. Finally, the sensors were evaluated during caster operation. The sensors successfully detected the load increase associated with the presence of a liquid core under each instrumented roll location. Incidents of bulging and roll eccentricity were also detected using frequency analysis of the optical strain signal. The liquid core position measurements were compared using predictions from computer models (digital twins) in use at the production site. The measurements were in good agreement with the model predictions, with a few exceptions. Under certain transient caster operating conditions, such as spraying practice changes and SEN exchanges, the model predictions deviated slightly from the liquid core position determined from strain measurements.

47 OTHER INSTRUMENTATION↗
DOE OSTI2014

THE BLUE STRAGGLER STAR POPULATION IN NGC 1261: EVIDENCE FOR A POST-CORE-COLLAPSE BOUNCE STATE

We present a multi-passband photometric study of the Blue Straggler Star (BSS) population in the Galactic globular cluster (GC) NGC 1261, using available space- and ground-based survey data. The inner BSS population is found to have two distinct sequences in the color-magnitude diagram (CMD), similar to double BSS sequences detected in other GCs. These well defined sequences are presumably linked to single short-lived events such as core collapse, which are expected to boost the formation of BSSs. In agreement with this, we find a BSS sequence in NGC 1261 which can be well reproduced individually by a theoretical model prediction of a 2 Gyr old population of stellar collision products, which are expected to form in the denser inner regions during short-lived core contraction phases. Additionally, we report the occurrence of a group of BSSs with unusually blue colors in the CMD, which are consistent with a corresponding model of a 200 Myr old population of stellar collision products. The properties of the NGC 1261 BSS populations, including their spatial distributions, suggest an advanced dynamical evolutionary state of the cluster, but the core of this GC does not show the classical signatures of core collapse. We argue that these apparent contradictions provide evidence for a post-core-collapse bounce state seen in dynamical simulations of old GCs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Molecular dynamics study of hydrogen Cottrell atmosphere in aluminum: Influence of solute-solute interactions in the dislocation core

Cottrell atmospheres form when solute atoms segregate to dislocations, and is one of the most basic processes by which solutes alter mechanical properties of materials. In the case of hydrogen, Cottrell atmospheres are believed to contribute to hydrogen embrittlement. Here, in this study, we use direct molecular dynamics simulations to study the formation of hydrogen Cottrell atmospheres around edge dislocations in aluminum. Using coarse-graining techniques, we resolve the temporal and spatial evolution of the hydrogen concentration field. By comparing our results with theories of Cottrell atmosphere formation and thermodynamics, we show that a repulsive, concentration-dependent H-H interaction reduces the peak atmosphere concentration at the dislocation core. Furthermore, this interaction differs significantly from the H-H interaction that is experienced in the bulk away from a dislocation.

Aluminum↗
DOE OSTI2021

Abscisic acid supports colonization of Eucalyptus grandis roots by the mutualistic ectomycorrhizal fungus Pisolithus microcarpus

Summary The pathways regulated in ectomycorrhizal (EcM) plant hosts during the establishment of symbiosis are not as well understood when compared to the functional stages of this mutualistic interaction. Our study used the EcM host Eucalyptus grandis to elucidate symbiosis‐regulated pathways across the three phases of this interaction. Using a combination of RNA sequencing and metabolomics we studied both stage‐specific and core responses of E. grandis during colonization by Pisolithus microcarpus . Using exogenous manipulation of the abscisic acid (ABA), we studied the role of this pathway during symbiosis establishment. Despite the mutualistic nature of this symbiosis, a large number of disease signalling TIR‐NBS‐LRR genes were induced. The transcriptional regulation in E. grandis was found to be dynamic across colonization with a small core of genes consistently regulated at all stages. Genes associated to the carotenoid/ABA pathway were found within this core and ABA concentrations increased during fungal integration into the root. Supplementation of ABA led to improved accommodation of P. microcarpus into E. grandis roots. The carotenoid pathway is a core response of an EcM host to its symbiont and highlights the need to understand the role of the stress hormone ABA in controlling host–EcM fungal interactions.

Hill, Richard A.↗
DOE OSTI2018

Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation

Tight oil formations such as those in the Bakken petroleum system are known to hold hundreds of billions of barrels of oil in place; however, the primary recovery factor for these plays is typically less than 10%. Tight oil formations, including the Bakken Formation, therefore, may be attractive candidates for enhanced oil recovery (EOR) using CO 2 . Multiphase fluid behavior and flow in fluid-rich shales can vary substantially depending on the size of pore throats, and properties such as fluid viscosity and density are much different in nanoscale pores than in macroscale pores. Thus it is critical to understand the nature and distribution of nano-, micro-, and macroscale pores and fracture networks. To address these issues, the Energy & Environmental Research Center (EERC) has been conducting a research program entitled “Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation.” The objectives of the project are 1) the use of advanced characterization methods to better understand and quantify the petrophysical and geomechanical factors that control CO 2 and oil mobility within tight oil formation samples, 2) the determination of CO 2 permeation and oil extraction rates in tight reservoir rocks and organic-rich shales of the Bakken, and 3) the integration of the laboratory-based CO 2 permeation and oil extraction data and the characterization data into geologic models and dynamic simulations to develop predictions of CO 2 storage resource and EOR in the Bakken tight oil formation. A combination of standard and advanced petrophysical characterization techniques were applied to characterize samples of Bakken Formation tight reservoir rock and shales from multiple wells. Techniques included advanced computer tomography (CT) imaging, scanning electron microscopy (SEM) techniques, whole-core and micro x-ray CT imaging, field emission (FE) SEM, and focused ion beam (FIB) SEM. Selected samples were also analyzed for geomechanical properties. X-ray CT imaging yielded information on the occurrence of fractures, bedding planes, fossils, and bioturbation in core, as well as data on bulk density and photoelectric factor logs, which were used to interpret porosity, organic content, and mineralogy. FESEM was used for characterization of nano- and microscale features, including nanoscale pore visualization and micropore and pore throat mineralogy. FIBSEM yielded micro- to nanoscale visualization of fracture networks, porosity and pore-size distribution, connected versus isolated porosity, and distribution of organics. Results from the characterization activities provide insight on nanoscale fracture properties, pore throat mineralogy and connectivity, rock matrix characteristics, mineralogy, and organic content. Laboratory experiments demonstrated that CO 2 can permeate the tight matrix of Bakken shale and nonshale reservoir samples and mobilize oil from those samples. Geologic models were created at scales ranging from the core plug to the reservoir, and dynamic simulations were conducted. The data from the characterization and laboratory-based activities were integrated into modeling research activities to determine the fundamental mechanisms controlling fluid transport in the Bakken, which support EOR scheme design and estimation of CO 2 storage potential in tight oil formations. Simulation results suggest a CO 2 storage resource estimate range of 169 million to 1.5 billion tonnes for the Bakken in North Dakota, possibly resulting in 1.8 billion to 16 billion barrels of incremental oil.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2016

Comparison of HFIR Core Hydraulic Models with Operational Data

The core flow channels of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL) have been modeled with three different computer codes and calculated results have been compared with available plant operational data. The three codes employed are RELAP5-3D, TRACE, and COMSOL. This work was performed to develop consistent boundary conditions for three-dimensional (3-D) computational fluid dynamics (CFD) models of the HFIR core. The variables of interest for the CFD models include: (1) the inlet velocity and (2) outlet pressure for each coolant channel, and (3) the axial pressure difference between narrow and wide adjacent channels. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Triggering Collapse of the Presolar Dense Cloud Core and Injecting Short-lived Radioisotopes with a Shock Wave. V. Nonisothermal Collapse Regime

Recent meteoritical analyses support an initial abundance of the short-lived radioisotope (SLRI) {sup 60}Fe that may be high enough to require nucleosynthesis in a core-collapse supernova, followed by rapid incorporation into primitive meteoritical components, rather than a scenario where such isotopes were inherited from a well-mixed region of a giant molecular cloud polluted by a variety of supernovae remnants and massive star winds. This paper continues to explore the former scenario, by calculating three-dimensional, adaptive mesh refinement, hydrodynamical code (FLASH 2.5) models of the self-gravitational, dynamical collapse of a molecular cloud core that has been struck by a thin shock front with a speed of 40 km s{sup −1}, leading to the injection of shock front matter into the collapsing cloud through the formation of Rayleigh–Taylor fingers at the shock–cloud intersection. These models extend the previous work into the nonisothermal collapse regime using a polytropic approximation to represent compressional heating in the optically thick protostar. The models show that the injection efficiencies of shock front materials are enhanced compared to previous models, which were not carried into the nonisothermal regime, and so did not reach such high densities. The new models, combined with the recent estimates of initial {sup 60}Fe abundances, imply that the supernova triggering and injection scenario remains a plausible explanation for the origin of the SLRIs involved in the formation of our solar system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Code modernization strategies for short-range non-bonded molecular dynamics simulations

Modern HPC systems are increasingly relying on greater core counts and wider vector registers. Thus, applications need to be adapted to fully utilize these hardware capabilities. One class of applications that can benefit from this increase in parallelism are molecular dynamics simulations. In this paper, we describe our efforts at modernizing the ESPResSo++ simulation package for molecular dynamics by restructuring its particle data layout for efficient memory accesses and applying vectorization techniques to benefit the calculation of short-range non-bonded forces, which results in an overall three times speedup and serves as a baseline for further optimizations. We also implement fine-grained parallelism for multi-core CPUs through HPX, a C++ runtime system which uses lightweight threads and an asynchronous many-task approach to maximize concurrency. Our goal is to evaluate the performance of an HPX-based approach compared to the bulk-synchronous MPI-based implementation. This requires the introduction of an additional layer to the domain decomposition scheme that defines the task granularity. On spatially inhomogeneous systems, which impose a corresponding load-imbalance in traditional MPI-based approaches, we demonstrate that by choosing an optimal task size, the efficient work-stealing mechanisms of HPX can overcome the overhead of communication resulting in an overall 1.4 times speedup compared to the baseline MPI version.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2018

Performing a Large-Scale Modal Test on the B2 Stand Crane at NASA's Stennis Space Center

A modal test of NASA's Space Launch System (SLS) Core Stage is scheduled to occur at the Stennis Space Center B2 test stand. A derrick crane with a 150-ft long boom, located at the top of the stand, will be used to suspend the Core Stage in order to achieve defined boundary conditions. During this suspended modal test, it is expected that dynamic coupling will occur between the crane and the Core Stage. Therefore, a separate modal test was performed on the B2 crane itself, in order to evaluate the varying dynamic characteristics and correlate math models of the crane. Performing a modal test on such a massive structure was challenging and required creative test setup and procedures, including implementing both AC and DC accelerometers, and performing both classical hammer and operational modal analysis. This paper describes the logistics required to perform this large-scale test, as well as details of the test setup, the modal test methods used, and an overview and application of the results.

Stasiunas, Eric C.↗
NASA NTRS2018

Performing a Large-Scale Modal Test on the B2 Stand Crane at NASA's Stennis Space Center

A modal test of NASA’s Space Launch System (SLS) Core Stage is scheduled to occur prior to propulsion system verification testing at the Stennis Space Center B2 test stand. A derrick crane with a 180-ft long boom, located at the top of the stand, will be used to suspend the Core Stage in order to achieve defined boundary conditions. During this suspended modal test, it is expected that dynamic coupling will occur between the crane and the Core Stage. Therefore, a separate modal test was performed on the B2 crane itself, in order to evaluate the varying dynamic characteristics and correlate math models of the crane. Performing a modal test on such a massive structure was challenging and required creative test setup and procedures, including implementing both AC and DC accelerometers, and performing both classical hammer and operational modal analysis. This paper describes the logistics required to perform this large-scale test, as well as details of the test setup, the modal test methods used, and an overview of the results.

Stasiunas, Eric C.↗
DOE OSTI2015

BROWN DWARF BINARIES FROM DISINTEGRATING TRIPLE SYSTEMS

Binaries in which both components are brown dwarfs (BDs) are being discovered at an increasing rate, and their properties may hold clues to their origin. We have carried out 200,000 N-body simulations of three identical stellar embryos with masses drawn from a Chabrier IMF and embedded in a molecular core. The bodies are initially non-hierarchical and undergo chaotic motions within the cloud core, while accreting using Bondi–Hoyle accretion. The coupling of dynamics and accretion often leads to one or two dominant bodies controlling the center of the cloud core, while banishing the other(s) to the lower-density outskirts, leading to stunted growth. Eventually each system transforms either to a bound hierarchical configuration or breaks apart into separate single and binary components. The orbital motion is followed for 100 Myr. In order to illustrate 200,000 end-states of such dynamical evolution with accretion, we introduce the “triple diagnostic diagram,” which plots two dimensionless numbers against each other, representing the binary mass ratio and the mass ratio of the third body to the total system mass. Numerous freefloating BD binaries are formed in these simulations, and statistical properties are derived. The separation distribution function is in good correspondence with observations, showing a steep rise at close separations, peaking around 13 AU and declining more gently, reaching zero at separations greater than 200 AU. Unresolved BD triple systems may appear as wider BD binaries. Mass ratios are strongly peaked toward unity, as observed, but this is partially due to the initial assumptions. Eccentricities gradually increase toward higher values, due to the lack of viscous interactions in the simulations, which would both shrink the orbits and decrease their eccentricities. Most newborn triple systems are unstable and while there are 9209 ejected BD binaries at 1 Myr, corresponding to about 4% of the 200,000 simulations, this number has grown to 15,894 at 100 Myr (∼8%). The total binary fraction among freefloating BDs is 0.43, higher than indicated by current observations, which, however, are still incomplete. Also, the gradual breakup of higher-order multiples leads to many more singles, thus lowering the binary fraction. The main threat to newly born triple systems is internal instabilities, not external perturbations. At 1 Myr there are 1325 BD binaries still bound to a star, corresponding to 0.66% of the simulations, but only 253 (0.13%) are stable on timescales >100 Myr. These simulations indicate that dynamical interactions in newborn triple systems of stellar embryos embedded in and accreting from a cloud core naturally form a population of freefloating BD binaries, and this mechanism may constitute a significant pathway for the formation of BD binaries.

79 ASTRONOMY AND ASTROPHYSICS↗