Sharp Shock Research Articles

An improved efficient adaptive method for large-scale multi-explosives explosion simulations

Shock wave caused by a sudden release of high-energy, such as explosion and blast, usually affects a significant range of areas. The utilization of a uniform fine mesh to capture sharp shock wave and to obtain precise results is inefficient in terms of computational resource. This is particularly evident when large-scale fluid field simulations are conducted with significant differences in computational domain size. In this work, a variable-domain-size adaptive mesh enlargement (vAME) method is developed based on the proposed adaptive mesh enlargement (AME) method for modeling multi-explosives explosion problems. The vAME method reduces the division of numerous empty areas or unnecessary computational domains by adaptively suspending enlargement operation in one or two directions, rather than in all directions as in AME method. A series of numerical tests via AME and vAME with varying nonintegral enlargement ratios and different mesh numbers are simulated to verify the efficiency and order of accuracy. An estimate of speedup ratio is analyzed for further efficiency comparison. Several large-scale near-ground explosion experiments with single/multiple explosives are performed to analyze the shock wave superposition formed by the incident wave, reflected wave, and Mach wave. Additionally, the vAME method is employed to validate the accuracy, as well as to investigate the performance of the fluid field and shock wave propagation, considering explosive quantities ranging from 1 to 5 while maintaining a constant total mass. The results show a satisfactory correlation between the overpressure versus time curves for experiments and numerical simulations. The vAME method yields a competitive efficiency, increasing the computational speed to 3.0 and approximately 120,000 times in comparison to AME and the fully fine mesh method, respectively. It indicates that the vAME method reduces the computational cost with minimal impact on the results for such large-scale high-energy release problems with significant differences in computational domain size.

Decomposition of oil price shock, changes in market conditions and systemic risk transmission – an analysis from the Chinese perspective

ABSTRACT This research aims to investigate the heterogeneous responses of each specific risk variable across different states in the context of varying levels of crude oil price shocks. Using data period from 27 December 2007 to 28 October 2022, we construct proxy variables for systemic financial risks based on the considerations of four specific perspectives. The results of the empirical analysis indicate that oil price fluctuation observed in short, medium, and long-term time periods has an obvious impact on each risk variable. From the evaluation of conditional value-at-risk difference as well as the marginal expected shortfall, both short-term and long-term shocks have an apparent amplifying effect on market risk when it is at a high level. When evaluated from the angle of liquidity, oil price shocks generally exhibit positive effects on market liquidity, where this effect is more significant in the long term rather than terms. When considering the effect of volatility specifically, oil price shocks are found to cause sharp shocks to the market, whether the market itself is in a stable or volatile state.

Are shocks to entrepreneurship persistence? Case of a Resource-based economy

PurposeThis paper aims to investigate the stationary properties of entrepreneurship capital in the Iranian business sector. The investigation is conducted based on firm size (micro, small and medium, and large), sector (ISIC classification), and location (23 provinces) over the period 1981–2021.Design/methodology/approachTo achieve the purpose of the paper, we apply the Lagrange Multiplier (LM) unit root test with structural breaks.FindingsThe results of our study confirm the stationarity of entrepreneurship for the majority of our categorizations. Therefore, entrepreneurs' willingness to start a business is temporarily affected by sharp shocks, suggesting that the entrepreneurship trends in these categories will eventually revert to their long-run equilibrium. However, the time series of entrepreneurship in provinces such as Zanjan and Lorestan, as well as in the electrical machines and devices (code 31) sector, remain permanent. These findings can assist policymakers in each sector and location in designing effective policies to promote entrepreneurial activities.Research limitations/implicationsIn this study, the time period utilized was restricted, and there were no data accessible for an extended duration. Another limitation of this research is the absence of access to firm-level data on a shorter time scale, such as weekly or monthly.Originality/valueThe economic literature reveals that empirical studies on the persistence of entrepreneurship have received relatively less attention in the context of emerging and resource-based economies, compared to the increasing focus on them in developed countries. Therefore, to address this gap, this paper aims to extend the current empirical literature by presenting new evidence for the case of Iran, which has an emerging and resource-based economy.

Mental Health Impacts of COVID-19: Does Prepandemic Cognition and Dementia Status Matter?

