Search NASA⌕ Search

SEARCH · Search NASA

Results for “Tidal Modeling”

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 163 records · Page 9

Tidal disruption of dissipative planetesimals

A self-consistent numerical model is developed for the tidal disruption of a solid planetesimal. The planetesimal is treated as a highly viscous, slightly compressible fluid whose disturbed parts are an inviscid, pressureless fluid undergoing distortion and disruption. The distortions were constrained to being symmetrical above and below the equatorial plane. The tidal potential is expanded in terms of Legendre polynomials, which eliminates the center of mass acceleration effects, permitting definition of equations of motion in a noninertial frame. Consideration is given to viscous dissipation and to characteristics of the solid-atmosphere boundary. The model is applied to sample cases in one, two and three dimensions.

Mizuno, H.↗

Orbital, Rotational and Climatic Interactions: Energy Dissipation and Angular Momentum Exchange in the Earth-Moon System

A numerical ocean tide model has been developed and tested using highly accurate TOPEX/Poseidon (T/P) tidal solutions. The hydrodynamic model is based on time stepping a finite difference approximation to the non-linear shallow water equations. Two novel features of our implementation are a rigorous treatment of self attraction and loading (SAL), and a physically based parameterization for internal tide (IT) radiation drag. The model was run for a range of grid resolutions, and with variations in model parameters and bathymetry. For a rational treatment of SAL and IT drag, the model run at high resolution (1/12 degree) fits the T/P solutions to within 5 cm RMS in the open ocean. Both the rigorous SAL treatment and the IT drag parameterization are required to obtain solutions of this quality. The sensitivity of the solution to perturbations in bathymetry suggest that the fit to T/P is probably now limited by errors in this critical input. Since the model is not constrained by any data, we can test the effect of dropping sea-level to match estimated bathymetry from the last glacial maximum (LGM). Our results suggest that the 100 m drop in sea-level in the LGM would have significantly increased tidal amplitudes in the North Atlantic, and increased overall tidal dissipation by about 40%. However, details in tidal solutions for the past 20 ka are sensitive to the assumed stratification. IT drag accounts for a significant fraction of dissipation, especially in the LGM when large areas of present day shallow sea were exposed, and this parameter is poorly constrained at present.

Egbert, Gary D.↗

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.↗

Shallow Lunar Seismic Activity and the Current Stress State of the Moon

A vast, global network of more than 3200 lobate thrust fault scarps has been revealed in high resolution Lunar Reconnaissance Orbiter Camera (LROC) images. The fault scarps very young, less than 50 Ma, based on their small scale and crisp appearance, crosscutting relations with small-diameter impact craters, and rates of infilling of associated small, shallow graben and may be actively forming today. The population of young thrust fault scarps provides a window into the recent stress state of the Moon and offers insight into the origin of global lunar stresses. The distribution of orientations of the fault scarps is non-random, inconsistent with isotropic stresses from late-stage global contraction as the sole source of stress Modeling shows that tidal stresses contribute significantly to the current stress state of the lunar crust. Tidal stresses (orbital recession and diurnal tides) superimposed on stresses from global contraction result in non-isotropic compressional stress and thrust faults consistent with lobate scarp orientations. Stresses due to orbital recession do not change with orbital position, thus it is with the addition of diurnal stresses that peak stresses are reached. At apogee, diurnal and recession stresses are most compressive near the tidal axis, while at perigee they are most compressive 90 degrees away from the tidal axis. Coseismic slip events on currently active thrust faults are expected to be triggered when peak stresses are reached. Analysis of the timing of the 28 the shallow moonquakes recorded by the Apollo seismic network shows that 19 indeed occur when the Moon is closer to apogee, while only 9 shallow events occur when the Moon is closer to perigee. Here we show the results of relocating the shallow moonquake using an algorithm designed for sparse networks to better constrain their epicentral locations in order to compare them with stress models. The model for the current stress state of the Moon is refined by investigating the contribution of polar wander.

siesmology↗

KELT-22Ab: A Massive, Short-Period Hot Jupiter Transiting a Near-solar Twin

