SEARCH · Engineering Papers

Results for “monotonicity”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,405 records · Page 78

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Viscoplastic Model Development with an Eye Toward Characterization

A viscoplastic theory is developed that reduces analytically to creep theory under steady-state conditions. A viscoplastic model is constructed within this theoretical framework by defining material functions that have close ties to the physics of inelasticity. As a consequence, this model is easily characterized-only steady-state creep data, monotonic stress-strain curves, and saturated stress-strain hysteresis loops are required.

Freed, Alan D.↗

NASA NTRS1995

Unoccupied electronic resonances of Sc adsorbed on W(001) by k-resolved inverse photoemission

Scandium adsorbed on the (001) face of tungsten has been studied using Auger-electron spectroscopy, low-energy electron diffraction, k-resolved inverse photoelectron spectroscopy, work-function measurements, and relativistic-electronic-structure calculations. We find that the work function of W(001) does not decrease monotonically as a function of scandium adsorption but reaches a minimum value of Phi = 3.25 eV at 50% coverage. For a complete monolayer, an order 1 X 1 Sc overlayer is formed. Inverse photoemission in the isochromat mode was used to map the unoccupied energy levels of this overlayer along the bar-Gamma(bar-Mu) symmetry direction of the surface Brillouin zone. Local-density-functional calculations using the muffin-tin orbital method were performed for a 1 X 1 Sc overlayer on W and are compared with the experimental two-dimensional band structure.

Lamouri, A.↗

NASA NTRS1994

Temperature Field During Flame Spread over Alcohol Pools: Measurements and Modelling

A principal difference between flame spread over solid fuels and over liquid fuels is, in the latter case, the presence of liquid-phase convection ahead of the leading edge of the flame. The details of the fluid dynamics and heat transfer mechanisms in both the pulsating and uniform flame spread regimes were heavily debated, without resolution, in the 1960s and 1970s; recently, research on flame spread over pools was reinvigorated by the advent of enhanced diagnostic techniques and computational power. Temperature fields in the liquid, which enable determination of the extent of preheating ahead of the flame, were determined previously by the use of thermocouples and repetitive tests, and suggested that the surface temperature does not decrease monotonically ahead of the pulsating flame front, but that there exists a surface temperature valley. Recent predictions support this suggestion. However, others' thermocouple measurements and the recent field measurements using Holographic Interferometry (HI) did not find a similar valley. In this work we examine the temperature field using Rainbow Schlieren Deflectometry (RSD), with a measurement threshold exceeding that of conventional interferometry by a factor of 20:1, for uniform and pulsating flame spread using propanol and butanol as fuels. This technique was not applied before to flame spread over liquid pools, except in some preliminary measurements reported earlier. Noting that HI is sensitive to the refractive index while RSD responds to refractive index gradients, and that these two techniques might therefore be difficult to compare, we utilized a numerical simulation, described below, to predict and compare both types of field for the uniform and pulsating spread regimes. The experimental data also allows a validation of the model at a level of detail greater than has been attempted before.

Miller, Fletcher J.↗

NASA NTRS1995

Experimental Verification of a Progressive Damage Model for IM7/5260 Laminates Subjected to Tension-Tension Fatigue

The durability and damage tolerance of laminated composites are critical design considerations for airframe composite structures. Therefore, the ability to model damage initiation and growth and predict the life of laminated composites is necessary to achieve structurally efficient and economical designs. The purpose of this research is to experimentally verify the application of a continuum damage model to predict progressive damage development in a toughened material system. Damage due to monotonic and tension-tension fatigue was documented for IM7/5260 graphite/bismaleimide laminates. Crack density and delamination surface area were used to calculate matrix cracking and delamination internal state variables to predict stiffness loss in unnotched laminates. A damage dependent finite element code predicted the stiffness loss for notched laminates with good agreement to experimental data. It was concluded that the continuum damage model can adequately predict matrix damage progression in notched and unnotched laminates as a function of loading history and laminate stacking sequence.

Coats, Timothy W.↗

NASA NTRS1995

Cyclic Axial-Torsional Deformation Behavior of a Cobalt-Base Superalloy