The coronavirus disease 2019 (COVID-19) pandemic disrupted daily life and led to sharp shocks in trends for various health outcomes. Although substantial evidence exists linking the pandemic and mental health outcomes and linking dementia and mental health outcomes, little evidence exists on how cognitive status may alter the impact of COVID-19 on mental health. We used prepandemic data from the Longitudinal Aging Study in India-Diagnostic Assessment of Dementia study and 9 waves of data from the Real-Time Insights of COVID-19 in India study (N = 1182). We estimated associations between measures of prepandemic cognition (continuous cognition based on 22 cognitive tests, dementia status) and mental health measures during the pandemic (Patient Health Questionnaire [PHQ]-4 [9 time points], PHQ-9 [2 time points], Beck Anxiety Inventory [3 time points]), adjusting for age, gender, rural/urban residence, state, education, and prepandemic mental health. Summarizing across time points, PHQ-9 score was marginally or significantly associated with prepandemic cognition (PHQ-9 difference: -0.38 [-0.78 to 0.14] points per SD higher cognition; p = .06), and prepandemic dementia (PHQ-9 difference: 0.61 [0.11-1.13] points for those with dementia compared to no dementia; p = .02). Associations with BAI were null, whereas associations with PHQ-4 varied over time (p value for interaction = .02) and were strongest during the delta wave, when pandemic burden was highest. We present initial evidence that mental health impacts of COVID-19 or other acute stressors may be unequally distributed across strata of cognitive outcomes. In dynamically changing environments, those with cognitive impairment or dementia may be more vulnerable to adverse mental health outcomes.

The Coupling of Innovation in Traditional Chinese Cultural Diffusion and Smart Cultural Aging

Abstract Based on the coupled communication of smart elderly care and traditional Chinese culture, this paper proposes four paths of coupled communication, including in-depth communication, innovative cultural forms of communication, precise communication, and promotion of coupled development. Two systematic index systems have been constructed to disseminate traditional Chinese culture and develop smart culture and pension. The data is standardized using principal component analysis to establish the correlation coefficient matrix. The coordination index between the two subsystems is calculated, and the coordination model of the two is created based on the concepts of coordination coefficient and coefficient of variation. The subject area, Province S, is in the third stage of the coupled coordinated development during 2010-2023, with traditional culture dissemination and smart elderly care services coming to a high level, and is in a good development level during 2016-2023, with a coupled coordination degree of 8 or more. Using the VAR model, the effect of coupling development on the quality of life of elderly individuals is analyzed. Impulse response analysis of LNCCD and LNQOE: after LNCCD receives a positive shock, LNQOE shows a sharp shock effect. At this time, the response value is 0.07021, which plays a large positive effect, and then enters into a smooth development. The traditional Chinese culture dissemination and the wisdom culture of old age are in the stage of coupled development.

Effect of symmetrical and asymmetrical tilt grain boundaries with as the tilt axis on the shock response of bi-crystalline Ni

Aim of this article is to study the effect of tilt grain boundaries on the shock response of single and bi-crystalline nickel. Nickel is commonly used in many structural applications, but research on dynamic response is very limited or in an immature stage. In this article, molecular dynamics-based simulations were performed in conjunction with embedded atom method potential to quantify the response of shock in bi-crystalline Ni containing symmetrical or asymmetrical tilt grain boundaries. Simulations were performed under the influence of adiabatic conditions. Compressive ultra-shock pulse of width 2 ps was simulated in the simulation box; stress, particle velocity, dislocation emission and phase transformation were quantified for symmetrical and asymmetrical tilt grain boundaries in bi-crystalline Ni. It was concluded from the simulations that higher mis-orientation angle between the crystals of Ni help in retarding the shock front velocity and stress. The alternation in the shock front behavior was associated with the change in deformation governing mechanism. The emission of dislocations and formation of Lomer Cottrell lock governs deformation in symmetrical tilt grain boundary configurations under shock pulse loading. It was concluded from the set of simulations that plastic deformation blunts the shock front in bi-crystalline Ni, whereas the sharp shock front was observed in single crystal Ni. Radial distribution function in post shock atomistic configurations was performed to capture the stability of the crystal in post shock scenario. Asymmetrical tilt grain boundaries are not as efficient in restricting the shock front as symmetrical tilt grain boundary configurations.

A learned conservative semi-Lagrangian finite volume scheme for transport simulations

