Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model diagnostics”

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,657 records · Page 92

Simulations of future particle accelerators: issues and mitigations

The ever increasing demands placed upon machine performance have resulted in the need for more comprehensive particle accelerator modeling. Computer simulations are key to the success of particle accelerators. Many aspects of particle accelerators rely on computer modeling at some point, sometimes requiring complex simulation tools and massively parallel supercomputing. Examples include the modeling of beams at extreme intensities and densities (toward the quantum degeneracy limit), and with ultra-fine control (down to the level of individual particles). In the future, adaptively tuned models might also be relied upon to provide beam measurements beyond the resolution of existing diagnostics. Much time and effort has been put into creating accelerator software tools, some of which are highly successful. However, there are also shortcomings such as the general inability of existing software to be easily modified to meet changing simulation needs. In this paper possible mitigating strategies are discussed for issues faced by the accelerator community as it endeavors to produce better and more comprehensive modeling tools. This includes lack of coordination between code developers, lack of standards to make codes portable and/or reusable, lack of documentation, among others.

43 PARTICLE ACCELERATORS↗

MHD Modeling of Coronal Loops: the Transition Region Throat

Context. The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross-sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. Aims. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. Methods. We study the area response with a time-dependent 2D magnetohydrodynamic (MHD) loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction, and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 millikelvin. Results. We find that the area can change substantially with the quasi-steady heating rate, e.g., by approx. 40% at 0.5 millikelvin as the loop temperature varies between 1 millikelvin and 4 millikelvin, and, therefore, affects the interpretation of the differential emission measure vs. temperature (DEM(T)) curves.

magnetohydrodynamics (MHD)↗

Model-based real-time surface heat flux and temperature estimation for the DIII-D tokamak

A control-oriented model for monitoring of wall power flux densities on the DIII-D tokamak has been successfully implemented and validated experimentally. Future reactors will have to withstand severe steady state high heat flux loads on plasma-facing components (PFCs). Due to the difficulty of directly-measuring local heat fluxes on these components, monitoring and protection of PFCs during the plasma discharge can benefit from simplified physics-based real-time (RT) functional models to estimate and guide heat load control. As a first step into the development, a control-oriented model for monitoring of wall power flux densities and temperatures on DIII-D tokamak has been successfully implemented. The paper discusses the experimental demonstration and comparison of the 2-D model-based wall heat flux algorithm on the DIII-Dinner wall limiter (IWL) against infra-red camera heat flux measurements for limited plasma configurations. The paper also reports on the benchmarking of the field line tracing environment, SMITER, developed at ITER organization on DIII-D tokamak against experimental IR diagnostic data and the derivation of the component shaping weighting factors for the 2-D model-based approach. Here, the extension of the model-based approach for surface temperature estimation on the DIII-D IWL is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Formation of IRIS Diagnostics. IX. The Formation of the C i 135.58 NM Line in the Solar Atmosphere

The C i 135.58 nm line is located in the wavelength range of NASA’s Interface Region Imaging Spectrograph ( IRIS ) small explorer mission. We study the formation and diagnostic potential of this line by means of non local-thermodynamic-equilibrium modeling, employing both 1D and 3D radiation-magnetohydrodynamic models. The C i/C ii ionization balance is strongly influenced by photoionization by Ly α emission. The emission in the C i 135.58 nm line is dominated by a recombination cascade and the line forming region is optically thick. The Doppler shift of the line correlates strongly with the vertical velocity in its line forming region, which is typically located at 1.5 Mm height. With IRIS , the C i 135.58 nm line is usually observed together with the O i 135.56 nm line, and from the Doppler shift of both lines, we obtain the velocity difference between the line forming regions of the two lines. From the ratio of the C i/O i line core intensity, we can determine the distance between the C i and the O i forming layers. Combined with the velocity difference, the velocity gradient at mid-chromospheric heights can be derived. The C i/O i total intensity line ratio is correlated with the inverse of the electron density in the mid-chromosphere. We conclude that the C i 135.58 nm line is an excellent probe of the middle chromosphere by itself, and together with the O i 135.56 nm line the two lines provide even more information, which complements other powerful chromospheric diagnostics of IRIS such as the Mg ii h and k lines and the C ii lines around 133.5 nm.

