Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Design pattern”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Regiospecific N -alkyl substitution tunes the molecular packing of high-performance non-fullerene acceptors

The rapid development of non-fullerene acceptors (NFAs) with strong near-infrared absorption has led to remarkably enhanced short-circuit current density (J sc ) values in organic solar cells (OSCs). NFAs based on the benzotriazole (Bz) fused-ring π-core have great potential in delivering both high J sc and decent open-circuit voltage values due to their strong intramolecular charge transfer with reasonably low energy loss. Here in this work, we have designed and synthesized a series of Bz-based NFAs, PN6SBO-4F, AN6SBO-4F and EHN6SEH-4F, via regiospecific N-alkyl engineering based on the high-performance NFA mBzS-4F that was reported previously. The molecular packing of mBzS-4F, AN6SBO-4F, and EHN6SEH-4F single crystals was analyzed using X-ray crystallography in order to provide a comprehensive understanding of the correlation between the molecular structure, the charge-transporting properties, and the solar cell performance. Compared with the typical honeycomb single-crystal structure of Y6 derivatives, these NFAs exhibit distinctly different molecular packing patterns. The strong interactions of terminal indanone groups in mBzS-4F and the J-aggregate-like packing in EHN6SEH-4F lead to the formation of ordered 3D networks in single-crystals with channels for efficient charge transport. Consequently, OSCs based on mBzS-4F and EHN6SEH-4F show efficient photon-to-current conversions, achieving the highest power conversion efficiency of 17.48% with a Jsc of 28.83 mA cm -2 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Protecting Websites from Cross-Site Scripting (XSS) Attacks: A Novel Configuration using Pulse Secure © Pulse Connect Secure © and Virtual Web Application Firewall (vWAF)

Cross-site scripting (XSS), one of the most prevalent forms of client-side attacks, is when bad actors attempt to access sensitive information from the backend web server and other systems on the backend network. Some XSS attacks attempt to access client-side sensitive information, such as cookies. Web application firewalls (WAFs) are a first line of defense where common Uniform Resource Locator (URL) patterns are analyzed to detect and block known attacks. This paper describes a novel configuration using the Pulse Secure © Pulse Connect Secure © (PCS © ) Secure Socket Layer Virtual Private Network software and Virtual Web Application Firewall (vWAF) that protects a website from XSS attacks. This paper also presents novel aspects of the configuration that control the redirection of traffic through the vWAF and provide fine-grained behavioral control at the application level while decoupling the PCS and vWAF configurations. The intended audience for this paper comprises system and site administrators who are familiar with standard web server environments. These configuration details might prove useful during the design of a more secure infrastructure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
DOE OSTI2024

The Pixel Anomaly Detection Tool : a user-friendly GUI for classifying detector frames using machine-learning approaches

Data collection at X-ray free electron lasers has particular experimental challenges, such as continuous sample delivery or the use of novel ultrafast high-dynamic-range gain-switching X-ray detectors. This can result in a multitude of data artefacts, which can be detrimental to accurately determining structure-factor amplitudes for serial crystallography or single-particle imaging experiments. Here, a new data-classification tool is reported that offers a variety of machine-learning algorithms to sort data trained either on manual data sorting by the user or by profile fitting the intensity distribution on the detector based on the experiment. This is integrated into an easy-to-use graphical user interface, specifically designed to support the detectors, file formats and software available at most X-ray free electron laser facilities. The highly modular design makes the tool easily expandable to comply with other X-ray sources and detectors, and the supervised learning approach enables even the novice user to sort data containing unwanted artefacts or perform routine data-analysis tasks such as hit finding during an experiment, without needing to write code.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A Resonant Metamaterial Line Array for Ultrasound Compressive Imaging