Semi-Lagrangian (SL) schemes are known as a major numerical tool for solving transport equations with many advantages and have been widely deployed in the fields of computational fluid dynamics, plasma physics modeling, numerical weather prediction, among others. In this work, we develop a novel machine learning-assisted approach to accelerate the conventional SL finite volume (FV) schemes. The proposed scheme avoids the expensive tracking of upstream cells but attempts to learn the SL discretization from the data by incorporating specific inductive biases in the neural network, significantly simplifying the algorithm implementation and leading to improved efficiency. In addition, the method delivers sharp shock transitions and a level of accuracy that would typically require a much finer grid with traditional transport solvers. Numerical tests demonstrate the effectiveness and efficiency of the proposed method.

BI06 Short sharp shock or slow burn: clinician perspectives on treatment of actinic keratoses with 5-fluorouracil 5% cream monotherapy vs. 5-fluorouracil/calcipotriol ointment combination therapy

Abstract Actinic keratoses (AK) are common in immunocompromised patients and associated with increased risk of cutaneous squamous cell carcinoma (cSCC). The current standard of care is 5-fluorouracil cream 5% (5-FU) monotherapy for 4 weeks. In recent years, 5-FU cream and calcipotriol ointment (5-FU-Cal) combination AK therapy for 4–6 days has emerged as potentially more effective and better tolerated, although the original randomized controlled trial (RCT) excluded immunosuppressed patients. Previously, we shared immunosuppressed patients’ real-world perspectives of using 5-FU monotherapy vs. 5-FU-Cal combination therapy (Junejo MH, Demirel S, Matin R et al. Short shock or slow burn? Patient perspectives on treatment of actinic keratoses with 5-fluorouracil 5% cream monotherapy versus 5-fluorouracil plus calcipotriol combination therapy. Clin Exp Dermatol 2023; 48:31–2). The majority (81%) found 5-FU-Cal to be more effective than 5-FU monotherapy. Despite 47% reporting worse local skin reactions (LSRs), most preferred the 5-FU-Cal combination. We now report results from a clinician survey on the use of 5-FU monotherapy vs. 5-FU-Cal combination therapy for the treatment of AK. An online survey was sent to 200 UK consultant dermatologists. Fifty-seven (28.5%) completed the survey, all of whom reported using 5-FU monotherapy for treating AK in patients with prior cSCC; 43 (75%) reported this to be their preferred treatment. Most clinicians (n = 40; 70%) were aware of 5-FU-Cal combination therapy, although only 24 (42%) reported using it to treat AK in patients with prior cSCC. Of these, 18 (75%) used combination treatment for AK as part of field change, 14 (58%) used it for the treatment of patches of more than five AK and 13 (54%) used it for grouped AK. Very few (n = 4; 17%) used combination therapy for treatment of cSCC in situ (n = 4), single AK (n = 3) or superficial basal cell carcinoma (n = 1). Of 24 clinicians using 5-FU-Cal, the majority felt that patients developed LSRs often (n = 11; 46%) or always (n = 7; 29%); fewer felt that LSRs occurred sometimes (n = 5; 21%) or rarely (n = 1; 4%). Clinicians using combination therapy primarily felt the LSRs were moderate but not severe enough to require supplementary treatment or discontinuation (n = 17; 71%). Half (n = 12; 50%) of those using combination therapy felt patients did not signal a preference between 5-FU monotherapy vs. 5-FU-Cal combination therapy. Two-thirds (n = 38/57; 67%) of clinicians reported an interest in participating in a RCT comparing 5-FU vs. 5-FU-Cal vs. sunscreen alone for the treatment of AK in patients with prior cSCC. This is the first evaluation of clinician perspectives on the use of the newly emerging 5-FU-Cal therapy for AK. These data provide a rationale for the inclusion of combination 5-FU-Cal therapy in future studies using topical AK treatments for cSCC prevention in the UK.

Short sharp shock and long deep rumble: Natural disasters and private property loss in New Zealand

This article contributes to the nascent research agenda on accounting for natural disasters by addressing the question: What responsibility did New Zealand's central government assume in respect of private property loss caused by natural disasters, and did this change over time? It considers the source of financial support to assist New Zealanders who suffered a loss due to a natural disaster in the period from early settlement until 1990. Although change occurred, throughout the period the state did not assume responsibility for private property loss. This article adopts historical institutionalism to contribute to the research agenda of accounting for natural disasters by considering policy and response over the long term, rather than focusing on a single disaster event or specific disaster type.

Transport-Equation-Based Shock Sensor for Shock–Turbulent Boundary-Layer Interaction