79 ASTRONOMY AND ASTROPHYSICS↗

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗

Bridging the PSI Knowledge Gap: A Multi-Scale Approach

Plasma-surface interactions (PSI) pose an immense scientific hurdle in magnetic confinement fusion and our present understanding of PSI in confinement environments is highly inadequate; indeed, a recent Fusion Energy Sciences Advisory Committee report found that 4 out of the 5 top five fusion knowledge gaps were related to PSI. The time is appropriate to develop a concentrated and synergistic science effort that would expand, exploit and integrate the wealth of laboratory ion-beam and plasma research, as well as exciting new computational tools, towards the goal of bridging the PSI knowledge gap. This effort would broadly advance plasma and material sciences, while providing critical knowledge towards progress in fusion PSI. This project involves the development of a Science Center focused on a new approach to PSI science; an approach that both exploits access to state-of-the-art PSI experiments and modeling, as well as confinement devices. The organizing principle is to develop synergistic experimental and modeling tools that treat the truly coupled multi-scale aspect of the PSI issues in confinement devices. This is motivated by the simple observation that while typical lab experiments and models allow independent manipulation of controlling variables, the confinement PSI environment is essentially self-determined with few outside controls. This means that processes that may be treated independently in laboratory experiments, because they involve vastly different physical and time scales, will now affect one another in the confinement environment. Also, lab experiments cannot simultaneously match all exposure conditions found in confinement devices typically forcing a linear extrapolation of lab results. At the same time programmatic limitations prevent confinement experiments alone from answering many key PSI questions. The resolution to this problem is to usefully exploit access to PSI science in lab devices, while retooling our thinking from a linear and de-coupled extrapolation to a multi-scale, coupled approach. The PSI Plasma Center consisted of three equal co-centers; one located at the MIT Plasma Science and Fusion Center, one at UC San Diego Center for Energy Research and one at the UC Berkeley Department of Nuclear Engineering, which moved to the University of Tennessee, Knoxville (UTK) with Professor Brian Wirth in July 2010. The Center had three co-directors: Prof. Dennis Whyte led the MIT co-center, the UCSD co-center was led by Dr. Russell Doerner, and Prof. Brian Wirth led the UCB/UTK center. The directors have extensive experience in PSI and material research, and have been internationally recognized in the magnetic fusion, materials and plasma research fields. The co-centers feature keystone PSI experimental and modeling facilities dedicated to PSI science: the DIONISOS/CLASS facility at MIT, the PISCES facility at UCSD, and the state-of-the-art numerical modeling capabilities at UCB/UTK. A collaborative partner in the center is Sandia National Laboratory at Livermore (SNL/CA), which has extensive capabilities with low energy ion beams and surface diagnostics, as well as supporting plasma facilities, including the Tritium Plasma Experiment, all of which significantly augment the Center. Interpretive, continuum material models are available through SNL/CA, UCSD and MIT. The participating institutions of MIT, UCSD, UCB/UTK, SNL/CA and LLNL brought a formidable array of experimental tools and personnel abilities into the PSI Plasma Center. Our work has focused on modeling activities associated with plasma surface interactions that are involved in effects of He and H plasma bombardment on tungsten surfaces. This involved performing computational material modeling of the surface evolution during plasma bombardment using molecular dynamics modeling. The principal outcomes of the research efforts within the combined experimental – modeling PSI center are to provide a knowledgebase of the mechanisms of surface degradation, and the influence of the surface on plasma conditions.

36 MATERIALS SCIENCE↗

A seasonal radiative-dynamic model of Saturn's troposphere

The present steady-state, linearized and zonally-averaged model of the thermal response of Saturn's atmosphere to seasonal insolation is compared with extant data in order to extract diagnostic information on the planet's troposphere and stratosphere. The planet's rings are shown to constitute a strong modulator of the seasonal insolation function. A detailed comparison of the model with existing temperature measurements shows that the large thermal inertia of Saturn and uncertainties in the wind and thermal measurements overwhelm the rings' thermal effects. A meridional asymmetry in heating is noted which is consistent with an enhancement of aerosols and of photochemically generated ethane and acetylene in the southern hemisphere's summer.

Barnet, C. D.↗

Development of spatially and temporally resolved electron density measurements for the assessment of radiation hydrodynamics simulations of laboratory X-ray photoionized plasmas