Abstract Acoustic metamaterials have been proposed for numerous applications including subwavelength imaging, impedance matching, and lensing. Yet, their application in compressive sensing and imaging has not been fully investigated. When metamaterials are used as resonators at certain frequencies, they can generate random radiation patterns in the transmitted waves from the transducers and received waves from a target. Compressive sensing favors such randomness inasmuch as it can increase incoherence by decreasing the amount of mutual information between any two different measurements. This study aims at assessing whether the use of resonating metamaterial unit cells in a single-layered non-optimized array between a number of ultrasound transceivers and targets can improve the sensing capacity, point-spread function of the sensing array (their beam focusing ability), and imaging performance in point-like target detection. The theoretical results are promising and can open the way for more efficient metamaterial designs with the aim of enhancing ultrasound imaging with lower number of transceivers compared to the regular systems.

Acoustics↗
DOE OSTI2022

Overview of the MOSAiC expedition—Atmosphere

With the Arctic rapidly changing, the needs to observe, understand, and model the changes are essential. To support these needs, an annual cycle of observations of atmospheric properties, processes, and interactions were made while drifting with the sea ice across the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition from October 2019 to September 2020. An international team designed and implemented the comprehensive program to document and characterize all aspects of the Arctic atmospheric system in unprecedented detail, using a variety of approaches, and across multiple scales. These measurements were coordinated with other observational teams to explore cross-cutting and coupled interactions with the Arctic Ocean, sea ice, and ecosystem through a variety of physical and biogeochemical processes. This overview outlines the breadth and complexity of the atmospheric research program, which was organized into 4 subgroups: atmospheric state, clouds and precipitation, gases and aerosols, and energy budgets. Atmospheric variability over the annual cycle revealed important influences from a persistent large-scale winter circulation pattern, leading to some storms with pressure and winds that were outside the interquartile range of past conditions suggested by long-term reanalysis. Similarly, the MOSAiC location was warmer and wetter in summer than the reanalysis climatology, in part due to its close proximity to the sea ice edge. The comprehensiveness of the observational program for characterizing and analyzing atmospheric phenomena is demonstrated via a winter case study examining air mass transitions and a summer case study examining vertical atmospheric evolution. Overall, the MOSAiC atmospheric program successfully met its objectives and was the most comprehensive atmospheric measurement program to date conducted over the Arctic sea ice. The obtained data will support a broad range of coupled-system scientific research and provide an important foundation for advancing multiscale modeling capabilities in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The role of driver head pose dynamics and instantaneous driving in safety critical events: Application of computer vision in naturalistic driving

This paper investigates the role of driver behavior especially head pose dynamics in safety–critical events (SCEs). Using a large dataset collected in a naturalistic driving study, this paper analyzes the head pose dynamics and driving behavior in moments leading up to crashes or near-crashes. The study uses advanced computer vision and mixed logit modeling techniques to identify patterns and relationships between drivers’ head pose dynamics and crash involvement. The results suggest that driver-head pose dynamics, especially poses that indicate distraction and movement volatility, are important factors that can contribute to undesirable safety outcomes. Marginal effects show that angular deviation for head pose dynamics indicated by yaw, pitch and roll increase the likelihood of crash intensity by 4.56%, 4.92% and 8.26% respectively. Furthermore, traffic flow and lane changing also contribute to increase in likelihood of crash intensity. These findings provide new insights into pre-crash factors, especially human factors and safety–critical events. The study highlights the importance of considering human factors in designing driver assistance systems and developing safer vehicles. This research contributes by examining naturalistic driving data at the microscopic level with early detection of behaviors that lead to SCEs and provides a basis for future research on automation.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

A method for determining optimal electric range by considering electric vehicle lightweighting on perceived ownership cost

