Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Practical Effects of Integrating Temperature with Strang Split Reactions

For astrophysical reacting flows, operator splitting is commonly used to couple hydrodynamics and reactions. Each process operates independent of one another, but by staggering the updates in a symmetric fashion (via Strang splitting) second order accuracy in time can be achieved. However, approximations are often made to the reacting system, including the choice of whether or not to integrate temperature with the species. Here we demonstrate through a simple convergence test that integrating an energy equation together with reactions achieves the best convergence when modeling reactive flows with Strang splitting. Additionally, second order convergence cannot be achieved without integrating an energy or temperature equation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Accurate hydrodynamic force and torque correlations for prolate spheroids from Stokes regime to high Reynolds numbers

Detailed simulations of flow around various prolate spheroids are performed using multi-relaxation time lattice Boltzmann method (LBM). The simulations are performed in the Reynolds numbers range 0.1 ≤ Re ≤ 2000 at different incident angles 0°≤ φ ≤ 90° for various prolate spheroids of aspect ratios 1 ≤ λ ≤ 8. The Re is based on the volume equivalent sphere diameter. From the simulations, accurate correlations for average drag, lift and torque coefficients (C D , C L , and C T respectively) are proposed. The mean deviations between the correlations and the simulation data for C D , C L , and C T are 2%, 6.5%, and 4.7% respectively. Furthermore, the correlations are fitted with the physics in mind such that they can be extrapolated beyond the regimes simulated. The fits are designed such that correlations reduce to analytical (C D and C L ) solutions available for prolate spheroids in the limit Re ≈ 0. At high Re, the correlations mimic the asymptotic flattening of C D as for spheres and therefore, we expect the correlations to be valid until the critical Re (≈ 10 5 ). Furthermore, extensive comparison of our correlations with other high-fidelity simulations from literature demonstrate the accuracy of our work compared to other correlations from literature. This work enables researchers to perform accurate unresolved Euler-Lagrangian simulations for a wide range of elongated particles.

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

Bayesian analysis of (3 +1)⁢D relativistic nuclear dynamics with the RHIC beam energy scan data

This work presents a Bayesian inference study for relativistic heavy-ion collisions in the beam energy scan program at the BNL Relativistic Heavy-Ion Collider. The theoretical model simulates event-by-event (3+1)-dimensional [(3+1)⁢D] collision dynamics using hydrodynamics and hadronic transport theory. We analyze the model's 20-dimensional posterior distributions obtained using three model emulators with different accuracy and demonstrate the essential role of training an accurate model emulator in the Bayesian analysis. Our analysis provides robust constraints on the quark-gluon plasma's transport properties and various aspects of (3+1)⁢D relativistic nuclear dynamics. By running full model simulations with 100 parameter sets sampled from the posterior distribution, we make predictions for p T -differential observables and estimate their systematic theory uncertainty. Here, a sensitivity analysis is performed to elucidate how individual experimental observables respond to different model parameters, providing useful physics insights into the phenomenological model for heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes

The Kelvin–Helmholtz (KH) instability of a shear layer with an initially uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with hydrodynamic KH, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited, thus transporting less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear because of the stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

RAFT: Reconfigurable Array of High-Efficiency Ducted Turbines for Hydrokinetic Energy Harvesting

Diversifying the energy harvesting portfolio is crucial to achieving the ambitious goal of transitioning to clean energy by 2030. Marine hydrokinetic energy has garnered renewed interest due to its high harvesting potential in the U.S., and the resource's reliability and predictability—remaining relatively constant on a daily basis and available 24/7. However, there are currently few commercial devices capable of harnessing the energy from flowing water. This project aims to bridge that gap by designing and evaluating a novel hydrokinetic turbine concept that can efficiently harvest energy from both rivers and tidal streams. The RAFT (Reconfigurable Array of High-Efficiency Ducted Turbines) concept introduces a duct surrounding the turbine rotor and creates an array of small 5-kW units. The duct serves two primary purposes: (1) it enhances hydrodynamic efficiency by accelerating flow to the rotor, and (2) it functions as a structural component, facilitating the formation of modular arrays that lower costs. This project focuses on demonstrating this concept and validating these benefits through simulations and scaled prototype testing. The project team includes 8 faculty members and over 20 students from 3 universities, organized into three core areas: hydrodynamics, electrical systems, and structural analysis, with additional teams dedicated to system integration, environmental assessment and risk management, and tech-to-market strategy. The team successfully demonstrated the increased hydrodynamic efficiency of a ducted turbine compared to an unducted version using high-fidelity simulations and prototype tests. Moreover, design optimization efforts led to surpassing the SHARKS program's goal of 60% reduction of the levelized cost of energy with a significant margin.

