Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear dynamics in fluids”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

[?]Nonlinear Issues in the Aerothermochemistry of Gases and Materials and the Associated Physics and Dynamics of Interfaces

Our research and technology are focused on nonlinear issues in the aerothermochemistry of gases and materials and the associated physics and dynamics of interfaces. Our program is now organized to aggressively support the NASA Aeronautics Enterprise so as to: (a) develop a new generation of environmentally compatible, economic subsonic aircraft; (b) develop the technology base for an economically viable and environmentally compatible high-speed civil transport; (c) develop the technology options for new capabilities in high-performance aircraft; (d) develop hypersonic technologies for air-breathing flight; and (e) develop advanced concepts, understanding of physical phenomena, and theoretical, experimental, and computational tools for advanced aerospace systems. The implications from our research for aeronautical and aerospace technology have been both broad and deep. For example, using advanced computational techniques, we have determined exact solutions for the Schrodinger equation in electron-molecule scattering allowing us to evaluate atmospheric models important to reentry physics. We have also found a new class of exact solutions for the Navier Stokes equations. In experimental fluid dynamics, we have found explicit evidence of turbulence modification of droplet sizes in shock tube flow with condensation. We have developed a new diagnostic tool for the direct estimation of flow velocities at MHz sampling rates in quasi-one dimensional turbulent flow. This procedure suggests an unexpected confirmation of the possibility of 'natural' closure in Reynolds stresses with deep implications for the development of turbulent models. A transient increase is observed in both the spectral energy decay rate and the degree of chaotic complexity at the interface of a shock wave and a turbulent ionized gas. Even though the gas is apparently brought to rest by the shock wave, no evidence is found either of the expected relaminarization. A unique diamond-shaped nozzle has been designed for a detailed investigation of the effect of significant streamwise vorticity on the acoustic and IR characteristics of supersonic jets. Our results provide convincing evidence of the significant effect of vorticity on the far-field noise for the diamond jet as compared to the conventional round jets. We have found that the countercurrent shear layer mixes much more efficiently than conventional coflowing shear layers. We also developed the fluid thrust vectoring procedures which use counter flow to vector a jet. Our materials research has shown that the steep stress gradients at the fiber-matrix interface could be the primary cause of interface cracks after the processing of metallic and intermetallic matrix composites. New techniques have been evolved for: the microcharacterization of materials including microplastic strain and, point by point, the misorientation and plasticity for matrix composites; thermally induced stress measurements and load relaxation; the growth and characterization of metallic matrix composite interfaces; and for the growth of ferrite materials by pulsed laser deposition. The FAMU commitment to the HBCU Research Center also continues to be broad and deep.

Johnson, Joseph A., III↗
NASA NTRS1996

Nonlinear Bubble Interactions in Acoustic Pressure Fields

The systems consisting of a two-phase mixture, as clouds of bubbles or drops, have shown many common features in their responses to different external force fields. One of particular interest is the effect of an unsteady pressure field applied to these systems, case in which the coupling of the vibrations induced in two neighboring components (two drops or two bubbles) may result in an interaction force between them. This behavior was explained by Bjerknes by postulating that every body that is moving in an accelerating fluid is subjected to a 'kinetic buoyancy' equal with the product of the acceleration of the fluid multiplied by the mass of the fluid displaced by the body. The external sound wave applied to a system of drops/bubbles triggers secondary sound waves from each component of the system. These secondary pressure fields integrated over the surface of the neighboring drop/bubble may result in a force additional to the effect of the primary sound wave on each component of the system. In certain conditions, the magnitude of these secondary forces may result in significant changes in the dynamics of each component, thus in the behavior of the entire system. In a system containing bubbles, the sound wave radiated by one bubble at the location of a neighboring one is dominated by the volume oscillation mode and its effects can be important for a large range of frequencies. The interaction forces in a system consisting of drops are much smaller than those consisting of bubbles. Therefore, as a first step towards the understanding of the drop-drop interaction subject to external pressure fluctuations, it is more convenient to study the bubble interactions. This paper presents experimental results and theoretical predictions concerning the interaction and the motion of two levitated air bubbles in water in the presence of an acoustic field at high frequencies (22-23 KHz).

Barbat, Tiberiu↗
NASA NTRS1999

Simulation of Non-Linear Transonic Aeroelastic Behavior on the B-2