The limited driving range due to high costs and low energy densities of batteries constraints the battery electric vehicle (BEV) market growth. Lightweighting in theory can reduce energy consumption rate and extend the driving range. The knowledge gap is to quantitatively link the cost-effectiveness of light-weight technologies with range extension and consumer acceptance of BEVs. In this study, a physics-based energy consumption model of BEVs is constructed and associated with a statistics-based model on the basis of travel surveys. A perceived cost of ownership (PCO) is then developed by adding intangible costs to traditional total cost of ownership models. We estimate, at the disaggregate vehicle model and driver level and the aggregate market level, 1) the extended range due to lightweighting for a given battery size; and 2) the optimal electric range based on lightweighting decisions. The cost-effectiveness of lightweighting for BEV range extension is found to vary with income-dependent daily range limitation value, driving patterns and lightweighting technology costs. In general, adopting lightweighting in BEVs is more cost-effective for consumers with higher daily limitation value, as well as for those with higher driving intensity or suitable daily driving patterns. When the lightweighting involves a higher vehicle production cost, less lightweighting could reduce the overall PCO for BEV owners. 4 of the selected top ten BEV models are found to benefit from additional 2.09%–4.45% lightweighting. Finally, the method built in this study can guide automakers in planning R&D investments in battery and lightweighting technologies.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Tuning commensurability in twisted van der Waals bilayers

Moiré superlattices based on van der Waals bilayers created at small twist angles lead to a long wavelength pattern with approximate translational symmetry. At large twist angles ( θ t ), moiré patterns are, in general, incommensurate except for a few discrete angles. Here we show that large-angle twisted bilayers offer distinctly different platforms. More specifically, by using twisted tungsten diselenide bilayers, we create the incommensurate dodecagon quasicrystals at θ t = 30° and the commensurate moiré crystals at θ t = 21.8° and 38.2°. Further, valley-resolved scanning tunnelling spectroscopy shows disparate behaviours between moiré crystals (with translational symmetry) and quasicrystals (with broken translational symmetry). In particular, the K valley shows rich electronic structures exemplified by the formation of mini-gaps near the valence band maximum. These discoveries demonstrate that bilayers with large twist angles offer a design platform to explore moiré physics beyond those formed with small twist angles.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Oxygen redox chemistry in P2-Na 0.6 Li 0.11 Fe 0.27 Mn 0.62 O 2 cathode for high-energy Na-ion batteries

Owing to the abundance of raw material reserves and low cost, Na-ion batteries (NIBs) have successfully gained widespread attention from academic and industrial communities in the past few decades. However, the insufficient cathode energy density is still one of the critical bottlenecks restricting the development of NIBs. Following a strategy of introducing Li+ into the transition-metal (TM) layer to enhance the oxygen redox reaction, a novel layered cathode material P2-Na 0.6 Li 0.11 Fe 0.27 Mn 0.62 O 2 (NLFMO) was designed and successfully synthesized. This NLFMO cathode not only delivers a large initial reversible capacity of 207.3 mA h g -1 , but also shows a good cycling performance (104.2 mA h g -1 after 80 cycles) and rate capability (126.2 mA h g -1 at 1C). The ultrahigh capacity is contributed by both cationic (Fe 3+ /Fe 4+ and Mn 3+ /Mn 4+ ) and partially reversible anionic redox (O 2- /O n- ) reactions, revealed by in situ X-ray absorption spectroscopy (XAS) and X-ray photoelectron spectroscopy (XPS) techniques. Moreover, no detrimental P2–O2 phase transition was observed in ex situ X-ray diffraction (XRD) patterns, confirming the high structural stability during Na + deintercalation/intercalation processes. These results provide valuable information about the high-energy density layered cathode materials based on anionic redox reactions for NIBs.

25 ENERGY STORAGE↗
DOE OSTI2026

Data-Driven Modeling and Control of Systems with Plasma-Surface Interactions (Final Technical Report)