13 HYDRO ENERGY↗
DOE OSTI2023

Development of an accelerator-based neutron source to prototype Mo-99 production, part I: A liquid LBE windowless target

In this article, Molybdenum-99 (Mo-99)’s decay product, technetium-99 (Tc-99 m), is one of the most critical isotopes for medical diagnostics. To provide U.S. domestic supply of Mo-99 without using high-enriched uranium (HEU), a subcritical uranium target assembly (UTA) is irradiated by an accelerator-based neutron source to create Mo-99 through fission. This study discusses the development of the accelerator-based neutron source. The high-energy electrons from the accelerator irradiate a liquid lead-bismuth eutectic (LBE) target to produce neutrons. Part I of this work focuses on numerical and experimental analysis towards the development of a liquid LBE windowless target. Unlike the existing windowless targets in literature, the current design creates a vertical free surface for a beam to irradiate. First, a hydrodynamic analysis of the LBE windowless target is performed. Simplified analytical calculations are assisted by 2D computational fluid dynamics (CFD) simulations to design the target, with the focus on eliminating recirculation zones and avoiding cavitation. With the optimized geometry, the experimental study is performed to investigate the flow hydrodynamics using liquid LBE. The experiments (1) compare pressure drop in the system to correlation predictions; (2) visualize the free surface liquid LBE flow from the beam view; (3) validate the LBE flow profile using temperature sensitive paint from the side view; and (4) validate the liquid LBE film thickness using gamma densitometer measurements. Second, the power handling capability of the designed windowless target is investigated. The divider plate in the current design is susceptible to overheating due to the thin LBE film in front. As LBE erosion and corrosion is likely to occur at an LBE velocity of 2.0 m/s and temperature above 500 °C, a power limit of 10 kW of beam power was established to prevent this corrosion from occurring, which is calculated by a Nusselt number correlation. The divider plate surface temperature at 10 kW agrees well with the 3D CFD simulation results. Part I demonstrates the fundamental physics in liquid LBE windowless target design and associated testing. A companion paper, Part II will demonstrate how to couple this windowless target into the Mo-99 production system, including an accelerator system operating under an ultra-high vacuum and the UTA cooled by water at room temperature.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Examining the effects of soil entrainment during nuclear cloud rise on fallout predictions using a multiscale atmospheric modeling framework

Current operational models for nuclear cloud rise over land were developed and validated using observations from shallow-buried or surface detonations, where lofted soil quickly mixed with fission products from the detonation. These models poorly predict fallout from elevated detonations near the fallout-free height of burst (FFHOB), where interactions with the ground are limited and the mixing of fission products and lofted soil is incomplete. Fallout-free is a misnomer at this HOB, as fallout was observed in these cases, but was below the levels of concern, especially off-grounds of the nuclear test site. To correctly characterize and model fallout from detonations near the FFHOB, models must be developed which can capture the stratified nature of the particle and activity-size distributions within the cloud. Previously, it was shown that the Weather Research and Forecasting (WRF) model can accurately simulate nuclear cloud rise for airbursts with little to no ground interactions (Arthur et al., 2021). That work is expanded here by (1) using a radiation-hydrodynamics code to improve the fireball initialization in WRF, (2) further developing an aerosol package from WRF-Chem to simulate lofted soil, and (3) combining the WRF cloud rise simulations with the operational models used at the National Atmospheric Release Advisory Center (NARAC) for fallout modeling. Using this combination of codes, the Upshot-Knothole Grable detonation, which was just below the FFHOB, is simulated from seconds after detonation through cloud rise and fallout, and results are compared to historical test data. Here, the results show improved prediction of dose rate and highlight the need to correctly characterize the entrainment of material into the cloud and the subsequent mixing of fission products with entrained material.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2026

Cholla Galactic OutfLow Simulations (CGOLS)