At high subsonic flight speeds, large flexible aircraft begin to encounter unsteady airloads which are not predicted by most currently available aerodynamic analysis and design methods. With increasing speed and the development of transonic flow and shocks, viscous effects quickly become very important, and flow separation can occur. The Northrop Grumman USAF B-2 Bomber encountered a nonlinear aeroelastic Residual Pitch Oscillation (RPO) under these conditions. Simulation studies were performed with the Computational Aeroelasticity Program-Transonic Small Disturbance, Viscous (CAPTSDv) computer program to evaluate its ability to predict these nonlinear aeroelastic responses. Open and closed loop simulations were performed to assess the participation of the flight control system. Control, system actuator hysteresis characteristics were modeled and found to be a significant participant in the RPO phenomenon. Simulations were also performed for varying Mach numbers and altitudes to establish the stability boundaries and compare with flight test data. These CAPTSDv simulations compared well with flight data and revealed many potential further modeling enhancements.

Residual pitch oscillation↗
NASA NTRS2000

Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-capturing Methods

In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Yee, H. C.↗
NASA NTRS1986

Multiple grid problems on concurrent-processing computers

Three computer codes were studied which make use of concurrent processing computer architectures in computational fluid dynamics (CFD). The three parallel codes were tested on a two processor multiple-instruction/multiple-data (MIMD) facility at NASA Ames Research Center, and are suggested for efficient parallel computations. The first code is a well-known program which makes use of the Beam and Warming, implicit, approximate factored algorithm. This study demonstrates the parallelism found in a well-known scheme and it achieved speedups exceeding 1.9 on the two processor MIMD test facility. The second code studied made use of an embedded grid scheme which is used to solve problems having complex geometries. The particular application for this study considered an airfoil/flap geometry in an incompressible flow. The scheme eliminates some of the inherent difficulties found in adapting approximate factorization techniques onto MIMD machines and allows the use of chaotic relaxation and asynchronous iteration techniques. The third code studied is an application of overset grids to a supersonic blunt body problem. The code addresses the difficulties encountered when using embedded grids on a compressible, and therefore nonlinear, problem. The complex numerical boundary system associated with overset grids is discussed and several boundary schemes are suggested. A boundary scheme based on the method of characteristics achieved the best results.

Eberhardt, D. S.↗
NASA NTRS2023

Breaking Wave Field Statistics With A Multi-Layer Model

The statistics of breaking wave fields are characterised within a novel multi-layer framework, which generalises the single-layer Saint-Venant system into a multi-layer and non-hydrostatic formulation of the Navier–Stokes equations. We simulate an ensemble of phase-resolved surface wave fields in physical space, where strong nonlinearities, including directional wave breaking and the subsequent highly rotational flow motion, are modelled, without surface overturning. We extract the kinematics of wave breaking by identifying breaking fronts and their speed, for freely evolving wave fields initialised with typical wind wave spectra. The Λ(c) distribution, defined as the length of breaking fronts (per unit area) moving with speed c to c + dc following Phillips (J. Fluid Mech., vol. 156, 1985, pp. 505–531), is reported for a broad range of conditions. We recover the Λ(c) ∝ c−6 scaling without wind forcing for sufficiently steep wave fields. A scaling of Λ(c) based solely on the root-mean-square slope and peak wave phase speed is shown to describe the modelled breaking distributions well. The modelled breaking distributions are in good agreement with field measurements and the proposed scaling can be applied successfully to the observational data sets. The present work paves the way for simulations of the turbulent upper ocean directly coupled to a realistic breaking wave dynamics, including Langmuir turbulence, and other sub-mesoscale processes.

Jiarong Wu↗
NASA NTRS2003

Assessing blood brain barrier dynamics or identifying or measuring selected substances or toxins in a subject by analyzing Raman spectrum signals of selected regions in the eye

A non-invasive method for analyzing the blood-brain barrier includes obtaining a Raman spectrum of a selected portion of the eye and monitoring the Raman spectrum to ascertain a change to the dynamics of the blood brain barrier. Also, non-invasive methods for determining the brain or blood level of an analyte of interest, such as glucose, drugs, alcohol, poisons, and the like, comprises: generating an excitation laser beam (e.g., at a wavelength of 600 to 900 nanometers); focusing the excitation laser beam into the anterior chamber of an eye of the subject so that aqueous humor, vitreous humor, or one or more conjunctiva vessels in the eye is illuminated; detecting (preferably confocally detecting) a Raman spectrum from the illuminated portion of the eye; and then determining the blood level or brain level (intracranial or cerebral spinal fluid level) of an analyte of interest for the subject from the Raman spectrum. In certain embodiments, the detecting step may be followed by the step of subtracting a confounding fluorescence spectrum from the Raman spectrum to produce a difference spectrum; and determining the blood level and/or brain level of the analyte of interest for the subject from that difference spectrum, preferably using linear or nonlinear multivariate analysis such as partial least squares analysis. Apparatus for carrying out the foregoing methods are also disclosed.