The photoionized plasma gas cell experiment is an established platform we use to make at-parameter (ξ > > 1 ergs cm s –1 ) measurements of plasma properties with application to high-energy astrophysical systems. We model the experiments with 1D radiation hydrodynamics simulations using the HELIOS-CR code to inform our understanding and assist in the interpretation of results. The simulations predict that the bulk of the plasma is in a quasi-uniform and hydrodynamically unperturbed state throughout the duration of the experiment. To evaluate this prediction, we introduced a photonic Doppler velocimetry (PDV) diagnostic to measure spatially and temporally resolved plasma electron density. The initial measurements were successful but had limitations that made model-data comparisons challenging. To address this, we re-designed the gas cell PDV diagnostic and doubled the number of measurement locations to sample across two thirds of the depth of the cell. We also present a comparison of the results from the upgraded PDV diagnostic to the HELIOS-CR simulations for the first time. As a result, the experimental data confirms the prediction of an unperturbed region in the bulk of the plasma but reveals discrepancies in the time evolution and spatial distribution of the simulated electron density.

79 ASTRONOMY AND ASTROPHYSICS↗

Photochemistry and dynamics of the ozone layer

The paper presents a broad review of the photochemical and dynamic theories of the ozone layer. The two theories are combined into the MIT three-dimensional dynamic-chemical quasi-geostrophic model with 26 levels in the vertical spaced in logarithmic pressure coordinates between the ground and 72-km altitude. The chemical scheme incorporates the important odd nitrogen, odd hydrogen, and odd oxygen chemistry, but is simplified in the sense that it requires specification of the distributions of NO2, OH and HO2. The prognostic equations are the vorticity equation, the perturbation thermodynamic equation, and the global mean and perturbation continuity equations for ozone; diagnostic equations include the hydrostatic equation, the balance condition, and the mass continuity equation. The model is applied to the investigation of the impact of supersonic aircraft on the ozone layer.

Prinn, R. G.↗

Cratering on Titan: A Pre-Cassini Perspective

The NASA-ESA Cassini mission, comprising a formidably instrumented orbiter and parachute-borne probe to be launched this October, promises to reveal a crater population on Titan that has been heretofore hidden by atmospheric haze. This population on the largest remaining unexplored surface in the solar system will be invaluable in comparative planetological studies, since it introduces evidence of the atmospheric effects of cratering on an icy satellite. Here, I highlight some impact features we may hope to find and could devote some modeling effort toward. Titan in a Nutshell: Radius= 2575 km. Density= 1880 kg/cubic m consistent with rock-ice composition. Surface pressure = 1.5 bar. Surface gravity = 1.35 m/square s Atmosphere -94% N2 6% CH, Surface temperature = 94K Tropopause temperature = 70K at 40 km alt. Probable liquid hydrocarbon deposits exist on or near the surface.Titan in a Nutshell: Radius= 2575 km. Density= 1880 kg/cubic m consistent with rock-ice composition. Surface pressure = 1.5 bar. Surface gravity = 1.35 m/square s; Atmosphere about 94% N2 6% CH, Surface temperature = 94K Tropopause temperature = 70K at 40 km alt. Probable liquid hydrocarbon deposits exist on or near the surface. Titan is comparable to Callisto and Ganymede for strength/gravity, Mars/Earth/Venus for atmospheric interaction, and Hyperion, Rhea, and Iapetus for impactor distribution. The leading/trailing asymmetry of crater density from heliocentric impactors is expected to be about 5-6, in the absence of resurfacing. Any Saturnocentric impactor population is likely to alter this. In particular the impact disruption of Hyperion is noted; because of the 3:4 orbital resonance with Titan, fragments from the proto-Hyperion breakup would have rapidly accreted onto Titan. Titan's resurfacing history is of course unknown. The disruption of impactors into fragments that individually create small craters is expected to occur. A crude estimate suggests a maximum separation of about 2 km (compared with 4 km on Venus, or 0.5 km on Earth). Crater chains are unlikely on Titan, since impactors must pass close enough to Saturn to be tidally disrupted; as a result, they would suffer aerodynamic disruption. Crater counting on adjacent satellites gives densities of about 200 per 10 (exp 6) square km for 20-km-diameter craters. However, the presence of a thick atmosphere leads to atmospheric shielding, depleting the relative abundance of small craters. This has been evaluated by models, and the relative abundance of small craters may be due to a diagnostic atmospheric collapse. A number of radar-dark "splotches" have been detected on Venus; these have been attributed to the interaction of the surface with the atmospheric shockwave produced by the Tunguska-like explosion of a bolide in the atmosphere. Simple analogy suggests that similar features might occur on Titan, but the shocked mass density (which controls the momentum coupling between the surface and the shockwave) of Titan's cold N2 atmosphere is about 20x smaller than that of Venus's hot CO2 atmosphere. Unless ice is much more easily turned to rubble than is rock, such features seem less probable on Titan. When the energy deposited by an impact forms a fireball with an equilibrate greater than one scale height, the fireball expands upward and can distribute ejecta. on ballistic exoatmospheric trajectories. On Venus this process is believed to be responsible for the parabolic features; the interaction of various-sized particles falling through the atmosphere with the zonal wind field winnows the particles to form a parabolic deposit. Although such a process is possible on Titan, the large scale height at higher altitudes would make it more difficult. Comparison with craters on other icy satellites suggests that craters on Titan will be fairly shallow (depth/diameter about 0.1) and craters greater than 10 km in diameter will have central peaks or domed bases, perhaps with central pits. The formation of ejecta. blankets may involve the atmosphere in a significant way, both by restraining the expansion of the ejecta cloud and by influencing the thermal history of the ejecta. Compared with Venus, Titan's atmosphere will chill an impact melt somewhat quickly, so the long ejecta flows seen on Venus seem less likely, detailed modeling needs to be performed to determine the impact melt production. Crater topography on Titan may be highlighted by the influence of liquids forming crater lakes. Craters with central peaks will typically form ring-shaped lakes, although horseshoe-shaped takes may be common; domed craters with central pits may even form bullseye lakes with islands with central ponds. If liquids have covered a substantial part of Titan's surface for a substantial period, hydroblemes and tsunami deposits may be common.