We present the discovery of KELT-22Ab, a hot Jupiter from the KELT-South survey. KELT-22Ab transits the moderately bright (V11.1) Sun-like G2V star TYC 7518-468-1. The planet has an orbital period of P = 1.3866529 0.0000027 days, a radius of = -R 1.285+ R P 0.071 J 0.12 , and a relatively large mass of = - M 3.47+ M P 0.14 J 0.15 . The star has = - + R 1.099 0.046 R 0.079 , = - + M 1.092 0.041 M 0.045 , = - T 5767+ eff 49 50 K, = - log g 4.393+0.0600.039 (cgs), and [m/H]=+ -0.259+0.0830.085; thus other than its slightly super-solar metallicity, it appears to be a near-solar twin. Surprisingly, KELT-22A exhibits kinematics and a Galactic orbit that are somewhat atypical for thin-disk stars. Nevertheless, the star is rotating rapidly for its estimated age, and shows evidence of chromospheric activity. Imaging reveals a slightly fainter companion to KELT-22A that is likely bound, with a projected separation of 6″ (∼1400 au). In addition to the orbital motion caused by the transiting planet, we detect a possible linear trend in the radial velocity of KELT-22A, suggesting the presence of another relatively nearby body that is perhaps non-stellar. KELT-22Ab is highly irradiated (as a consequence of the small semimajor axis of a R = 4.97), and is mildly inflated. At such small separations, tidal forces become significant. The configuration of this system is optimal for measuring the rate of tidal dissipation within the host star. Our models predict that, due to tidal forces, the semimajor axis is decreasing rapidly, and KELT-22Ab is predicted to spiral into the star within the next Gyr.

Knicole Colon↗

Shallow Lunar Seismic Activity and the Current Stress State of the Moon

A vast, global network of more than 3200 lobate thrust fault scarps has been revealed in high resolution Lunar Reconnaissance Orbiter Camera (LROC) images. The fault scarps are very young, less than 50 Ma, based on their small scale and crisp appearance, crosscutting relations with small-diameter impact craters, and rates of infilling of associated small, shallow graben and may be actively forming today. The population of young thrust fault scarps provides a window into the recent stress state of the Moon and offers insight into the origin of global lunar stresses. The distribution of orientations of the fault scarps is non-random, inconsistent with isotropic stresses from late-stage global contraction as the sole source of stress. Modeling shows that tidal stresses contribute significantly to the current stress state of the lunar crust. Tidal stresses (orbital recession and diurnal tides) superimposed on stresses from global contraction result in non-isotropic compressional stress and may produce thrust faults consistent with lobate scarp orientations. At any particular point on the lunar surface, peak compressive stress will be reached at a certain time in the diurnal cycle. Coseismic slip events on currently active thrust faults are expected to be triggered when peak stresses are reached. Analysis of the timing of the 28 the shallow moonquakes recorded by the Apollo seismic network shows that 19 indeed occur when the Moon is closer to apogee, while only 9 shallow events occur when the Moon is closer to perigee. Here we report efforts to refine the model for the current stress state of the Moon by investigating the contribution of polar wander. Progress on relocating the epicentral locations of the shallow moonquakes using an algorithm designed for sparse networks is also reported.

structure↗

StreamGen: Connecting Populations of Streams and Shells to Their Host Galaxies

In this work, we study how the abundance and dynamics of populations of disrupting satellite galaxies change systematically as a function of host galaxy properties. We apply a theoretical model of the phase-mixing process to classify intact satellite galaxies and stellar streamlike and shell-like debris in ∼1500 Milky Way–mass systems generated by a semi-analytic galaxy formation code, SatGen. In particular, we test the effect of host galaxy halo mass, disk mass, ratio of disk scale height to length, and stellar feedback model on disrupting satellite populations. We find that the counts of tidal debris are consistent across all host galaxy models, within a given host mass range, and that all models can have streamlike debris on low-energy orbits, consistent with that observed around the Milky Way. However, we find a preference for streamlike debris on lower-energy orbits in models with a thicker (lower-density) host disk or on higher-energy orbits in models with a more massive host disk. Importantly, we observe significant halo-to-halo variance across all models. These results highlight the importance of simulating and observing large samples of Milky Way–mass galaxies and accounting for variations in host properties when using disrupting satellites in studies of near-field cosmology.