The cyclic, high-temperature deformation behavior of a wrought cobalt-base super-alloy, Haynes 188, is investigated under combined axial and torsional loads. This is accomplished through the examination of hysteresis loops generated from a biaxial fatigue test program. A high-temperature axial, torsional, and combined axial-torsional fatigue database has been generated on Haynes 188 at 760 C. Cyclic loading tests have been conducted on uniform gage section tubular specimens in a servohydraulic axial-torsional test rig. Test control and data acquisition were accomplished with a minicomputer. The fatigue behavior of Haynes 188 at 760 C under axial, torsional, and combined axial-torsional loads and the monotonic and cyclic deformation behaviors under axial and torsional loads have been previously reported. In this paper, the cyclic hardening characteristics and typical hysteresis loops in the axial stress versus axial strain, shear stress ,versus engineering shear strain, axial strain versus engineering shear strain. and axial stress versus shear stress spaces are presented for cyclic in-phase and out-of-phase axial-torsional tests. For in-phase tests, three different values of the proportionality constant lambda (the ratio of engineering shear strain amplitude to axial strain amplitude, are examined, viz. 0.86, 1.73, and 3.46. In the out-of-phase tests, three different values of the phase angle, phi (between the axial and engineering shear strain waveforms), are studied, viz., 30, 60, and 90 degrees with lambda equals 1.73. The cyclic hardening behaviors of all the tests conducted on Haynes 188 at 760 C are evaluated using the von Mises equivalent stress-strain and the maximum shear stress-maximum engineering shear strain (Tresca) curves. Comparisons are also made between the hardening behaviors of cyclic axial, torsional, and combined in-phase (lambda = 1.73 and phi = 0) and out-of-phase (lambda = 1.73 and phi = 90') axial-torsional fatigue tests. These comparisons are accomplished through simple Ramberg-Osgood type stress-strain functions for cyclic, axial stress-strain and shear stress-engineering shear strain curves.

Bonacuse, Peter J.↗

NASA NTRS1996

Multigrid Method for Modeling Multi-Dimensional Combustion with Detailed Chemistry

A highly accurate and efficient numerical method is developed for modeling 3-D reacting flows with detailed chemistry. A contravariant velocity-based governing system is developed for general curvilinear coordinates to maintain simplicity of the continuity equation and compactness of the discretization stencil. A fully-implicit backward Euler technique and a third-order monotone upwind-biased scheme on a staggered grid are used for the respective temporal and spatial terms. An efficient semi-coarsening multigrid method based on line-distributive relaxation is used as the flow solver. The species equations are solved in a fully coupled way and the chemical reaction source terms are treated implicitly. Example results are shown for a 3-D gas turbine combustor with strong swirling inflows.

Zheng, Xiaoqing↗

NASA NTRS1995

MHD Simulations of Earth's Bow Shock at low Mach Numbers: Standoff Distances

Global, three-dimensional, ideal MHD simulations of Earth's bow shock are reported for low Alfven Mach numbers M(sub A) and quasi-perpendicular magnetic field orientations. The simulations use a hard, infinitely conducting magnetopauause obstacle, with axisymmetric three-dimensional location given by scaled standard model, to directly address previous gasdynamic (GD) and field-aligned MHD (FA-MHD) work. Tests of the simulated shocks' density jumps X for 1.4 approx. less than MA approx. less than 10 and the high M(sub A) shock location, and reproduction of the GD relation between magnetosheath thickness and X for quasi-gasdynamic MHD runs with M(sub A) much greater than M(sub s), confirm that the MHD code is working correctly. The MHD simulations show the standoff distance a(sub s), increasing monotonically with decreasing M(sub A). Significantly larger a(sub s), are found at low M(sub A) than predicted by GD and phenomenological MHD models and FA-MHD simulations, as required qualitatively by observations. The GD and FA-MHD predictions err qualitatively, predicting either constant or decreasing a(sub s), with decreasing M(sub A). This qualitative difference between quasi- perpendicular MHD and FA-MHD simulations is direct evidence for a(sub s), depending on the magnetic field orientation Theta. The enhancement factor over the phenomenological MHD predictions at MA approx. 2.4 agrees quantitatively with one observatiorial estimate. A linear relationship is found between the magnetosheath thickness and X, modified both quantitatively and intrinsically by MHD effects from the GD result. The MHD and GD results agree in the high M(sub A) limit. An MHD theory is developed for a(sub s), restricted to sufficiently perpendicular Theta and high sonic Mach numbers M(sub s). It explains the simulation results with excellent accuracy. Observational and further simulation testing of this MHD theory, and of its predicted M(sub A), Theta, and M(sub s) effects, is desirable.

Cairns, Iver H.↗

NASA NTRS1997

Prediction and Verification of Ductile Crack Growth from Simulated Defects in Strength Overmatched Butt Welds

Defects that develop in welds during the fabrication process are frequently manifested as embedded flaws from lack of fusion or lack of penetration. Fracture analyses of welded structures must be able to assess the effect of such defects on the structural integrity of weldments; however, the transferability of R-curves measured in laboratory specimens to defective structural welds has not been fully examined. In the current study, the fracture behavior of an overmatched butt weld containing a simulated buried, lack-of-penetration defect is studied. A specimen designed to simulate pressure vessel butt welds is considered; namely, a center crack panel specimen, of 1.25 inch by 1.25 inch cross section, loaded in tension. The stress-relieved double-V weld has a yield strength 50% higher than that of the plate material, and displays upper shelf fracture behavior at room temperature. Specimens are precracked, loaded monotonically while load-CMOD measurements are made, then stopped and heat tinted to mark the extent of ductile crack growth. These measurements are compared to predictions made using finite element analysis of the specimens using the fracture mechanics code Warp3D, which models void growth using the Gurson-Tvergaard dilitant plasticity formulation within fixed sized computational cells ahead of the crack front. Calibrating data for the finite element analyses, namely cell size and initial material porosities are obtained by matching computational predictions to experimental results from tests of welded compact tension specimens. The R-curves measured in compact tension specimens are compared to those obtained from multi-specimen weld tests, and conclusions as to the transferability of R-curves is discussed.

Nishioka, Owen S.↗

NASA NTRS1994

Molecular Dynamics of a Water-Lipid Bilayer Interface

We present results of molecular dynamics simulations of a glycerol 1-monooleate bilayer in water. The total length of analyzed trajectories is 5ns. The calculated width of the bilayer agrees well with the experimentally measured value. The interior of the membrane is in a highly disordered fluid state. Atomic density profile, orientational and conformational distribution functions, and order parameters indicate that disorder increases toward the center of the bilayer. Analysis of out-of-plane thermal fluctuations of the bilayer surfaces occurring at the time scale of the present calculations reveals that the distribution of modes agrees with predictions of the capillary wave model. Fluctuations of both bilayer surfaces are uncorrelated, yielding Gaussian distribution of instantaneous widths of the membrane. Fluctuations of the width produce transient thinning defects in the bilayer which occasionally span almost half of the membrane. The leading mechanism of these fluctuations is the orientational and conformational motion of head groups rather than vertical motion of the whole molecules. Water considerably penetrates the head group region of the bilayer but not its hydrocarbon core. The total net excess dipole moment of the interfacial water points toward the aqueous phase, but the water polarization profile is non-monotonic. Both water and head groups significantly contribute to the surface potential across the interface. The calculated sign of the surface potential is in agreement with that from experimental measurements, but the value is markedly overestimated. The structural and electrical properties of the water-bilayer system are discussed in relation to membrane functions, in particular transport of ions and nonelectrolytes across membranes.

Wilson, Michael A.↗

NASA NTRS1997

Mechanical Characterization and Micromechanical Modeling of Woven Carbon/Copper Composites

The present investigation examines the in-plane mechanical behavior of a particular woven metal matrix composite (MMC); 8-harness (8H) satin carbon/copper (C/Cu). This is accomplished via mechanical testing as well as micromechanical modeling. While the literature is replete with experimental and modeling efforts for woven and braided polymer matrix composites, little work has been done on woven and braided MMC's. Thus, the development and understanding of woven MMC's is at an early stage. 8H satin C/Cu owes its existence to the high thermal conductivity of copper and low density and thermal expansion of carbon fibers. It is a candidate material for high heat flux applications, such as space power radiator panels. The experimental portion of this investigation consists of monotonic and cyclic tension, compression, and Iosipescu shear tests, as well as combined tension-compression tests. Tests were performed on composite specimens with three copper matrix alloy types: pure Cu, Cu-0.5 weight percent Ti (Cu-Ti), and Cu-0.7 weight percent Cr (Cu-Cr). The small alloying additions are present to promote fiber/matrix interfacial bonding. The analytical modeling effort utilizes an approach in which a local micromechanical model is embedded in a global micromechanical model. This approach differs from previously developed analytical models for woven composites in that a true repeating unit cell is analyzed. However, unlike finite element modeling of woven composites, the geometry is sufficiently idealized to allow efficient geometric discretization and efficient execution.

Bednarcyk, Brett A.↗

NASA NTRS1997

Nonlinear Performance Seeking Control using Fuzzy Model Reference Learning Control and the Method of Steepest Descent

Performance Seeking Control (PSC) attempts to find and control the process at the operating condition that will generate maximum performance. In this paper a nonlinear multivariable PSC methodology will be developed, utilizing the Fuzzy Model Reference Learning Control (FMRLC) and the method of Steepest Descent or Gradient (SDG). This PSC control methodology employs the SDG method to find the operating condition that will generate maximum performance. This operating condition is in turn passed to the FMRLC controller as a set point for the control of the process. The conventional SDG algorithm is modified in this paper in order for convergence to occur monotonically. For the FMRLC control, the conventional fuzzy model reference learning control methodology is utilized, with guidelines generated here for effective tuning of the FMRLC controller.

Kopasakis, George↗

NASA NTRS1997

Quasi-Static Viscoelasticity Loading Measurements of an Aircraft Tire

Stair-step loading, cyclic loading, and long-term relaxation tests were performed on an aircraft tire to observe the quasi-static viscoelastic response of the tire. The data indicate that the tire continues to respond viscoelastically even after it has been softened by deformation. Load relaxation data from the stair-step test at the 15,000-lb loading was fit to a monotonically decreasing Prony series.

Mason, Angela J.↗

NASA NTRS1995

Modified Interior Distance Functions (Theory and Methods)

In this paper we introduced and developed the theory of Modified Interior Distance Functions (MIDF's). The MIDF is a Classical Lagrangian (CL) for a constrained optimization problem which is equivalent to the initial one and can be obtained from the latter by monotone transformation both the objective function and constraints. In contrast to the Interior Distance Functions (IDF's), which played a fundamental role in Interior Point Methods (IPM's), the MIDF's are defined on an extended feasible set and along with center, have two extra tools, which control the computational process: the barrier parameter and the vector of Lagrange multipliers. The extra tools allow to attach to the MEDF's very important properties of Augmented Lagrangeans. One can consider the MIDFs as Interior Augmented Lagrangeans. It makes MIDF's similar in spirit to Modified Barrier Functions (MBF's), although there is a fundamental difference between them both in theory and methods. Based on MIDF's theory, Modified Center Methods (MCM's) have been developed and analyzed. The MCM's find an unconstrained minimizer in primal space and update the Lagrange multipliers, while both the center and the barrier parameter can be fixed or updated at each step. The MCM's convergence was investigated, and their rate of convergence was estimated. The extension of the feasible set and the special role of the Lagrange multipliers allow to develop MCM's, which produce, in case of nondegenerate constrained optimization, a primal and dual sequences that converge to the primal-dual solutions with linear rate, even when both the center and the barrier parameter are fixed. Moreover, every Lagrange multipliers update shrinks the distance to the primal dual solution by a factor 0 less than gamma less than 1 which can be made as small as one wants by choosing a fixed interior point as a 'center' and a fixed but large enough barrier parameter. The numericai realization of MCM leads to the Newton MCM (NMCM). The approximation for the primal minimizer one finds by Newton Method followed by the Lagrange multipliers update. Due to the MCM convergence, when both the center and the barrier parameter are fixed, the condition of the MDF Hessism and the neighborhood of the primal ninimizer where Newton method is 'well' defined remains stable. It contributes to both the complexity and the numerical stability of the NMCM.

Polyak, Roman A.↗

NASA NTRS1997

Micromechanical Modeling of Woven Metal Matrix Composites

This report presents the results of an extensive micromechanical modeling effort for woven metal matrix composites. The model is employed to predict the mechanical response of 8-harness (8H) satin weave carbon/copper (C/Cu) composites. Experimental mechanical results for this novel high thermal conductivity material were recently reported by Bednarcyk et al. along with preliminary model results. The micromechanics model developed herein is based on an embedded approach. A micromechanics model for the local (micro-scale) behavior of the woven composite, the original method of cells (Aboudi), is embedded in a global (macro-scale) micromechanics model (the three-dimensional generalized method of cells (GMC-3D) (Aboudi). This approach allows representation of true repeating unit cells for woven metal matrix composites via GMC-3D, and representation of local effects, such as matrix plasticity, yarn porosity, and imperfect fiber-matrix bonding. In addition, the equations of GMC-3D were reformulated to significantly reduce the number of unknown quantities that characterize the deformation fields at the microlevel in order to make possible the analysis of actual microstructures of woven composites. The resulting micromechanical model (WCGMC) provides an intermediate level of geometric representation, versatility, and computational efficiency with respect to previous analytical and numerical models for woven composites, but surpasses all previous modeling work by allowing the mechanical response of a woven metal matrix composite, with an elastoplastic matrix, to be examined for the first time. WCGMC is employed to examine the effects of composite microstructure, porosity, residual stresses, and imperfect fiber-matrix bonding on the predicted mechanical response of 8H satin C/Cu. The previously reported experimental results are summarized, and the model predictions are compared to monotonic and cyclic tensile and shear test data. By considering appropriate levels of porosity, residual stresses, and imperfect fiber-matrix debonding, reasonably good qualitative and quantitative correlation is achieved between model and experiment.

Bednarcyk, Brett A.↗

NASA NTRS1995

Effect of Processing Route on Strain Controlled Low Cycle Fatigue Behavior of Polycrystalline NiAl

The present investigation examines the effects of manufacturing process on the total axial strain controlled low cycle fatigue behavior of polycrystalline NiAl at 1000 K, a temperature above the monotonic Brittle-to-Ductile Transition Temperature (BDTT). The nickel aluminide samples were produced by three different processing routes: hot isostatic pressing of pre- alloyed powders, extrusion of prealloyed powders, and extrusion of vacuum induction melted ingots. The LCF behavior of the cast plus extruded material was also determined at room temperature (below the BD77) for comparison to the high temperature data. The cyclic stress response, cyclic stress-strain behavior, and strain-life relationships were influenced by the alloy preparation technique and the testing temperature. Detailed characterization of the LCF tested samples was conducted by optical and electron microscopy to determine the variations in fracture and deformation modes and to determine any microstructural changes that occurred during LCF testing. The dependence of LCF properties on processing route was rationalized on the basis of starting microstructure, brittle-to-ductile transition temperature, deformation induced changes in the basic microstructure, deformation substructure, and synergistic interaction between the damage modes.

Rao, K. Bhanu Sankara↗

NASA NTRS1997

Interchange for Joint Research Entitled: Measurement of Stable Nitrogen and Sulfur Isotopes

Viking measurements of the Martian atmosphere indicate a value of N-15/N-14 which is markedly greater than that found in Earth's atmosphere. These isotopic measurements provide a powerful diagnostic tool which may be used to derive valuable information regarding the past history of Mars and they have been used to place important constraints on the evolution of Mars' atmosphere. Initial partial pressures of nitrogen, outgassing rates, and integrated deposition of nitrogen into minerals have been calculated from this important atmospheric data (McElroy et al., 1976 and 1977; Fox and Dalgarno, 1983). The greater precision obtained in laser spectrometer isotopic measurements compared to the Viking data will greatly improve these calculated values. It has also been proposed that the N-15/N-14 value in Mars' atmosphere has increased monotonically over time (McElroy et al., 1977; Fox and Dalgarno, 1983; Wallis, 1989) owing to preferential escape of atmospheric 14N to space. Nitrogen isotopic ratios might be used to identify relatively ancient crustal rocks (R. Mancinelli, personal communication), and perhaps determine relative aces of surface samples. As a first step in successfully measuring nitrogen isotopes optically we have demonstrated the measurement of 15NI14N to a precision of 0.1% (See Figures 1-4) using a tunable diode laser and an available gas (N-,O) with spectral lines in the 2188 cm-1 region. The sample and reference gas cells contained gases of identical isotopic composition so that the 15NI14N absorption ratio determined from the sample cell, when divided by the 15NI14N absorption ratio determined from the reference cell, should yield an ideal value of unity. The average measured value of this "ratio of ratios" was 0.9983 with a standard deviation (20 values) of 0.0010. This corresponds to a precision of 0.1% (1 per mil) for nitrogen isotopes, a value sufficiently precise to provide important isotopic data of interest to exobiologists. The precision presently attainable in gases is sufficient to permit the instrument to be used in the measurement of isotopic ratios of interest to exobiologists as well as geologists and planetary scientists.

Becker, Joseph F.↗

NASA NTRS1996

Chemical Characteristics of Continental Outflow Over the Tropical South Atlantic Ocean from Brazil and Africa

The chemical characteristics of air parcels over the tropical South Atlantic during September - October 1992 are summarized by analysis of aged marine and continental outflow classifications. Positive correlations between CO and CH3CL and minimal enhancements of C2CL40, and various ChloroFluoroCarbon (CFC) species in air parcels recently advected over the South Atlantic basin strongly suggest an impact on tropospheric chemistry from biomass burning on adjacent continental areas of Brazil and Africa. Comparison of the composition of aged Pacific air with aged marine air over the South Atlantic basin from 0.3 to 12.5 km altitude indicates potential accumulation of long-lived species during the local dry season. This may amount to enhancements of up to two-fold for C2H6, 30% for CO, and 10% for CH3Cl. Nitric oxide and NO(x) were significantly enhanced (up to approx. 1 part per billion by volume (ppbv)) above 10 km altitude and poorly correlated with CO and CH3Cl. In addition, median mixing ratios of NO and NO(x) were essentially identical in aged marine and continental outflow air masses. It appears that in addition to biomass burning, lightning or recycled reactive nitrogen may be an important source of NO(x) to the upper troposphere. Methane exhibited a monotonic increase with altitude from approx. 1690 to 1720 ppbv in both aged marine and continental outflow air masses. The largest mixing ratios in the upper troposphere were often anticorrelated with CO, CH3Cl, and CO2, suggesting CH, contributions from natural sources. We also argue, based on CH4/CO ratios and relationships with various hydrocarbon and CFC species, that inputs from biomass burning and the northern hemisphere are unlikely to be the dominant sources of CO, CH4 and C2H6 in aged marine air. Emissions from urban areas would seem to be necessary to account for the distribution of at least CH4 and C2H6. Over the African and South American continents an efficient mechanism of convective vertical transport coupled with large-scale circulations conveys biomass burning, urban, and natural emissions to the upper troposphere over the South Atlantic basin. Slow subsidence over the eastern South Atlantic basin may play an important role in establishing and maintaining the rather uniform vertical distribution of long-lived species over this region. The common occurrence of values greater than 1 for the ratio CH3OOH/H2O2 in the upper troposphere suggests that precipitation scavenging effectively removed highly water soluble gases (H2O2, HNO3, HCOOH, and CH3COOH) and aerosols during vertical convective transport over the continents. However, horizontal injection of biomass burning products over the South Atlantic, particularly water soluble species and aerosol particles, was frequent below 6 km altitude.

Talbot, R. W.↗

NASA NTRS1998

Statistical Properties of Line Centroid Velocity Increments in the rho Ophiuchi Cloud

We present a comparison of histograms of CO (2-1) line centroid velocity increments in the rho Ophiuchi molecular cloud with those computed for spectra synthesized from a three-dimensional, compressible, but non-starforming and non-gravitating hydrodynamic simulation. Histograms of centroid velocity increments in the rho Ophiuchi cloud show clearly non-Gaussian wings, similar to those found in histograms of velocity increments and derivatives in experimental studies of laboratory and atmospheric flows, as well as numerical simulations of turbulence. The magnitude of these wings increases monotonically with decreasing separation, down to the angular resolution of the data. This behavior is consistent with that found in the phase of the simulation which has most of the properties of incompressible turbulence. The time evolution of the magnitude of the non-Gaussian wings in the histograms of centroid velocity increments in the simulation is consistent with the evolution of the vorticity in the flow. However, we cannot exclude the possibility that the wings are associated with the shock interaction regions. Moreover, in an active starforming region like the rho Ophiuchi cloud, the effects of shocks may be more important than in the simulation. However, being able to identify shock interaction regions in the interstellar medium is also important, since numerical simulations show that vorticity is generated in shock interactions.

Lis, D. C.↗