Lambert, James L.↗
DOE OSTI2021

Effects of Strain Rate and Temperature on the Mechanical Properties of Simulated Silica Ionogels

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Additionally, owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2008

Assessing blood brain barrier dynamics or identifying or measuring selected substances, including ethanol or toxins, in a subject by analyzing Raman spectrum signals

A non-invasive method for analyzing the blood-brain barrier includes obtaining a Raman spectrum of a selected portion of the eye and monitoring the Raman spectrum to ascertain a change to the dynamics of the blood brain barrier.Also, non-invasive methods for determining the brain or blood level of an analyte of interest, such as glucose, drugs, alcohol, poisons, and the like, comprises: generating an excitation laser beam at a selected wavelength (e.g., at a wavelength of about 400 to 900 nanometers); focusing the excitation laser beam into the anterior chamber of an eye of the subject so that aqueous humor, vitreous humor, or one or more conjunctiva vessels in the eye is illuminated; detecting (preferably confocally detecting) a Raman spectrum from the illuminated portion of the eye; and then determining the blood level or brain level (intracranial or cerebral spinal fluid level) of an analyte of interest for the subject from the Raman spectrum. In certain embodiments, the detecting step may be followed by the step of subtracting a confounding fluorescence spectrum from the Raman spectrum to produce a difference spectrum; and determining the blood level and/or brain level of the analyte of interest for the subject from that difference spectrum, preferably using linear or nonlinear multivariate analysis such as partial least squares analysis. Apparatus for carrying out the foregoing methods are also disclosed.

Lambert, James L.↗
DOE OSTI2024

Microphysics of shock-grain interaction for inertial confinement fusion ablators in a fluid approach

Ablator materials used for inertial confinement fusion, such as high-density carbon (HDC) and beryllium, have grain structure which may lead to small-scale density nonuniformity and the generation of perturbations when the materials are shocked and compressed. Here, we use a combination of a linear theory of shock interaction with density nonuniformity [Velikovich et al., Phys. Plasmas 14, 072706 (2007)] and numerical simulations to study shock interaction with a model representation of HDC grains. While the shock-grain interaction is nonlinear, the linear theory shows some key features of the shock-grain interaction, which also hold for the (nonlinear) simulations. The postshock perturbations are made up of sonic reflections off of grain boundaries and vorticity deposition along them, with the latter dominating the perturbed energy content. The mean (per mass) postshock perturbed kinetic energy decreases with increasing grain size, but energy will be deposited at increasing spatial scale. From the perspective of the postshock perturbed energy, the detailed linear theory largely supports a proposed method [S. Davidovits et al., Phys. Plasmas 29, 112708 (2022)] for deresolving the grains (in a similar grains model) that treats the grains statistically. Finally, our simulation results highlight the influence of thermal conduction on the perturbation dynamics at grain scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Halo formation from Yukawa forces in the very early Universe

If long-range attractive forces exist and are stronger than gravity then cosmic halo formation can begin in the radiation-dominated era. We study a simple realization of this effect in a system where dark matter fermions have Yukawa interactions mediated by scalar particles, analogous to the Higgs boson in the standard model. We develop a self-consistent description of the system including exact background dynamics of the scalar field, and precise modeling of the fermion density fluctuations. For the latter, we provide accurate approximations for the linear growth as well as quantitative modeling of the nonlinear evolution using N-body simulations. We find that halo formation occurs exponentially fast and on scales substantially larger than simple estimates predict. Here, the final fate of these halos remains uncertain, but could be annihilation, dark stars, primordial black holes, or even the existence of galaxy-sized halos at matter-radiation equality. More generally, our results demonstrate the importance of mapping scalar-mediated interactions onto structure formation outcomes and constraints for beyond the standard model theories.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Using adaptive grid in modeling rocket nozzle flow

The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which cannot be solved analytically. However, this system of equations called the Navier-Stokes equations can be solved numerically. The accuracy and the convergence of the solution of the system of equations will depend largely on how precisely the sharp gradients in the domain of interest can be resolved. With the advances in computer technology, more sophisticated algorithms are available to improve the accuracy and convergence of the solutions. An adaptive grid generation is one of the schemes which can be incorporated into the algorithm to enhance the capability of numerical modeling. It is equivalent to putting intelligence into the algorithm to optimize the use of computer memory. With this scheme, the finite difference domain of the flow field called the grid does neither have to be very fine nor strategically placed at the location of sharp gradients. The grid is self adapting as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzles by taking the refinement part of grid generation out of the hands of computational fluid dynamics (CFD) specialists and place it into the computer algorithm itself.