Lorenz, R. D.↗

Orbit tomography of energetic particle distribution functions

Both fast ions and runaway electrons are described by distribution functions, the understanding of which are of critical importance for the success of future fusion devices such as ITER. Typically, energetic particle diagnostics are only sensitive to a limited subsection of the energetic particle phase-space which is often insufficient for model validation. However, previous publications show that multiple measurements of a single spatially localized volume can be used to reconstruct a distribution function of the energetic particle velocity-space by using the diagnostics' velocity-space weight functions, i.e. velocity-space Tomography. In this work we use the recently formulated orbit weight functions to remove the restriction of spatially localized measurements and present orbit tomography, which is used to reconstruct the 3D phase-space distribution of all energetic particle orbits in the plasma. Through a transformation of the orbit distribution, the full energetic particle distribution function can be determined in the standard {energy, pitch, r, z}-space. We benchmark the technique by reconstructing the fast-ion distribution function of an MHD-quiescent DIII-D discharge using synthetic and experimental FIDA measurements. We also use the method to study the redistribution of fast ions during a sawtooth crash at ASDEX upgrade using FIDA measurements. Lastly, a comparison between the orbit tomography and velocity-space tomography is shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing the laser wakefield in underdense plasmas by induced terahertz emission

Terahertz (THz) radiation can be produced from a laser wakefield driven in underdense plasmas in the presence of a transverse DC magnetic field. It is shown that the radiation usually contains a component at the electron plasma frequency and its harmonics when the wakefield is excited at high amplitudes. In the highly nonlinear bubble/blowout regime, the radiation contains a smooth component peaked at the reduced electron plasma frequency and an irregular spectrum extending to tens of the electron plasma frequency. The latter is due to the broken-wave structure behind the bubble. A theoretical model is presented and validated via two-dimensional particle-in-cell simulations. The measurement of such THz emission may provide a diagnostic of the laser wakefield structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Friction Stir Welding at MSFC: Kinematics