High-speed flows are often associated with shock–boundary-layer interaction (SBLI), which causes excessive vehicle surface heating and total pressure loss. Therefore, precise computational prediction of shock interactions becomes essential. Commonly used simulation techniques, including numerical schemes and turbulence models, give erroneous results at shock discontinuities. Local modifications are applied at shocks based on their location and strength. Most shock sensors in the literature predict shock location but do not give shock properties. We propose a shock function to locate shocks and estimate their strength in terms of density ratio by solving a transport equation, along with the governing equations of fluid flow. The sensitivity of the shock function to the transport equation’s parameters and grid resolution are explored. The shock function is used to define a shock sensor that can identify thin and sharp shocks. The shock function and shock sensor are then applied to an oblique shock impinging on a turbulent boundary layer. We use the shock unsteadiness modified model, which uses the estimated shock strength to predict the turbulence levels at shocks accurately. We study the effect of nonadiabatic wall temperature on the SBLI region, and the results are compared with available Direct Numerical Simulation (DNS) data.

Chemically bonded multi-nanolayer inorganic aerogel with a record-low thermal conductivity in a vacuum.

Inorganic aerogels have exhibited many superior characteristics with extensive applications, but are still plagued by a nearly century-old tradeoff between their mechanical and thermal properties. When reducing thermal conductivity by ultralow density, inorganic aerogels generally suffer from large fragility due to their brittle nature or weak joint crosslinking, while enhancing the mechanical robustness by material design and structural engineering, they easily sacrifice thermal insulation and stability. Here, we report a chemically bonded multi-nanolayer design and synthesis of a graphene/amorphous boron nitride aerogel to address this typical tradeoff to further enhance mechanical and thermal properties. Attributed to the chemically bonded interface and coupled toughening effect, our aerogels display a low density of 0.8mg cm-3 with ultrahigh flexibility (elastic compressive strain up to 99% and bending strain up to 90%), and exceptional thermostability (strength degradation <3% after sharp thermal shocks), as well as the lowest thermal conductivities in a vacuum (only 1.57mW m-1 K-1 at room temperature and 10.39mW m-1 K-1 at 500°C) among solid materials to date. This unique combination of mechanical and thermal properties offers an attractive material system for thermal superinsulation at extreme conditions.

Polyimide dynamically compressed to decomposition pressures: Two-wave structures captured by velocimetry and modeling

We performed a series of six plate impact experiments on polyimide and modeled them using new reactant and products equations of state combined with an Arrhenius rate model. The first experiment was diagnosed with embedded electromagnetic velocity gauges through which we directly observed attenuation of the lead shock to an approximately constant state over a propagation distance of roughly 4 mm. Simulated gauge profiles were in excellent qualitative agreement with experiment and suggested a sluggish chemical reaction that did not proceed to completion. The remaining five experiments were conducted in a transmission geometry and diagnosed velocimetrically at the sample/window interface. All five of these yielded profiles with a sharp shock followed by a more gradual approach to maximum interface velocity that was “rounded” to varying degree. These profiles proved difficult to interpret unambiguously due to the convolution of the reactive wave upon first shock with reflection of the lead wave and reshock or release by the window. Comparison with thermochemical calculations strongly suggests that the point of maximum interface velocity corresponds to the equilibrium reshock or release locus. We discuss the implications of this point for the practice of impedance matching based on the reflected Hugoniot of reactive materials such as polymers. The reactant and thermochemical products equations of state are developmental SESAME tables 97710 and 97720, respectively.

Single-Atom Metallophilic Sites for Liquid NaK Alloy Confinement toward Stable Alkali-Metal Anodes.

Room temperature liquid NaK alloy is a promising candidate for high performance metal batteries, due to its dendrite-free property and high energy density. However, its practical application is hindered by the high surface tension of liquid NaK, which causes difficulties in maintaining a stable contact with a current collector. Here, the authors demonstrate the extraordinary stable confinement of NaK alloy at room temperature by constructing a super-wetting substrate, which is based on highly dispersed cobalt-single-atom carbon nanoarrays. The developed liquid anode electrode prevented successfully the leakage of NaK alloy even in harsh stress (>5MPa) or sharp shock conditions. The symmetric cells achieved ultra-long reversible plating/stripping cycling life in both Na-ion (>1010hrs) and K-ion electrolytes (>4000hrs) at 10mA cm-2 /10mAh cm-2 . Upon fitting with Na3 V2 (PO4 )3 , the NaK assembled full battery provided high energy density (332.6kWh kg-1 ) and power density (11.05kW kg-1 ) with excellent stability after >21600cycles, which is the best value reported so far. The prepared pouch cell was able to drive a four-axis aircraft, demonstrating a great prospect in practical application. This work offers a new approach in the preparation of advanced dendrite-free liquid metal anodes with promising applications in electrochemical energy storage.