Chow, Alan S.↗
NASA NTRS1995

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗
DOE OSTI2025

Enhanced shear stabilization of turbulence in NSTX

In studying a particular non-stationary NSTX L-mode plasma, we observed unexpectedly high levels of flux—first with quasilinear (TGLF) modeling and subsequently with nonlinear gyrokinetic simulations. Upon more detailed analysis, a novel confinement regime was discovered in which a modest increase in E x B shear (beyond baseline experimental estimates) rapidly reduced turbulent transport to levels consistent with power balance. This modest increase is plausible given the errors inherent to the estimation of shearing rates, and the added complexity of the non-stationary (time-dependent) power balance. Remarkably, an additional small increase in shear yields the familiar ion-neoclassical transport level with what appears to be the onset of high-k electron transport only. Although analyses using the TGLF-SAT2 model successfully capture numerous parametric dependencies of this plasma, TGLF does not reproduce the rapid E x B stabilization seen in CGYRO. We believe the results presented should help to better characterize the nonlinear physics of spherical tokamak confinement regimes, provide useful ST datasets for reduced model development, and motivate more accurate experimental diagnosis of E x B shearing rates.

Atomic and molecular collisions↗
DOE OSTI2020

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2025

Assimilating partial observation to enhance feedback control of stochastic dynamical systems

Here, in this paper, we present a novel methodology to tackle feedback optimal control problems in scenarios where the exact state of the controlled process is unknown. It integrates data assimilation techniques and optimal control solvers to manage partial observation of the state process, a common occurrence in practical scenarios. Traditional stochastic optimal control methods assume full state observation, which is often not feasible in real-world fluid dynamics control problems. Our approach underscores the significance of utilizing observational data to inform control policy design. Specifically, we introduce a kernel learning backward stochastic differential equation (SDE) filter to enhance data assimilation efficiency and propose a sample-wise stochastic optimization method within the stochastic maximum principle framework. We demonstrate the efficacy and accuracy of our method in the control of advection-diffusion-reaction flow problem and the Dubins airplane maneuvering problem with model uncertainty.

data driven↗
DOE OSTI2026

Effects of external magnetic fields on the multi-mode Rayleigh–Taylor instability

The magneto-Rayleigh–Taylor instability (mRTI) is a key process in inertial confinement fusion and is thought to be widespread in the interstellar medium, where it can concentrate plasma into discrete structures. We present resistive MHD simulations of the nonlinear evolution of multi-mode mRTI in both two and three dimensions, examining the effects of uniform external magnetic fields oriented either parallel or perpendicular to the initial interface. In both 2D and 3D, weak parallel fields enhance mixing-zone growth, whereas stronger fields suppress it. For perpendicular fields, growth is initially inhibited but becomes enhanced at later times. These behaviors arise from magnetic tension, which modifies flow anisotropy, buoyancy, drag, and vortex dynamics. The interplay of these mechanisms governs the distinct ways in which magnetic fields influence mRTI evolution.

Flow instabilities↗
NASA NTRSDate not supplied

A Moving Discontinuous Galerkin Method with Interface Condition Enforcement for Reacting Hypersonic Flows

The necessity of enforcing conservation in computational elements or cells (element conservation) for discontinuous solutions is well understood and respected for solving conservation laws in computational fluid dynamics (CFD). In contrast, interface conservation, where the conservation across cell interfaces is enforced, is long ignored, and yet is also ruled and required by the underlying physics just like element conservation. Violation of the interface conservation across discontinuities is the root cause why an exact discontinuous solution can never be achieved in shock capturing methods. The interface conservation is examined and explored in this talk. Moving discontinuous Galerkin (MDG) finite element method with interface conservation enforcement (MDG-ICE) [1],[2] are then presented for solving compressible flow problems with discontinuities based on the observation that the interface conservation can only be satisfied, only when mesh interfaces are aligned with discontinuities. In the MDG-ICE formulation, both conservative quantities and grid geometry are considered as independent variables. A space-time DG formulation is used to solve the multi-material compressible Euler equations in the standard discontinuous solution space and the discrete grid geometry is solved using a variational formulation in a continuous space. A self-adaptive Levenberg-Marquardt method is utilized to solve the resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation. A number of numerical experiments for a variety of flow problems are conducted to assess the accuracy and performance of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to deliver the designed order of both h- and p-convergence even for discontinuous solutions, and detect all types of interfaces, via interface condition enforcement and satisfy, via grid movement, the compressible Euler equations and the associated interface condition.

Discontinuous Galekin↗