dark matter↗

Evolution of satellite resonances by tidal dissipation.

Analysis of a realistic model shows how satellites' gravitational interaction can halt their differential tidal evolution when resonant commensurabilities of their orbital periods are reached. The success of this study lends support to the hypothesis that orbit-orbit resonances among satellites in the solar system, including the Titan-Hyperion case, did evolve as a result of tidal energy dissipation. Consideration of the time scale for this evolution process, possible now that the capture mechanism has been revealed, can offer more sophisticated constraints on the tidal dissipation function, Q, and on past orbital conditions.

Greenberg, R.↗

From Existing and New Nuclear and Astrophysical Constraints to Stringent Limits on the Equation of State of Neutron-Rich Dense Matter

Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star (NS) observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. To constrain the EOS across its entire density range, this information needs to be combined consistently. However, the impact and model dependency of individual observations vary. Given their growing number, assessing the various methods is crucial to compare the respective effects on the EOS and discover potential biases. For this purpose, we present a broad compendium of different constraints and apply them individually to a large set of EOS candidates within a Bayesian framework. Specifically, we explore different ways of how chiral effective field theory and perturbative quantum chromodynamics can be used to place a likelihood on EOS candidates. We also investigate the impact of nuclear experimental constraints, as well as different radio and x-ray observations of NS masses and radii. This is augmented by reanalyses of the existing data from binary neutron star coalescences, in particular of GW170817, with improved models for the tidal waveform and kilonova light curves, which we also utilize to construct a tight upper limit of 2.39 M ⊙ on the TOV mass based on GW170817’s remnant. Our diverse set of constraints is eventually combined to obtain stringent limits on NS properties. We organize the combination in a way to distinguish between constraints where the systematic uncertainties are deemed small and those that rely on less conservative assumptions. For the former, we find the radius of the canonical 1.4 M ⊙ neutron star to be R 1.4 = 12.2 6 − 0.91 + 0.80 km and the TOV mass at M TOV = 2.2 5 − 0.22 + 0.42 M ⊙ (95% credibility). Including all the presented constraints yields R 1.4 = 12.2 0 − 0.48 + 0.50 km and M TOV = 2.3 0 − 0.20 + 0.07 M ⊙ . When comparing these limits to individual data points, we find that the quoted radius of HESS J1731-347 displays noticeable tension with other constraints. Constraining microphysical properties of the EOS proves more challenging. For instance, the symmetry energy slope is restricted to L sym = 48 − 25 + 21 MeV , where this constraint is mainly dominated by our reanalysis of the PREX-II and CREX experiment. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Libration in the earth's rotation

External luni-solar torque exerted on the difference (B-A) of the earth's two equatorial principal moments of inertia gives rise to two types of librational motions in the earth's rotation: the semidiurnal libration in spin and the prograde diurnal libration in polar motion. Formulas for the librations considering a realistic earth model and their tidal decompositions are derived and evaluated. The spin libration has a maximum peal-to-peak amplitude of 0.90 milliarcseconds, that of the polar libration is 0.06 milliarcseconds. Implications concerning their detectability and role in the tidal variation of earth rotation are discussed.

Chao, B. F.↗

Mars Global Reference Atmospheric Model (Mars-GRAM): Release No. 2 - Overview and applications

The Mars Global Reference Atmospheric Model (Mars-GRAM), a science and engineering model for empirically parameterizing the temperature, pressure, density, and wind structure of the Martian atmosphere, is described with particular attention to the model's newest version, Mars-GRAM, Release No. 2 and to the improvements incorporated into the Release No. 2 model as compared with the Release No. 1 version. These improvements include (1) an addition of a new capability to simulate local-scale Martian dust storms and the growth and decay of these storms; (2) an addition of the Zurek and Haberle (1988) wave perturbation model, for simulating tidal perturbation effects; and (3) a new modular version of Mars-GRAM, for incorporation as a subroutine into other codes.