Pengenalan Luka Dan Macam-Macam Luka Dalam Kegiatan Tri Bakti PMR Di SDN Sukabumi I

If the continuity of the cell structure is cut off due to imperfect or missing cell elements, it can cause damage to the protective function of the skin and can be accompanied by damage to other cell elements, which can be called a wound. Wounds can also be accompanied by other tissue damage. Accidents involving moving vehicles, falls, sudden or sharp shocks, and surgical procedures are all potential causes of injury. There may be abrasions, lacerations, bruises, cuts, knife wounds, and gunshot wounds among the several types of injuries that may occur. The purpose of the Tri Bakti activities that researchers carried out at SDN SUKABUMI 1 was to introduce injuries and their various types in PPMR Tri Bakti activities. In various wound problems, wound care generally uses procedures. To ensure the wound healing process goes as it should, wound care must adapt to the circumstances and challenges of the injuries that arise.

Multi-label multi-class COVID-19 Arabic Twitter dataset with fine-grained misinformation and situational information annotations.

Since the inception of the current COVID-19 pandemic, related misleading information has spread at a remarkable rate on social media, leading to serious implications for individuals and societies. Although COVID-19 looks to be ending for most places after the sharp shock of Omicron, severe new variants can emerge and cause new waves, especially if the variants can evade the insufficient immunity provided by prior infection and incomplete vaccination. Fighting the fake news that promotes vaccine hesitancy, for instance, is crucial for the success of the global vaccination programs and thus achieving herd immunity. To combat the proliferation of COVID-19-related misinformation, considerable research efforts have been and are still being dedicated to building and sharing COVID-19 misinformation detection datasets and models for Arabic and other languages. However, most of these datasets provide binary (true/false) misinformation classifications. Besides, the few studies that support multi-class misinformation classification deal with a small set of misinformation classes or mix them with situational information classes. False news stories about COVID-19 are not equal; some tend to have more sinister effects than others (e.g., fake cures and false vaccine info). This suggests that identifying the sub-type of misinformation is critical for choosing the suitable action based on their level of seriousness, ranging from assigning warning labels to the susceptible post to removing the misleading post instantly. We develop comprehensive annotation guidelines in this work that define 19 fine-grained misinformation classes. Then, we release the first Arabic COVID-19-related misinformation dataset comprising about 6.7K tweets with multi-class and multi-label misinformation annotations. In addition, we release a version of the dataset to be the first Twitter Arabic dataset annotated exclusively with six different situational information classes. Identifying situational information (e.g., caution, help-seeking) helps authorities or individuals understand the situation during emergencies. To confirm the validity of the collected data, we define three classification tasks and experiment with various machine learning and transformer-based classifiers to offer baseline results for future research. The experimental results indicate the quality and validity of the data and its suitability for constructing misinformation and situational information classification models. The results also demonstrate the superiority of AraBERT-COV19, a transformer-based model pretrained on COVID-19-related tweets, with micro-averaged F-scores of 81.6% and 78.8% for the multi-class misinformation and situational information classification tasks, respectively. Label Powerset with linear SVC achieved the best performance among the presented methods for multi-label misinformation classification with micro-averaged F-scores of 76.69%.

A High-Order Residual-Based Viscosity Finite Element Method for the Ideal MHD Equations

We present a high order, robust, and stable shock-capturing technique for finite element approximations of ideal MHD. The method uses continuous Lagrange polynomials in space and explicit Runge-Kutta schemes in time. The shock-capturing term is based on the residual of MHD which tracks the shock and discontinuity positions, and adds sufficient amount of viscosity to stabilize them. The method is tested up to third order polynomial spaces and an expected fourth-order convergence rate is obtained for smooth problems. Several discontinuous benchmarks such as Orszag-Tang, MHD rotor, Brio-Wu problems are solved in one, two, and three spacial dimensions. Sharp shocks and discontinuity resolutions are obtained.

Hypocrystalline ceramic aerogels for thermal insulation at extreme conditions