In 1991 The Welding Institute of the United Kingdom patented the Friction Stir Welding (FSW) process. In FSW a rotating pin-tool is inserted into a weld seam and literally stirs the faying surfaces together as it moves up the seam. By April 2000 the American Welding Society International Welding and Fabricating Exposition featured several exhibits of commercial FSW processes and the 81st Annual Convention devoted a technical session to the process. The FSW process is of interest to Marshall Space Flight Center (MSFC) as a means of avoiding hot-cracking problems presented by the 2195 aluminum-lithium alloy, which is the primary constituent of the Lightweight Space Shuttle External Tank. The process has been under development at MSFC for External Tank applications since the early 1990's. Early development of the FSW process proceeded by cut-and-try empirical methods. A substantial and complex body of data resulted. A theoretical model was wanted to deal with the complexity and reduce the data to concepts serviceable for process diagnostics, optimization, parameter selection, etc. A first step in understanding the FSW process is to determine the kinematics, i.e., the flow field in the metal in the vicinity of the pin-tool. Given the kinematics, the dynamics, i.e., the forces, can be targeted. Given a completed model of the FSW process, attempts at rational design of tools and selection of process parameters can be made.

Nunes, A. C., Jr.↗

Synergy of coupled kink and tearing modes in fast ion transport

One aspect of the interaction between fast ions and magnetohydrodynamic (MHD) instabilities is the fast ion transport. Coupled kink and tearing MHD instabilities have also been reported to cause fast ion transport. Recently, the 'kick' model has been developed to compute the evolution of the fast ion distribution from neutral beam injection using instabilities as phase-space resonance sources. Overall, the goal of this paper is to utilize the kick model to understand the physics of fast ion transport caused by the coupled kink and tearing modes. Soft x-ray diagnostics are used to identify the mode parameters in the National Spherical Torus Experiment. The comparison of neutron rates measured and computed from time-dependent TRANSP simulation with the kick model shows that the coupling of kink and tearing mode is important in determination of the fast ion transport. The numerical scan of the mode parameters shows that the relative phase of the kink and tearing modes and the overlapping of kink and tearing mode resonances in the phase space can affect the fast ion transport, suggesting that the synergy of the coupled modes may be causing the fast ion transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mode and Intermediate Waters in Earth System Models

This report describes work done as part of a joint Princeton-Johns Hopkins project to look at the impact of mode and intermediate waters in Earth System Models. The Johns Hopkins portion of this work focussed on the role of lateral mixing in ventilating such waters, with important implications for hypoxia, the uptake of anthropogenic carbon, the dynamics of El Nino and carbon pumps. The Johns Hopkins group also collaborated with the Princeton Group to help develop a watermass diagnostics framework.

54 ENVIRONMENTAL SCIENCES↗

Hybrid time dependent/independent solution for the He I line ratio temperature and density diagnostic for a thermal helium beam with applications in the scrape-off layer-edge regions in tokamaks

Spectroscopic studies of line emission intensities and ratios offer an attractive option in the development of non-invasive plasma diagnostics. Evaluating ratios of selected He I line emission profiles from the singlet and triplet neutral helium spin systems allows for simultaneous measurement of electron density (n{sub e}) and temperature (T{sub e}) profiles. Typically, this powerful diagnostic tool is limited by the relatively long relaxation times of the {sup 3}S metastable term of helium that populates the triplet spin system, and on which electron temperature sensitive lines are based. By developing a time dependent analytical solution, we model the time evolution of the two spin systems. We present a hybrid time dependent/independent line ratio solution that improves the range of application of this diagnostic technique in the scrape-off layer (SOL) and edge plasma regions when comparing it against the current equilibrium line ratio helium model used at TEXTOR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accelerators for critical experiments involving single-particle upset in solid-state microcircuits

Charged-particle interactions in microelectronic circuit chips (integrated circuits) present a particularly insidious problem for solid-state electronic systems due to the generation of soft errors or single-particle event upset (SEU) by either cosmic rays or other radiation sources. Particle accelerators are used to provide both light and heavy ions in order to assess the propensity of integrated circuit chips for SEU. Critical aspects of this assessment involve the ability to analytically model SEU for the prediction of error rates in known radiation environments. In order to accurately model SEU, the measurement and prediction of energy deposition in the form of an electron-hole plasma generated along an ion track is of paramount importance. This requires the use of accelerators which allow for ease in both energy control (change of energy) and change of ion species. This and other aspects of ion-beam control and diagnostics (e.g., uniformity and flux) are of critical concern for the experimental verification of theoretical SEU models.

Zoutendyk, J. A.↗

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