James, B.↗

Simulations of the Surface Pressure Tides on Mars Incorporating New MOLA Topography Data

Simulations of surface pressure tides on Mars have been conducted previously with the NASA Ames Mars General Circulation Model (MGCM), as well as by Wilson and Hamilton. In the case of the MGCM, annual simulations with realistic (Viking-derived) dust loading have been reasonably successful at reproducing observed tides at the Viking Lander sites. One obvious discrepancy between observed and simulated tides is the amplification found in the model's diurnal tidal amplitude around L(sub s) 90. The observations show no such amplification. We report here on recent simulations of the surface pressure tide on Mars conducted with the MGCM and with the newly acquired MOLA topography data. The data is provided at 1deg x 1deg resolution, and is used by the MGCM at 7.5deg x 9deg resolution. The vertical domain of the MGCM extends to around 80 km. In the simulations reported in, the vertical domain extended to around 50 km, and the Consortium topography dataset was used. Additional information is contained in the original extended abstract.

Bridger, A. F. C.↗

The Mountains of Io: Global and Geological Perspectives from Voyager and Galileo

To search for local and global scale geologic associations that may be related to the internal dynamics of Io, we have completed a global catalog of all mountains and volcanic centers. We have identified 115 mountain structures (covering approx. 3% of the surface) and 541 volcanic centers, including paterae (calderas and dark spots) and shield volcanoes. The average length of an Ionian mountain is 157 km, with the longest being 570 km. The mean height of Ionian mountains is 6.3 km, and the highest known mountain is Boosaule Montes (17.5 +/- 3 km). Five basic morphologic types of mountains have been identified; mesa, plateau peak, ridge, and massif. Very few mountains bear any physical similarity. to classic volcanic landforms, but many resemble flatiron mountains on Earth and are interpreted as tilted crustal blocks. This would be consistent with the hypothesis that most mountains are thrust blocks formed as a result of compressive stresses built up in the lower crust due to the global subsidence of volcanic layers as they are buried over time. More than one mechanism may be responsible for all Ionian mountains, however. The proximity of some mountains to paterae may indicate a direct link between some mountains and volcanism, although it is not always clear which came first. In contrast to earlier studies, a pronounced bimodal pattern is observed in the global distribution of both mountains and volcanic centers. The regions of highest areal densities of volcanic centers are near the sub- and anti-Jovian regions, but are offset roughly 90deg in longitude from the two, regions of greatest concentration of mountains. This anticorrelation may indicate the overprinting of a second stress field on the global compressive stresses due to subsidence. The bimodal distribution of volcanic centers and mountains is consistent with models of asthenospheric tidal heating and internal convection developed by Tackley et al.Over regions of mantle upwelling, compressive stresses in the lower crust induced by global subsidence might be reduced, encouraging volcanism and discouraging mountain building. In regions of mantle downwelling, these compressive stresses in the lower crust might be increased, discouraging volcanism and encouraging mountain building. Alternatively, the global pattern may be related to possible (but undocumented) nonsynchronous rotation of lo, which would produce two regions each of compression and extension in the crust. Evidence of layering and of mass wasting, including landslides, block sliding, debris aprons and downslope creep, on Ionian mountains suggests that the crust of Io is essentially a layered stack of partially consolidated volcanic lavas and plume deposits, becoming more consolidated with depth. The lower crust especially may also be ductily deformed, punctuated by volcanic intrusions and faulting at paterae, and broken into blocks, some of which have been uplifted to form mountains.

Schenk, Paul↗

Saturn's Rings

