Dynamical behavior of ejecta from the moon. Part I - Initial conditions
Initial orbital elements of moon-earth trajectory for lunar dust or chunks ejected from surface
SEARCH · Search NASA
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.
Initial orbital elements of moon-earth trajectory for lunar dust or chunks ejected from surface
The effect of the initial condition at the jet exit on the downstream evolution, particularly within the potential core length, were numerically investigated as well as with available experimental data. In order to select the most dependable computational model for the present numerical experiment, a comparative study has been performed with different turbulence models at k-epsilon level, and it was found that the k-epsilon-gammma model yields superior prediction accuracy over other conventional models. The calculated results show that the potential core length and the spreading rate the initial mixing layer are dependent on the initial length scale as well as the turbulent kinetic energy at the jet exit. Such effect of the initial length scale increases with higher initial turbulence level. An empirical parameter has been devised to collapse the calculated data of the potential core length and the spreading rate with various initial conditions onto a single curve.
The initial conditions in heavy-ion collisions are calculated in many different frameworks. The importance of nucleon position fluctuations within the nucleus and subnucleon structure has been established when modeling initial conditions for input to hydrodynamic calculations. However, there remain outstanding puzzles regarding these initial conditions, including the measurement of the near equivalence of the elliptical v 2 and triangular v 3 flow coefficients in ultracentral 0–1% Pb + Pb collisions at the CERN Large Hadron Collider. Recently a calculation termed magma incorporating gluonic hot spots via two-point correlators in the color glass condensate framework, and no nucleons, provided a simultaneous match to these flow coefficients measured by the ATLAS experiment, including in ultracentral 0–1% collisions. Our calculations reveal that the magma initial conditions do not describe the experimental data when run through full hydrodynamic sonic simulations or when the hot spots from one nucleus resolve hot spots from the other nucleus, as predicted in the color glass condensate framework. Finally, we also explore alternative initial condition calculations and discuss their implications.
A theory of initial conditions is necessary for a complete explanation of the presently observed large scale structural features of the universe, and a quantum theory of cosmology is probably needed for its formulation. The kinematics of quantum cosmology are reviewed, and some candidates for a law of initial conditions are discussed. The proposal that the quantum state of a closed universe is the natural analog of the ground state for closed cosmologies and is specified by a Euclidean sum over histories is sketched. When implemented in simple models, this proposal is consistent with the most important large-scale observations.
An experimental study was conducted concerning the influence of small changes in initial conditions on the near- and far-field evolution of the three-dimensional structure of a plan mixing layer. A two-stream mixing layer with a velocity ratio of 0.6 was generated with the initial boundary layers on the splitter plate laminar and was nominally two-dimensional. The initial conditions were changed slightly by interchanging the high- and low-speed sides of the wind tunnel, while maintaining the same velocities, and hence velocity ratio. This resulted in small changes in the initial boundary layer properties, and the perturbations present in the boundary layers were interchanged between the high- and low-speed sides for the two cases. The results indicate that, even with this relatively minor change in initial conditions, the near-field regions of the two cases differ significantly. The peak Reynolds stress levels in the near-field differ by up to 100 percent, and this is attributed to a difference in the location of the initial spanwise vortex roll-up. In addition, the positions and shapes of the individual streamwise vortical structures differ for the two cases, although the overall structures differ for the two cases, although the overall qualitative description of these structures is comparable. The subsequent reorganization and decay of the streamwise vortical structures is very similar for the two cases. As a result, in the far field, both mixing layers achieve similar structure, yielding comparable growth rates, Reynolds stress, distribution, and spectral content.
The previous study of Young and Villere concerning growth of planetary scale waves forced by wave-wave interactions of amplifying intermediate scale baroclinic eddies is extended to investigate effects of different eddy initial conditions. A global, spectral, primitive equation model is used for the calculations. For every set of eddy initial conditions considered, growth rates of planetary modes are considerably greater than growth rates computed from linear instability theory for a fixed zonally independent basic state. However, values of growth rates ranged over a factor of 3 depending on the particular set of eddy initial conditions used. Nonlinear forcing of planetary modes via wave-wave coupling becomes more important than baroclinic growth on the basic state at small values of the intermediate-scale modal amplitudes. The relative importance of direct transfer of kinetic energy from intermediate scales of motion to a planetary mode, compared to baroclinic conversion of available potential energy to kinetic energy within that planetary mode, depends on the individual case. In all cases, however, the transfer of either kinetic or available potential energy to the planetary modes was accomplished principally by wave-wave transfer from intermediate scale eddies, rather than from the zonally averaged state. The zonal wavenumber 2 planetary mode was prominent in all solutions, even in those for which eddy initial conditions were such that a different planetary mode was selectively forced at the start. General characteristics of the structural evolution of the planetary wave components of total heat and momentum flux, and modal structures themselves, were relatively insensitive to variations in eddy initial conditions, even though quantitative details varied from case to case.
Computational models of plasma technologies often solve for the system operating conditions by time-stepping an initial value problem to a quasi-steady solution. However, the strongly nonlinear and multi-timescale nature of plasma dynamics often necessitate millions, or even hundreds of millions, of steps to reach convergence, reducing the effectiveness of these simulations for computer-aided engineering. We consider acceleration of kinetic plasma simulations via data-driven machine-learning-generated initial conditions, which initialize the simulations close to their final quasi-steady-state, thereby reducing the number of steps to reach convergence. Three machine-learning models are developed to predict the density and ion kinetic profiles of capacitively coupled plasma discharges relevant to the microelectronics industry. The models are trained on kinetic simulations over a range of device operating frequencies and pressures. Best performance was observed when simulations were initialized with ion kinetic profiles generated by a convolutional neural network, reducing the mean number of steps to reach convergence by 17.1× when compared to initialization with a zero-dimensional global model. We also outline a workflow for continuous data-driven model improvement and simulation speedup, with the aim of generating sufficient data for full device digital twins.
Numerical experimens are presented which were intended to verify the independence of bar formation in rotating galaxy models from special initial conditions. Three-dimensional n-body programs were run starting from aspherical initial conditions characterized by particle rotational and peculiar velocities and offsets from sphericity for a set of 100,000 particles. All the initial configurations examined are observed to collapse to a thin form, overshoot, collapse again and again overshoot, with diameters and central densities similar to those observed under spheric initial conditions. Every aspheric collapse is also found to result in bar formation within two rotation periods, which is even sooner than in initially spherical conditions. In the intermediate stages of collapse, transient sheet patterns are found to be destroyed at asymmetries of 2-3%, while transient rings survived asymmetries of 6-7% but not 15-18%. It is noted that the results may be applicable to galaxy formation in the early universe, and possible mechanisms for the arresting of protogalactic collapse are indicated.
The collapse and fragmentation of rotating protostellar clouds is explored, starting from nonuniform density and nonuniform rotation initial conditions. Whether binary fragmentation occurs during the first dynamic collapse phase depends strongly on the initial density profile. Exponential clouds are only somewhat more resistant to fragmentation than uniform-density clouds, but power-law clouds do not undergo fragmentation for likely values of a relevant parameter. Because binary fragments start from profiles intermediate between uniform density and exponential clouds, minimum protostellar mass for population I stars should be increased to approximately 0.02 solar mass. The axisymmetric Terey et al. (1984) model should be stable with respect to nonaxisymmetric perturbations. Considering the observed binary frequency, collapse from power-law initial conditions appears to be less common than collapse from more uniform initial conditions.
A pair of microphones (separated axially by 5.08 cm and laterally by 1.3 cm) are placed on either side of the jet centerline to investigate coherent pressure fluctuations in an axisymmetric jet at Strouhal numbers less than unity. Auto-spectra, transfer-function, and coherence measurements are made for a tripped and untripped boundary layer initial condition. It was found that coherent acoustic pressure waves originating in the upstream plenum chamber propagate a greater distance downstream for the tripped initial condition than for the untripped initial condition. In addition, for the untripped initial condition the development of the coherent hydrodynamic pressure waves shifts downstream.
A pair of microphones (separated axially by 5.08 cm and laterally by 1.3 cm) are placed on either side of the jet centerline to investigate coherent pressure fluctuations in an axisymmetric jet at Strouhal numbers less than unity. Auto-spectra, transfer-function, and coherence measurements are made for a tripped and untripped boundary layer initial condition. It was found that coherent acoustic pressure waves originating in the upstream plenum chamber propagate a greater distance downstream for the tripped initial condition than for the untripped initial condition. In addition, for the untripped initial condition the developmet of the coherent hydrodynamic pressure waves shifts downstream.
The effects of relatively small changes in the initial conditions of the development of the 3D structure of a plane mixing layer originating from laminar boundary layers were experimentally studied. It is found that, while the exact shapes and positions of the streamwise vortex structures are not the same for the two initial conditions, their overall distribution, reorganization, and decay are very similar. The results likely that while some of the specific details of the streamwise vortex structure may be facility-dependent, a relatively strong structure which produces significant three-dimensionality should form in all mixing layers originating from laminar boundary layers. After some initial readjustments, the structure will appear in the form of counterrotating pairs of streamwise vortices which in the mean, grow with the mixing layer and decay in strength. The results also serve to warn that small changes in initial conditions may significantly affect the Reynolds stress distributions in the near field.
High-energy nuclear collisions encompass three key stages: the structure of the colliding nuclei, informed by low-energy nuclear physics, the initial condition , leading to the formation of quark–gluon plasma (QGP), and the hydrodynamic expansion and hadronization of the QGP, leading to final-state hadron distributions that are observed experimentally. Recent advances in both experimental and theoretical methods have ushered in a precision era of heavy-ion collisions, enabling an increasingly accurate understanding of these stages. However, most approaches involve simultaneously determining both QGP properties and initial conditions from a single collision system, creating complexity due to the coupled contributions of these stages to the final-state observables. To avoid this, we propose leveraging established knowledge of low-energy nuclear structures and hydrodynamic observables to independently constrain the QGP’s initial condition. By conducting comparative studies of collisions involving isobar-like nuclei—species with similar mass numbers but different ground-state geometries—we can disentangle the initial condition’s impacts from the QGP properties. This approach not only refines our understanding of the initial stages of the collisions but also turns high-energy nuclear experiments into a precision tool for imaging nuclear structures, offering insights that complement traditional low-energy approaches. Opportunities for carrying out such comparative experiments at the Large Hadron Collider and other facilities could significantly advance both high-energy and low-energy nuclear physics. Additionally, this approach has implications for the future electron-ion collider. While the possibilities are extensive, we focus on selected proposals that could benefit both the high-energy and low-energy nuclear physics communities. Originally prepared as input for the long-range plan of U.S. nuclear physics, this white paper reflects the status as of September 2022, with a brief update on developments since then.
We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.
The effects of rectangular initial data on the evolution of solitons governed by the Korteweg-de Vries equation is studied. Both isolated and separated disturbances are considered, providing some general insight into how the nature of the initial condition influences the appearance of solitons in the asymptotic state. The analytic approach is based on the inverse scattering transform which relates the initial condition to Schroedinger's equation. The results are used to model the initial shallow-water disturbances, and the results can be summarized by stating that the number of evolved solitons depends on the strength of each rectangular disturbance, the relative amplitudes of the rectangular disturbances, and the relative proximity of the disturbances.
Nonlinear response of a flexible curved panel exhibiting bifurcation to fully developed chaos is demonstrated along with the sensitivity to small perturbation from the initial conditions. The response is determined from the measured time series at two fixed points. The panel is forced by an external nonharmonic multifrequency and monofrequency sound field. Using a low power time-continuous feedback control, carefully tuned at each initial condition, produces large long-term effects on the dynamics toward taming chaos. Without the knowledge of the initial conditions, control may be achieved by destructive interference. In this case, the control power is proportional to the loading power. Calculation of the correlation dimension and the estimation of positive Lyapunov exponents, in practice, are the proof of chaotic response.
Relationships between diverse thermodynamic quantities appropriate to weakly compressible turbulence are derived. It is shown that for turbulence of a finite turbulent Mach number there is a finite element of compressibility. A methodology for generating initial conditions for the fluctuating pressure, density and dilatational velocity is given which is consistent with finite Mach number effects. Use of these initial conditions gives rise to a smooth development of the flow, in contrast to cases in which these fields are specified arbitrarily or set to zero. Comparisons of the effect of different types of initial conditions are made using direct numerical simulation of decaying isotropic turbulence.
Heavy quarks produced in relativistic heavy-ion collisions are known to be sensitive probes of the hot and dense QCD matter they traverse. In this paper we study how their dynamics is affected by the nature of the bulk evolution of the QCD matter, the initial condition of the system, and the treatment of elementary interactions between heavy quarks and the surrounding medium. For the same initial condition and the same quark-gluon plasma (QGP) expansion scenario we discuss the consequences of the assumption of a local equilibrium by comparing the consequences for the nuclear modification factor $R_{AA}$ and the elliptic flows of charm quarks, scrutinizing the different components of the final distribution of charm quarks. For this purpose we employ the parton-hadron-string dynamics (PHSD) model, which is an off-shell microscopic transport approach, as well as the linearized-Boltzmann (LB) scheme obtained by coarse graining the PHSD bulk and assuming local equilibrium for the interactions of the charm quarks with the bulk. The $R_{AA}$ of charm quarks stemming from the later LB approach is also compared to a genuine fluid dynamics evolution initiated by the coarse grained PHSD, which allows us to further assess the consequences of reducing the full $\textit{n}$-body dynamics. We then proceed to a systematic comparison of PHSD (in its LB approximation) with MC@HQ, another transport model for heavy flavors which also relies on the LB approach. In particular, we investigate the consequences for the nuclear modification factor of charm quarks if we vary separately the initial heavy quark distribution function in matter, the expansion dynamics of the QGP, and the elementary interactions of heavy quarks of these models. We find that the results for both models vary significantly depending on the details of the calculation. However, both models achieve very similar predictions for key heavy quark observables for certain combinations of initial condition, bulk evolution, and interactions. In this work, we conclude that this ambiguity limits our ability to determine the different properties of the system based on the current set of observables.