Computational modelling of multi-material energetic materials and systems

ABSTRACT

In this paper, we present a systematic roadmap for developing a robust and parallel multi-material reactive hydrodynamic solver that integrates historically stable algorithms with new and current modern methods to solve explosive system design problems. The Ghost Fluid Method and Riemann solvers were used to enforce appropriate interface boundary conditions. Improved performance in terms of computational work and convergence properties was achieved by modifying a local node sorting strategy that decouples ghost nodes, allowing us to set material boundary conditions via an explicit procedure, removing the need to solve a coupled system of equations numerically. The locality and explicit nature of the node sorting concept allows for greater levels of parallelism and lower computational cost when populating ghost nodes. Non-linear numerical issues endemic to the use of real Equations of State in hydro-codes were resolved by using more thermodynamically consistent forms allowing us to accurately resolve large density gradients associated with high energy detonation problems at material interfaces. Pre-computed volume tables were implemented adding to the robustness of the solver base.     

Contingent derivatives and regularization for noncoercive inverse problems

Abstract

We study the inverse problem of parameter identification in noncoercive variational problems that commonly appear in applied models. We examine the differentiability of the set-valued parameter-to-solution map using the first-order and the second-order contingent derivatives. We explore the inverse problem using the output least-squares and the modified output least-squares objectives. By regularizing the noncoercive variational problem, we obtain a single-valued regularized parameter-to-solution map and investigate its smoothness and boundedness. We also consider optimization problems using the output least-squares and the modified output least-squares objectives for the regularized variational problem. We give a complete convergence analysis showing that for the output least-squares and the modified output least-squares, the regularized minimization problems approximate the original optimization problems suitably. We also provide the first-order and the second-order adjoint method for the computation of the first-order and the second-order derivatives of the output least-squares objective. We provide discrete formulas for the gradient and the Hessian calculation and present numerical results.     

Crystal structure and DFT calculations of methanol {E-N′-(2-hydroxybenzlidene)benzohydrazido}dioxidomolybdenum (VI) complex

The crystal and molecular structure of methanol {E-N′-(2-hydroxybenzlidene)benzohydrazido}dioxidomolybdenum (VI) was determined by single crystal X-ray diffraction. B3LYP/DZP basis set theoretical calculations nicely reproduce the X-ray experimental geometry. The obtained results are discussed in connection with the electronic and structural properties of the compound.     

Extremal A-statistical limit points via ideals

In this paper, following the line of recent work of Savaş et al. [20] we apply the notion of ideals to A-statistical limit superior and inferior for a sequence of real numbers.     

QSAR Modeling of Large Heterogeneous Sets of Molecules

Abstract

In aquatic toxicology, QSAR models are generally designed for chemicals presenting the same mode of toxic action. Their proper use provides good simulation results. Problems arise when the mechanism of toxicity of a chemical is not clearly identified. Indeed, in that case, the inappropriate application of a specific QSAR model can lead to a dramatic error in the toxicity estimation. With the advent of powerful computers and easy access to them, and the introduction of soft modeling and artificial intelligence in SAR and QSAR, radically different models, designed from large non-congeneric sets of chemicals have been proposed. Some of these new QSAR models are reviewed and their originality, advantages, and limitations are stressed.     

Gaussian basis sets for molecular applications

The choice of basis set in quantum chemical calculations can have a huge impact on the quality of the results, especially for correlated ab initio methods. This article provides an overview of the development of Gaussian basis sets for molecular calculations, with a focus on four popular families of modern atom‐centered, energy‐optimized bases: atomic natural orbital, correlation consistent, polarization consistent, and def2. The terminology used for describing basis sets is briefly covered, along with an overview of the auxiliary basis sets used in a number of integral approximation techniques and an outlook on possible future directions of basis set design. © 2012 Wiley Periodicals, Inc.     

Polarization functions for the modified m6‐31G basis sets for atoms Ga through Kr

The 2df polarization functions for the modified m6‐31G basis sets of the third‐row atoms Ga through Kr (Int J Quantum Chem, 2007, 107, 3028; Int J. Quantum Chem, 2009, 109, 1158) are proposed. The performances of the m6‐31G, m6‐31G(d,p), and m6‐31G(2df,p) basis sets were examined in molecular calculations carried out by the density functional theory (DFT) method with B3LYP hybrid functional, Møller‐Plesset perturbation theory of the second order (MP2), quadratic configuration interaction method with single and double substitutions and were compared with those for the known 6‐31G basis sets as well as with the other similar 641 and 6‐311G basis sets with and without polarization functions. Obtained results have shown that the performances of the m6‐31G, m6‐31G(d,p), and m6‐31G(2df,p) basis sets are better in comparison with the performances of the known 6‐31G, 6‐31G(d,p) and 6‐31G(2df,p) basis sets. These improvements are mainly reached due to better approximations of different electrons belonging to the different atomic shells in the modified basis sets. Applicability of the modified basis sets in thermochemical calculations is also discussed. © 2013 Wiley Periodicals, Inc.     

New basis sets for the evaluation of interaction‐induced electric properties in hydrogen‐bonded complexes

Interaction‐induced static electric properties, that is, dipole moment, polarizability, and first hyperpolarizability, of the CO? (HF)_n and N₂? (HF)_n, n = 1–9 hydrogen‐bonded complexes are evaluated within the finite field approach using the Hartree–Fock, density functional theory, Møller–Plesset second‐order perturbation theory, and coupled cluster methods, and the LPol‐n (n = ds, dl, fs, fl) basis sets. To compare the performance of the different methods with respect to the increase of the complex size, we consider as model systems linear chains of the complexes. We analyze the results in terms of the many‐body and cooperative effects. © 2012 Wiley Periodicals, Inc.     

Resonant equations with the Neumann p-Laplacian plus an indefinite potential

We consider a nonlinear Neumann problem driven by the p  -Laplacian plus an indefinite potential and a Carathéodory reaction which at ±∞ is resonant with respect to any nonprincipal variational eigenvalue of the differential operator. Using critical point theory and Morse theory (critical groups), we show that the problem has at least three nontrivial smooth solutions, two of which have constant sign. In the process we prove some results of independent interest concerning the unique continuation property of eigenfunctions and the critical groups at infinity of a C¹$C^1$-functionals.     

Reachable set bounding for nonlinear perturbed time-delay systems: The smallest bound

In this letter, we propose a new approach to obtain the smallest box which bounds all reachable sets of a class of nonlinear time-delay systems with bounded disturbances. A numerical example is studied to illustrate the obtained result.