The rings are changing before our eyes; structure varies on all timescales and unexpected things have been discovered. Many questions have been answered, but some answers remain elusive. Here we highlight the major ring science progress over the mission to date, and describe new observations planned for Cassini’s final three years. Ring Composition and particle sizes: The rings are nearly all water ice with no other ices – so why are they reddish? The C Ring and Cassini Division are “dirtier” than the more massive B and A Rings, as shown by near- IR and, recently, microwave observations. Particle sizes, from stellar and radio occultation's, vary from place to place. Ring structure, micro and macro: numerous spiral density waves and ubiquitous “self-gravity wakes” reveal processes which fostered planet formation in the solar system and elsewhere. However, big puzzles remain regarding the main ring divisions, the C Ring plateau structures, and the B Ring irregular structure. Moonlets, inside and out, seen and unseen: Two gaps contain sizeable moonlets, but more gaps seem to contain none; even smaller embedded “propeller” objects wander, systematically or randomly, through the A ring. Rubble pile ring moons just outside the rings may escaped from the rings, and the recently discovered “Peggy” may be trying this as we watch. Impact bombardment of the rings: Comet fragments set the rings to rippling on century-timescales, and boulders crash through hourly; meanwhile, the constant hail of in falling Kuiper belt material has a lower mass flux than previously thought. Origin and Age of the Rings: The ring mass and bombardment play key roles. The ring mass is well known everywhere but in the B Ring (where most of it is). New models suggest how tidal breakup of evolving moons may have formed massive ancient rings, of which the current ring is just a shadow. During its last three years, the Cassini tour profile will allow entirely new observations: direct measurement of the still-unknown ring mass; direct in-situ sampling of ring particle composition (targeting the iron- or carbon based red non-icy component); and radar backscattering observations.

Cassini↗

Dataset for scientific paper "Simulated plant‑mediated oxygen input has strong impacts on fine‑scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands", a modeling study based on field observation at the tidal salt marshes of the Parker River Estuary, Massachusetts, United States