These datasets contain full hydro-field snapshots from the galactic outflow simulations in the CGOLS suite, models I-V. The datasets were generated using the Cholla hydrodynamics code (https://github.com/cholla-hydro/cholla); descriptions of the models are in the associated publications (Schneider & Robertson 2018, ApJ; Schneider et al. 2018, ApJ; Schneider et al. 2020, ApJ; and Schneider & Mao, 2024, ApJ). Each hdf5 dataset is numbered according to the simulation time of the snapshot, in Myr. Fields include density, x momentum, y momentum, z momentum, total energy, and thermal energy (for models I - III), as well as a passive scalar field (models IV and V). 2 dimensional density and temperature projections, as well as slices along each midplane are also included if they exist.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Direct numerical simulation of open-channel flow over a heterogeneous particle bed at low relative submergence

In this article, this study investigates turbulent open-channel flows over beds of irregularly arranged particles, using direct numerical simulations at a friction Reynolds number of Re T = 300⁠. Two distinct cases are examined: a polydisperse bed (P800) composed of multiple layers of randomly distributed spheres of varying sizes, and a monodisperse bed (M1015) formed by a random distribution of uniform sized spheres, with a bottommost single layer of varied-sized particles to introduce realistic randomness. Our investigation unveils a rich network of low- and high-speed streaks within the flow field, exhibiting distinctive behaviors in different bed configurations. The P800 case presents a poorly organized flow pattern induced by the varied particle sizes and arrangements, while the M1015 case shows a more regular flow pattern, marked by larger streaks. We also observe that total wall shear stress is substantially influenced by surface roughness-induced drag, extending beyond the effects documented in existing studies of open-channel flows. The present study reveals intricate secondary flow patterns over irregular particle beds. Large-scale circulations are discerned around particle crests in the P800 case and localized circulations with increased turbulence in the M1015 case. Furthermore, analysis of Reynolds stress tensor components indicates that roughness disrupts coherent turbulent eddies, consequently mitigating peak stress. We quantify correlations between drag force and local fluid velocity fluctuations. Notably, a larger deviation in drag is observed in the P800 case compared to M1015, accentuating the influence of particle size and distribution on fluid–particle interactions.

42 ENGINEERING↗
DOE OSTI2021

The collapse and three-dimensional explosion of three-dimensional massive-star supernova progenitor models

ABSTRACT The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M⊙ model, evolved in three dimensions from collapse to bounce through explosion with the radiation-hydrodynamic code fornax. We compare the results using those starting from three-dimensional progenitors to three-dimensional simulations of spherically symmetric, one-dimensional progenitors of the same mass. We find that the models evolved in three dimensions during the final stages of massive star evolution are more prone to explosion. The turbulence arising in these multidimensional initial models serves as seed turbulence that promotes shock revival. Detection of gravitational waves and neutrinos signals could reveal signatures of pre-bounce turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Realistic Full-Scale 3D Modeling of Turning Using Coupled Smoothed Particle Hydrodynamics and Finite Element Method for Predicting Cutting Forces

Computational modelling is an effective technique for understanding the complex physics of machining. Large deformations, material separation, and high computational requirements are the key challenges faced while simulating machining. This work introduces a full-scale three-dimensional model of turning operations using a combined approach based on the Smoothed Particle Hydrodynamics (SPH) and Finite Element (FE) methods. By exploiting the advantages of each method, this approach leads to high-fidelity coupled SPH-FE machining models. Cutting forces and chip morphology are the primary results of interest. The machining models are validated with the results of turning experiments. Two-dimensional machining model underpredicts the cutting force and feed force by approximately 49% and 70%, respectively. Moreover, passive force cannot be predicted using the two-dimensional model. On the other hand, with the three-dimensional models developed in this manuscript, the difference between the total simulated force and experimentally measured force is ∼17%. The chip morphologies correlate with experiments in terms of the direction of the chip movement and the “long” continuous chips observed while turning Al 6061. This work expands the realm of machining simulations from two-dimensional orthogonal machining or sectional three-dimensional model to a full-scale realistic simulation. The encouraging simulation results show the potential to study more complex phenomena, such as machining stability and tool path modulation.

42 ENGINEERING↗
DOE OSTI2024

Project DarkStar: Vision for LLNL in 2030

DarkStar was a Strategic Initiative (FY2021-FY2024) to investigate applications of Artificial Intelligence (AI) and Machine Learning (ML) to scientific problems of complex hydrodynamics, shockwave physics and energetic materials. The research focused on physics and engineering design as a process that can be tremendously accelerated through merging AI with advanced physics simulation on exascale-class platforms, and to experimentally validate this revolutionary new approach through dynamic materials campaigns. A central thread of scientific inquiry was in the application of AI to enable human understanding of how to control hydrodynamic instability (which has impacts to areas such as inertial confinement fusion) via engineering features and time-dependent sources. Motivated by an unfinished line of research started by Dr. Johnny von Neumann, AI-enabled simulation approaches were developed that allowed DarkStar researchers to uncover several ground-breaking discoveries regarding hydrodynamic instability, including how to completely suppress Richtmyer-Meshkov instability (RMI). These S&T discoveries, along with other advances, have shown the way for an entirely new approach to time-dependent problems known as inverse design – the idea that complex systems can be developed directly from a final state that is to be achieved and resolve the initial design via satisfying several constraints simultaneously via AI/ML. Through experimental campaigns conducted across a wide range of facilities in the NNSA complex (the High Explosive Application Facility at LLNL, the Dynamic Compression Sector/Advanced Photon Source at Argonne National Lab, and Special Technologies Laboratory at MSTS) the radical new AI/ML approach to engineering complex material dynamics was verified, establishing a new field of study within the realm of shock physics. As advanced manufacturing capabilities continue to develop, the great importance of inverse design as a means to apply that technology effectively for NNSA missions will feature prominently over this decade. DarkStar has positioned NNSA as a world-leader in this newly emerging cross-disciplinary area of AI methods for advanced physics simulation and pioneered multiple novel approaches that have enabled the broader scientific community. By allowing us to see past the horizon, to 2030 and beyond, DarkStar has illuminated the vast potential of AI/ML to impact a wide range of new national security missions and, consequently, multiple areas of further research have already emerged across the NNSA and DOD complex.

42 ENGINEERING↗
DOE OSTI2021

Tidal disruption discs formed and fed by stream–stream and stream–disc interactions in global GRHD simulations

When a star passes close to a supermassive black hole (BH), the BH’s tidal forces rip it apart into a thin stream, leading to a tidal disruption event (TDE). In this work, we study the post-disruption phase of TDEs in general relativistic hydrodynamics (GRHD) using our GPU-accelerated code h-amr. We carry out the first grid-based simulation of a deep-penetration TDE (β = 7) with realistic system parameters: a black hole-to-star mass ratio of 10 6 , a parabolic stellar trajectory, and a non-zero BH spin. We also carry out a simulation of a tilted TDE whose stellar orbit is inclined relative to the BH midplane. We show that for our aligned TDE, an accretion disc forms due to the dissipation of orbital energy with ~20 percent of the infalling material reaching the BH. The dissipation is initially dominated by violent self-intersections and later by stream–disc interactions near the pericentre. The self-intersections completely disrupt the incoming stream, resulting in five distinct self-intersection events separated by approximately 12 h and a flaring in the accretion rate. We also find that the disc is eccentric with mean eccentricity e ≈ 0.88. For our tilted TDE, we find only partial self-intersections due to nodal precession near pericentre. Although these partial intersections eject gas out of the orbital plane, an accretion disc still forms with a similar accreted fraction of the material to the aligned case. These results have important implications for disc formation in realistic tidal disruptions. For instance, the periodicity in accretion rate induced by the complete stream disruption may explain the flaring events from Swift J1644+57.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Spatially and temporally resolved plasma parameter estimations of laser heated MagLIF relevant gas pipes at NIF

The ability to control laser pre-heat is an integral part of the inertial confinement fusion concept known as Magnetized Liner Inertial Fusion. This process is studied at the National Ignition Facility (NIF) where 4 of the 192 laser beams are propagated through a 1-cm long gas cell where they deposit >20 kJ of energy into the gaseous fuel via inverse bremsstrahlung absorption. This process ionizes the gas, producing a plasma that follows behind the laser front and expands over the radius of the cell. Emission from this plasma, as viewed by a gated x-ray detector, can be used to build spatially and temporally resolved estimations of the pre-heat plasma's density and temperature profiles. This can then be used to estimate the plasma pressure, internal energy, and radiation losses. Estimations show the evolution of the plasma in magnetized and unmagnetized gas cells filled with ambient temperature neopentane (C5H12) +1% Ar, as well as unmagnetized cryogenically cooled (32 K) deuterium +1% Ne filled targets. This analysis shows the effects of initial gas-fill density, composition, and axial magnetization on the time-dependent plasma parameters. Previously, these parameters at the NIF had not been experimentally characterized, and these estimations provided a potential new means of testing radiation magneto-hydrodynamic predictive capability models. Results in unmagnetized targets have strong agreement with simulations. However, in targets with a 19 T applied axial magnetic field, this method yields electron temperatures up to 100% hotter than those predicted by HYDRA codes.

Bremsstrahlung↗
DOE OSTI2021

Light Regulation of Phytoplankton Growth in San Francisco Bay Studied Using a 3D Sediment Transport Model

In San Francisco Bay (SFB), light availability is largely determined by the concentration of suspended particulate matter (SPM) in the water column. SPM exhibits substantial variation with time, depth, and location. To study how SPM influences light and phytoplankton growth, we coupled a sediment transport model with a hydrodynamic model and a biogeochemical model. The coupled models were used to simulate conditions for the year of 2011 with a focus on northern SFB. For comparison, two simulations were conducted with ecosystem processes driven by SPM concentrations supplied by the sediment transport model and by applying a constant SPM concentration of 20 mg l –1 . The sediment transport model successfully reproduced the general pattern of SPM variation in northern SFB, which improved the chlorophyll-a simulation resulting from the biogeochemical model, with vertically integrated primary productivity varying greatly, from 40 g[C] m –2 year –1 over shoals to 160 g[C] m –2 year –1 in the deep channel. Primary productivity in northern SFB is influenced by euphotic zone depth ( Ze ). Our results show that Ze in shallow water regions (<2 m) is mainly determined by water depth, while Ze in deep water regions is controlled by SPM concentration. As a result, Ze has low (high) values in shallow (deep) water regions. Large (small) differences in primary productivity exist between the two simulations in deep (shallow) water regions. Furthermore, we defined a new parameter F light for “averaged light limitation” in the euphotic zone. The averaged chlorophyll-a concentration in the euphotic zone and F light share a similar distribution such that both have high (low) values in shallow (deep) water regions. Our study demonstrates that light is a critical factor in regulating the phytoplankton growth in northern SFB, and a sediment transport model improves simulation of light availability in the water column.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Geometry Optimization of Cable-Based Actuation for Small-Scale Model Testing of a Floating Marine Turbine: Preprint

This research aims to apply combined wave and tidal current loads to a small-scale floating marine current turbine in a wave tank, where an actuation system applies hydrodynamic and mooring forces on the hardware based on results from a simulation. We use a real-time hybrid test setup with physical wave forcing from the wave tank and simulated current and mooring forces implemented through a tensioned cable array. For this paper, we developed an optimization algorithm that adjusts the hardware geometry of the cable array to achieve more efficient tension allocation across the cables for the loads that need to be actuated. By adjusting the points where the cables are attached to the floating platform and the angles between the platform and the winches, an optimal cable geometry can be found to minimize tension variations in the lines, maintain the desired pretension, and prevent excessive tensions. We present the optimization problem formulation, the actuation system evaluation approach, and optimization results that show effective cable actuation setups that are being considered for implementation in the wave tank tests.

cable-based actuation↗
DOE OSTI2024

A Gaussian process based surrogate approach for the optimization of cylindrical targets

Simulating direct-drive inertial confinement experiments presents significant computational challenges, both due to the complexity of the codes required for such simulations and the substantial computational expense associated with target design studies. Machine learning models, and in particular, surrogate models, offer a solution by replacing simulation results with a simplified approximation. In this study, we apply surrogate modeling and optimization techniques that are well established in the existing literature to one-dimensional simulation data of a new cylindrical target design containing deuterium–tritium fuel. These models predict yields without the need for expensive simulations. We find that Bayesian optimization with Gaussian process surrogates enhances sampling efficiency in low-dimensional design spaces but becomes less efficient as dimensionality increases. Nonetheless, optimization routines within two-dimensional and five-dimensional design spaces can identify designs that maximize yield, while also aligning with established physical intuition. Optimization routines, which ignore constraints on hydrodynamic instability growth, are shown to lead to unstable designs in 2D, resulting in yield loss. However, routines that utilize 1D simulations and impose constraints on the in-flight aspect ratio converge on novel cylindrical target designs that are stable against hydrodynamic instability growth in 2D and achieve high yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Anisotropic fluid dynamical simulations of heavy-ion collisions

We present VAH, a (3+1)–dimensional simulation that evolves the far-from-equilibrium quark-gluon plasma produced in ultrarelativistic heavy-ion collisions with anisotropic fluid dynamics. We solve the hydrodynamic equations on an Eulerian grid using the Kurganov–Tadmor algorithm in combination with a new adaptive Runge–Kutta method. Our numerical scheme allows us to start the simulation soon after the nuclear collision, largely avoiding the need to integrate it with a separate pre-equilibrium dynamics module. We test the code's performance by simulating on the Eulerian grid conformal and non-conformal Bjorken flow as well as conformal Gubser flow, whose (0+1)–dimensional solutions are precisely known. Finally, we compare non-conformal anisotropic hydrodynamics to second-order viscous hydrodynamics in central Pb+Pb collisions and find that the former's longitudinal flow profile responds more consistently to the fluid's gradients along the spacetime rapidity direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