This final technical report summarizes the activities and accomplishments in the period from February 2023 thru January 2026. The objective of the proposed research is to investigate the physical mechanisms and processes underlying the formation of structures and patterns in systems with plasma-surface interactions. In the past decades, there have been extensive studies on the interaction of glow discharges, dielectric barrier discharges, and arc discharges with confining or intervening surfaces. The advancement of the understanding of these phenomena is not only of fundamental scientific interest and relevance to the knowledge of the plasma state, but also with profound implications in various technological applications. The research will integrate theoretical, computational, and experimental work within an innovative framework of data assimilation, i.e., optimally combining model predictions with measurements. The scientific merit of this research has three aspects. Firstly, it extends the studies of plasma-surface interactions to systems with insulator surfaces and multi-layer systems, while existing studies are predominantly on electrode surfaces. Secondly, it expects to develop a novel data-driven modeling approach based on data assimilation to enhance the predictive and control capabilities, which could make transformative contributions to basic plasma research. Thirdly, it will shed new light on outstanding problems related to formation of patterns interfacing plasmas. This project also aims to launch an education and outreach initiative at Texas A&M University-Kingsville, a non-R1, minority-serving institution in South Texas. The initiative is structured as a four-tier pyramid. Tier one will be a webinar series for culture and capacity building to inform broader audience in the region about the research fields of plasma science and engineering. Tier two will be the creation and offering of an upper-level undergraduate course on introductory plasma physics, which will help with the recruitment for the upper tiers. On tier three, we will engage and mentor senior design students to conduct work toward the research goal of this project. There will also be a certificate program on general plasma science for undergrad and graduate students, part of which will be lab training at Princeton University. Tier four will be the supervision and mentoring of Ph.D. students. Therefore, this project will systematically expand the talent pipeline, broaden participation from communities historically and geographically underrepresented in DOE SC research portfolio, significantly improve the research and education capacity at the PI’s institution, and contribute to developing a diverse workforce in plasma science and engineering.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Efficient Parallelization of Irregular Applications on GPU Architectures

With the enlarging computation capacity of general Graphics Processing Units (GPUs), leveraging GPUs to accelerate parallel applications has become a critical topic in academia and industry. However, a wide range of irregular applications with the computation-/memory-intensive nature cannot easily achieve high GPU utilization. The challenges mainly involve the following aspects: first, data dependence leads to coarse-grained kernel and inefficient parallelism; second, heavy GPU memory usage may cause frequent memory evictions and extra overhead of I/O; third, specific computation patterns produce memory redundancies; last, workload balance and data reusability conjunctly benefit the overall performance, but there may exist a dynamic trade-off between them. Targeting these challenges, this dissertation proposes multiple optimizations to accelerate two real-world applications: many-body correlation functions to simulate nuclear physics in a large-scale scientific system; the other is the eALS-based matrix factorization recommendation system. To accelerate the calculations of many-body correlation functions, this dissertation presents three frameworks in GPU memory management and multi-GPU scheduling. Firstly, an optimized systematic GPU memory management framework, MemHC, utilizes a series of new memory reduction designs in GPU memory allocation, CPU/GPU communications, and GPU memory oversubscription. Secondly, an enhanced multi-GPU scheduling framework, MICCO, particularly by taking both data dimension (e.g., data reuse and data eviction) and computation dimension into account. MICCO designs a heuristic scheduling algorithm and a machine learning-based regression model to generate the optimal settings of a proposed new concept to manage the trade-off. Thirdly, a locality-aware multi-GPU scheduling framework. This scheduler leverages pipeline batch generation with a looking-ahead strategy by building local dependency graphs for memory transfer reduction and better data reuse, achieving up to 79.92% memory cost reduction and 1.67x speedup. To parallelize the eALS-based recommendation system, this dissertation proposes an efficient CPU/GPU heterogeneous recommendation system, HEALS. HEALS employs newly designed architecture-adaptive data formats to achieve load balance and good data locality on CPU and GPU. To mitigate the data dependence, HEALS presents a CPU/GPU collaboration model for both task parallelism and data parallelism with multiple kernel computation optimizations. In summary, this dissertation efficiently accelerates two typical irregular applications on GPUs by building four frameworks, including CPU/GPU collaboration, GPU memory management, and multi-GPU scheduling.

Wang, Qihan↗
DOE OSTI2022

Unveiling the interaction of nanopatterned void superlattices with irradiation cascades