This dataset is the raw and processed data for the paper "Simulated plant ‑ mediated oxygen input has strong impacts on fine ‑ scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands". This study investigated how plant-mediated oxygen input affects subsurface biogeochemical reactions of organic carbon degradation and the resulting methane emissions of coastal wetlands by model simulation. We used the subsurface geochemical simulator PFLOTRAN for the modeling, which produced the simulated changes in porewater chemical substances and methane emissions over 10 days under different scenarios of plant-mediated oxygen input.Specifically, this dataset contains: 1) the input files for PFLOTRAN of all simulation runs conducted in this study. Those files are with an extension of ".in", containing information of the biogeochemical reaction network (stoichiometry, reaction rate, Monod constants, etc), fluid flow rate and oxygen concentration in the fluid which together simulated the plant-mediated oxygen input, the configuration of artificial reactions that simulated the methane fluxes, etc. The PFLOTRAN input files are text files, which can be opened by NotePad, but running these input files will require proper installation of PFLOTRAN (instruction: https://documentation.pflotran.org/user_guide/how_to/installation/installation.html). 2) the raw and processed model output from PFLOTRAN of all simulation runs, and 3) the python scripts used to process the raw model output, including random allocation of root cells, converting raw data into organized formats, calculating the methane fluxes based on the model output, data visualization, etc. The raw and processed model output from PFLOTRAN are in .spydata format, which can be viewed with Python. and 3) the python scripts for data processing and analysis are programming scripts, which can be opened with Python.This modeling work, in particular the model parameterization of root density and initial conditions of porewater concentrations of biogeochemical substances, was based on field measurements at the salt marsh of the Upper Parker River Estuary, Massachusetts, United States.

54 ENVIRONMENTAL SCIENCES↗

TOPEX/POSEIDON tides estimated using a global inverse model

Altimetric data from the TOPEX/POSEIDON mission will be used for studies of global ocean circulation and marine geophysics. However, it is first necessary to remove the ocean tides, which are aliased in the raw data. The tides are constrained by the two distinct types of information: the hydrodynamic equations which the tidal fields of elevations and velocities must satisfy, and direct observational data from tide gauges and satellite altimetry. Here we develop and apply a generalized inverse method, which allows us to combine rationally all of this information into global tidal fields best fitting both the data and the dynamics, in a least squares sense. The resulting inverse solution is a sum of the direct solution to the astronomically forced Laplace tidal equations and a linear combination of the representers for the data functionals. The representer functions (one for each datum) are determined by the dynamical equations, and by our prior estimates of the statistics or errors in these equations. Our major task is a direct numerical calculation of these representers. This task is computationally intensive, but well suited to massively parallel processing. By calculating the representers we reduce the full (infinite dimensional) problem to a relatively low-dimensional problem at the outset, allowing full control over the conditioning and hence the stability of the inverse solution. With the representers calculated we can easily update our model as additional TOPEX/POSEIDON data become available. As an initial illustration we invert harmonic constants from a set of 80 open-ocean tide gauges. We then present a practical scheme for direct inversion of TOPEX/POSEIDON crossover data. We apply this method to 38 cycles of geophysical data records (GDR) data, computing preliminary global estimates of the four principal tidal constituents, M(sub 2), S(sub 2), K(sub 1) and O(sub 1). The inverse solution yields tidal fields which are simultaneously smoother, and in better agreement with altimetric and ground truth data, than previously proposed tidal models. Relative to the 'default' tidal corrections provided with the TOPEX/POSEIDON GDR, the inverse solution reduces crossover difference variances significantly (approximately 20-30%), even though only a small number of free parameters (approximately equal to 1000) are actually fit to the crossover data.

Egbert, Gary D.↗

Topographic enhancement of tidal motion in the western Barents Sea

A high-resolution numerical lattice is used to study a topographically trapped motion around islands and shallow banks of the western Barents Sea caused both by the semidiurnal and diurnal tidal waves. Observations and model computations in the vicinity of Bear Island show well-developed trapped motion with distinctive tidal oscillatory motion. Numerical investigations demonstrate that one source of the trapped motion is tidal current rectification over shallow topgraphy. Tidal motion supports residual currents of the order of 8 cm/s around Bear Island and shallow Spitsbergenbanken. The structures of enhanced tidal currents for the semidiurnal components are generated in the shallow areas due to topographic amplification. In the diurnal band of oscillations the maximum current is associated with the shelf wave occurrence. Residual currents due to diurnal tides occur at both the shallow areas and the shelf slope in regions of maximum topographic gradients. Surface manifestation of the diurnal current enhancement is the local maximum of tidal amplitude at the shelf break of the order of 5 to 10 cm. Tidal current enhancement and tidally generated residual currents in the Bear Island and Spitsbergenabanken regions cause an increased generation of ice leads, ridges and, trapped motion of the ice floes.

Kowalik, Z.↗

Numerical simulations of interacting disk galaxies

Galaxy-galaxy interactions have long attracted many extragalactic astronomers in various aspects. A number of computer simulations performed in the 1970s have successfully reproduced the peculiar morphologies observed in interacting disk galaxies and clarified that tidal deformation explains most of the observed global peculiarities. However, most of these simulations have used test particles in modelling the disk component. Tidal response of a self-gravitating disk remains to be further clarified. Another topic which is intensely discussed at present is the relation between galaxy-galaxy interactions and activity. Many observations suggest that interactions trigger strong starbursts and possibly active galactic nuclei (AGN). However, the detailed mechanism of triggering is not yet clear. It is vital here to understand the dynamics of interstellar gas. In order to understand various phenomena related to galaxy-galaxy interactions (mainly for disk galaxies), the author performed a series of numerical simulations on close galaxy encounters which includes both interstellar gas and self-gravitating disk components. In these simulations, the galaxy model to be perturbed (target galaxy) consists of a halo and a disk. The halo was treated as a rigid spherical gravitational field which is assumed to remain fixed during the interaction. The disk is composed of stars and gas. The stellar disk was constructed by 20000 collisionless particles of the same mass. Those particles move in the halo gravitational field, interacting with each other and with the perturber. Therefore, the self-graviy of the disk is properly taken into account. Stellar particles were initially given circular velocities with small random motions required to stabilize the disk against local axisymmetric disturbances. The gravitational field of the stellar disk was calculated by the particle-mesh scheme (e.g. Hockney and Eastwood 1981). The gaseous component was modelled by the cloud-particle scheme (e.g. Roberts and Hausman 1984). Here, the authors represent the gas as an ensemble of small spheres (i.e. clouds) and include the creation of an OB star in a cloud-cloud collision and subsequent velocity push on nearby clouds due to a supernova explosion.

Noguchi, Masafumi↗