Thermal insulation under extreme conditions requires materials that can withstand complex thermomechanical stress and retain excellent thermal insulation properties at temperatures exceeding 1,000 degrees Celsius1–3. Ceramic aerogels are attractive thermal insulating materials; however, at very high temperatures, they often show considerably increased thermal conductivity and limited thermomechanical stability that can lead to catastrophic failure4–6. Here we report a multiscale design of hypocrystalline zircon nanofibrous aerogels with a zig-zag architecture that leads to exceptional thermomechanical stability and ultralow thermal conductivity at high temperatures. The aerogels show a near-zero Poisson’s ratio (3.3 × 10−4) and a near-zero thermal expansion coefficient (1.2 × 10−7 per degree Celsius), which ensures excellent structural flexibility and thermomechanical properties. They show high thermal stability with ultralow strength degradation (less than 1 per cent) after sharp thermal shocks, and a high working temperature (up to 1,300 degrees Celsius). By deliberately entrapping residue carbon species in the constituent hypocrystalline zircon fibres, we substantially reduce the thermal radiation heat transfer and achieve one of the lowest high-temperature thermal conductivities among ceramic aerogels so far—104 milliwatts per metre per kelvin at 1,000 degrees Celsius. The combined thermomechanical and thermal insulating properties offer an attractive material system for robust thermal insulation under extreme conditions.

“Finger coral-like” ceramic fiber aerogels with enhanced high-temperature thermal insulation, anti-oxidation, and mechanical performance

Ceramic fiber aerogels are widely used in high temperature insulation. Thermal insulation under extreme conditions requires materials that can withstand sharp thermal shocks and extended high-temperature exposure. However, the oxidation resistance of traditional ceramic aerogels under a high temperature environment(above 1000 °C) is significantly reduced. Moreover, the strength degradation and structural collapse easily occurred, restricting their application. In this study, “Finger Coral” ceramic fibers aerogels were obtained using precursor infiltration pyrolysis (PIP) and electrospinning technology. The chemical composition, morphology, thermal insulation, mechanical properties and isothermal oxidation properties of SiCBN/SiC/disilicates nanoparticles (IP) ceramic aerogels were tested. The experimental results showed that the aerogel exhibited excellent mechanical and antioxidant properties over a wide range of temperatures. SiCBN/SiC/IP aerogel with a thickness of 2cm can achieve the heat insulation effect of the hot surface temperature from 1600 to 75 °C, and the thermal conductivity at 1600 °C was only 0.036 W·m−1·K−1. After 12 h of cyclic isothermal oxidation, the weight gain was only 5.17%. The material is expected to be used as a thermal insulation material in the hypersonic vehicles.

Corporate diversification and downsizing decisions: International evidence from sharp and sudden performance shocks

While firms regularly reduce workforce following sharp performance decline, diversified firms may abstain from employment downsizing by transferring capital and labor between segments (the allocative flexibility effect). However, downsizing may be more likely if a performance shock leads to efforts to reduce inefficiency in resource allocation (the inefficient internal market effect). Using a large cross-country dataset, our results provide strong support for the inefficient internal market effect. We find that diversified firms are more likely to downsize and the national employment protection and union power laws moderate this link. We also find that diversified firms with more excess employment are more likely to downsize and that downsizing following major adverse performance shocks is associated with lower level of diversification and excess employment.

Investigating the migration of immiscible contaminant fluid flow in homogeneous and heterogeneous aquifers with high-precision numerical simulations.

Numerical modeling of the migration of three-phase immiscible fluid flow in variably saturated zones is challenging due to the different behavior of the system between unsaturated and saturated zones. This behavior results in the use of different numerical methods for the numerical simulation of the fluid flow depending on whether it is in the unsaturated or saturated zones. This paper shows that using a high-resolution shock-capturing conservative method to resolve the nonlinear governing coupled partial differential equations of a three-phase immiscible fluid flow allows the numerical simulation of the system through both zones providing a unitary vision (and resolution) of the migration of an immiscible contaminant problem within a porous medium. In particular, using different initial scenarios (including impermeable “lenses” in heterogeneous aquifers), three-dimensional numerical simulation results are presented on the temporal evolution of the contaminant migration following the saturation profiles of the three-phases fluids flow in variably saturated zones. It is considered either light nonaqueous phase liquid with a density less than the water, or dense nonaqueous phase liquid, which has densities greater than the water initially released in unsaturated dry soil. Our study shows that the fate of the migration of immiscible contaminants in variably saturated zones can be accurately described, using a unique mathematical conservative model, with different evolution depending on the value of the system’s physical parameters, including the contaminant density, and accurately tracking the evolution of the sharp (shock) contaminant front.