Nanopatterned microstructures in materials can have a profound impact on materials’ physical and chemical properties. While voids are typically considered as detrimental defects in irradiated materials, the patterning of nanoscale voids causes the formation of void superlattices and provides a highly efficient mechanism for gas storage. Despite the important applications of nanopatterned defect superlattices, how they degrade under irradiation remains unclear. Here we provide direct observation of the evolution of void superlattices under irradiation and elucidate the interaction of void superlattices with irradiation cascades. We reveal that the instability of void superlattices under irradiation is caused by heterogenous void shrinkage and demonstrate the imperative role of mixed 1D/3D diffusion of self-interstitial atoms and injected inert gas atoms on void shrinkage and void superlattice instability. Understanding the degradation mechanisms of nanopatterned microstructures is essential to designing damage-tolerant materials and broadening their applications in extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Atomic-scale manipulation of buried graphene–silicon carbide interface by local electric field

Abstract Precision of scanning tunneling microscopy (STM) enables control of matter at scales of single atoms. However, transition from atomic-scale manipulation strategies to practical devices encounters fundamental problems in protection of the designer structures formed atop the surface. In this context, STM manipulation of subsurface structures on technologically relevant materials is encouraging. Here, we propose a material platform and protocols for precise manipulation of a buried graphene interface. We show that an electric field from the STM tip reversibly controls breaking and restoring of covalent bonds between the graphene buffer layer and the SiC substrate. The process involves charge redistribution at the atomically sharp interface plane under the epitaxial graphene layer(s). This buried manipulation platform is laterally defined by unit cells from the corresponding (6×6) SiC moiré lattice of the epitaxial graphene. Local and reversible electric-field-induced patterning of graphene heterostructures from the bottom interface creates an alternative architecture concept for their applications.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Resilient stellarator divertor characteristics in the Helically Symmetric eXperiment

Resilient divertor features connected to open chaotic edge structures in the Helically Symmetric eXperiment are investigated. For the first time, an expanded vessel wall was considered that would give space for implementation of a physical divertor target structure. The analysis was done for four different magnetic configurations with very different chaotic plasma edges. A resilient plasma wall interaction pattern was identified across all configurations. This manifests as qualitatively very similar footprint behavior across the different plasma equilibria. Overall, the resilient field lines of interest with high connection length L C lie within a helical band along the wall for all configurations. This resiliency can be used to identify the best location of a divertor. The details of the magnetic footprint’s resilient helical band is subject to specific field line structures which are linked to the penetration depth of field lines into the plasma and directly influence the heat and particle flux patterns. The differences arising from these details are characterized by introducing a new metric, the minimum radial connection min $(\delta_N)$ of a field line from the last closed flux surface. The relationship, namely the deviation from a scaling law, between min $(\delta_N)$ and L C of the field lines in the plasma edge field line behavior suggests that the field lines are associated with structures such as resonant islands, cantori, and turnstiles. This helps determine the relevant magnetic flux channels based on the radial location of these chaotic edge structures and the divertor target footprint. These details will need to be taken into account for resilient divertor design.

divertor↗
DOE OSTI2021

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Experimental and theoretical rationalization for the base pairing abilities of inosine, guanosine, adenosine, and their corresponding 8-oxo -7,8-dihydropurine, and 8-bromopurine analogues within A-form duplexes of RNA

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (Tm) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

On differentiable local bounds preserving stabilization for Euler equations

This work presents the design of nonlinear stabilization techniques for the finite element discretization of Euler equations in both steady and transient form. Implicit time integration is used in the case of the transient form. A differentiable local bounds preserving method has been developed, which combines a Rusanov artificial diffusion operator and a differentiable shock detector. Nonlinear stabilization schemes are usually stiff and highly nonlinear. This issue is mitigated by the differentiability properties of the proposed method. Moreover, in order to further improve the nonlinear convergence, we also propose a continuation method for a subset of the stabilization parameters. The resulting method has been successfully applied to steady and transient problems with complex shock patterns. Numerical experiments show that it is able to provide sharp and well resolved shocks. Furthermore, the importance of the differentiability is assessed by comparing the new scheme with its non-differentiable counterpart. Numerical experiments suggest that, for up to moderate nonlinear tolerances, the method exhibits improved robustness and nonlinear convergence behavior for steady problems. Additionally, in the case of transient problem, we also observe a reduction in the computational cost.

42 ENGINEERING↗