A Hybrid Numerical Analysis Method for Structural Health Monitoring
NASA Technical Reports Server (NTRS)
Forth, Scott C.; Staroselsky, Alexander
2001-01-01
A new hybrid surface-integral-finite-element numerical scheme has been developed to model a three-dimensional crack propagating through a thin, multi-layered coating. The finite element method was used to model the physical state of the coating (far field), and the surface integral method was used to model the fatigue crack growth. The two formulations are coupled through the need to satisfy boundary conditions on the crack surface and the external boundary. The coupling is sufficiently weak that the surface integral mesh of the crack surface and the finite element mesh of the uncracked volume can be set up independently. Thus when modeling crack growth, the finite element mesh can remain fixed for the duration of the simulation as the crack mesh is advanced. This method was implemented to evaluate the feasibility of fabricating a structural health monitoring system for real-time detection of surface cracks propagating in engine components. In this work, the authors formulate the hybrid surface-integral-finite-element method and discuss the mechanical issues of implementing a structural health monitoring system in an aircraft engine environment.
Asymmetric MRI magnet design using a hybrid numerical method.
Zhao, H; Crozier, S; Doddrell, D M
1999-12-01
This paper describes a hybrid numerical method for the design of asymmetric magnetic resonance imaging magnet systems. The problem is formulated as a field synthesis and the desired current density on the surface of a cylinder is first calculated by solving a Fredholm equation of the first kind. Nonlinear optimization methods are then invoked to fit practical magnet coils to the desired current density. The field calculations are performed using a semi-analytical method. A new type of asymmetric magnet is proposed in this work. The asymmetric MRI magnet allows the diameter spherical imaging volume to be positioned close to one end of the magnet. The main advantages of making the magnet asymmetric include the potential to reduce the perception of claustrophobia for the patient, better access to the patient by attending physicians, and the potential for reduced peripheral nerve stimulation due to the gradient coil configuration. The results highlight that the method can be used to obtain an asymmetric MRI magnet structure and a very homogeneous magnetic field over the central imaging volume in clinical systems of approximately 1.2 m in length. Unshielded designs are the focus of this work. This method is flexible and may be applied to magnets of other geometries. PMID:10579958
Trigonometrically fitted two step hybrid method for the numerical integration of second order IVPs
NASA Astrophysics Data System (ADS)
Monovasilis, Th.; Kalogiratou, Z.; Simos, T. E.
2016-06-01
In this work we consider the numerical integration of second order ODEs where the first derivative is missing. We construct trigonometrically fitted two step hybrid methods. We apply the new methods on the numerical integration of several test problems.
Hybrid Numerical Methods for Multiscale Simulations of Subsurface Biogeochemical Processes
Scheibe, Timothy D.; Tartakovsky, Alexandre M.; Tartakovsky, Daniel M.; Redden, George D.; Meakin, Paul
2007-08-01
Many subsurface flow and transport problems of importance today involve coupled non-linear flow, transport, and reaction in media exhibiting complex heterogeneity. In particular, problems involving biological mediation of reactions fall into this class of problems. Recent experimental research has revealed important details about the physical, chemical, and biological mechanisms involved in these processes at a variety of scales ranging from molecular to laboratory scales. However, it has not been practical or possible to translate detailed knowledge at small scales into reliable predictions of field-scale phenomena important for environmental management applications. A large assortment of numerical simulation tools have been developed, each with its own characteristic scale including molecular (e.g., molecular dynamics), microbial (e.g., cellular automata or particle individual-based models), pore (e.g., lattice-Boltzmann, pore network models, and discrete particle methods such as smoothed particle hydrodynamics) and continuum scales (e.g., traditional partial differential equations solved by finite difference or finite element methods). While many problems can be effectively addressed by one of these models at a single scale, some problems may require explicit integration of models across multiple scales. We are developing a hybrid multi-scale subsurface reactive transport modeling framework that integrates models with diverse representations of physics, chemistry and biology at different scales (sub-pore, pore and continuum). The modeling framework is being designed to take advantage of advanced computational technologies including parallel code components using the Common Component Architecture, parallel solvers, gridding, data and workflow management, and visualization. This paper describes the specific methods/codes being used at each scale, techniques used to directly and adaptively couple across model scales, and preliminary results of application to a
A Hybrid Numerical Method for Turbulent Mixing Layers. Degree awarded by Case Western Reserve Univ.
NASA Technical Reports Server (NTRS)
Georgiadis, Nicholas J.
2001-01-01
A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modern day aircraft and also those of hypersonic vehicles currently under development. The method configurations in which a dominant structural feature provides an unsteady mechanism to drive the turbulent development in the mixing layer. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS-LES method on stretched, non-Cartesian grids. Closure for the RANS equations was obtained using the Cebeci-Smith algebraic turbulence model in conjunction with the wall-function approach of Ota and Goldberg. The wall-function approach enabled a continuous computational grid from the RANS regions to the LES region. The LES equations were closed using the Smagorinsky subgrid scale model. The hybrid RANS-LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Vortex shedding from the base region of a splitter plate separating the upstream flows was observed to eventually transition to turbulence. The location of the transition, however, was much further downstream than indicated by experiments. Actual LES calculations, performed in three spatial directions, also indicated vortex shedding, but the transition to turbulence was found to occur much closer to the beginning of the mixing section. which is in agreement with experimental observations. These calculations demonstrated that LES simulations must be performed in three dimensions. Comparisons of time-averaged axial velocities and turbulence intensities indicated reasonable agreement with experimental
Two step hybrid methods of 7th and 8th order for the numerical integration of second order IVPs
NASA Astrophysics Data System (ADS)
Kalogiratou, Z.; Monovasilis, Th.; Simos, T. E.
2016-06-01
In this work we consider the numerical integration of second order ODEs where the first derivative is missing. We construct two step hybrid methods with six and seven stages and seventh and eighth algebraic order. We apply the new methods on the numerical integration of several test problems.
Comparison of symbolic and numerical integration methods for an assumed-stress hybrid shell element
NASA Technical Reports Server (NTRS)
Rengarajan, Govind; Knight, Norman F., Jr.; Aminpour, Mohammad A.
1993-01-01
Hybrid shell elements have long been regarded with reserve by the commercial finite element developers despite the high degree of reliability and accuracy associated with such formulations. The fundamental reason is the inherent higher computational cost of the hybrid approach as compared to the displacement-based formulations. However, a noteworthy factor in favor of hybrid elements is that numerical integration to generate element matrices can be entirely avoided by the use of symbolic integration. In this paper, the use of the symbolic computational approach is presented for an assumed-stress hybrid shell element with drilling degrees of freedom and the significant time savings achieved is demonstrated through an example.
Hybrid undulator numerical optimization
Hairetdinov, A.H.; Zukov, A.A.
1995-12-31
3D properties of the hybrid undulator scheme arc studied numerically using PANDIRA code. It is shown that there exist two well defined sets of undulator parameters which provide either maximum on-axis field amplitude or minimal higher harmonics amplitude of the basic undulator field. Thus the alternative between higher field amplitude or pure sinusoidal field exists. The behavior of the undulator field amplitude and harmonics structure for a large set of (undulator gap)/(undulator wavelength) values is demonstrated.
Stress analysis and damage evaluation of flawed composite laminates by hybrid-numerical methods
NASA Technical Reports Server (NTRS)
Yang, Yii-Ching
1992-01-01
Structural components in flight vehicles is often inherited flaws, such as microcracks, voids, holes, and delamination. These defects will degrade structures the same as that due to damages in service, such as impact, corrosion, and erosion. It is very important to know how a structural component can be useful and survive after these flaws and damages. To understand the behavior and limitation of these structural components researchers usually do experimental tests or theoretical analyses on structures with simulated flaws. However, neither approach has been completely successful. As Durelli states that 'Seldom does one method give a complete solution, with the most efficiency'. Examples of this principle is seen in photomechanics which additional strain-gage testing can only average stresses at locations of high concentration. On the other hand, theoretical analyses including numerical analyses are implemented with simplified assumptions which may not reflect actual boundary conditions. Hybrid-Numerical methods which combine photomechanics and numerical analysis have been used to correct this inefficiency since 1950's. But its application is limited until 1970's when modern computer codes became available. In recent years, researchers have enhanced the data obtained from photoelasticity, laser speckle, holography and moire' interferometry for input of finite element analysis on metals. Nevertheless, there is only few of literature being done on composite laminates. Therefore, this research is dedicated to this highly anisotropic material.
A numerical study of hybrid optimization methods for the molecular conformation problems
Meza, J.C.; Martinez, M.L.
1993-05-01
An important area of research in computational biochemistry is the design of molecules for specific applications. The design of these molecules depends on the accurate determination of their three-dimensional structure or conformation. Under the assumption that molecules will settle into a configuration for which their energy is at a minimum, this design problem can be formulated as a global optimization problem. The solution of the molecular conformation problem can then be obtained, at least in principle, through any number of optimization algorithms. Unfortunately, it can easily be shown that there exist a large number of local minima for most molecules which makes this an extremely difficult problem for any standard optimization method. In this study, we present results for various optimization algorithms applied to a molecular conformation problem. We include results for genetic algorithms, simulated annealing, direct search methods, and several gradient methods. The major result of this study is that none of these standard methods can be used in isolation to efficiently generate minimum energy configurations. We propose instead several hybrid methods that combine properties of several local optimization algorithms. These hybrid methods have yielded better results on representative test problems than single methods.
A hybrid method for the numerical prediction of enthalpy transport in fluid flow
Stevanovic, V.D.; Jovanovic, Z.L.
2000-01-01
The solution of the transient thermal energy transport by convection in forced fluid flow is a necessary step in thermal design, simulation and analyses regarding the operational conditions of various energy and chemical plants. It is of primary importance for the prediction of cooling or heating rate changes in the thermal equipment flow channels, as well as for the prediction of the boiling boundary location in a boiling channel. For instance, prediction of water enthalpy front propagation in complex pipeline networks of a district heating system is necessary for the setup of operating procedures, which will enable punctual and optimal heat supply to customers. Here, hybrid method (HM) is proposed for the prediction f convective enthalpy transport and the boiling boundary location. The momentum equation is solved with the SIMPLE procedure, while the energy equation is solved with the Method of Characteristics (MOC) and with the application of the Lagrange interpolation polynomial (LIP). The MOC method is applied on an equidistant grid. The initial values at the starting points of characteristic paths are calculated with the LIP. The accuracy of the method is demonstrated on tests of the transient boiling boundary prediction and the well known problem of a propagating discontinuity. The dependence of the hybrid method on the spatial integration step size and the degree of the LIP is analyzed. The LIP of the third degree is recommended. by which practically exact solutions with acceptable number of nodes are obtained.
NASA Astrophysics Data System (ADS)
Briones-Torres, J. A.; Pernas-Salomón, R.; Pérez-Álvarez, R.; Rodríguez-Vargas, I.
2016-05-01
Gapless bilayer graphene (GBG), like monolayer graphene, is a material system with unique properties, such as anti-Klein tunneling and intrinsic Fano resonances. These properties rely on the gapless parabolic dispersion relation and the chiral nature of bilayer graphene electrons. In addition, propagating and evanescent electron states coexist inherently in this material, giving rise to these exotic properties. In this sense, bilayer graphene is unique, since in most material systems in which Fano resonance phenomena are manifested an external source that provides extended states is required. However, from a numerical standpoint, the presence of evanescent-divergent states in the eigenfunctions linear superposition representing the Dirac spinors, leads to a numerical degradation (the so called Ωd problem) in the practical applications of the standard Coefficient Transfer Matrix (K) method used to study charge transport properties in Bilayer Graphene based multi-barrier systems. We present here a straightforward procedure based in the hybrid compliance-stiffness matrix method (H) that can overcome this numerical degradation. Our results show that in contrast to standard matrix method, the proposed H method is suitable to study the transmission and transport properties of electrons in GBG superlattice since it remains numerically stable regardless the size of the superlattice and the range of values taken by the input parameters: the energy and angle of the incident electrons, the barrier height and the thickness and number of barriers. We show that the matrix determinant can be used as a test of the numerical accuracy in real calculations.
NASA Astrophysics Data System (ADS)
Burago, N. G.; Nikitin, I. S.; Yakushev, V. L.
2016-06-01
Techniques that improve the accuracy of numerical solutions and reduce their computational costs are discussed as applied to continuum mechanics problems with complex time-varying geometry. The approach combines shock-capturing computations with the following methods: (1) overlapping meshes for specifying complex geometry; (2) elastic arbitrarily moving adaptive meshes for minimizing the approximation errors near shock waves, boundary layers, contact discontinuities, and moving boundaries; (3) matrix-free implementation of efficient iterative and explicit-implicit finite element schemes; (4) balancing viscosity (version of the stabilized Petrov-Galerkin method); (5) exponential adjustment of physical viscosity coefficients; and (6) stepwise correction of solutions for providing their monotonicity and conservativeness.
Malik, Suheel Abdullah; Qureshi, Ijaz Mansoor; Amir, Muhammad; Malik, Aqdas Naveed; Haq, Ihsanul
2015-01-01
In this paper, a new heuristic scheme for the approximate solution of the generalized Burgers'-Fisher equation is proposed. The scheme is based on the hybridization of Exp-function method with nature inspired algorithm. The given nonlinear partial differential equation (NPDE) through substitution is converted into a nonlinear ordinary differential equation (NODE). The travelling wave solution is approximated by the Exp-function method with unknown parameters. The unknown parameters are estimated by transforming the NODE into an equivalent global error minimization problem by using a fitness function. The popular genetic algorithm (GA) is used to solve the minimization problem, and to achieve the unknown parameters. The proposed scheme is successfully implemented to solve the generalized Burgers'-Fisher equation. The comparison of numerical results with the exact solutions, and the solutions obtained using some traditional methods, including adomian decomposition method (ADM), homotopy perturbation method (HPM), and optimal homotopy asymptotic method (OHAM), show that the suggested scheme is fairly accurate and viable for solving such problems. PMID:25811858
NASA Astrophysics Data System (ADS)
Sentyabov, A. V.; Gavrilov, A. A.; Dekterev, A. A.; Minakov, A. V.
2014-12-01
Numerical modeling of the unsteady flow in the draft tube of the test bench hydro turbine is conducted. The hybrid RANS-LES methods for modeling turbulent flows are compared. The intensity and frequency of pressure fluctuations, which are induced by the vortex core precession under the runner, and the integral characteristics are considered. An analysis of the synchronous and asynchronous parts of pressure fluctuations is done; the generating and influence of the synchronous component of fluctuations are considered. The vortex core interaction with the draft tube elbow is considered.
A hybrid numerical scheme for the numerical solution of the Burgers' equation
NASA Astrophysics Data System (ADS)
Jiwari, Ram
2015-03-01
In this article, a hybrid numerical scheme based on Euler implicit method, quasilinearization and uniform Haar wavelets has been developed for the numerical solutions of Burgers' equation. Most of the numerical methods available in the literature fail to capture the physical behavior of the equations when viscosity ν → 0. In Jiwari (2012), the author presented the numerical results up to ν = 0.003 and the scheme failed for values smaller than ν = 0.003. The main aim in the development of the present scheme is to overcome the drawback of the scheme developed in Jiwari (2012). Lastly, three test problems are chosen to check the accuracy of the proposed scheme. The approximated results are compared with existing numerical and exact solutions found in literature. The use of uniform Haar wavelet is found to be accurate, simple, fast, flexible, convenient and at small computation costs.
Numerically controlled oscillators with hybrid function generators.
Jainiszewski, Ireneusz; Hoppe, Bernhard; Meuth, Hermann
2002-07-01
Numerically controlled oscillators (NCOs), with a hybrid scheme of both look-up tables (LUT) and coordinate transformation digital computer (CORDIC) algorithms for a hardware efficient, high performance sine/cosine function generation are investigated. This scheme combines fast access and power efficiency of reasonably sized LUTs, and arbitrary precision obtainable from a rigorous iteration algorithm. Systematic studies using hardware description language (HDL) models and synthesis lead to optimum LUT/CORDIC ratios, which minimize power consumption and silicon area for a given operating clock frequency. First order error models are presented as guidelines for choosing internal NCO parameters. The NCO accuracy is tested with HDL simulations for all algorithmic states to limit output errors to 1 least significant bit (LSB) and by spectra derived from discrete Fourier transform (DFT) for typical frequency inputs f, resulting in a signal to noise ratio (SNR) of better than 100 dB for an amplitude word length AW of 16 Bit. Two benchmark designs were adopted for the two clock frequencies 200 MHz and 20 MHz, as "high" and "moderate" performance, respectively. The NCO models are synthesized in a 0.35 microm CMOS standard cell target technology and optimized to actually achieve after layout maximum clock frequencies exceeding 310 MHz, i.e., signal frequencies of up to 100 MHz. PMID:12152954
Hybrid codes: Methods and applications
Winske, D. ); Omidi, N. )
1991-01-01
In this chapter we discuss hybrid'' algorithms used in the study of low frequency electromagnetic phenomena, where one or more ion species are treated kinetically via standard PIC methods used in particle codes and the electrons are treated as a single charge neutralizing massless fluid. Other types of hybrid models are possible, as discussed in Winske and Quest, but hybrid codes with particle ions and massless fluid electrons have become the most common for simulating space plasma physics phenomena in the last decade, as we discuss in this paper.
Hybrid particles and associated methods
Fox, Robert V; Rodriguez, Rene; Pak, Joshua J; Sun, Chivin
2015-02-10
Hybrid particles that comprise a coating surrounding a chalcopyrite material, the coating comprising a metal, a semiconductive material, or a polymer; a core comprising a chalcopyrite material and a shell comprising a functionalized chalcopyrite material, the shell enveloping the core; or a reaction product of a chalcopyrite material and at least one of a reagent, heat, and radiation. Methods of forming the hybrid particles are also disclosed.
Numerical methods in structural mechanics
NASA Astrophysics Data System (ADS)
Obraztsov, I. F.
The papers contained in this volume focus on numerical, numerical-analytical, and theoretical methods for dealing with strength, stability, and dynamics problems in the design of the structural elements of flight vehicles. Topics discussed include the solution of homogeneous boundary value problems for systems of ordinary differential equations modified by a difference factorization method, a study of the rupture strength of a welded joint between plates, singular solutions in mixed problems for a wedge and a half-strip, and a thermoelasticity problem for an open-profile cylindrical shell with a localized temperature field.
Numerical simulations of plasma brush behavior in hybrid armatures
NASA Astrophysics Data System (ADS)
Hawke, R. S.; Pincosy, P. A.
1993-01-01
Hybrid armatures used to accelerate projectiles in railguns are often the consequence of using a solid armature or in some cases the preferred armature type. Although hybrid armatures are often used, their design has been empirical and their performance sporadic. As a first step towards understanding hybrid design and performance, we have begun a combined numerical simulation and experimental verification effort. This paper will describe numerical simulations performed with a quasi 1-D MHD code (CONFUSE) which has been applied to simulate the behavior of plasma brushes used in hybrid armatures. The simulations have provided estimates of the plasma brush length, resistive voltage drop and temperatures corresponding to a range of; 1) brush gap size, 2) fuse thickness, and 3) magnetic pressure. The results of these simulations is presented and discussed.
Numerical methods for turbulent flow
NASA Astrophysics Data System (ADS)
Turner, James C., Jr.
1988-09-01
It has generally become accepted that the Navier-Strokes equations predict the dynamic behavior of turbulent as well as laminar flows of a fluid at a point in space away form a discontinuity such as a shock wave. Turbulence is also closely related to the phenomena of non-uniqueness of solutions of the Navier-Strokes equations. These second order, nonlinear partial differential equations can be solved analytically for only a few simple flows. Turbulent flow fields are much to complex to lend themselves to these few analytical methods. Numerical methods, therefore, offer the only possibility of achieving a solution of turbulent flow equations. In spite of recent advances in computer technology, the direct solution, by discrete methods, of the Navier-Strokes equations for turbulent flow fields is today, and in the foreseeable future, impossible. Thus the only economically feasible way to solve practical turbulent flow problems numerically is to use statistically averaged equations governing mean-flow quantities. The objective is to study some recent developments relating to the use of numerical methods to study turbulent flow.
Numerical methods for molecular dynamics
Skeel, R.D.
1991-01-01
This report summarizes our research progress to date on the use of multigrid methods for three-dimensional elliptic partial differential equations, with particular emphasis on application to the Poisson-Boltzmann equation of molecular biophysics. This research is motivated by the need for fast and accurate numerical solution techniques for three-dimensional problems arising in physics and engineering. In many applications these problems must be solved repeatedly, and the extremely large number of discrete unknowns required to accurately approximate solutions to partial differential equations in three-dimensional regions necessitates the use of efficient solution methods. This situation makes clear the importance of developing methods which are of optimal order (or nearly so), meaning that the number of operations required to solve the discrete problem is on the order of the number of discrete unknowns. Multigrid methods are generally regarded as being in this class of methods, and are in fact provably optimal order for an increasingly large class of problems. The fundamental goal of this research is to develop a fast and accurate numerical technique, based on multi-level principles, for the solutions of the Poisson-Boltzmann equation of molecular biophysics and similar equations occurring in other applications. An outline of the report is as follows. We first present some background material, followed by a survey of the literature on the use of multigrid methods for solving problems similar to the Poisson-Boltzmann equation. A short description of the software we have developed so far is then given, and numerical results are discussed. Finally, our research plans for the coming year are presented.
Hybrid natural element method for viscoelasticity problems
NASA Astrophysics Data System (ADS)
Zhou, Yan-Kai; Ma, Yong-Qi; Dong, Yi; Feng, Wei
2015-01-01
A hybrid natural element method (HNEM) for two-dimensional viscoelasticity problems under the creep condition is proposed. The natural neighbor interpolation is used as the test function, and the discrete equation system of the HNEM for viscoelasticity problems is obtained using the Hellinger-Reissner variational principle. In contrast to the natural element method (NEM), the HNEM can directly obtain the nodal stresses, which have higher precisions than those obtained using the moving least-square (MLS) approximation. Some numerical examples are given to demonstrate the validity and superiority of this HNEM. Project supported by the Natural Science Foundation of Shanghai, China (Grant No.13ZR1415900).
Numerical solution of distributed order fractional differential equations by hybrid functions
NASA Astrophysics Data System (ADS)
Mashayekhi, S.; Razzaghi, M.
2016-06-01
In this paper, a new numerical method for solving the distributed fractional differential equations is presented. The method is based upon hybrid functions approximation. The properties of hybrid functions consisting of block-pulse functions and Bernoulli polynomials are presented. The Riemann-Liouville fractional integral operator for hybrid functions is introduced. This operator is then utilized to reduce the solution of the distributed fractional differential equations to a system of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique.
Numerical methods for multibody systems
NASA Technical Reports Server (NTRS)
Glowinski, Roland; Nasser, Mahmoud G.
1994-01-01
This article gives a brief summary of some results obtained by Nasser on modeling and simulation of inequality problems in multibody dynamics. In particular, the augmented Lagrangian method discussed here is applied to a constrained motion problem with impulsive inequality constraints. A fundamental characteristic of the multibody dynamics problem is the lack of global convexity of its Lagrangian. The problem is transformed into a convex analysis problem by localization (piecewise linearization), where the augmented Lagrangian has been successfully used. A model test problem is considered and a set of numerical experiments is presented.
NASA Astrophysics Data System (ADS)
Kuo, K. A.; Verbraken, H.; Degrande, G.; Lombaert, G.
2016-07-01
Along with the rapid expansion of urban rail networks comes the need for accurate predictions of railway induced vibration levels at grade and in buildings. Current computational methods for making predictions of railway induced ground vibration rely on simplifying modelling assumptions and require detailed parameter inputs, which lead to high levels of uncertainty. It is possible to mitigate against these issues using a combination of field measurements and state-of-the-art numerical methods, known as a hybrid model. In this paper, two hybrid models are developed, based on the use of separate source and propagation terms that are quantified using in situ measurements or modelling results. These models are implemented using term definitions proposed by the Federal Railroad Administration and assessed using the specific illustration of a surface railway. It is shown that the limitations of numerical and empirical methods can be addressed in a hybrid procedure without compromising prediction accuracy.
Hybrid power management system and method
NASA Technical Reports Server (NTRS)
Eichenberg, Dennis J. (Inventor)
2007-01-01
A system and method for hybrid power management. The system includes photovoltaic cells, ultracapacitors, and pulse generators. In one embodiment, the hybrid power management system is used to provide power for a highway safety flasher.
Hybrid Power Management System and Method
NASA Technical Reports Server (NTRS)
Eichenberg, Dennis J. (Inventor)
2008-01-01
A system and method for hybrid power management. The system includes photovoltaic cells, ultracapacitors, and pulse generators. In one embodiment, the hybrid power management system is used to provide power for a highway safety flasher.
Numerical methods used in fusion science numerical modeling
NASA Astrophysics Data System (ADS)
Yagi, M.
2015-04-01
The dynamics of burning plasma is very complicated physics, which is dominated by multi-scale and multi-physics phenomena. To understand such phenomena, numerical simulations are indispensable. Fundamentals of numerical methods used in fusion science numerical modeling are briefly discussed in this paper. In addition, the parallelization technique such as open multi processing (OpenMP) and message passing interface (MPI) parallel programing are introduced and the loop-level parallelization is shown as an example.
NASA Astrophysics Data System (ADS)
Terrana, S.; Vilotte, J. P.; Guillot, L.
2015-12-01
New seismological monitoring networks combine broadband seismic receivers, hydrophones and micro-barometers antenna, providing complementary observation of source-radiated waves. Exploiting these observations requires accurate and multi-media - elastic, hydro-acoustic, infrasound - wave simulation methods, in order to improve our physical understanding of energy exchanges at material interfaces.We present here a new development of a high-order Hybridized Discontinuous Galerkin (HDG) method, for the simulation of coupled seismic and acoustic wave propagation, within a unified framework ([1],[2]) allowing for continuous and discontinuous Spectral Element Methods (SEM) to be used in the same simulation, with conforming and non-conforming meshes. The HDG-SEM approximation leads to differential - algebraic equations, which can be solved implicitly using energy-preserving time-schemes.The proposed HDG-SEM is computationally attractive, when compared with classical Discontinuous Galerkin methods, involving only the approximation of the single-valued traces of the velocity field along the element interfaces as globally coupled unknowns. The formulation is based on a variational approximation of the physical fluxes, which are shown to be the explicit solution of an exact Riemann problem at each element boundaries. This leads to a highly parallel and efficient unstructured and high-order accurate method, which can be space-and-time adaptive.A numerical study of the accuracy and convergence of the HDG-SEM is performed through a number of case studies involving elastic-acoustic (infrasound) coupling with geometries of increasing complexity. Finally, the performance of the method is illustrated through realistic case studies involving ground wave propagation associated to topography effects.In conclusion, we outline some on-going extensions of the method.References:[1] Cockburn, B., Gopalakrishnan, J., Lazarov, R., Unified hybridization of discontinuous Galerkin, mixed and
A hybrid neurocomputing/numerical strategy for nonlinear structural analysis
NASA Technical Reports Server (NTRS)
Szewczyk, Z. Peter; Noor, Ahmed K.
1995-01-01
A hybrid neurocomputing/numerical strategy is presented for geometrically nonlinear analysis of structures. The strategy combines model-free data processing capabilities of computational neural networks with a Pade approximants-based perturbation technique to predict partial information about the nonlinear response of structures. In the hybrid strategy, multilayer feedforward neural networks are used to extend the validity of solutions by using training samples produced by Pade approximations to the Taylor series expansion of the response function. The range of validity of the training samples is taken to be the radius of convergence of Pade approximants and is estimated by setting a tolerance on the diverging approximants. The norm of residual vector of unbalanced forces in a given element is used as a measure to assess the quality of network predictions. To further increase the accuracy and the range of network predictions, additional training data are generated by either applying linear regression to weight matrices or expanding the training data by using predicted coefficients in a Taylor series. The effectiveness of the hybrid strategy is assessed by performing large-deflection analysis of a doubly-curved composite panel with a circular cutout, and postbuckling analyses of stiffened composite panels subjected to an in-plane edge shear load. In all the problems considered, the hybrid strategy is used to predict selective information about the structural response, namely the total strain energy and the maximum displacement components only.
High numerical aperture hybrid optics for two-photon polymerization.
Burmeister, Frank; Zeitner, Uwe D; Nolte, Stefan; Tünnermann, Andreas
2012-03-26
We report on an immersion hybrid optics specially designed for focusing ultrashort laser pulses into a polymer for direct laser writing via two-photon polymerization. The hybrid optics allows for well-corrected focusing over a large working distance range of 577 μm with a numerical aperture (NA) of 1.33 and low internal dispersion. We combine the concepts of an aplanatic solid immersion lens (ASIL) for achieving a high NA with a diffractive optical element (DOE) for correction of aberrations. To demonstrate the improvements for volume structuring of the polymer, we compare the achievable structure sizes of our optics with a commercially available oil-immersion objective (100x, NA=1.4). PMID:22453471
On Hybrid and mixed finite element methods
NASA Technical Reports Server (NTRS)
Pian, T. H. H.
1981-01-01
Three versions of the assumed stress hybrid model in finite element methods and the corresponding variational principles for the formulation are presented. Examples of rank deficiency for stiffness matrices by the hybrid stress model are given and their corresponding kinematic deformation modes are identified. A discussion of the derivation of general semi-Loof elements for plates and shells by the hybrid stress method is given. It is shown that the equilibrium model by Fraeijs de Veubeke can be derived by the approach of the hybrid stress model as a special case of semi-Loof elements.
A fast multipole hybrid boundary node method for composite materials
NASA Astrophysics Data System (ADS)
Wang, Qiao; Miao, Yu; Zhu, Hongping
2013-06-01
This article presents a multi-domain fast multipole hybrid boundary node method for composite materials in 3D elasticity. The hybrid boundary node method (hybrid BNM) is a meshless method which only requires nodes constructed on the surface of a domain. The method is applied to 3D simulation of composite materials by a multi-domain solver and accelerated by the fast multipole method (FMM) in this paper. The preconditioned GMRES is employed to solve the final system equation and precondition techniques are discussed. The matrix-vector multiplication in each iteration is divided into smaller scale ones at the sub-domain level and then accelerated by FMM within individual sub-domains. The computed matrix-vector products at the sub-domain level are then combined according to the continuity conditions on the interfaces. The algorithm is implemented on a computer code written in C + +. Numerical results show that the technique is accurate and efficient.
Numerical Methods of Computational Electromagnetics for Complex Inhomogeneous Systems
Cai, Wei
2014-05-15
Understanding electromagnetic phenomena is the key in many scientific investigation and engineering designs such as solar cell designs, studying biological ion channels for diseases, and creating clean fusion energies, among other things. The objectives of the project are to develop high order numerical methods to simulate evanescent electromagnetic waves occurring in plasmon solar cells and biological ion-channels, where local field enhancement within random media in the former and long range electrostatic interactions in the latter are of major challenges for accurate and efficient numerical computations. We have accomplished these objectives by developing high order numerical methods for solving Maxwell equations such as high order finite element basis for discontinuous Galerkin methods, well-conditioned Nedelec edge element method, divergence free finite element basis for MHD, and fast integral equation methods for layered media. These methods can be used to model the complex local field enhancement in plasmon solar cells. On the other hand, to treat long range electrostatic interaction in ion channels, we have developed image charge based method for a hybrid model in combining atomistic electrostatics and continuum Poisson-Boltzmann electrostatics. Such a hybrid model will speed up the molecular dynamics simulation of transport in biological ion-channels.
Efficient hybrid-symbolic methods for quantum mechanical calculations
NASA Astrophysics Data System (ADS)
Scott, T. C.; Zhang, Wenxing
2015-06-01
We present hybrid symbolic-numerical tools to generate optimized numerical code for rapid prototyping and fast numerical computation starting from a computer algebra system (CAS) and tailored to any given quantum mechanical problem. Although a major focus concerns the quantum chemistry methods of H. Nakatsuji which has yielded successful and very accurate eigensolutions for small atoms and molecules, the tools are general and may be applied to any basis set calculation with a variational principle applied to its linear and non-linear parameters.
Hybrid methods for cybersecurity analysis :
Davis, Warren Leon,; Dunlavy, Daniel M.
2014-01-01
Early 2010 saw a signi cant change in adversarial techniques aimed at network intrusion: a shift from malware delivered via email attachments toward the use of hidden, embedded hyperlinks to initiate sequences of downloads and interactions with web sites and network servers containing malicious software. Enterprise security groups were well poised and experienced in defending the former attacks, but the new types of attacks were larger in number, more challenging to detect, dynamic in nature, and required the development of new technologies and analytic capabilities. The Hybrid LDRD project was aimed at delivering new capabilities in large-scale data modeling and analysis to enterprise security operators and analysts and understanding the challenges of detection and prevention of emerging cybersecurity threats. Leveraging previous LDRD research e orts and capabilities in large-scale relational data analysis, large-scale discrete data analysis and visualization, and streaming data analysis, new modeling and analysis capabilities were quickly brought to bear on the problems in email phishing and spear phishing attacks in the Sandia enterprise security operational groups at the onset of the Hybrid project. As part of this project, a software development and deployment framework was created within the security analyst work ow tool sets to facilitate the delivery and testing of new capabilities as they became available, and machine learning algorithms were developed to address the challenge of dynamic threats. Furthermore, researchers from the Hybrid project were embedded in the security analyst groups for almost a full year, engaged in daily operational activities and routines, creating an atmosphere of trust and collaboration between the researchers and security personnel. The Hybrid project has altered the way that research ideas can be incorporated into the production environments of Sandias enterprise security groups, reducing time to deployment from months and
NASA Astrophysics Data System (ADS)
Li, Chengen; Cai, Guobiao; Tian, Hui
2016-06-01
This paper is aimed to analyse the combustion characteristics of hybrid rocket motor with multi-section swirl injection by simulating the combustion flow field. Numerical combustion flow field and combustion performance parameters are obtained through three-dimensional numerical simulations based on a steady numerical model proposed in this paper. The hybrid rocket motor adopts 98% hydrogen peroxide and polyethylene as the propellants. Multiple injection sections are set along the axis of the solid fuel grain, and the oxidizer enters the combustion chamber by means of tangential injection via the injector ports in the injection sections. Simulation results indicate that the combustion flow field structure of the hybrid rocket motor could be improved by multi-section swirl injection method. The transformation of the combustion flow field can greatly increase the fuel regression rate and the combustion efficiency. The average fuel regression rate of the motor with multi-section swirl injection is improved by 8.37 times compared with that of the motor with conventional head-end irrotational injection. The combustion efficiency is increased to 95.73%. Besides, the simulation results also indicate that (1) the additional injection sections can increase the fuel regression rate and the combustion efficiency; (2) the upstream offset of the injection sections reduces the combustion efficiency; and (3) the fuel regression rate and the combustion efficiency decrease with the reduction of the number of injector ports in each injection section.
A hybrid method with deviational particles for spatial inhomogeneous plasma
NASA Astrophysics Data System (ADS)
Yan, Bokai
2016-03-01
In this work we propose a Hybrid method with Deviational Particles (HDP) for a plasma modeled by the inhomogeneous Vlasov-Poisson-Landau system. We split the distribution into a Maxwellian part evolved by a grid based fluid solver and a deviation part simulated by numerical particles. These particles, named deviational particles, could be both positive and negative. We combine the Monte Carlo method proposed in [31], a Particle in Cell method and a Macro-Micro decomposition method [3] to design an efficient hybrid method. Furthermore, coarse particles are employed to accelerate the simulation. A particle resampling technique on both deviational particles and coarse particles is also investigated and improved. This method is applicable in all regimes and significantly more efficient compared to a PIC-DSMC method near the fluid regime.
A hybrid artificial bee colony algorithm for numerical function optimization
NASA Astrophysics Data System (ADS)
Alqattan, Zakaria N.; Abdullah, Rosni
2015-02-01
Artificial Bee Colony (ABC) algorithm is one of the swarm intelligence algorithms; it has been introduced by Karaboga in 2005. It is a meta-heuristic optimization search algorithm inspired from the intelligent foraging behavior of the honey bees in nature. Its unique search process made it as one of the most competitive algorithm with some other search algorithms in the area of optimization, such as Genetic algorithm (GA) and Particle Swarm Optimization (PSO). However, the ABC performance of the local search process and the bee movement or the solution improvement equation still has some weaknesses. The ABC is good in avoiding trapping at the local optimum but it spends its time searching around unpromising random selected solutions. Inspired by the PSO, we propose a Hybrid Particle-movement ABC algorithm called HPABC, which adapts the particle movement process to improve the exploration of the original ABC algorithm. Numerical benchmark functions were used in order to experimentally test the HPABC algorithm. The results illustrate that the HPABC algorithm can outperform the ABC algorithm in most of the experiments (75% better in accuracy and over 3 times faster).
Numerical Methods for Radiation Magnetohydrodynamics in Astrophysics
Klein, R I; Stone, J M
2007-11-20
We describe numerical methods for solving the equations of radiation magnetohydrodynamics (MHD) for astrophysical fluid flow. Such methods are essential for the investigation of the time-dependent and multidimensional dynamics of a variety of astrophysical systems, although our particular interest is motivated by problems in star formation. Over the past few years, the authors have been members of two parallel code development efforts, and this review reflects that organization. In particular, we discuss numerical methods for MHD as implemented in the Athena code, and numerical methods for radiation hydrodynamics as implemented in the Orion code. We discuss the challenges introduced by the use of adaptive mesh refinement in both codes, as well as the most promising directions for future developments.
A hybrid Godunov method for radiation hydrodynamics
NASA Astrophysics Data System (ADS)
Sekora, Michael D.; Stone, James M.
2010-09-01
From a mathematical perspective, radiation hydrodynamics can be thought of as a system of hyperbolic balance laws with dual multiscale behavior (multiscale behavior associated with the hyperbolic wave speeds as well as multiscale behavior associated with source term relaxation). With this outlook in mind, this paper presents a hybrid Godunov method for one-dimensional radiation hydrodynamics that is uniformly well behaved from the photon free streaming (hyperbolic) limit through the weak equilibrium diffusion (parabolic) limit and to the strong equilibrium diffusion (hyperbolic) limit. Moreover, one finds that the technique preserves certain asymptotic limits. The method incorporates a backward Euler upwinding scheme for the radiation energy density Er and flux Fr as well as a modified Godunov scheme for the material density ρ, momentum density m, and energy density E. The backward Euler upwinding scheme is first-order accurate and uses an implicit HLLE flux function to temporally advance the radiation components according to the material flow scale. The modified Godunov scheme is second-order accurate and directly couples stiff source term effects to the hyperbolic structure of the system of balance laws. This Godunov technique is composed of a predictor step that is based on Duhamel's principle and a corrector step that is based on Picard iteration. The Godunov scheme is explicit on the material flow scale but is unsplit and fully couples matter and radiation without invoking a diffusion-type approximation for radiation hydrodynamics. This technique derives from earlier work by Miniati and Colella (2007) [41]. Numerical tests demonstrate that the method is stable, robust, and accurate across various parameter regimes.
On Numerical Methods For Hypersonic Turbulent Flows
NASA Astrophysics Data System (ADS)
Yee, H. C.; Sjogreen, B.; Shu, C. W.; Wang, W.; Magin, T.; Hadjadj, A.
2011-05-01
Proper control of numerical dissipation in numerical methods beyond the standard shock-capturing dissipation at discontinuities is an essential element for accurate and stable simulation of hypersonic turbulent flows, including combustion, and thermal and chemical nonequilibrium flows. Unlike rapidly developing shock interaction flows, turbulence computations involve long time integrations. Improper control of numerical dissipation from one time step to another would be compounded over time, resulting in the smearing of turbulent fluctuations to an unrecognizable form. Hypersonic turbulent flows around re- entry space vehicles involve mixed steady strong shocks and turbulence with unsteady shocklets that pose added computational challenges. Stiffness of the source terms and material mixing in combustion pose yet other types of numerical challenges. A low dissipative high order well- balanced scheme, which can preserve certain non-trivial steady solutions of the governing equations exactly, may help minimize some of these difficulties. For stiff reactions it is well known that the wrong propagation speed of discontinuities occurs due to the under-resolved numerical solutions in both space and time. Schemes to improve the wrong propagation speed of discontinuities for systems of stiff reacting flows remain a challenge for algorithm development. Some of the recent algorithm developments for direct numerical simulations (DNS) and large eddy simulations (LES) for the subject physics, including the aforementioned numerical challenges, will be discussed.
A Hybrid Method for Abstracting Newspaper Articles.
ERIC Educational Resources Information Center
Liu, James; Wu, Yan; Zhou, Lina
1999-01-01
Introduces a hybrid method for abstracting Chinese text that integrates the statistical approach with language understandings, incorporating some linguistics heuristics and segmentation into the abstracting process. Initial responses from application to Chinese newspaper articles show that the method contributes much to the flexibility and…
Hybrid immersed interface-immersed boundary methods for AC dielectrophoresis
Hossan, Mohammad Robiul; Dillon, Robert; Dutta, Prashanta
2014-08-01
Dielectrophoresis, a nonlinear electrokinetic transport mechanism, has become popular in many engineering applications including manipulation, characterization and actuation of biomaterials, particles and biological cells. In this paper, we present a hybrid immersed interface–immersed boundary method to study AC dielectrophoresis where an algorithm is developed to solve the complex Poisson equation using a real variable formulation. An immersed interface method is employed to obtain the AC electric field in a fluid media with suspended particles and an immersed boundary method is used for the fluid equations and particle transport. The convergence of the proposed algorithm as well as validation of the hybrid scheme with experimental results is presented. In this paper, the Maxwell stress tensor is used to calculate the dielectrophoretic force acting on particles by considering the physical effect of particles in the computational domain. Thus, this study eliminates the approximations used in point dipole methods for calculating dielectrophoretic force. A comparative study between Maxwell stress tensor and point dipole methods for computing dielectrophoretic forces are presented. The hybrid method is used to investigate the physics of dielectrophoresis in microfluidic devices using an AC electric field. The numerical results show that with proper design and appropriate selection of applied potential and frequency, global electric field minima can be obtained to facilitate multiple particle trapping by exploiting the mechanism of negative dielectrophoresis. Our numerical results also show that electrically neutral particles form a chain parallel to the applied electric field irrespective of their initial orientation when an AC electric field is applied. This proposed hybrid numerical scheme will help to better understand dielectrophoresis and to design and optimize microfluidic devices.
Hybrid immersed interface-immersed boundary methods for AC dielectrophoresis
NASA Astrophysics Data System (ADS)
Hossan, Mohammad Robiul; Dillon, Robert; Dutta, Prashanta
2014-08-01
Dielectrophoresis, a nonlinear electrokinetic transport mechanism, has become popular in many engineering applications including manipulation, characterization and actuation of biomaterials, particles and biological cells. In this paper, we present a hybrid immersed interface-immersed boundary method to study AC dielectrophoresis where an algorithm is developed to solve the complex Poisson equation using a real variable formulation. An immersed interface method is employed to obtain the AC electric field in a fluid media with suspended particles and an immersed boundary method is used for the fluid equations and particle transport. The convergence of the proposed algorithm as well as validation of the hybrid scheme with experimental results is presented. In this paper, the Maxwell stress tensor is used to calculate the dielectrophoretic force acting on particles by considering the physical effect of particles in the computational domain. Thus, this study eliminates the approximations used in point dipole methods for calculating dielectrophoretic force. A comparative study between Maxwell stress tensor and point dipole methods for computing dielectrophoretic forces are presented. The hybrid method is used to investigate the physics of dielectrophoresis in microfluidic devices using an AC electric field. The numerical results show that with proper design and appropriate selection of applied potential and frequency, global electric field minima can be obtained to facilitate multiple particle trapping by exploiting the mechanism of negative dielectrophoresis. Our numerical results also show that electrically neutral particles form a chain parallel to the applied electric field irrespective of their initial orientation when an AC electric field is applied. This proposed hybrid numerical scheme will help to better understand dielectrophoresis and to design and optimize microfluidic devices.
A multigrid solution method for mixed hybrid finite elements
Schmid, W.
1996-12-31
We consider the multigrid solution of linear equations arising within the discretization of elliptic second order boundary value problems of the form by mixed hybrid finite elements. Using the equivalence of mixed hybrid finite elements and non-conforming nodal finite elements, we construct a multigrid scheme for the corresponding non-conforming finite elements, and, by this equivalence, for the mixed hybrid finite elements, following guidelines from Arbogast/Chen. For a rectangular triangulation of the computational domain, this non-conforming schemes are the so-called nodal finite elements. We explicitly construct prolongation and restriction operators for this type of non-conforming finite elements. We discuss the use of plain multigrid and the multilevel-preconditioned cg-method and compare their efficiency in numerical tests.
Fast numerical algorithms for fitting multiresolution hybrid shape models to brain MRI.
Vemuri, B C; Guo, Y; Lai, S H; Leonard, C M
1997-09-01
In this paper, we present new and fast numerical algorithms for shape recovery from brain MRI using multiresolution hybrid shape models. In this modeling framework, shapes are represented by a core rigid shape characterized by a superquadric function and a superimposed displacement function which is characterized by a membrane spline discretized using the finite-element method. Fitting the model to brain MRI data is cast as an energy minimization problem which is solved numerically. We present three new computational methods for model fitting to data. These methods involve novel mathematical derivations that lead to efficient numerical solutions of the model fitting problem. The first method involves using the nonlinear conjugate gradient technique with a diagonal Hessian preconditioner. The second method involves the nonlinear conjugate gradient in the outer loop for solving global parameters of the model and a preconditioned conjugate gradient scheme for solving the local parameters of the model. The third method involves the nonlinear conjugate gradient in the outer loop for solving the global parameters and a combination of the Schur complement formula and the alternating direction-implicit method for solving the local parameters of the model. We demonstrate the efficiency of our model fitting methods via experiments on several MR brain scans. PMID:9873915
Ostermann, Lars; Seidel, Christian
2015-03-10
The numerical analysis of hydro power stations is an important method of the hydraulic design and is used for the development and optimisation of hydro power stations in addition to the experiments with the physical submodel of a full model in the hydraulic laboratory. For the numerical analysis, 2D and 3D models are appropriate and commonly used.The 2D models refer mainly to the shallow water equations (SWE), since for this flow model a large experience on a wide field of applications for the flow analysis of numerous problems in hydraulic engineering already exists. Often, the flow model is verified by in situ measurements. In order to consider 3D flow phenomena close to singularities like weirs, hydro power stations etc. the development of a hybrid fluid model is advantageous to improve the quality and significance of the global model. Here, an extended hybrid flow model based on the principle of the SWE is presented. The hybrid flow model directly links the numerical model with the experimental data, which may originate from physical full models, physical submodels and in-situ measurements. Hence a wide field of application of the hybrid model emerges including the improvement of numerical models and the strong coupling of numerical and experimental analysis.
A numerical method of detecting singularity
NASA Technical Reports Server (NTRS)
Laporte, M.; Vignes, J.
1978-01-01
A numerical method is reported which determines a value C for the degree of conditioning of a matrix. This value is C = 0 for a singular matrix and has progressively larger values for matrices which are increasingly well-conditioned. This value is C sub = C max sub max (C defined by the precision of the computer) when the matrix is perfectly well conditioned.
A numerical method for predicting hypersonic flowfields
NASA Technical Reports Server (NTRS)
Maccormack, Robert W.; Candler, Graham V.
1989-01-01
The flow about a body traveling at hypersonic speed is energetic enough to cause the atmospheric gases to chemically react and reach states in thermal nonequilibrium. The prediction of hypersonic flowfields requires a numerical method capable of solving the conservation equations of fluid flow, the chemical rate equations for specie formation and dissociation, and the transfer of energy relations between translational and vibrational temperature states. Because the number of equations to be solved is large, the numerical method should also be as efficient as possible. The proposed paper presents a fully implicit method that fully couples the solution of the fluid flow equations with the gas physics and chemistry relations. The method flux splits the inviscid flow terms, central differences of the viscous terms, preserves element conservation in the strong chemistry source terms, and solves the resulting block matrix equation by Gauss Seidel line relaxation.
The hybrid Eulerian Lagrangian numerical scheme tested with Chemistry
NASA Astrophysics Data System (ADS)
Hansen, A. B.; Sørensen, B.; Tarning-Andersen, P.; Christensen, J. H.; Brandt, J.; Kaas, E.
2012-12-01
A newly developed transport scheme, the Hybrid Eulerian Lagrangian (HEL) scheme, has been tested using a module for atmospheric chemistry, including 58 chemical species, and compared to two other traditional advection schemes; a classical pseudospectral Eulerian method the Accurate Space Derivative (ASD) scheme and the bi-cubic semi-Lagrangian (SL) scheme using classical rotation tests. The rotation tests have been designed to test and compare the advection schemes for different spatial and temporal resolutions in different chemical conditions (rural and urban) and for different shapes (cone and slotted cylinder). This gives the advection schemes different challenges with respect to relatively slow or fast chemistry and smooth or sharp gradients. In every test, error measures have been calculated and used for ranking of the advection schemes with respect to performance, i.e. lowest overall errors for all chemical species. The results presented show that the new transport scheme, HEL, by far outperforms both the Eulerian and semi-Lagrangian schemes with very low error estimates compared to the two other schemes.
The hybrid Eulerian Lagrangian numerical scheme tested with Chemistry
NASA Astrophysics Data System (ADS)
Hansen, A. B.; Sørensen, B.; Tarning-Andersen, P.; Christensen, J. H.; Brandt, J.; Kaas, E.
2012-11-01
A newly developed advection scheme, the Hybrid Eulerian Lagrangian (HEL) scheme, has been tested, including a module for atmospheric chemistry, including 58 chemical species, and compared to two other traditional advection schemes; a classical pseudospectral Eulerian method the Accurate Space Derivative (ASD) scheme and the bi-cubic semi-Lagrangian (SL) scheme using classical rotation tests. The rotation tests have been designed to test and compare the advection schemes for different spatial and temporal resolutions in different chemical conditions (rural and urban) and for different shapes (cone and slotted cylinder) giving the advection schemes different challenges with respect to relatively slow or fast chemistry and smooth or sharp gradients, respectively. In every test, error measures have been calculated and used for ranking of the advection schemes with respect to performance, i.e. lowest overall errors for all chemical species. Furthermore, the HEL and SL schemes have been compared in a shallow water model, demonstrating the performance in a more realistic non-linear deformation flow. The results in this paper show that the new advection scheme, HEL, by far outperforms both the Eulerian and semi-Lagrangian schemes with very low error estimates compared to the two other schemes. Although no analytic solution can be obtained for the performance in the non-linear shallow water model flow, the tracer distribution appears realistic as compared to LMCSL when a mixing between local parcel concentrations is introduced in HEL.
Numerical Analysis of the Symmetric Methods
NASA Astrophysics Data System (ADS)
Xu, Ji-Hong; Zhang, A.-Li
1995-03-01
Aimed at the initial value problem of the particular second-order ordinary differential equations,y ″=f(x, y), the symmetric methods (Quinlan and Tremaine, 1990) and our methods (Xu and Zhang, 1994) have been compared in detail by integrating the artificial earth satellite orbits in this paper. In the end, we point out clearly that the integral accuracy of numerical integration of the satellite orbits by applying our methods is obviously higher than that by applying the same order formula of the symmetric methods when the integration time-interval is not greater than 12000 periods.
A hybrid numerical fluid dynamics code for resistive magnetohydrodynamics
Energy Science and Technology Software Center (ESTSC)
2006-04-01
Spasmos is a computational fluid dynamics code that uses two numerical methods to solve the equations of resistive magnetohydrodynamic (MHD) flows in compressible, inviscid, conducting media[1]. The code is implemented as a set of libraries for the Python programming language[2]. It represents conducting and non-conducting gases and materials with uncomplicated (analytic) equations of state. It supports calculations in 1D, 2D, and 3D geometry, though only the 1D configuation has received significant testing to date. Becausemore » it uses the Python interpreter as a front end, users can easily write test programs to model systems with a variety of different numerical and physical parameters. Currently, the code includes 1D test programs for hydrodynamics (linear acoustic waves, the Sod weak shock[3], the Noh strong shock[4], the Sedov explosion[5], magnetic diffusion (decay of a magnetic pulse[6], a driven oscillatory "wine-cellar" problem[7], magnetic equilibrium), and magnetohydrodynamics (an advected magnetic pulse[8], linear MHD waves, a magnetized shock tube[9]). Spasmos current runs only in a serial configuration. In the future, it will use MPI for parallel computation.« less
Numerical Methods for Stochastic Partial Differential Equations
Sharp, D.H.; Habib, S.; Mineev, M.B.
1999-07-08
This is the final report of a Laboratory Directed Research and Development (LDRD) project at the Los Alamos National laboratory (LANL). The objectives of this proposal were (1) the development of methods for understanding and control of spacetime discretization errors in nonlinear stochastic partial differential equations, and (2) the development of new and improved practical numerical methods for the solutions of these equations. The authors have succeeded in establishing two methods for error control: the functional Fokker-Planck equation for calculating the time discretization error and the transfer integral method for calculating the spatial discretization error. In addition they have developed a new second-order stochastic algorithm for multiplicative noise applicable to the case of colored noises, and which requires only a single random sequence generation per time step. All of these results have been verified via high-resolution numerical simulations and have been successfully applied to physical test cases. They have also made substantial progress on a longstanding problem in the dynamics of unstable fluid interfaces in porous media. This work has lead to highly accurate quasi-analytic solutions of idealized versions of this problem. These may be of use in benchmarking numerical solutions of the full stochastic PDEs that govern real-world problems.
High order hybrid numerical simulations of two dimensional detonation waves
NASA Technical Reports Server (NTRS)
Cai, Wei
1993-01-01
In order to study multi-dimensional unstable detonation waves, a high order numerical scheme suitable for calculating the detailed transverse wave structures of multidimensional detonation waves was developed. The numerical algorithm uses a multi-domain approach so different numerical techniques can be applied for different components of detonation waves. The detonation waves are assumed to undergo an irreversible, unimolecular reaction A yields B. Several cases of unstable two dimensional detonation waves are simulated and detailed transverse wave interactions are documented. The numerical results show the importance of resolving the detonation front without excessive numerical viscosity in order to obtain the correct cellular patterns.
Numerical methods for supersonic astrophysical jets
NASA Astrophysics Data System (ADS)
Ha, Youngsoo
2003-09-01
The Euler equations of gas dynamics are used for the simulation of general astrophysical fluid flows including high Mach number astrophysical jets with radiative cooling. To accurately compute supersonic jet solutions with sharp resolution of shock waves, three modern numerical methods for gas dynamics were used: (1)a second-order Godunov method in LeVeque's software package CLAWPACK, (2)the Nessyahu-Tadmor-Kurganov (NTK) central hyperbolic scheme, and (3)the WENO-LF (Weighted Essentially Non-Oscillatory Lax-Friedrichs) scheme. Then simulations of supersonic astrophysical jets were compared, first without and then with radiative cooling. CLAWPACK consists of routines for solving time-dependent nonlinear hyperbolic conservation laws based on higher order Godunov methods and approximate Riemann problem solutions; the NTK scheme solves conservation laws using a modified Lax-Friedrichs central difference method without appealing to Riemann problem solutions; and the WENO-LF finite difference scheme is based on the Essentially Non-Oscillatory (ENO) idea by using Lax- Friedrichs flux splitting. The ENO method constructs a solution using the smoothness of the interpolating polynomial on given stencils; on the other hand, the WENO scheme uses a convex combination of the interpolate functions on all candidate stencils. The third-order and fifth-order WENO-LF methods were used to simulate the high Mach number jets. Appropriate numerical methods for incorporating radiative cooling in these numerical methods are also discussed. Interactions of supersonic jets with their environments (jet-“blob” interactions) are shown after modifying the codes to handle high Mach numbers and radiative cooling.
Jahantigh, Nabi; Keshavarz, Ali; Mirzaei, Masoud
2015-01-01
The aim of this study is to determine optimum hybrid heating systems parameters, such as temperature, surface area of a radiant heater and vent area to have thermal comfort conditions. DOE, Factorial design method is used to determine the optimum values for input parameters. A 3D model of a virtual standing thermal manikin with real dimensions is considered in this study. Continuity, momentum, energy, species equations for turbulent flow and physiological equation for thermal comfort are numerically solved to study heat, moisture and flow field. K - ɛRNG Model is used for turbulence modeling and DO method is used for radiation effects. Numerical results have a good agreement with the experimental data reported in the literature. The effect of various combinations of inlet parameters on thermal comfort is considered. According to Pareto graph, some of these combinations that have significant effect on the thermal comfort require no more energy can be used as useful tools. A better symmetrical velocity distribution around the manikin is also presented in the hybrid system. PMID:26052442
Numerical optimization of actuator trajectories for ITER hybrid scenario profile evolution
NASA Astrophysics Data System (ADS)
van Dongen, J.; Felici, F.; Hogeweij, G. M. D.; Geelen, P.; Maljaars, E.
2014-12-01
Optimal actuator trajectories for an ITER hybrid scenario ramp-up are computed using a numerical optimization method. For both L-mode and H-mode scenarios, the time trajectory of plasma current, EC heating and current drive distribution is determined that minimizes a chosen cost function, while satisfying constraints. The cost function is formulated to reflect two desired properties of the plasma q profile at the end of the ramp-up. The first objective is to maximize the ITG turbulence threshold by maximizing the volume-averaged s/q ratio. The second objective is to achieve a stationary q profile by having a flat loop voltage profile. Actuator and physics-derived constraints are included, imposing limits on plasma current, ramp rates, internal inductance and q profile. This numerical method uses the fast control-oriented plasma profile evolution code RAPTOR, which is successfully benchmarked against more complete CRONOS simulations for L-mode and H-mode mode ITER hybrid scenarios. It is shown that the optimized trajectories computed using RAPTOR also result in an improved ramp-up scenario for CRONOS simulations using the same input trajectories. Furthermore, the optimal trajectories are shown to vary depending on the precise timing of the L-H transition.
Hybrid Lanczos-type product methods
Ressel, K.J.
1996-12-31
A general framework is proposed to construct hybrid iterative methods for the solution of large nonsymmetric systems of linear equations. This framework is based on Lanczos-type product methods, whose iteration polynomial consists of the Lanczos polynomial multiplied by some other arbitrary, {open_quotes}shadow{close_quotes} polynomial. By using for the shadow polynomial Chebyshev (more general Faber) polynomials or L{sup 2}-optimal polynomials, hybrid (Chebyshev-like) methods are incorporated into Lanczos-type product methods. In addition, to acquire spectral information on the system matrix, which is required for such a choice of shadow polynomials, the Lanczos-process can be employed either directly or in an QMR-like approach. The QMR like approach allows the cheap computation of the roots of the B-orthogonal polynomials and the residual polynomials associated with the QMR iteration. These roots can be used as a good approximation for the spectrum of the system matrix. Different choices for the shadow polynomials and their construction are analyzed. The resulting hybrid methods are compared with standard Lanczos-type product methods, like BiOStab, BiOStab({ell}) and BiOS.
Hyperbolic conservation laws and numerical methods
NASA Technical Reports Server (NTRS)
Leveque, Randall J.
1990-01-01
The mathematical structure of hyperbolic systems and the scalar equation case of conservation laws are discussed. Linear, nonlinear systems and the Riemann problem for the Euler equations are also studied. The numerical methods for conservation laws are presented in a nonstandard manner which leads to large time steps generalizations and computations on irregular grids. The solution of conservation laws with stiff source terms is examined.
Hybrid method for the chemical master equation
Hellander, Andreas Loetstedt, Per
2007-11-10
The chemical master equation is solved by a hybrid method coupling a macroscopic, deterministic description with a mesoscopic, stochastic model. The molecular species are divided into one subset where the expected values of the number of molecules are computed and one subset with species with a stochastic variation in the number of molecules. The macroscopic equations resemble the reaction rate equations and the probability distribution for the stochastic variables satisfy a master equation. The probability distribution is obtained by the Stochastic Simulation Algorithm due to Gillespie. The equations are coupled via a summation over the mesoscale variables. This summation is approximated by Quasi-Monte Carlo methods. The error in the approximations is analyzed. The hybrid method is applied to three chemical systems from molecular cell biology.
Hybrid least squares multivariate spectral analysis methods
Haaland, David M.
2004-03-23
A set of hybrid least squares multivariate spectral analysis methods in which spectral shapes of components or effects not present in the original calibration step are added in a following prediction or calibration step to improve the accuracy of the estimation of the amount of the original components in the sampled mixture. The hybrid method herein means a combination of an initial calibration step with subsequent analysis by an inverse multivariate analysis method. A spectral shape herein means normally the spectral shape of a non-calibrated chemical component in the sample mixture but can also mean the spectral shapes of other sources of spectral variation, including temperature drift, shifts between spectrometers, spectrometer drift, etc. The shape can be continuous, discontinuous, or even discrete points illustrative of the particular effect.
Hybrid least squares multivariate spectral analysis methods
Haaland, David M.
2002-01-01
A set of hybrid least squares multivariate spectral analysis methods in which spectral shapes of components or effects not present in the original calibration step are added in a following estimation or calibration step to improve the accuracy of the estimation of the amount of the original components in the sampled mixture. The "hybrid" method herein means a combination of an initial classical least squares analysis calibration step with subsequent analysis by an inverse multivariate analysis method. A "spectral shape" herein means normally the spectral shape of a non-calibrated chemical component in the sample mixture but can also mean the spectral shapes of other sources of spectral variation, including temperature drift, shifts between spectrometers, spectrometer drift, etc. The "shape" can be continuous, discontinuous, or even discrete points illustrative of the particular effect.
Numerical methods for molecular dynamics. Progress report
Skeel, R.D.
1991-12-31
This report summarizes our research progress to date on the use of multigrid methods for three-dimensional elliptic partial differential equations, with particular emphasis on application to the Poisson-Boltzmann equation of molecular biophysics. This research is motivated by the need for fast and accurate numerical solution techniques for three-dimensional problems arising in physics and engineering. In many applications these problems must be solved repeatedly, and the extremely large number of discrete unknowns required to accurately approximate solutions to partial differential equations in three-dimensional regions necessitates the use of efficient solution methods. This situation makes clear the importance of developing methods which are of optimal order (or nearly so), meaning that the number of operations required to solve the discrete problem is on the order of the number of discrete unknowns. Multigrid methods are generally regarded as being in this class of methods, and are in fact provably optimal order for an increasingly large class of problems. The fundamental goal of this research is to develop a fast and accurate numerical technique, based on multi-level principles, for the solutions of the Poisson-Boltzmann equation of molecular biophysics and similar equations occurring in other applications. An outline of the report is as follows. We first present some background material, followed by a survey of the literature on the use of multigrid methods for solving problems similar to the Poisson-Boltzmann equation. A short description of the software we have developed so far is then given, and numerical results are discussed. Finally, our research plans for the coming year are presented.
RELAP-7 Numerical Stabilization: Entropy Viscosity Method
R. A. Berry; M. O. Delchini; J. Ragusa
2014-06-01
The RELAP-7 code is the next generation nuclear reactor system safety analysis code being developed at the Idaho National Laboratory (INL). The code is based on the INL's modern scientific software development framework, MOOSE (Multi-Physics Object Oriented Simulation Environment). The overall design goal of RELAP-7 is to take advantage of the previous thirty years of advancements in computer architecture, software design, numerical integration methods, and physical models. The end result will be a reactor systems analysis capability that retains and improves upon RELAP5's capability and extends the analysis capability for all reactor system simulation scenarios. RELAP-7 utilizes a single phase and a novel seven-equation two-phase flow models as described in the RELAP-7 Theory Manual (INL/EXT-14-31366). The basic equation systems are hyperbolic, which generally require some type of stabilization (or artificial viscosity) to capture nonlinear discontinuities and to suppress advection-caused oscillations. This report documents one of the available options for this stabilization in RELAP-7 -- a new and novel approach known as the entropy viscosity method. Because the code is an ongoing development effort in which the physical sub models, numerics, and coding are evolving, so too must the specific details of the entropy viscosity stabilization method. Here the fundamentals of the method in their current state are presented.
A new hybrid electro-numerical model of the left ventricle.
Kozarski, Maciej; Ferrari, Gianfranco; Zieliński, Krzysztof; Górczyńska, Krystyna; Pałko, Krzysztof J; Tokarz, Arkadiusz; Darowski, Marek
2008-09-01
The paper presents a new project of a hybrid numerical-physical model of the left ventricle. A physical part of the model can be based on electrical or hydraulic structures. Four variants of the model with numerical and physical heart valves have been designed to investigate an effect of a heart assistance connected in series and in parallel to the natural heart. The LabVIEW real time environment has been used in the model to increase its accuracy and reliability. A prototype of the hybrid electro-numerical model of the left ventricle has been tested in an open loop and closed loop configuration. PMID:18762290
Numerical methods for engine-airframe integration
Murthy, S.N.B.; Paynter, G.C.
1986-01-01
Various papers on numerical methods for engine-airframe integration are presented. The individual topics considered include: scientific computing environment for the 1980s, overview of prediction of complex turbulent flows, numerical solutions of the compressible Navier-Stokes equations, elements of computational engine/airframe integrations, computational requirements for efficient engine installation, application of CAE and CFD techniques to complete tactical missile design, CFD applications to engine/airframe integration, and application of a second-generation low-order panel methods to powerplant installation studies. Also addressed are: three-dimensional flow analysis of turboprop inlet and nacelle configurations, application of computational methods to the design of large turbofan engine nacelles, comparison of full potential and Euler solution algorithms for aeropropulsive flow field computations, subsonic/transonic, supersonic nozzle flows and nozzle integration, subsonic/transonic prediction capabilities for nozzle/afterbody configurations, three-dimensional viscous design methodology of supersonic inlet systems for advanced technology aircraft, and a user's technology assessment.
Numerical analysis method for linear induction machines.
NASA Technical Reports Server (NTRS)
Elliott, D. G.
1972-01-01
A numerical analysis method has been developed for linear induction machines such as liquid metal MHD pumps and generators and linear motors. Arbitrary phase currents or voltages can be specified and the moving conductor can have arbitrary velocity and conductivity variations from point to point. The moving conductor is divided into a mesh and coefficients are calculated for the voltage induced at each mesh point by unit current at every other mesh point. Combining the coefficients with the mesh resistances yields a set of simultaneous equations which are solved for the unknown currents.
Numerical methods for finding stationary gravitational solutions
NASA Astrophysics Data System (ADS)
Dias, Óscar J. C.; Santos, Jorge E.; Way, Benson
2016-07-01
The wide applications of higher dimensional gravity and gauge/gravity duality have fuelled the search for new stationary solutions of the Einstein equation (possibly coupled to matter). In this topical review, we explain the mathematical foundations and give a practical guide for the numerical solution of gravitational boundary value problems. We present these methods by way of example: resolving asymptotically flat black rings, singly spinning lumpy black holes in anti-de Sitter (AdS), and the Gregory–Laflamme zero modes of small rotating black holes in AdS{}5× {S}5. We also include several tools and tricks that have been useful throughout the literature.
Cilfone, Nicholas A.; Kirschner, Denise E.; Linderman, Jennifer J.
2015-01-01
Biologically related processes operate across multiple spatiotemporal scales. For computational modeling methodologies to mimic this biological complexity, individual scale models must be linked in ways that allow for dynamic exchange of information across scales. A powerful methodology is to combine a discrete modeling approach, agent-based models (ABMs), with continuum models to form hybrid models. Hybrid multi-scale ABMs have been used to simulate emergent responses of biological systems. Here, we review two aspects of hybrid multi-scale ABMs: linking individual scale models and efficiently solving the resulting model. We discuss the computational choices associated with aspects of linking individual scale models while simultaneously maintaining model tractability. We demonstrate implementations of existing numerical methods in the context of hybrid multi-scale ABMs. Using an example model describing Mycobacterium tuberculosis infection, we show relative computational speeds of various combinations of numerical methods. Efficient linking and solution of hybrid multi-scale ABMs is key to model portability, modularity, and their use in understanding biological phenomena at a systems level. PMID:26366228
Tsiambas, E; Karameris, A; Lygeros, M; Athanasiou, A E; Salemis, N S; Gourgiotis, S; Ragkos, V; Metaxas, G E; Vilaras, G; Patsouris, E
2012-01-01
Design and development of novel targeted therapeutic strategies is an innovation in handling patients with solid malignancies including breast, colon, lung, head & neck or even pancreatic and hepatocellular carcinoma. For a long time, immunohistocytochemistry (IHC/ICC) has been performed as a routine method in almost all labs for evaluating protein expression. Modern molecular approaches show that identification of specific structural and numerical imbalances regarding genes involved in signal transduction pathways provide important data to the oncologists. Alterations in molecules such as epidermal growth factor receptor (EGFR), HER2/neu, PTEN or Topoisomerase IIa affect the response rates to specific chemotherapeutic agents modifying also patients' prognostic rates. In situ hybridization (ISH) techniques based on fluorescence and chromogenic variants (FISH/CISH) or silver in situ hybridization (SISH) are applicable in both tissue and cell substrates. Concerning cytological specimens, FISH/CISH analysis appears to be a fast and very accurate method in estimating gene/chromosome ratios. In this paper, we sought to evaluate the usefulness of FISH/ CISH analysis in cytological specimens, describing also the advantages and disadvantages of these methods from the technical point of view. PMID:23033306
Application of numerical methods to elasticity imaging.
Castaneda, Benjamin; Ormachea, Juvenal; Rodríguez, Paul; Parker, Kevin J
2013-03-01
Elasticity imaging can be understood as the intersection of the study of biomechanical properties, imaging sciences, and physics. It was mainly motivated by the fact that pathological tissue presents an increased stiffness when compared to surrounding normal tissue. In the last two decades, research on elasticity imaging has been an international and interdisciplinary pursuit aiming to map the viscoelastic properties of tissue in order to provide clinically useful information. As a result, several modalities of elasticity imaging, mostly based on ultrasound but also on magnetic resonance imaging and optical coherence tomography, have been proposed and applied to a number of clinical applications: cancer diagnosis (prostate, breast, liver), hepatic cirrhosis, renal disease, thyroiditis, arterial plaque evaluation, wall stiffness in arteries, evaluation of thrombosis in veins, and many others. In this context, numerical methods are applied to solve forward and inverse problems implicit in the algorithms in order to estimate viscoelastic linear and nonlinear parameters, especially for quantitative elasticity imaging modalities. In this work, an introduction to elasticity imaging modalities is presented. The working principle of qualitative modalities (sonoelasticity, strain elastography, acoustic radiation force impulse) and quantitative modalities (Crawling Waves Sonoelastography, Spatially Modulated Ultrasound Radiation Force (SMURF), Supersonic Imaging) will be explained. Subsequently, the areas in which numerical methods can be applied to elasticity imaging are highlighted and discussed. Finally, we present a detailed example of applying total variation and AM-FM techniques to the estimation of elasticity. PMID:24010245
Mathematica with a Numerical Methods Course
NASA Astrophysics Data System (ADS)
Varley, Rodney
2003-04-01
An interdisciplinary "Numerical Methods" course has been shared between physics, mathematics and computer science since 1992 at Hunter C. Recently, the lectures and workshops for this course have become formalized and placed on the internet at http://www.ph.hunter.cuny.edu (follow the links "Course Listings and Websites" >> "PHYS385 (Numerical Methods)". Mathematica notebooks for the lectures are available for automatic download (by "double clicking" the lecture icon) for student use in the classroom or at home. AOL (or Netscape/Explorer) can be used provided Mathematica (or the "free" MathReader) has been made a "helper application". Using Mathematica has the virtue that mathematical equations (no LaTex required) can easily be included with the text and Mathematica's graphing is easy to use. Computational cells can be included within the notebook and students may easily modify the calculation to see the result of "what if..." questions. Homework is sent as Mathematica notebooks to the instructor via the internet and the corrected workshops are returned in the same manner. Most exam questions require computational solutions.
Lee, H.; Lee, D.
2013-07-01
This paper presents a new hybrid method of continuous energy Monte Carlo (MC) and multi-group Method of Characteristics (MOC). For a continuous energy neutron transport analysis, the hybrid method employs a continuous energy MC for resonance energy range to treat the resonances accurately and a multi-group MOC for high and low energy ranges for efficiency. Numerical test with a model problem confirms that the hybrid method can produce consistent results with the reference continuous energy MC-only calculation as well as multi-group MOC-only calculation. (authors)
Optimized Vertex Method and Hybrid Reliability
NASA Technical Reports Server (NTRS)
Smith, Steven A.; Krishnamurthy, T.; Mason, B. H.
2002-01-01
A method of calculating the fuzzy response of a system is presented. This method, called the Optimized Vertex Method (OVM), is based upon the vertex method but requires considerably fewer function evaluations. The method is demonstrated by calculating the response membership function of strain-energy release rate for a bonded joint with a crack. The possibility of failure of the bonded joint was determined over a range of loads. After completing the possibilistic analysis, the possibilistic (fuzzy) membership functions were transformed to probability density functions and the probability of failure of the bonded joint was calculated. This approach is called a possibility-based hybrid reliability assessment. The possibility and probability of failure are presented and compared to a Monte Carlo Simulation (MCS) of the bonded joint.
A hybrid method for image interpolation
NASA Astrophysics Data System (ADS)
Qi, Chun; Huang, Hua; Wang, Wen-Bing; Zhang, Jing; Di, Shuangliang
2003-09-01
A hybrid method for image interpolation is proposed. The method consists of three different approaches: Circular arc or B-spline interpolation, linear interpolatino and human visual sensitivity based on interpolation. The image can be divided into three regions: linear smooth region, sharp edge region and human visual insensitive region. The method uses local variance and mean value to find different regions adaptively. The linear interpolation is used for linear smooth region. The human visual sensitivity based interpolation is used for human visual insensitive region and the circular arc or B-spline interpolation is used for sharp edge region. Experiments show that proposed method produces results that are more visually realistic than standard function-fitting methods.
NASA Technical Reports Server (NTRS)
Homicz, G. F.; Moselle, J. R.
1985-01-01
A hybrid numerical procedure is presented for the prediction of the aerodynamic and acoustic performance of advanced turboprops. A hybrid scheme is proposed which in principle leads to a consistent simultaneous prediction of both fields. In the inner flow a finite difference method, the Approximate-Factorization Alternating-Direction-Implicit (ADI) scheme, is used to solve the nonlinear Euler equations. In the outer flow the linearized acoustic equations are solved via a Boundary-Integral Equation (BIE) method. The two solutions are iteratively matched across a fictitious interface in the flow so as to maintain continuity. At convergence the resulting aerodynamic load prediction will automatically satisfy the appropriate free-field boundary conditions at the edge of the finite difference grid, while the acoustic predictions will reflect the back-reaction of the radiated field on the magnitude of the loading source terms, as well as refractive effects in the inner flow. The equations and logic needed to match the two solutions are developed and the computer program implementing the procedure is described. Unfortunately, no converged solutions were obtained, due to unexpectedly large running times. The reasons for this are discussed and several means to alleviate the situation are suggested.
Numerical methods for problems in computational aeroacoustics
NASA Astrophysics Data System (ADS)
Mead, Jodi Lorraine
1998-12-01
A goal of computational aeroacoustics is the accurate calculation of noise from a jet in the far field. This work concerns the numerical aspects of accurately calculating acoustic waves over large distances and long time. More specifically, the stability, efficiency, accuracy, dispersion and dissipation in spatial discretizations, time stepping schemes, and absorbing boundaries for the direct solution of wave propagation problems are determined. Efficient finite difference methods developed by Tam and Webb, which minimize dispersion and dissipation, are commonly used for the spatial and temporal discretization. Alternatively, high order pseudospectral methods can be made more efficient by using the grid transformation introduced by Kosloff and Tal-Ezer. Work in this dissertation confirms that the grid transformation introduced by Kosloff and Tal-Ezer is not spectrally accurate because, in the limit, the grid transformation forces zero derivatives at the boundaries. If a small number of grid points are used, it is shown that approximations with the Chebyshev pseudospectral method with the Kosloff and Tal-Ezer grid transformation are as accurate as with the Chebyshev pseudospectral method. This result is based on the analysis of the phase and amplitude errors of these methods, and their use for the solution of a benchmark problem in computational aeroacoustics. For the grid transformed Chebyshev method with a small number of grid points it is, however, more appropriate to compare its accuracy with that of high- order finite difference methods. This comparison, for an order of accuracy 10-3 for a benchmark problem in computational aeroacoustics, is performed for the grid transformed Chebyshev method and the fourth order finite difference method of Tam. Solutions with the finite difference method are as accurate. and the finite difference method is more efficient than, the Chebyshev pseudospectral method with the grid transformation. The efficiency of the Chebyshev
Small RNA Detection by in Situ Hybridization Methods
Urbanek, Martyna O.; Nawrocka, Anna U.; Krzyzosiak, Wlodzimierz J.
2015-01-01
Small noncoding RNAs perform multiple regulatory functions in cells, and their exogenous mimics are widely used in research and experimental therapies to interfere with target gene expression. MicroRNAs (miRNAs) are the most thoroughly investigated representatives of the small RNA family, which includes short interfering RNAs (siRNAs), PIWI-associated RNA (piRNAs), and others. Numerous methods have been adopted for the detection and characterization of small RNAs, which is challenging due to their short length and low level of expression. These include molecular biology methods such as real-time RT-PCR, northern blotting, hybridization to microarrays, cloning and sequencing, as well as single cell miRNA detection by microscopy with in situ hybridization (ISH). In this review, we focus on the ISH method, including its fluorescent version (FISH), and we present recent methodological advances that facilitated its successful adaptation for small RNA detection. We discuss relevant technical aspects as well as the advantages and limitations of ISH. We also refer to numerous applications of small RNA ISH in basic research and molecular diagnostics. PMID:26068454
Small RNA Detection by in Situ Hybridization Methods.
Urbanek, Martyna O; Nawrocka, Anna U; Krzyzosiak, Wlodzimierz J
2015-01-01
Small noncoding RNAs perform multiple regulatory functions in cells, and their exogenous mimics are widely used in research and experimental therapies to interfere with target gene expression. MicroRNAs (miRNAs) are the most thoroughly investigated representatives of the small RNA family, which includes short interfering RNAs (siRNAs), PIWI-associated RNA (piRNAs), and others. Numerous methods have been adopted for the detection and characterization of small RNAs, which is challenging due to their short length and low level of expression. These include molecular biology methods such as real-time RT-PCR, northern blotting, hybridization to microarrays, cloning and sequencing, as well as single cell miRNA detection by microscopy with in situ hybridization (ISH). In this review, we focus on the ISH method, including its fluorescent version (FISH), and we present recent methodological advances that facilitated its successful adaptation for small RNA detection. We discuss relevant technical aspects as well as the advantages and limitations of ISH. We also refer to numerous applications of small RNA ISH in basic research and molecular diagnostics. PMID:26068454
Hybrid perturbation methods based on statistical time series models
NASA Astrophysics Data System (ADS)
San-Juan, Juan Félix; San-Martín, Montserrat; Pérez, Iván; López, Rosario
2016-04-01
In this work we present a new methodology for orbit propagation, the hybrid perturbation theory, based on the combination of an integration method and a prediction technique. The former, which can be a numerical, analytical or semianalytical theory, generates an initial approximation that contains some inaccuracies derived from the fact that, in order to simplify the expressions and subsequent computations, not all the involved forces are taken into account and only low-order terms are considered, not to mention the fact that mathematical models of perturbations not always reproduce physical phenomena with absolute precision. The prediction technique, which can be based on either statistical time series models or computational intelligence methods, is aimed at modelling and reproducing missing dynamics in the previously integrated approximation. This combination results in the precision improvement of conventional numerical, analytical and semianalytical theories for determining the position and velocity of any artificial satellite or space debris object. In order to validate this methodology, we present a family of three hybrid orbit propagators formed by the combination of three different orders of approximation of an analytical theory and a statistical time series model, and analyse their capability to process the effect produced by the flattening of the Earth. The three considered analytical components are the integration of the Kepler problem, a first-order and a second-order analytical theories, whereas the prediction technique is the same in the three cases, namely an additive Holt-Winters method.
Nunn, D. )
1993-09-01
This paper reports a simple novel technique for the numerical simulation of hot collision-free plasmas. The method is termed Vlasov hybrid simulation (VHS). A time varying phase space simulation box and grid are defined, and the phase fluid within the box is filled with simulation particles. The distribution function F (or [sigma]F) is defined on the phase trajectory of each particle. At each timestep F (or [sigma]F) is interpolated from the simulation particles onto the phase space grid. Particles are followed continuously until exiting from the phase box and are not constantly recreated at phase space grid points. The algorithm is very efficient, stable, and has low noise levels. Distribution function fine structure is tolerated and the formalism does not require diffusion of the distribution function. The VHS method is particularly valuable when the flux of phase fluid across the phase box boundary is significant. In this case VHS codes have a dynamic population of particles-giving great efficiency gains over PIC codes with fixed particle populations. The VHS method has been applied to the numerical simulation of triggered VLF emissions in the magnetosphere and gives results in close agreement with observations. 27 refs., 13 figs., 1 tab.
Hybrid optimization methods for Full Waveform Inversion
NASA Astrophysics Data System (ADS)
Datta, D.; Sen, M. K.
2014-12-01
FWI is slowly becoming the mainstream method to estimate velocity models of the subsurface from seismic data. Typically it makes use of a gradient descent approach in which a model update is computed by back propagating the residual seismograms and cross correlating with the forward propagating wavefields at each grid point in the subsurface model. FWI is a local optimization technique, which requires the starting model to be very close to the true model. Because the objective function is multimodal with many local minima, the requirement of good starting model becomes essential. A starting model is generated using travel time tomography. We propose two hybrid FWI algorithms one of which generates a very good starting model for a conventional FWI and the other, which works with a population of models uses gradient information from multiple starting locations in guiding the search. The first approach uses a sparse parameterization of model space using non-oscillatory splines, whose coeffiencts are estimated using an optimization algorithm like very fast simulated annealing (VFSA) by minimizing the misfit between the observed and synthetic data. The estimated velocity model is then used as a starting model for gradient-based FWI. This is done in the shot domain by converting the end-on marine geometry to a split spread geometry using the principle of reciprocity. The second approach is to uses an alternate global optimization algorithm called particle swarm optimization (PSO) where PSO update rules are applied. However, we employ a new gradient guided PSO that exploits the gradient information as well. This approach avoids the local minima and converges faster than a conventional PSO. We demonstrate our methods with application to 2D marine data sets from offshore India. Each line comprises over 1000 shots; our hybrid methods produce geologically meaningful velocity models fairly rapidly on a GPU cluster. We show that starting with the hybrid model gives a much
NASA Astrophysics Data System (ADS)
Reverdy, Frédéric; Mahaut, Steve; Dominguez, Nicolas; Dubois, Philippe
2015-03-01
Carbon Fiber reinforced composites are increasingly used in structural parts in the aeronautics industry, as they allow to reduce the weight of aircrafts while maintaining high mechanical performances. However, such structures can be complicated to inspect due to their complex geometries and complex composite properties, leading to highly heterogeneous and anisotropic materials. Different potential damages and manufacturing flaws related to these parts are to be detected: porosities, ply waviness, delaminations after impact. Ultrasonic inspection, which is commonly used to test the full volume of composite panels, thus has to cope with both complex wave propagation (within anisotropic parts whose crystallographic orientation varies according to the layers structure) and flaw interaction (local distortion of plies such as ply waviness, small pores, structural noise due to periodicity patterns…). Developing NDT procedures for those parts therefore requires simulation tools to help for understanding those phenomena, and to optimize probes and techniques. Within the CIVA multi-techniques platform, CEA-LIST has developed semi-analytical tools for ultrasonic techniques, which have the advantages of high computational efficiency (fast calculations), but with limited range of application due to some hypothesis (for instance, homogenization approaches which don't allow to take account of structural noise). On the other hand, numerical methods such as finite element (FEM) or finite difference in time domain (FDTD) are more suitable to compute ultrasonic wave propagation and defect scattering in complex materials such as composite but require more computational efforts. Hybrid methods couple semi-analytical solutions and numerical computations in limited spatial domains to handle complex cases with high computation performances. In CIVA we have integrated a hybrid model that combines the semi-analytical methods developed at CEA to FDTD codes developed at Airbus Group
Numerical modeling of lower hybrid heating and current drive
Valeo, E.J.; Eder, D.C.
1986-03-01
The generation of currents in toroidal plasma by application of waves in the lower hybrid frequency range involves the interplay of several physical phenomena which include: wave propagation in toroidal geometry, absorption via wave-particle resonances, the quasilinear generation of strongly nonequilibrium electron and ion distribution functions, and the self-consistent evolution of the current density in such a nonequilibrium plasma. We describe a code, LHMOD, which we have developed to treat these aspects of current drive and heating in tokamaks. We present results obtained by applying the code to a computation of current ramp-up and to an investigation of the possible importance of minority hydrogen absorption in a deuterium plasma as the ''density limit'' to current drive is approached.
Transport Test Problems for Hybrid Methods Development
Shaver, Mark W.; Miller, Erin A.; Wittman, Richard S.; McDonald, Benjamin S.
2011-12-28
This report presents 9 test problems to guide testing and development of hybrid calculations for the ADVANTG code at ORNL. These test cases can be used for comparing different types of radiation transport calculations, as well as for guiding the development of variance reduction methods. Cases are drawn primarily from existing or previous calculations with a preference for cases which include experimental data, or otherwise have results with a high level of confidence, are non-sensitive, and represent problem sets of interest to NA-22.
A separable shadow Hamiltonian hybrid Monte Carlo method
NASA Astrophysics Data System (ADS)
Sweet, Christopher R.; Hampton, Scott S.; Skeel, Robert D.; Izaguirre, Jesús A.
2009-11-01
Hybrid Monte Carlo (HMC) is a rigorous sampling method that uses molecular dynamics (MD) as a global Monte Carlo move. The acceptance rate of HMC decays exponentially with system size. The shadow hybrid Monte Carlo (SHMC) was previously introduced to reduce this performance degradation by sampling instead from the shadow Hamiltonian defined for MD when using a symplectic integrator. SHMC's performance is limited by the need to generate momenta for the MD step from a nonseparable shadow Hamiltonian. We introduce the separable shadow Hamiltonian hybrid Monte Carlo (S2HMC) method based on a formulation of the leapfrog/Verlet integrator that corresponds to a separable shadow Hamiltonian, which allows efficient generation of momenta. S2HMC gives the acceptance rate of a fourth order integrator at the cost of a second-order integrator. Through numerical experiments we show that S2HMC consistently gives a speedup greater than two over HMC for systems with more than 4000 atoms for the same variance. By comparison, SHMC gave a maximum speedup of only 1.6 over HMC. S2HMC has the additional advantage of not requiring any user parameters beyond those of HMC. S2HMC is available in the program PROTOMOL 2.1. A Python version, adequate for didactic purposes, is also in MDL (http://mdlab.sourceforge.net/s2hmc).
Novel Parallel Numerical Methods for Radiation& Neutron Transport
Brown, P N
2001-03-06
In many of the multiphysics simulations performed at LLNL, transport calculations can take up 30 to 50% of the total run time. If Monte Carlo methods are used, the percentage can be as high as 80%. Thus, a significant core competence in the formulation, software implementation, and solution of the numerical problems arising in transport modeling is essential to Laboratory and DOE research. In this project, we worked on developing scalable solution methods for the equations that model the transport of photons and neutrons through materials. Our goal was to reduce the transport solve time in these simulations by means of more advanced numerical methods and their parallel implementations. These methods must be scalable, that is, the time to solution must remain constant as the problem size grows and additional computer resources are used. For iterative methods, scalability requires that (1) the number of iterations to reach convergence is independent of problem size, and (2) that the computational cost grows linearly with problem size. We focused on deterministic approaches to transport, building on our earlier work in which we performed a new, detailed analysis of some existing transport methods and developed new approaches. The Boltzmann equation (the underlying equation to be solved) and various solution methods have been developed over many years. Consequently, many laboratory codes are based on these methods, which are in some cases decades old. For the transport of x-rays through partially ionized plasmas in local thermodynamic equilibrium, the transport equation is coupled to nonlinear diffusion equations for the electron and ion temperatures via the highly nonlinear Planck function. We investigated the suitability of traditional-solution approaches to transport on terascale architectures and also designed new scalable algorithms; in some cases, we investigated hybrid approaches that combined both.
NASA Astrophysics Data System (ADS)
Schuster, Jonathan
Infrared (IR) detectors are well established as a vital sensor technology for military, defense and commercial applications. Due to the expense and effort required to fabricate pixel arrays, it is imperative to develop numerical simulation models to perform predictive device simulations which assess device characteristics and design considerations. Towards this end, we have developed a robust three-dimensional (3D) numerical simulation model for IR detector pixel arrays. We used the finite-difference time-domain technique to compute the optical characteristics including the reflectance and the carrier generation rate in the device. Subsequently, we employ the finite element method to solve the drift-diffusion equations to compute the electrical characteristics including the I(V) characteristics, quantum efficiency, crosstalk and modulation transfer function. We use our 3D numerical model to study a new class of detector based on the nBn-architecture. This detector is a unipolar unity-gain barrier device consisting of a narrow-gap absorber layer, a wide-gap barrier layer, and a narrow-gap collector layer. We use our model to study the underlying physics of these devices and to explain the anomalously long lateral collection lengths for photocarriers measured experimentally. Next, we investigate the crosstalk in HgCdTe photovoltaic pixel arrays employing a photon-trapping (PT) structure realized with a periodic array of pillars intended to provide broadband operation. The PT region drastically reduces the crosstalk; making the use of the PT structures not only useful to obtain broadband operation, but also desirable for reducing crosstalk, especially in small pitch detector arrays. Then, the power and flexibility of the nBn architecture is coupled with a PT structure to engineer spectrally filtering detectors. Last, we developed a technique to reduce the cost of large-format, high performance HgCdTe detectors by nondestructively screen-testing detector arrays prior
Optimization methods applied to hybrid vehicle design
NASA Technical Reports Server (NTRS)
Donoghue, J. F.; Burghart, J. H.
1983-01-01
The use of optimization methods as an effective design tool in the design of hybrid vehicle propulsion systems is demonstrated. Optimization techniques were used to select values for three design parameters (battery weight, heat engine power rating and power split between the two on-board energy sources) such that various measures of vehicle performance (acquisition cost, life cycle cost and petroleum consumption) were optimized. The apporach produced designs which were often significant improvements over hybrid designs already reported on in the literature. The principal conclusions are as follows. First, it was found that the strategy used to split the required power between the two on-board energy sources can have a significant effect on life cycle cost and petroleum consumption. Second, the optimization program should be constructed so that performance measures and design variables can be easily changed. Third, the vehicle simulation program has a significant effect on the computer run time of the overall optimization program; run time can be significantly reduced by proper design of the types of trips the vehicle takes in a one year period. Fourth, care must be taken in designing the cost and constraint expressions which are used in the optimization so that they are relatively smooth functions of the design variables. Fifth, proper handling of constraints on battery weight and heat engine rating, variables which must be large enough to meet power demands, is particularly important for the success of an optimization study. Finally, the principal conclusion is that optimization methods provide a practical tool for carrying out the design of a hybrid vehicle propulsion system.
Validation of a numerical method for unsteady flow calculations
Giles, M.; Haimes, R. . Dept. of Aeronautics and Astronautics)
1993-01-01
This paper describes and validates a numerical method for the calculation of unsteady inviscid and viscous flows. A companion paper compares experimental measurements of unsteady heat transfer on a transonic rotor with the corresponding computational results. The mathematical model is the Reynolds-averaged unsteady Navier-Stokes equations for a compressible ideal gas. Quasi-three-dimensionality is included through the use of a variable streamtube thickness. The numerical algorithm is unusual in two respects: (a) For reasons of efficiency and flexibility, it uses a hybrid Navier-Stokes/Euler method, and (b) to allow for the computation of stator/rotor combinations with arbitrary pitch ratio, a novel space-time coordinate transformation is used. Several test cases are presented to validate the performance of the computer program, UNSFLO. These include: (a) unsteady, inviscid flat plate cascade flows (b) steady and unsteady, viscous flat plate cascade flows, (c) steady turbine heat transfer and loss prediction. In the first two sets of cases comparisons are made with theory, and in the third the comparison is with experimental data.
NASA Astrophysics Data System (ADS)
Niknam, A.; Rajabi, A. A.; Solaimani, M.
2016-03-01
Solution of the radial Schrodinger equation for the Woods-Saxon potential together with spin-orbit interaction, coulomb and centrifugal terms by using usual Nikiforov-Uvarov (NU) method is not possible. Here, we have presented a new NU procedure with which we are able to solve this Schrodinger equation and any other one-dimensional ones with any shape of the potential profile. For this purpose, we have combined the NU method with numerical fitting schema. The energy eigenvalues and corresponding eigenfunctions for various values of n, l, and j quantum numbers have been obtained. Good agreement with experimental values is also achieved. We have calculated the 1/2+ state energy with more accuracy (our absolute error = 0.023 MeV and Hagen et al. absolute error = 0.0918 MeV), while Hagen et al. have calculated the 5/2+ state energy with higher accuracy (our absolute error = 0.71 MeV and Hagen et al. absolute error = 0.0337 MeV). Our wave functions are in agreement with Kim et al.'s work, too.
NASA Technical Reports Server (NTRS)
Banyukevich, A.; Ziolkovski, K.
1975-01-01
A number of hybrid methods for solving Cauchy problems are described on the basis of an evaluation of advantages of single and multiple-point numerical integration methods. The selection criterion is the principle of minimizing computer time. The methods discussed include the Nordsieck method, the Bulirsch-Stoer extrapolation method, and the method of recursive Taylor-Steffensen power series.
Numerical methods for analyzing electromagnetic scattering
NASA Technical Reports Server (NTRS)
Lee, S. W.; Lo, Y. T.; Chuang, S. L.; Lee, C. S.
1985-01-01
Numerical methods to analyze electromagnetic scattering are presented. The dispersions and attenuations of the normal modes in a circular waveguide coated with lossy material were completely analyzed. The radar cross section (RCS) from a circular waveguide coated with lossy material was calculated. The following is observed: (1) the interior irradiation contributes to the RCS much more than does the rim diffraction; (2) at low frequency, the RCS from the circular waveguide terminated by a perfect electric conductor (PEC) can be reduced more than 13 dB down with a coating thickness less than 1% of the radius using the best lossy material available in a 6 radius-long cylinder; (3) at high frequency, a modal separation between the highly attenuated and the lowly attenuated modes is evident if the coating material is too lossy, however, a large RCS reduction can be achieved for a small incident angle with a thin layer of coating. It is found that the waveguide coated with a lossy magnetic material can be used as a substitute for a corrugated waveguide to produce a circularly polarized radiation yield.
A numerical method for cardiac mechanoelectric simulations.
Pathmanathan, Pras; Whiteley, Jonathan P
2009-05-01
Much effort has been devoted to developing numerical techniques for solving the equations that describe cardiac electrophysiology, namely the monodomain equations and bidomain equations. Only a limited selection of publications, however, address the development of numerical techniques for mechanoelectric simulations where cardiac electrophysiology is coupled with deformation of cardiac tissue. One problem commonly encountered in mechanoelectric simulations is instability of the coupled numerical scheme. In this study, we develop a stable numerical scheme for mechanoelectric simulations. A number of convergence tests are carried out using this stable technique for simulations where deformations are of the magnitude typically observed in a beating heart. These convergence tests demonstrate that accurate computation of tissue deformation requires a nodal spacing of around 1 mm in the mesh used to calculate tissue deformation. This is a much finer computational grid than has previously been acknowledged, and has implications for the computational efficiency of the resulting numerical scheme. PMID:19263223
Huang, Xuehui; Yang, Shihua; Gong, Junyi; Zhao, Yan; Feng, Qi; Gong, Hao; Li, Wenjun; Zhan, Qilin; Cheng, Benyi; Xia, Junhui; Chen, Neng; Hao, Zhongna; Liu, Kunyan; Zhu, Chuanrang; Huang, Tao; Zhao, Qiang; Zhang, Lei; Fan, Danlin; Zhou, Congcong; Lu, Yiqi; Weng, Qijun; Wang, Zi-Xuan; Li, Jiayang; Han, Bin
2015-01-01
Exploitation of heterosis is one of the most important applications of genetics in agriculture. However, the genetic mechanisms of heterosis are only partly understood, and a global view of heterosis from a representative number of hybrid combinations is lacking. Here we develop an integrated genomic approach to construct a genome map for 1,495 elite hybrid rice varieties and their inbred parental lines. We investigate 38 agronomic traits and identify 130 associated loci. In-depth analyses of the effects of heterozygous genotypes reveal that there are only a few loci with strong overdominance effects in hybrids, but a strong correlation is observed between the yield and the number of superior alleles. While most parental inbred lines have only a small number of superior alleles, high-yielding hybrid varieties have several. We conclude that the accumulation of numerous rare superior alleles with positive dominance is an important contributor to the heterotic phenomena. PMID:25651972
Hybrid natural element method for large deformation elastoplasticity problems
NASA Astrophysics Data System (ADS)
Ma, Yong-Qi; Zhou, Yan-Kai
2015-03-01
We present the hybrid natural element method (HNEM) for two-dimensional elastoplastic large deformation problems. Sibson interpolation is adopted to construct the shape functions of nodal incremental displacements and incremental stresses. The incremental form of Hellinger-Reissner variational principle for elastoplastic large deformation problems is deduced to obtain the equation system. The total Lagrangian formulation is used to describe the discrete equation system. Compared with the natural element method (NEM), the HNEM has higher computational precision and efficiency in solving elastoplastic large deformation problems. Some numerical examples are selected to demonstrate the advantage of the HNEM for large deformation elastoplasticity problems. Project supported by the Natural Science Foundation of Shanghai, China (Grant No. 13ZR1415900).
Development of a Hybrid RANS/LES Method for Compressible Mixing Layer Simulations
NASA Technical Reports Server (NTRS)
Georgiadis, Nicholas J.; Alexander, J. Iwan D.; Reshotko, Eli
2001-01-01
A hybrid method has been developed for simulations of compressible turbulent mixing layers. Such mixing layers dominate the flows in exhaust systems of modem day aircraft and also those of hypersonic vehicles currently under development. The hybrid method uses a Reynolds-averaged Navier-Stokes (RANS) procedure to calculate wall bounded regions entering a mixing section, and a Large Eddy Simulation (LES) procedure to calculate the mixing dominated regions. A numerical technique was developed to enable the use of the hybrid RANS/LES method on stretched, non-Cartesian grids. The hybrid RANS/LES method is applied to a benchmark compressible mixing layer experiment. Preliminary two-dimensional calculations are used to investigate the effects of axial grid density and boundary conditions. Actual LES calculations, performed in three spatial directions, indicated an initial vortex shedding followed by rapid transition to turbulence, which is in agreement with experimental observations.
C deg continuity elements by Hybrid Stress method. M.S. Thesis, 1982 Final Report
NASA Technical Reports Server (NTRS)
Kang, David Sung-Soo
1991-01-01
An intensive study of the assumed variable distribution necessary for the Assumed Displacement Formulation, the Hellinger-Reissner Formulation, and the Hu-Washizu Formulation is made in a unified manner. With emphasis on physical explanation, a systematic method for the Hybrid Stress element construction is outlined. The numerical examples use four and eight node plane stress elements and eight and twenty node solid elements. Computation cost study indicates that the hybrid stress element derived using recently developed Uncoupled Stress Formulation is comparable in CPU time to the Assumed Displacement element. Overall, main emphasis is placed on providing a broader understanding of the Hybrid Stress Formulation.
Hybrid selection as a method of increasing mapping power for radiation hybrids.
Jones, H B
1996-08-01
Radiation hybrids have become a widely used tool for physical mapping. A drawback of the technique is that large numbers of hybrids are required to construct robust, high-resolution maps. The information contained within a panel of radiation hybrids is limited by the frequency of retention of chromosomal fragments from the donor cell line. In almost all experiments to date, the retention frequency has been below the optimal level; therefore, many hybrids are needed to produce high-quality maps. Because of the labor-intensive nature of large-scale mapping projects, it is important to make panels as small as possible. One method that has been adopted is to produce initially a large number of hybrids that are all typed with a few loci. Those hybrids showing satisfactorily high retention are admitted to the final panel and the rest are discarded. In this way, a panel of radiation hybrids with higher than expected retention can be created. Methods for conducting such a selection regime are discussed. To investigate the potential advantages of selecting hybrids based on their retention frequency, simulations were run under a variety of conditions. As expected panels with high retention (40%) provided better mapping resources than panels with lower (20%) retention. Beginning with an initial panel of 200 hybrids, comparisons of a random sample of 100 hybrids and the set of those 100 hybrids showing the highest marker retention demonstrated that selection may not be always the best strategy despite the increase in mean retention it yields. The selection of hybrids containing large numbers of fragments leads to an overestimation of the frequency of radiation-induced breaks. When breaks occur with high frequency (for example, when high radiation doses are used), the selection of hybrids leads to a loss of linkage and hence an inability to order the markers. As such, the merits of screening hybrids depends on both the radiation dose and the desired map resolution. PMID:8858351
NASA Astrophysics Data System (ADS)
Bendaoud, Issam; Matteï, Simone; Cicala, Eugen; Tomashchuk, Iryna; Andrzejewski, Henri; Sallamand, Pierre; Mathieu, Alexandre; Bouchaud, Fréderic
2014-03-01
The present study is dedicated to the numerical simulation of an industrial case of hybrid laser-MIG welding of high thickness duplex steel UR2507Cu with Y-shaped chamfer geometry. It consists in simulation of heat transfer phenomena using heat equivalent source approach and implementing in finite element software COMSOL Multiphysics. A numerical exploratory designs method is used to identify the heat sources parameters in order to obtain a minimal required difference between the numerical results and the experiment which are the shape of the welded zone and the temperature evolution in different locations. The obtained results were found in good correspondence with experiment, both for melted zone shape and thermal history.
Evaluation of the Performance of the Hybrid Lattice Boltzmann Based Numerical Flux
NASA Astrophysics Data System (ADS)
Zheng, H. W.; Shu, C.
2016-06-01
It is well known that the numerical scheme is a key factor to the stability and accuracy of a Navier-Stokes solver. Recently, a new hybrid lattice Boltzmann numerical flux (HLBFS) is developed by Shu's group. It combines two different LBFS schemes by a switch function. It solves the Boltzmann equation instead of the Euler equation. In this article, the main object is to evaluate the ability of this HLBFS scheme by our in-house cell centered hybrid mesh based Navier-Stokes code. Its performance is examined by several widely-used bench-mark test cases. The comparisons on results between calculation and experiment are conducted. They show that the scheme can capture the shock wave as well as the resolving of boundary layer.
Experimental and Numerical Characterization of a Hybrid Fabry-Pérot Cavity for Temperature Sensing
Lopez-Aldaba, Aitor; Pinto, Ana Margarida Rodrigues; Lopez-Amo, Manuel; Frazão, Orlando; Santos, José Luís; Baptista, José Manuel; Baierl, Hardy; Auguste, Jean-Louis; Jamier, Raphael; Roy, Philippe
2015-01-01
A hybrid Fabry-Pérot cavity sensing head based on a four-bridge microstructured fiber is characterized for temperature sensing. The characterization of this cavity is performed numerically and experimentally in the L-band. The sensing head output signal presents a linear variation with temperature changes, showing a sensitivity of 12.5 pm/°C. Moreover, this Fabry-Pérot cavity exhibits good sensitivity to polarization changes and high stability over time. PMID:25853404
Finite element methods in numerical relativity.
NASA Astrophysics Data System (ADS)
Mann, P. J.
The finite element method is very successful in Newtonian fluid simulations, and can be extended to relativitstic fluid flows. This paper describes the general method, and then outlines some preliminary results for spherically symmetric geometries. The mixed finite element - finite difference scheme is introduced, and used for the description of spherically symmetric collapse. Baker's (Newtonian) shock modelling method and Miller's moving finite element method are also mentioned. Collapse in double-null coordinates requires non-constant time slicing, so the full finite element method in space and time is described.
Numerical matrix method for quantum periodic potentials
NASA Astrophysics Data System (ADS)
Le Vot, Felipe; Meléndez, Juan J.; Yuste, Santos B.
2016-06-01
A numerical matrix methodology is applied to quantum problems with periodic potentials. The procedure consists essentially in replacing the true potential by an alternative one, restricted by an infinite square well, and in expressing the wave functions as finite superpositions of eigenfunctions of the infinite well. A matrix eigenvalue equation then yields the energy levels of the periodic potential within an acceptable accuracy. The methodology has been successfully used to deal with problems based on the well-known Kronig-Penney (KP) model. Besides the original model, these problems are a dimerized KP solid, a KP solid containing a surface, and a KP solid under an external field. A short list of additional problems that can be solved with this procedure is presented.
Method for numerical simulations of metastable states
Heller, U.M.; Seiberg, N.
1983-06-15
We present a numerical simulation of metastable states near a first-order phase transition in the example of a U(1) lattice gauge theory with a generalized action. In order to make measurements in these states possible their decay has to be prevented. We achieve this by using a microcanonical simulation for a finite system. We then obtain the coupling constant (inverse temperature) as a function of the action density. It turns out to be nonmonotonic and hence not uniquely invertible. From it we derive the effective potential for the action density. This effective potential is not always convex, a property that seems to be in contradiction with the standard lore about its convexity. This apparent ''paradox'' is resolved in a discussion about different definitions of the effective potential.
Interpolation Method Needed for Numerical Uncertainty
NASA Technical Reports Server (NTRS)
Groves, Curtis E.; Ilie, Marcel; Schallhorn, Paul A.
2014-01-01
Using Computational Fluid Dynamics (CFD) to predict a flow field is an approximation to the exact problem and uncertainties exist. There is a method to approximate the errors in CFD via Richardson's Extrapolation. This method is based off of progressive grid refinement. To estimate the errors, the analyst must interpolate between at least three grids. This paper describes a study to find an appropriate interpolation scheme that can be used in Richardson's extrapolation or other uncertainty method to approximate errors.
Numerical methods in Markov chain modeling
NASA Technical Reports Server (NTRS)
Philippe, Bernard; Saad, Youcef; Stewart, William J.
1989-01-01
Several methods for computing stationary probability distributions of Markov chains are described and compared. The main linear algebra problem consists of computing an eigenvector of a sparse, usually nonsymmetric, matrix associated with a known eigenvalue. It can also be cast as a problem of solving a homogeneous singular linear system. Several methods based on combinations of Krylov subspace techniques are presented. The performance of these methods on some realistic problems are compared.
Hybrid N-order Lagrangian interpolation Eulerian-Lagrangian method for salinity calculation
NASA Astrophysics Data System (ADS)
Wu, Yan-cheng; Zhu, Shou-xian; Zhou, Lin; You, Xiao-bao; Zhang, Wen-jing
2016-04-01
The Eulerian-Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order ( N=3-8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.
New developments in the multiscale hybrid energy density computational method
NASA Astrophysics Data System (ADS)
Min, Sun; Shanying, Wang; Dianwu, Wang; Chongyu, Wang
2016-01-01
Further developments in the hybrid multiscale energy density method are proposed on the basis of our previous papers. The key points are as follows. (i) The theoretical method for the determination of the weight parameter in the energy coupling equation of transition region in multiscale model is given via constructing underdetermined equations. (ii) By applying the developed mathematical method, the weight parameters have been given and used to treat some problems in homogeneous charge density systems, which are directly related with multiscale science. (iii) A theoretical algorithm has also been presented for treating non-homogeneous systems of charge density. The key to the theoretical computational methods is the decomposition of the electrostatic energy in the total energy of density functional theory for probing the spanning characteristic at atomic scale, layer by layer, by which the choice of chemical elements and the defect complex effect can be understood deeply. (iv) The numerical computational program and design have also been presented. Project supported by the National Basic Research Program of China (Grant No. 2011CB606402) and the National Natural Science Foundation of China (Grant No. 51071091).
An Effective Hybrid Firefly Algorithm with Harmony Search for Global Numerical Optimization
Guo, Lihong; Wang, Gai-Ge; Wang, Heqi; Wang, Dinan
2013-01-01
A hybrid metaheuristic approach by hybridizing harmony search (HS) and firefly algorithm (FA), namely, HS/FA, is proposed to solve function optimization. In HS/FA, the exploration of HS and the exploitation of FA are fully exerted, so HS/FA has a faster convergence speed than HS and FA. Also, top fireflies scheme is introduced to reduce running time, and HS is utilized to mutate between fireflies when updating fireflies. The HS/FA method is verified by various benchmarks. From the experiments, the implementation of HS/FA is better than the standard FA and other eight optimization methods. PMID:24348137
Chandekar, Gautam S.; Kelkar, Ajit D.
2014-01-01
In the present study experimental and numerical investigations were carried out to predict the low velocity impact response of four symmetric configurations: 10 ply E Glass, 10 ply AS4 Carbon, and two Hybrid combinations with 1 and 2 outer plies of E Glass and 8 and 6 inner plies of Carbon. All numerical investigations were performed using commercial finite element software, LS-DYNA. The test coupons were manufactured using the low cost Heated Vacuum Assisted Resin Transfer Molding (H-VARTM©) technique. Low velocity impact testing was carried out using an Instron Dynatup 8250 impact testing machine. Standard 6 × 6 Boeing fixture was used for all impact experiments. Impact experiments were performed over progressive damage, that is, from incipient damage till complete failure of the laminate in six successive impact energy levels for each configuration. The simulation results for the impact loading were compared with the experimental results. For both nonhybrid configurations, it was observed that the simulated results were in good agreement with the experimental results, whereas, for hybrid configurations, the simulated impact response was softer than the experimental response. Maximum impact load carrying capacity was also compared for all four configurations based on their areal density. It was observed that Hybrid262 configuration has superior impact load to areal density ratio. PMID:24719573
Chandekar, Gautam S; Kelkar, Ajit D
2014-01-01
In the present study experimental and numerical investigations were carried out to predict the low velocity impact response of four symmetric configurations: 10 ply E Glass, 10 ply AS4 Carbon, and two Hybrid combinations with 1 and 2 outer plies of E Glass and 8 and 6 inner plies of Carbon. All numerical investigations were performed using commercial finite element software, LS-DYNA. The test coupons were manufactured using the low cost Heated Vacuum Assisted Resin Transfer Molding (H-VARTM©) technique. Low velocity impact testing was carried out using an Instron Dynatup 8250 impact testing machine. Standard 6 × 6 Boeing fixture was used for all impact experiments. Impact experiments were performed over progressive damage, that is, from incipient damage till complete failure of the laminate in six successive impact energy levels for each configuration. The simulation results for the impact loading were compared with the experimental results. For both nonhybrid configurations, it was observed that the simulated results were in good agreement with the experimental results, whereas, for hybrid configurations, the simulated impact response was softer than the experimental response. Maximum impact load carrying capacity was also compared for all four configurations based on their areal density. It was observed that Hybrid262 configuration has superior impact load to areal density ratio. PMID:24719573
Modelling asteroid brightness variations. I - Numerical methods
NASA Technical Reports Server (NTRS)
Karttunen, H.
1989-01-01
A method for generating lightcurves of asteroid models is presented. The effects of the shape of the asteroid and the scattering law of a surface element are distinctly separable, being described by chosen functions that can easily be changed. The shape is specified by means of two functions that yield the length of the radius vector and the normal vector of the surface at a given point. The general shape must be convex, but spherical concavities producing macroscopic shadowing can also be modeled.
Thermal analysis of 3D composites by a new fast multipole hybrid boundary node method
NASA Astrophysics Data System (ADS)
Miao, Yu; Wang, Qiao; Zhu, Hongping; Li, Yinping
2014-01-01
This paper applies the hybrid boundary node method (Hybrid BNM) for the thermal analysis of 3D composites. A new formulation is derived for the inclusion-based composites. In the new formulation, the unknowns of the interfaces are assembled only once in the final system equation, which can reduce nearly one half of degrees of freedom (DOFs) compared with the conventional multi-domain solver when there are lots of inclusions. A new version of the fast multipole method (FMM) is also coupled with the new formulation and the technique is applied to thermal analysis of composites with many inclusions. In the new fast multipole hybrid boundary node method (FM-HBNM), a diagonal form for translation operators is used and the method presented can be applied to the computation of more than 1,000,000 DOFs on a personal computer. Numerical examples are presented to analyze the thermal behavior of composites with many inclusions.
A numerical method for power plant simulations
Carcasci, C.; Facchini, B.
1996-03-01
This paper describes a highly flexible computerized method of calculating operating data in a power cycle. The computerized method presented here permits the study of steam, gas and combined plants. Its flexibility is not restricted by any defined cycle scheme. A power plant consists of simple elements (turbine, compressor, combustor chamber, pump, etc.). Each power plant component is represented by its typical equations relating to fundamental mechanical and thermodynamic laws, so a power plant system is represented by algebraic equations, which are the typical equations of components, continuity equations, and data concerning plant conditions. This equation system is not linear, but can be reduced to a linear equation system with variable coefficients. The solution is simultaneous for each component and it is determined by an iterative process. An example of a simple gas turbine cycle demonstrates the applied technique. This paper also presents the user interface based on MS-Windows. The input data, the results, and any characteristic parameters of a complex cycle scheme are also shown.
Petrov-Galerkin's method hybrid with finite element into the Helmholtz equation solution. Part II
NASA Astrophysics Data System (ADS)
Rabadan Malda, Itzala; Salazar Cordero, Emigdio; Ortega Herrera, Jose Angel
2002-11-01
This work proposes a hybridization between Petrov-Galerkins numeric method and finite element method (FEM) to resolve Helmholtz equation when dominion is an open or semiopen tube-shaped configuration and with determinate number of holes over cylindrical surface. It's pretended to solve these kind of cavities, thereby it allows us to obtain very important design parameters like: cavity length, quantity, size and distance between toneholes, form and size of mouthpiece or outlet. These parameters are design basis into acoustic and musical instrumentation: baffles outlet pipes, diffusers, silencers, flutes, oboes, saxophones, trumpets, quenas, and many more. In this way it's expected to determine advantages of this numeric method above another using actually.
A hybrid method for solving time-domain integral equations in transient scattering
NASA Astrophysics Data System (ADS)
Tijhuis, A. G.; Wiemans, R.; Kuester, E. F.
A new hybrid method is proposed for the numerical solution of integral equations describing transient scattering problems. The basic idea behind the method is to first discretize in space, and then solve the resulting system of linear time-domain equations by carrying out a temporal Laplace transformation. This approach combines the efficiency of the marching-on-in-time method with the stability and the accuracy of frequency-domain techniques. Numerical results are presented for the scattering of E-polarized, pulsed waves by a one-dimensionally inhomogeneous, lossy dielectric slab located in between two homogeneous, lossless dielectric half-spaces, and by a radially inhomogeneous, lossy dielectric circular cylinder embedded in vacuum. For these problems, the hybrid method turns out to be more powerful than existing solution techniques.
Method of hybrid plume plasma propulsion
NASA Technical Reports Server (NTRS)
Chang, Franklin R. (Inventor)
1990-01-01
A technique for producing thrust by generating a hybrid plume plasma exhaust is disclosed. A plasma flow is generated and introduced into a nozzle which features one or more inlets positioned to direct a flow of neutral gas about the interior of the nozzle. When such a neutral gas flow is combined with the plasma flow within the nozzle, a hybrid plume is constructed including a flow of hot plasma along the center of the nozzle surrounded by a generally annular flow of neutral gas, with an annular transition region between the pure plasma and the neutral gas. The temperature of the outer gas layer is below that of the pure plasma and generally separates the pure plasma from the interior surfaces of the nozzle. The neutral gas flow both insulates the nozzle walls from the high temperatures of the plasma flow and adds to the mass flow rate of the hybrid exhaust. The rate of flow of neutral gas into the interior of the nozzle may be selectively adjusted to control the thrust and specific impulse of the device.
Numerical methods for analyzing electromagnetic scattering
NASA Technical Reports Server (NTRS)
Lee, S. W.; Lo, Y. T.; Chuang, S. L.; Lee, C. S.
1985-01-01
Attenuation properties of the normal modes in an overmoded waveguide coated with a lossy material were analyzed. It is found that the low-order modes, can be significantly attenuated even with a thin layer of coating if the coating material is not too lossy. A thinner layer of coating is required for large attenuation of the low-order modes if the coating material is magnetic rather than dielectric. The Radar Cross Section (RCS) from an uncoated circular guide terminated by a perfect electric conductor was calculated and compared with available experimental data. It is confirmed that the interior irradiation contributes to the RCS. The equivalent-current method based on the geometrical theory of diffraction (GTD) was chosen for the calculation of the contribution from the rim diffraction. The RCS reduction from a coated circular guide terminated by a PEC are planned schemes for the experiments are included. The waveguide coated with a lossy magnetic material is suggested as a substitute for the corrugated waveguide.
Application of Hybrid Method for Aerodynamic Noise Prediction
NASA Astrophysics Data System (ADS)
Yu, L.; Song, W. P.
2011-09-01
A hybrid prediction method for aerodynamic noise is performed using high order accuracy method in this paper. The method combines a two-dimensional Unsteady Reynolds-Averaged Navier-Stokes(URANS) solver with the acoustic analogy method using Ffowcs Williams-Hawkings equation with penetrable data surface (FW-Hpds). Tandem cylinders are chosen to validate the prediction method. The computations are conducted at a Reynolds number of 1.66 × 105 based on the cylinder diameter. Both the aerodynamic and acoustic results show good agreement with the experimental data, showing a successful application of the hybrid prediction method using two-dimensional URANS simulation.
Hybrid intelligent optimization methods for engineering problems
NASA Astrophysics Data System (ADS)
Pehlivanoglu, Yasin Volkan
The purpose of optimization is to obtain the best solution under certain conditions. There are numerous optimization methods because different problems need different solution methodologies; therefore, it is difficult to construct patterns. Also mathematical modeling of a natural phenomenon is almost based on differentials. Differential equations are constructed with relative increments among the factors related to yield. Therefore, the gradients of these increments are essential to search the yield space. However, the landscape of yield is not a simple one and mostly multi-modal. Another issue is differentiability. Engineering design problems are usually nonlinear and they sometimes exhibit discontinuous derivatives for the objective and constraint functions. Due to these difficulties, non-gradient-based algorithms have become more popular in recent decades. Genetic algorithms (GA) and particle swarm optimization (PSO) algorithms are popular, non-gradient based algorithms. Both are population-based search algorithms and have multiple points for initiation. A significant difference from a gradient-based method is the nature of the search methodologies. For example, randomness is essential for the search in GA or PSO. Hence, they are also called stochastic optimization methods. These algorithms are simple, robust, and have high fidelity. However, they suffer from similar defects, such as, premature convergence, less accuracy, or large computational time. The premature convergence is sometimes inevitable due to the lack of diversity. As the generations of particles or individuals in the population evolve, they may lose their diversity and become similar to each other. To overcome this issue, we studied the diversity concept in GA and PSO algorithms. Diversity is essential for a healthy search, and mutations are the basic operators to provide the necessary variety within a population. After having a close scrutiny of the diversity concept based on qualification and
Methods for assessing DNA hybridization of PNA-TiO2 nanoconjugates
Brown, Eric M. B.; Paunesku, Tatjana; Wu, AiGuo; Thurn, K. Ted; Haley, Benjamin; Clark, Jimmy; Priester, Taisa; Woloschak, Gayle E.
2008-01-01
We describe the synthesis of peptide nucleic acid (PNA)-titanium dioxide (TiO2) nanoconjugates and the several novel methods developed to investigate the DNA hybridization behaviors of these constructs. PNAs are synthetic DNA analogs resistant to degradation by cellular enzymes, which hybridize to single strand DNA (ssDNA) with higher affinity than DNA oligonucleotides, invade double strand DNA (dsDNA), and form different PNA-DNA complexes. Previously, we developed a DNA-TiO2 nanoconjugate capable of hybridizing to target DNA intracellularly in a sequence-specific manner, with the ability to cleave DNA when excited by electromagnetic radiation, but susceptible to degradation which may lower its intracellular targeting efficiency and retention time. PNA-TiO2 nanoconjugates described herein hybridize to target ssDNA, oligonucleotide dsDNA, and supercoiled plasmid DNA under physiological-like ionic and temperature conditions, enabling rapid and inexpensive, sequence-specific precipitation of nucleic acids in vitro. When modified by the addition of imaging agents or peptides, hybridization capabilities of PNA-TiO2 nanoconjugates are enhanced which provides essential benefits for numerous in vitro and in vivo applications. The series of experiments shown here could not be done with either TiO2-DNA nanoconjugates or PNAs alone, and the novel methods developed will benefit studies of numerous other nanoconjugate systems. PMID:18786502
Cleft Lip Repair: The Hybrid Subunit Method.
Tollefson, Travis T
2016-04-01
The unilateral cleft lip repair is one of the most rewarding and challenging of plastic surgery procedures. Surgeons have introduced a variety of straight line, geometric, and rotation-advancement designs, while in practice the majority of North American surgeons have been using hybrids of the rotation-advancement techniques. The anatomic subunit approach was introduced in 2005 by Fisher and has gained popularity, with early adopters of the design touting its simplicity and effectiveness. The objectives of this article are to summarize the basic tenets of respecting the philtral subunit, accurate measurement and planning, and tips for transitioning to this subunit approach. PMID:27097136
Validation with experiments on simplified numerical prediction of hybrid rocket internal ballistics
NASA Astrophysics Data System (ADS)
Funami, Yuki; Shimada, Toru
2012-11-01
In order to design hybrid rocket engines, we have developed a numerical prediction approach to the internal ballistics. The key point is its cost performance. Therefore simple but efficient models are required. Fluid phenomenon and thermal conduction phenomenon in a solid fuel should be treated time-dependently, because characteristic times of these phenomena are longer than those of other phenomena. Besides, they are solved with the energy-flux balance equation at the solid fuel surface to determine the regression rate. It is confirmed that numerical evaluation of time- and space-averaged regression rate is the same order of magnitude as that in experiments. However, the factors n in ṙ¯ = aG¯oxn differ between calculations and experiments.
Method of producing a hybrid matrix fiber composite
Deteresa, Steven J.; Lyon, Richard E.; Groves, Scott E.
2006-03-28
Hybrid matrix fiber composites having enhanced compressive performance as well as enhanced stiffness, toughness and durability suitable for compression-critical applications. The methods for producing the fiber composites using matrix hybridization. The hybrid matrix fiber composites comprised of two chemically or physically bonded matrix materials, whereas the first matrix materials are used to impregnate multi-filament fibers formed into ribbons and the second matrix material is placed around and between the fiber ribbons that are impregnated with the first matrix material and both matrix materials are cured and solidified.
Numerical Modelling of Staged Combustion Aft-Injected Hybrid Rocket Motors
NASA Astrophysics Data System (ADS)
Nijsse, Jeff
The staged combustion aft-injected hybrid (SCAIH) rocket motor is a promising design for the future of hybrid rocket propulsion. Advances in computational fluid dynamics and scientific computing have made computational modelling an effective tool in hybrid rocket motor design and development. The focus of this thesis is the numerical modelling of the SCAIH rocket motor in a turbulent combustion, high-speed, reactive flow framework accounting for solid soot transport and radiative heat transfer. The SCAIH motor is modelled with a shear coaxial injector with liquid oxygen injected in the center at sub-critical conditions: 150 K and 150 m/s (Mach ≈ 0.9), and a gas-generator gas-solid mixture of one-third carbon soot by mass injected in the annual opening at 1175 K and 460 m/s (Mach ≈ 0.6). Flow conditions in the near injector region and the flame anchoring mechanism are of particular interest. Overall, the flow is shown to exhibit instabilities and the flame is shown to anchor directly on the injector faceplate with temperatures in excess of 2700 K.
Hybrid methods for multisource information fusion and decision support
NASA Astrophysics Data System (ADS)
Braun, Jerome J.; Glina, Yan
2006-04-01
This paper presents the progress of an ongoing research effort in multisource information fusion for biodefense decision support. The effort concentrates on a novel machine-intelligence hybrid-of-hybrids decision support architecture termed FLASH (Fusion, Learning, Adaptive Super-Hybrid) we proposed. The highlights of FLASH discussed in the paper include its cognitive-processing orientation and the hybrid nature involving heterogeneous multiclassifier machine learning and approximate reasoning paradigms. Selected specifics of the FLASH internals, such as its feature selection techniques, supervised learning, clustering, recognition and reasoning methods, and their integration, are discussed. The results to date are presented, including the background type determination and bioattack detection computational experiments using data obtained with a multisensor fusion testbed we have also developed. The processing of imprecise information originating from sources other than sensors is considered. Finally, the paper discusses applicability of FLASH and its methods to complex battlespace management problems such as course-of-action decision support.
Hybrid architecture active wavefront sensing and control system, and method
NASA Technical Reports Server (NTRS)
Feinberg, Lee D. (Inventor); Dean, Bruce H. (Inventor); Hyde, Tristram T. (Inventor)
2011-01-01
According to various embodiments, provided herein is an optical system and method that can be configured to perform image analysis. The optical system can comprise a telescope assembly and one or more hybrid instruments. The one or more hybrid instruments can be configured to receive image data from the telescope assembly and perform a fine guidance operation and a wavefront sensing operation, simultaneously, on the image data received from the telescope assembly.
Hybrid materials and methods for producing the same
Luzzi, David E.; Smith, Brian W.
2008-02-19
A hybrid material is provided which comprises a first single-walled nanotube having a lumen, and a fill molecule contained within the lumen of the single-walled nanotube. A method for producing the hybrid material is also provided wherein a single-walled nanotube is contacted with a fill molecule to cause the fill molecule to enter the lumen of the single-walled nanotube.
Hybrid materials and methods for producing the same
Luzzi, David E.; Smith, Brian W.
2003-04-08
A hybrid material is provided which comprises a first single-walled nanotube having a lumen, and a fill molecule contained within the lumen of the single-walled nanotube. A method for producing the hybrid material is also provided wherein a single-walled nanotube is contacted with a fill molecule to cause the fill molecule to enter the lumen of the single-walled nanotube.
A Novel Numerical Method for Fuzzy Boundary Value Problems
NASA Astrophysics Data System (ADS)
Can, E.; Bayrak, M. A.; Hicdurmaz
2016-05-01
In the present paper, a new numerical method is proposed for solving fuzzy differential equations which are utilized for the modeling problems in science and engineering. Fuzzy approach is selected due to its important applications on processing uncertainty or subjective information for mathematical models of physical problems. A second-order fuzzy linear boundary value problem is considered in particular due to its important applications in physics. Moreover, numerical experiments are presented to show the effectiveness of the proposed numerical method on specific physical problems such as heat conduction in an infinite plate and a fin.
Battery control system for hybrid vehicle and method for controlling a hybrid vehicle battery
Bockelmann, Thomas R.; Hope, Mark E.; Zou, Zhanjiang; Kang, Xiaosong
2009-02-10
A battery control system for hybrid vehicle includes a hybrid powertrain battery, a vehicle accessory battery, and a prime mover driven generator adapted to charge the vehicle accessory battery. A detecting arrangement is configured to monitor the vehicle accessory battery's state of charge. A controller is configured to activate the prime mover to drive the generator and recharge the vehicle accessory battery in response to the vehicle accessory battery's state of charge falling below a first predetermined level, or transfer electrical power from the hybrid powertrain battery to the vehicle accessory battery in response to the vehicle accessory battery's state of charge falling below a second predetermined level. The invention further includes a method for controlling a hybrid vehicle powertrain system.
[A hybrid volume rendering method using general hardware].
Li, Bin; Tian, Lianfang; Chen, Ping; Mao, Zongyuan
2008-06-01
In order to improve the effect and efficiency of the reconstructed image after hybrid volume rendering of different kinds of volume data from medical sequential slices or polygonal models, we propose a hybrid volume rendering method based on Shear-Warp with economical hardware. First, the hybrid volume data are pre-processed by Z-Buffer method and RLE (Run-Length Encoded) data structure. Then, during the process of compositing intermediate image, a resampling method based on the dual-interpolation and the intermediate slice interpolation methods is used to improve the efficiency and the effect. Finally, the reconstructed image is rendered by the texture-mapping technology of OpenGL. Experiments demonstrate the good performance of the proposed method. PMID:18693424
Asymptotic-induced numerical methods for conservation laws
NASA Technical Reports Server (NTRS)
Garbey, Marc; Scroggs, Jeffrey S.
1990-01-01
Asymptotic-induced methods are presented for the numerical solution of hyperbolic conservation laws with or without viscosity. The methods consist of multiple stages. The first stage is to obtain a first approximation by using a first-order method, such as the Godunov scheme. Subsequent stages of the method involve solving internal-layer problems identified by using techniques derived via asymptotics. Finally, a residual correction increases the accuracy of the scheme. The method is derived and justified with singular perturbation techniques.
Method for production of sorghum hybrids with selected flowering times
Mullet, John E.; Rooney, William L.
2016-08-30
Methods and composition for the production of sorghum hybrids with selected and different flowering times are provided. In accordance with the invention, a substantially continual and high-yield harvest of sorghum is provided. Improved methods of seed production are also provided.
Parallel processing numerical method for confined vortex dynamics and applications
NASA Astrophysics Data System (ADS)
Bistrian, Diana Alina
2013-10-01
This paper explores a combined analytical and numerical technique to investigate the hydrodynamic instability of confined swirling flows, with application to vortex rope dynamics in a Francis turbine diffuser, in condition of sophisticated boundary constraints. We present a new approach based on the method of orthogonal decomposition in the Hilbert space, implemented with a spectral descriptor scheme in discrete space. A parallel implementation of the numerical scheme is conducted reducing the computational time compared to other techniques.
A hybrid formulation of a component mode synthesis method
NASA Technical Reports Server (NTRS)
Farhat, Charbel; Geradin, Michel
1992-01-01
Component mode synthesis is a substructuring technique frequently employed in structural dynamics. In this method, a given structure is subdivided into components or substructures, each of which is analyzed independently for natural frequencies and for mode shapes. The substructure mode shapes are then assembled to give displacement shapes or load patterns of the original structure. An analytical justification of the basic concept is presented using spectral decompositions, and a variant substructuring approach where intersubstructure continuity is enforced in a weak form is derived. This leads to a hybrid formulation of the basic method which is particularly suitable for assembling heterogeneous substructures and analyzing nonconforming and incompatible finite element substructure models. For problems where both the basic and hybrid methods are applicable, the hybrid variant can be computationally more advantageous.
Numerical and experimental studies of the hybrid rocket motor with multi-port fuel grain
NASA Astrophysics Data System (ADS)
Tian, Hui; Li, Xintian; Zeng, Peng; Yu, Nanjia; Cai, Guobiao
2014-03-01
This paper presents three-dimensional numerical simulations and experimental studies of the hybrid rocket motor with multi-port fuel grain. The numerical model is established based on the Navier-Stokes equations with turbulence, chemical reactions, fuel pyrolysis, and solid-gas boundary interactions. The simulation is performed based on the 98% hydrogen peroxide (HP) and hydroxyl terminated polybutadiene (HTPB) propellant combination. The results indicate that the flow field and fuel regression rate distributions present apparent three-dimensional characteristics. The fuel regression rates decrease first and then gradually increase with the axial location increasing. At a certain cross section, the fuel regression rates are lower in the points on arcs with smaller radius of curvature when the fuel port is a derivable convex figure. Two experiments are carried out on a full scale motor with the simulation one. The working process of the motor is steady and no evident oscillatory combustion is observed. The fuel port profiles before and after tests indicate that the fuel regression rate distributions at the cross section match well with the numerical simulation results.
NASA Astrophysics Data System (ADS)
Robertson, Eric D.
A comprehensive survey of available numerical methods and models was performed on the open source computational fluid dynamics solver OpenFOAM version 2.0 for incompressible turbulent bluff body flows. Numerical methods are illuminated using source code for side-by-side comparison. For validation, the accuracy of flow predictions over a sphere in the subcritical regime and delta wing with sharp leading edge is assessed. Solutions show mostly good agreement with experimental data and data obtained from commercial software. A demonstration of the numerical implementation of a dynamic hybrid RANS/LES framework is also presented, including results from test studies.
Collocation Method for Numerical Solution of Coupled Nonlinear Schroedinger Equation
Ismail, M. S.
2010-09-30
The coupled nonlinear Schroedinger equation models several interesting physical phenomena presents a model equation for optical fiber with linear birefringence. In this paper we use collocation method to solve this equation, we test this method for stability and accuracy. Numerical tests using single soliton and interaction of three solitons are used to test the resulting scheme.
Investigating Convergence Patterns for Numerical Methods Using Data Analysis
ERIC Educational Resources Information Center
Gordon, Sheldon P.
2013-01-01
The article investigates the patterns that arise in the convergence of numerical methods, particularly those in the errors involved in successive iterations, using data analysis and curve fitting methods. In particular, the results obtained are used to convey a deeper level of understanding of the concepts of linear, quadratic, and cubic…
A numerical method for solving singular De`s
Mahaver, W.T.
1996-12-31
A numerical method is developed for solving singular differential equations using steepest descent based on weighted Sobolev gradients. The method is demonstrated on a variety of first and second order problems, including linear constrained, unconstrained, and partially constrained first order problems, a nonlinear first order problem with irregular singularity, and two second order variational problems.
NASA Technical Reports Server (NTRS)
Rosen, I. G.
1984-01-01
A cubic spline based Galerkin-like method is developed for the identification of a class of hybrid systems which describe the transverse vibration to flexible beams with attached tip bodies. The identification problem is formulated as a least squares fit to data subject to the system dynamics given by a coupled system of ordnary and partial differential equations recast as an abstract evolution equation (AEE) in an appropriate infinite dimensional Hilbert space. Projecting the AEE into spline-based subspaces leads naturally to a sequence of approximating finite dimensional identification problems. The solutions to these problems are shown to exist, are relatively easily computed, and are shown to, in some sense, converge to solutions to the original identification problem. Numerical results for a variety of examples are discussed.
NASA Technical Reports Server (NTRS)
Graf, W.; Chang, T. Y.; Saleeb, A. F.
1986-01-01
Three-dimensional thick shell elements with 8, 16, and 18 nodes are formulated by using the hybrid/mixed method. In bending applications, these elements are free from locking effect and give improved stress predictions. Finite element equations are derived from the Hellinger-Reissner variational principle in which both the displacement and stress fields are approximated by independent interpolation functions. For the assumption of stress parameters, three guidelines are followed: (1) suppression of kinematic deformation modes, (2) invariant element property, and (3) the constraint index exhibited by the element, when applied to constrained-media problems, must be greater than or equal to one. Numerical results are presented to show the element's behavior characteristics regarding sensitivity to locking, distortion effect (patch tests), mesh convergence and the accuracy of stress evaluation.
Multiple-time-stepping generalized hybrid Monte Carlo methods
NASA Astrophysics Data System (ADS)
Escribano, Bruno; Akhmatskaya, Elena; Reich, Sebastian; Azpiroz, Jon M.
2015-01-01
Performance of the generalized shadow hybrid Monte Carlo (GSHMC) method [1], which proved to be superior in sampling efficiency over its predecessors [2-4], molecular dynamics and hybrid Monte Carlo, can be further improved by combining it with multi-time-stepping (MTS) and mollification of slow forces. We demonstrate that the comparatively simple modifications of the method not only lead to better performance of GSHMC itself but also allow for beating the best performed methods, which use the similar force splitting schemes. In addition we show that the same ideas can be successfully applied to the conventional generalized hybrid Monte Carlo method (GHMC). The resulting methods, MTS-GHMC and MTS-GSHMC, provide accurate reproduction of thermodynamic and dynamical properties, exact temperature control during simulation and computational robustness and efficiency. MTS-GHMC uses a generalized momentum update to achieve weak stochastic stabilization to the molecular dynamics (MD) integrator. MTS-GSHMC adds the use of a shadow (modified) Hamiltonian to filter the MD trajectories in the HMC scheme. We introduce a new shadow Hamiltonian formulation adapted to force-splitting methods. The use of such Hamiltonians improves the acceptance rate of trajectories and has a strong impact on the sampling efficiency of the method. Both methods were implemented in the open-source MD package ProtoMol and were tested on a water and a protein systems. Results were compared to those obtained using a Langevin Molly (LM) method [5] on the same systems. The test results demonstrate the superiority of the new methods over LM in terms of stability, accuracy and sampling efficiency. This suggests that putting the MTS approach in the framework of hybrid Monte Carlo and using the natural stochasticity offered by the generalized hybrid Monte Carlo lead to improving stability of MTS and allow for achieving larger step sizes in the simulation of complex systems.
Multiple-time-stepping generalized hybrid Monte Carlo methods
Escribano, Bruno; Akhmatskaya, Elena; Reich, Sebastian; Azpiroz, Jon M.
2015-01-01
Performance of the generalized shadow hybrid Monte Carlo (GSHMC) method [1], which proved to be superior in sampling efficiency over its predecessors [2–4], molecular dynamics and hybrid Monte Carlo, can be further improved by combining it with multi-time-stepping (MTS) and mollification of slow forces. We demonstrate that the comparatively simple modifications of the method not only lead to better performance of GSHMC itself but also allow for beating the best performed methods, which use the similar force splitting schemes. In addition we show that the same ideas can be successfully applied to the conventional generalized hybrid Monte Carlo method (GHMC). The resulting methods, MTS-GHMC and MTS-GSHMC, provide accurate reproduction of thermodynamic and dynamical properties, exact temperature control during simulation and computational robustness and efficiency. MTS-GHMC uses a generalized momentum update to achieve weak stochastic stabilization to the molecular dynamics (MD) integrator. MTS-GSHMC adds the use of a shadow (modified) Hamiltonian to filter the MD trajectories in the HMC scheme. We introduce a new shadow Hamiltonian formulation adapted to force-splitting methods. The use of such Hamiltonians improves the acceptance rate of trajectories and has a strong impact on the sampling efficiency of the method. Both methods were implemented in the open-source MD package ProtoMol and were tested on a water and a protein systems. Results were compared to those obtained using a Langevin Molly (LM) method [5] on the same systems. The test results demonstrate the superiority of the new methods over LM in terms of stability, accuracy and sampling efficiency. This suggests that putting the MTS approach in the framework of hybrid Monte Carlo and using the natural stochasticity offered by the generalized hybrid Monte Carlo lead to improving stability of MTS and allow for achieving larger step sizes in the simulation of complex systems.
A new numerical method of total solar eclipse photography processing
NASA Astrophysics Data System (ADS)
Druckmüller, M.; Rušin, V.; Minarovjech, M.
2006-10-01
A new numerical method of image processing suitable for visualization of corona images taken during total solar eclipses is presented. This method allows us to study both small- and large-scale coronal structures that remain invisible on original images because of their very high dynamic range of the coronal brightness. The method is based on the use of adaptive filters inspired by human vision and the sensitivity of resulting images is thus very close to that of the human eye during an eclipse. A high precision alignment method for white-light corona images is also discussed. The proposed method highly improves a widely used unsharp masking method employing a radially blurred mask. The results of these numerical image processing techniques are illustrated by a series of images taken during eclipses of the last decade. The method minimizes the risk of processing artifacts.
Robust Hybrid Finite Element Methods for Antennas and Microwave Circuits
NASA Technical Reports Server (NTRS)
Gong, J.; Volakis, John L.
1996-01-01
One of the primary goals in this dissertation is concerned with the development of robust hybrid finite element-boundary integral (FE-BI) techniques for modeling and design of conformal antennas of arbitrary shape. Both the finite element and integral equation methods will be first overviewed in this chapter with an emphasis on recently developed hybrid FE-BI methodologies for antennas, microwave and millimeter wave applications. The structure of the dissertation is then outlined. We conclude the chapter with discussions of certain fundamental concepts and methods in electromagnetics, which are important to this study.
25 Years of Self-organized Criticality: Numerical Detection Methods
NASA Astrophysics Data System (ADS)
McAteer, R. T. James; Aschwanden, Markus J.; Dimitropoulou, Michaila; Georgoulis, Manolis K.; Pruessner, Gunnar; Morales, Laura; Ireland, Jack; Abramenko, Valentyna
2016-01-01
The detection and characterization of self-organized criticality (SOC), in both real and simulated data, has undergone many significant revisions over the past 25 years. The explosive advances in the many numerical methods available for detecting, discriminating, and ultimately testing, SOC have played a critical role in developing our understanding of how systems experience and exhibit SOC. In this article, methods of detecting SOC are reviewed; from correlations to complexity to critical quantities. A description of the basic autocorrelation method leads into a detailed analysis of application-oriented methods developed in the last 25 years. In the second half of this manuscript space-based, time-based and spatial-temporal methods are reviewed and the prevalence of power laws in nature is described, with an emphasis on event detection and characterization. The search for numerical methods to clearly and unambiguously detect SOC in data often leads us outside the comfort zone of our own disciplines—the answers to these questions are often obtained by studying the advances made in other fields of study. In addition, numerical detection methods often provide the optimum link between simulations and experiments in scientific research. We seek to explore this boundary where the rubber meets the road, to review this expanding field of research of numerical detection of SOC systems over the past 25 years, and to iterate forwards so as to provide some foresight and guidance into developing breakthroughs in this subject over the next quarter of a century.
A hybrid perturbation-Galerkin method for differential equations containing a parameter
NASA Technical Reports Server (NTRS)
Geer, James F.; Andersen, Carl M.
1989-01-01
A two-step hybrid perturbation-Galerkin method to solve a variety of differential equations which involve a parameter is presented and discussed. The method consists of: (1) the use of a perturbation method to determine the asymptotic expansion of the solution about one or more values of the parameter; and (2) the use of some of the perturbation coefficient functions as trial functions in the classical Bubnov-Galerkin method. This hybrid method has the potential of overcoming some of the drawbacks of the perturbation method and the Bubnov-Galerkin method when they are applied by themselves, while combining some of the good features of both. The proposed method is illustrated first with a simple linear two-point boundary value problem and is then applied to a nonlinear two-point boundary value problem in lubrication theory. The results obtained from the hybrid method are compared with approximate solutions obtained by purely numerical methods. Some general features of the method, as well as some special tips for its implementation, are discussed. A survey of some current research application areas is presented and its degree of applicability to broader problem areas is discussed.
Numerical Simulation of Cavitation Phenomena for Hybrid Contra-Rotating Shaft Propellers
NASA Astrophysics Data System (ADS)
Kim, Sung-Hoon; Choo, Sung-Han; Park, Jeong-Yong; Choi, Gil-Hwan
2015-12-01
This paper deals with a numerical simulation of cavitation flow around a hybrid contra-rotating shaft propeller operating in wake field. The simulation for the cavitating flow is performed for straight operating and turning condition of podded propeller located behind the main propeller using unsteady Reynolds-Averaged Navier-Stokes. The behavior of the main propeller is almost similar regardless of the turning angle. In contrast, the cavitation behavior of the podded propeller depending on the turning angle appears to be entirely different due to the change of the load distribution on the podded propeller. At the large angle of the turning condition, the unstable cavity flow due to the large amount of cavitation and the hub vortex separated from the forward propeller as well as face cavitation is observed. Thus, a great caution on the cavitation phenomena is needed when designing and operating the HCRSP.
Comparison of methods for numerical calculation of continuum damping
Bowden, G. W.; Hole, M. J.; Dennis, G. R.; Könies, A.; Gorelenkov, N. N.
2014-05-15
Continuum resonance damping is an important factor in determining the stability of certain global modes in fusion plasmas. A number of analytic and numerical approaches have been developed to compute this damping, particularly, in the case of the toroidicity-induced shear Alfvén eigenmode. This paper compares results obtained using an analytical perturbative approach with those found using resistive and complex contour numerical approaches. It is found that the perturbative method does not provide accurate agreement with reliable numerical methods for the range of parameters examined. This discrepancy exists even in the limit where damping approaches zero. When the perturbative technique is implemented using a standard finite element method, the damping estimate fails to converge with radial grid resolution. The finite elements used cannot accurately represent the eigenmode in the region of the continuum resonance, regardless of the number of radial grid points used.
Hybrid Particle-Continuum Methods for Nonequilibrium Gas and Plasma Flows
Boyd, Iain D.
2011-05-20
Two different hybrid particle-continuum methods are described for simulation of nonequilibrium gas and plasma dynamics. The first technique, used for nonequilibrium hypersonic gas flows, uses either a continuum description or a particle method throughout a flow domain based on local conditions. This technique is successful in reproducing the results of full particle simulations at a small fraction of the cost. The second method uses a continuum model of the electrons combined with a particle description of the ions and atoms for simulating plasma jets. The physical accuracy of the method is assessed through comparisons with plasma plume measurements obtained in space. These examples illustrate that the complex physical phenomena associated with nonequilibrium conditions can be simulated with physical accuracy and numerical efficiency using such hybrid approaches.
Non-intrusive hybrid interval method for uncertain nonlinear systems using derivative information
NASA Astrophysics Data System (ADS)
Liu, Zhuang-Zhuang; Wang, Tian-Shu; Li, Jun-Feng
2016-02-01
This paper proposes a new non-intrusive hybrid interval method using derivative information for the dynamic response analysis of nonlinear systems with uncertain-but-bounded parameters and/or initial conditions. This method provides tighter solution ranges compared to the existing polynomial approximation interval methods. Interval arithmetic using the Chebyshev basis and interval arithmetic using the general form modified affine basis for polynomials are developed to obtain tighter bounds for interval computation. To further reduce the overestimation caused by the "wrapping effect" of interval arithmetic, the derivative information of dynamic responses is used to achieve exact solutions when the dynamic responses are monotonic with respect to all the uncertain variables. Finally, two typical numerical examples with nonlinearity are applied to demonstrate the effectiveness of the proposed hybrid interval method, in particular, its ability to effectively control the overestimation for specific timepoints.
A hybrid formulation for the numerical simulation of condensed phase explosives
NASA Astrophysics Data System (ADS)
Michael, L.; Nikiforakis, N.
2016-07-01
In this article we present a new formulation and an associated numerical algorithm, for the simulation of combustion and transition to detonation of condensed-phase commercial- and military-grade explosives, which are confined by (or in general interacting with one or more) compliant inert materials. Examples include confined rate-stick problems and interaction of shock waves with gas cavities or solid particles in explosives. This formulation is based on an augmented Euler approach to account for the mixture of the explosive and its products, and a multi-phase diffuse interface approach to solve for the immiscible interaction between the mixture and the inert materials, so it is in essence a hybrid (augmented Euler and multi-phase) model. As such, it has many of the desirable features of the two approaches and, critically for our applications of interest, it provides the accurate recovery of temperature fields across all components. Moreover, it conveys a lot more physical information than augmented Euler, without the complexity of full multi-phase Baer-Nunziato-type models or the lack of robustness of augmented Euler models in the presence of more than two components. The model can sustain large density differences across material interfaces without the presence of spurious oscillations in velocity and pressure, and it can accommodate realistic equations of state and arbitrary (pressure- or temperature-based) reaction-rate laws. Under certain conditions, we show that the formulation reduces to well-known augmented Euler or multi-phase models, which have been extensively validated and used in practice. The full hybrid model and its reduced forms are validated against problems with exact (or independently-verified numerical) solutions and evaluated for robustness for rate-stick and shock-induced cavity collapse case-studies.
NASA Astrophysics Data System (ADS)
Jones, Marvin Quenten, Jr.
The motion and behavior of quantum processes can be described by the Schrodinger equation using the wave function, Psi(x,t). The use of the Schrodinger equation to study quantum phenomena is known as Quantum Mechanics, akin to classical mechanics being the tool to study classical physics. This research focuses on the emphasis of numerical techniques: Finite-Difference, Fast Fourier Transform (spectral method), finite difference schemes such as the Leapfrog method and the Crank-Nicolson scheme and second quantization to solve and analyze the Schrodinger equation for the infinite square well problem, the free particle with periodic boundary conditions, the barrier problem, tight-binding hamiltonians and a potential wall problem. We discuss these techniques and the problems created to test how these different techniques draw both physical and numerical conclusions in a tabular summary. We observed both numerical stability and quantum stability (conservation of energy, probability, momentum, etc.). We found in our results that the Crank-Nicolson scheme is an unconditionally stable scheme and conserves probability (unitary), and momentum, though dissipative with energy. The time-independent problems conserved energy, momentum and were unitary, which is of interest, but we found when time-dependence was introduced, quantum stability (i.e. conservation of mass, momentum, etc.) was not implied by numerical stability. Hence, we observed schemes that were numerically stable, but not quantum stable as well as schemes that were quantum stable, but not numerically stable for all of time, t. We also observed that second quantization removed the issues with stability as the problem was transformed into a discrete problem. Moreover, all quantum information is conserved in second quantization. This method, however, does not work universally for all problems.
Hybrid recommendation methods in complex networks
NASA Astrophysics Data System (ADS)
Fiasconaro, A.; Tumminello, M.; Nicosia, V.; Latora, V.; Mantegna, R. N.
2015-07-01
We propose two recommendation methods, based on the appropriate normalization of already existing similarity measures, and on the convex combination of the recommendation scores derived from similarity between users and between objects. We validate the proposed measures on three data sets, and we compare the performance of our methods to other recommendation systems recently proposed in the literature. We show that the proposed similarity measures allow us to attain an improvement of performances of up to 20% with respect to existing nonparametric methods, and that the accuracy of a recommendation can vary widely from one specific bipartite network to another, which suggests that a careful choice of the most suitable method is highly relevant for an effective recommendation on a given system. Finally, we study how an increasing presence of random links in the network affects the recommendation scores, finding that one of the two recommendation algorithms introduced here can systematically outperform the others in noisy data sets.
Thomas, L.D.; Alexander, M.H.; Johnson, B.R.; Lester Jr., W. A.; Light, J.C.; McLenithan, K.D.; Parker, G.A.; Redmon, M.J.; Schmalz, T.G.; Secrest, D.; Walker, R.B.
1980-07-01
The numerical solution of coupled, second-order differential equations is a fundamental problem in theoretical physics and chemistry. There are presently over 20 commonly used methods. Unbiased comparisons of the methods are difficult to make and few have been attempted. Here we report a comparison of 11 different methods applied to 3 different test problems. The test problems have been constructed to approximate chemical systems of current research interest and to be representative of the state of the art in inelastic molecular collisions. All calculations were done on the same computer and the attempt was made to do all calculations to the same level of accuracy. The results of the initial tests indicated that an improved method might be obtained by using different methods in different integration regions. Such a hybrid program was developed and found to be at least 1.5 to 2.0 times faster than any individual method.
Numerical modeling of magnetic induction tomography using the impedance method.
Ramos, Airton; Wolff, Julia G B
2011-02-01
This article discusses the impedance method in the forward calculation in magnetic induction tomography (MIT). Magnetic field and eddy current distributions were obtained numerically for a sphere in the field of a coil and were compared with an analytical model. Additionally, numerical and experimental results for phase sensitivity in MIT were obtained and compared for a cylindrical object in a planar array of sensors. The results showed that the impedance method provides results that agree very well with reality in the frequency range from 100 kHz to 20 MHz and for low conductivity objects (10 S/m or less). This opens the possibility of using this numerical approach in image reconstruction in MIT. PMID:21229327
Extending a Hybrid Godunov Method for Radiation Hydrodynamics to Multiple Dimensions
NASA Astrophysics Data System (ADS)
Sekora, Michael D.
This paper presents a hybrid Godunov method for three-dimensional radiation hydrodynamics. The multidimensional technique outlined in this paper is an extension of the one-dimensional method that was developed by Sekora and Stone 2009, 2010. The earlier one-dimensional technique was shown to preserve certain asymptotic limits and be uniformly well behaved from the photon free streaming (hyperbolic) limit through the weak equilibrium diffusion (parabolic) limit and to the strong equilibrium diffusion (hyperbolic) limit. This paper gives the algorithmic details for constructing a multidimensional method. A future paper will present numerical tests that demonstrate the robustness of the computational technique across a wide-range of parameter space.
Wang, Sen; Feng, Qihong; Han, Xiaodong
2013-01-01
Due to the long-term fluid-solid interactions in waterflooding, the tremendous variation of oil reservoir formation parameters will lead to the widespread evolution of preferential flow paths, thereby preventing the further enhancement of recovery efficiency because of unstable fingering and premature breakthrough. To improve oil recovery, the characterization of preferential flow paths is essential and imperative. In efforts that have been previously documented, fluid flow characteristics within preferential paths are assumed to obey Darcy's equation. However, the occurrence of non-Darcy flow behavior has been increasingly suggested. To examine this conjecture, the Forchheimer number with the inertial coefficient estimated from different empirical formulas is applied as the criterion. Considering a 10% non-Darcy effect, the fluid flow in a preferential path may do experience non-Darcy behavior. With the objective of characterizing the preferential path with non-Darcy flow, a hybrid analytical/numerical model has been developed to investigate the pressure transient response, which dynamically couples a numerical model describing the non-Darcy effect of a preferential flow path with an analytical reservoir model. The characteristics of the pressure transient behavior and the sensitivities of corresponding parameters have also been discussed. In addition, an interpretation approach for pressure transient testing is also proposed, in which the Gravitational Search Algorithm is employed as a non-linear regression technology to match measured pressure with this hybrid model. Examples of applications from different oilfields are also presented to illustrate this method. This cost-effective approach provides more accurate characterization of a preferential flow path with non-Darcy flow, which will lay a solid foundation for the design and operation of conformance control treatments, as well as several other Enhanced Oil Recovery projects. PMID:24386224
Hybrid radical energy storage device and method of making
Gennett, Thomas; Ginley, David S.; Braunecker, Wade; Ban, Chunmei; Owczarczyk, Zbyslaw
2016-04-26
Hybrid radical energy storage devices, such as batteries or electrochemical devices, and methods of use and making are disclosed. Also described herein are electrodes and electrolytes useful in energy storage devices, for example, radical polymer cathode materials and electrolytes for use in organic radical batteries.
Hybrid radical energy storage device and method of making
Gennett, Thomas; Ginley, David S; Braunecker, Wade; Ban, Chunmei; Owczarczyk, Zbyslaw
2015-01-27
Hybrid radical energy storage devices, such as batteries or electrochemical devices, and methods of use and making are disclosed. Also described herein are electrodes and electrolytes useful in energy storage devices, for example, radical polymer cathode materials and electrolytes for use in organic radical batteries.
Teaching with the Case Method Online: Pure versus Hybrid Approaches
ERIC Educational Resources Information Center
Webb, Harold W.; Gill, Grandon; Poe, Gary
2005-01-01
The impact of hybrid classroom/distance education approaches is examined in the context of the case method. Four distinct semester-long treatments, which varied mixes of classroom and online discussion, were used to teach a graduate MIS survey course. Specific findings suggest that by using Web technology, college instructors may offer students…
3D magnetospheric parallel hybrid multi-grid method applied to planet-plasma interactions
NASA Astrophysics Data System (ADS)
Leclercq, L.; Modolo, R.; Leblanc, F.; Hess, S.; Mancini, M.
2016-03-01
We present a new method to exploit multiple refinement levels within a 3D parallel hybrid model, developed to study planet-plasma interactions. This model is based on the hybrid formalism: ions are kinetically treated whereas electrons are considered as a inertia-less fluid. Generally, ions are represented by numerical particles whose size equals the volume of the cells. Particles that leave a coarse grid subsequently entering a refined region are split into particles whose volume corresponds to the volume of the refined cells. The number of refined particles created from a coarse particle depends on the grid refinement rate. In order to conserve velocity distribution functions and to avoid calculations of average velocities, particles are not coalesced. Moreover, to ensure the constancy of particles' shape function sizes, the hybrid method is adapted to allow refined particles to move within a coarse region. Another innovation of this approach is the method developed to compute grid moments at interfaces between two refinement levels. Indeed, the hybrid method is adapted to accurately account for the special grid structure at the interfaces, avoiding any overlapping grid considerations. Some fundamental test runs were performed to validate our approach (e.g. quiet plasma flow, Alfven wave propagation). Lastly, we also show a planetary application of the model, simulating the interaction between Jupiter's moon Ganymede and the Jovian plasma.
NASA Astrophysics Data System (ADS)
Bellos, Vasilis; Tsakiris, George
2016-09-01
The study presents a new hybrid method for the simulation of flood events in small catchments. It combines a physically-based two-dimensional hydrodynamic model and the hydrological unit hydrograph theory. Unit hydrographs are derived using the FLOW-R2D model which is based on the full form of two-dimensional Shallow Water Equations, solved by a modified McCormack numerical scheme. The method is tested at a small catchment in a suburb of Athens-Greece for a storm event which occurred in February 2013. The catchment is divided into three friction zones and unit hydrographs of 15 and 30 min are produced. The infiltration process is simulated by the empirical Kostiakov equation and the Green-Ampt model. The results from the implementation of the proposed hybrid method are compared with recorded data at the hydrometric station at the outlet of the catchment and the results derived from the fully hydrodynamic model FLOW-R2D. It is concluded that for the case studied, the proposed hybrid method produces results close to those of the fully hydrodynamic simulation at substantially shorter computational time. This finding, if further verified in a variety of case studies, can be useful in devising effective hybrid tools for the two-dimensional flood simulations, which are lead to accurate and considerably faster results than those achieved by the fully hydrodynamic simulations.
Hybrid modeling of spatial continuity for application to numerical inverse problems
Friedel, Michael J.; Iwashita, Fabio
2013-01-01
A novel two-step modeling approach is presented to obtain optimal starting values and geostatistical constraints for numerical inverse problems otherwise characterized by spatially-limited field data. First, a type of unsupervised neural network, called the self-organizing map (SOM), is trained to recognize nonlinear relations among environmental variables (covariates) occurring at various scales. The values of these variables are then estimated at random locations across the model domain by iterative minimization of SOM topographic error vectors. Cross-validation is used to ensure unbiasedness and compute prediction uncertainty for select subsets of the data. Second, analytical functions are fit to experimental variograms derived from original plus resampled SOM estimates producing model variograms. Sequential Gaussian simulation is used to evaluate spatial uncertainty associated with the analytical functions and probable range for constraining variables. The hybrid modeling of spatial continuity is demonstrated using spatially-limited hydrologic measurements at different scales in Brazil: (1) physical soil properties (sand, silt, clay, hydraulic conductivity) in the 42 km2 Vargem de Caldas basin; (2) well yield and electrical conductivity of groundwater in the 132 km2 fractured crystalline aquifer; and (3) specific capacity, hydraulic head, and major ions in a 100,000 km2 transboundary fractured-basalt aquifer. These results illustrate the benefits of exploiting nonlinear relations among sparse and disparate data sets for modeling spatial continuity, but the actual application of these spatial data to improve numerical inverse modeling requires testing.
SRYTH: A New Yeast Two-Hybrid Method.
Mallick, Jaideep; Jansen, Gregor; Wu, Cunle; Whiteway, Malcolm
2016-01-01
Many biological processes are regulated by protein-protein interactions, and the analysis of these interactions has been a productive endeavor contributing to our understanding of cellular organization and function. The yeast two-hybrid technique is a widely used, powerful method of analyzing protein-protein interactions. The currently used formats, however, have inherent limitations, providing an opportunity to develop new alternatives that extend our ability to detect protein-protein interactions of biological relevance. Here we present a two-hybrid system named SRYTH (Ste11p/Ste50p related yeast two-hybrid) based on the Ste11p/Ste50p interaction that uses the activation of the HOG pathway of Saccharomyces cerevisiae as a reporter for interactions. The system is suitable for detecting cytoplasmic protein interactions in their natural subcellular environment, and has been successfully used to investigate protein-protein interactions, including transcription factor associations, in Candida albicans. PMID:26519063
Hybrid methods using genetic algorithms for global optimization.
Renders, J M; Flasse, S P
1996-01-01
This paper discusses the trade-off between accuracy, reliability and computing time in global optimization. Particular compromises provided by traditional methods (Quasi-Newton and Nelder-Mead's simplex methods) and genetic algorithms are addressed and illustrated by a particular application in the field of nonlinear system identification. Subsequently, new hybrid methods are designed, combining principles from genetic algorithms and "hill-climbing" methods in order to find a better compromise to the trade-off. Inspired by biology and especially by the manner in which living beings adapt themselves to their environment, these hybrid methods involve two interwoven levels of optimization, namely evolution (genetic algorithms) and individual learning (Quasi-Newton), which cooperate in a global process of optimization. One of these hybrid methods appears to join the group of state-of-the-art global optimization methods: it combines the reliability properties of the genetic algorithms with the accuracy of Quasi-Newton method, while requiring a computation time only slightly higher than the latter. PMID:18263027
COMPARING NUMERICAL METHODS FOR ISOTHERMAL MAGNETIZED SUPERSONIC TURBULENCE
Kritsuk, Alexei G.; Collins, David; Norman, Michael L.; Xu Hao E-mail: dccollins@lanl.gov
2011-08-10
Many astrophysical applications involve magnetized turbulent flows with shock waves. Ab initio star formation simulations require a robust representation of supersonic turbulence in molecular clouds on a wide range of scales imposing stringent demands on the quality of numerical algorithms. We employ simulations of supersonic super-Alfvenic turbulence decay as a benchmark test problem to assess and compare the performance of nine popular astrophysical MHD methods actively used to model star formation. The set of nine codes includes: ENZO, FLASH, KT-MHD, LL-MHD, PLUTO, PPML, RAMSES, STAGGER, and ZEUS. These applications employ a variety of numerical approaches, including both split and unsplit, finite difference and finite volume, divergence preserving and divergence cleaning, a variety of Riemann solvers, and a range of spatial reconstruction and time integration techniques. We present a comprehensive set of statistical measures designed to quantify the effects of numerical dissipation in these MHD solvers. We compare power spectra for basic fields to determine the effective spectral bandwidth of the methods and rank them based on their relative effective Reynolds numbers. We also compare numerical dissipation for solenoidal and dilatational velocity components to check for possible impacts of the numerics on small-scale density statistics. Finally, we discuss the convergence of various characteristics for the turbulence decay test and the impact of various components of numerical schemes on the accuracy of solutions. The nine codes gave qualitatively the same results, implying that they are all performing reasonably well and are useful for scientific applications. We show that the best performing codes employ a consistently high order of accuracy for spatial reconstruction of the evolved fields, transverse gradient interpolation, conservation law update step, and Lorentz force computation. The best results are achieved with divergence-free evolution of the
A hybrid encapsulation method for organic electronics
Kim, N.; Graham, S.; Potscavage, W. J. Jr.; Domercq, B.; Kippelen, B.
2009-04-20
We report a thin-film encapsulation method for organic electronics that combines the deposition of a layer of SiO{sub x} or SiN{sub x} (100 nm) by plasma enhanced chemical vapor deposition followed by a layer of Al{sub 2}O{sub 3} (10-50 nm) by atomic layer deposition and a 1-{mu}m-thick layer of parylene by chemical vapor deposition. The effective water vapor transmission rates of the encapsulation was (2{+-}1)x10{sup -5} g/m{sup 2} day at 20 deg. C and 50% relative humidity (RH). The encapsulation was integrated with pentacene/C{sub 60} solar cells, which showed no decrease in conversion efficiency after 5800 h of exposure to air demonstrating the effectiveness of the encapsulation methodology.
Atomistic hybrid DSMC/NEMD method for nonequilibrium multiscale simulations
Gu Kai; Watkins, Charles B. Koplik, Joel
2010-03-01
A multiscale hybrid method for coupling the direct simulation Monte Carlo (DSMC) method to the nonequilibrium molecular dynamics (NEMD) method is introduced. The method addresses Knudsen layer type gas flows within a few mean free paths of an interface or about an object with dimensions of the order of a few mean free paths. It employs the NEMD method to resolve nanoscale phenomena closest to the interface along with coupled DSMC simulation of the remainder of the Knudsen layer. The hybrid DSMC/NEMD method is a particle based algorithm without a buffer zone. It incorporates a new, modified generalized soft sphere (MGSS) molecular collision model to improve the poor computational efficiency of the traditional generalized soft sphere GSS model and to achieve DSMC compatibility with Lennard-Jones NEMD molecular interactions. An equilibrium gas, a Fourier thermal flow, and an oscillatory Couette flow, are simulated to validate the method. The method shows good agreement with Maxwell-Boltzmann theory for the equilibrium system, Chapman-Enskog theory for Fourier flow, and pure DSMC simulations for oscillatory Couette flow. Speedup in CPU time of the hybrid solver is benchmarked against a pure NEMD solver baseline for different system sizes and solver domain partitions. Finally, the hybrid method is applied to investigate interaction of argon gas with solid surface molecules in a parametric study of the influence of wetting effects and solid molecular mass on energy transfer and thermal accommodation coefficients. It is determined that wetting effect strength and solid molecular mass have a significant impact on the energy transfer between gas and solid phases and thermal accommodation coefficient.
A novel gas-droplet numerical method for spray combustion
NASA Technical Reports Server (NTRS)
Chen, C. P.; Shang, H. M.; Jiang, Y.
1991-01-01
This paper presents a non-iterative numerical technique for computing time-dependent gas-droplet flows. The method is a fully-interacting combination of Eulerian fluid and Lagrangian particle calculation. The interaction calculations between the two phases are formulated on a pressure-velocity coupling procedure based on the operator-splitting technique. This procedure eliminates the global iterations required in the conventional particle-source-in-cell (PSIC) procedure. Turbulent dispersion calculations are treated by a stochastic procedure. Numerical calculations and comparisons with available experimental data, as well as efficiency assessments are given for some sprays typical of spray combustion applications.
The TAB method for numerical calculation of spray droplet breakup
NASA Astrophysics Data System (ADS)
Orourke, P. J.; Amsden, A. A.
A short history is given of the major milestones in the development of the stochastic particle method for calculating liquid fuel sprays. The most recent advance has been the discovery of the importance of drop breakup in engine sprays. A new method, called TAB, for calculating drop breakup is presented. Some theoretical properties of the method are derived; its numerical implementation in the computer program KIVA is described; and comparisons are presented between TAB-method calculations and experiments and calculations using another breakup model.
Simple numerical method for predicting steady compressible flows
NASA Technical Reports Server (NTRS)
Vonlavante, Ernst; Nelson, N. Duane
1986-01-01
A numerical method for solving the isenthalpic form of the governing equations for compressible viscous and inviscid flows was developed. The method was based on the concept of flux vector splitting in its implicit form. The method was tested on several demanding inviscid and viscous configurations. Two different forms of the implicit operator were investigated. The time marching to steady state was accelerated by the implementation of the multigrid procedure. Its various forms very effectively increased the rate of convergence of the present scheme. High quality steady state results were obtained in most of the test cases; these required only short computational times due to the relative efficiency of the basic method.
Singularity Preserving Numerical Methods for Boundary Integral Equations
NASA Technical Reports Server (NTRS)
Kaneko, Hideaki (Principal Investigator)
1996-01-01
In the past twelve months (May 8, 1995 - May 8, 1996), under the cooperative agreement with Division of Multidisciplinary Optimization at NASA Langley, we have accomplished the following five projects: a note on the finite element method with singular basis functions; numerical quadrature for weakly singular integrals; superconvergence of degenerate kernel method; superconvergence of the iterated collocation method for Hammersteion equations; and singularity preserving Galerkin method for Hammerstein equations with logarithmic kernel. This final report consists of five papers describing these projects. Each project is preceeded by a brief abstract.
The development of accurate and efficient methods of numerical quadrature
NASA Technical Reports Server (NTRS)
Feagin, T.
1973-01-01
Some new methods for performing numerical quadrature of an integrable function over a finite interval are described. Each method provides a sequence of approximations of increasing order to the value of the integral. Each approximation makes use of all previously computed values of the integrand. The points at which new values of the integrand are computed are selected in such a way that the order of the approximation is maximized. The methods are compared with the quadrature methods of Clenshaw and Curtis, Gauss, Patterson, and Romberg using several examples.
NASA Technical Reports Server (NTRS)
Fahrenthold, Eric P.; Shivarama, Ravishankar
2004-01-01
The hybrid particle-finite element method of Fahrenthold and Horban, developed for the simulation of hypervelocity impact problems, has been extended to include new formulations of the particle-element kinematics, additional constitutive models, and an improved numerical implementation. The extended formulation has been validated in three dimensional simulations of published impact experiments. The test cases demonstrate good agreement with experiment, good parallel speedup, and numerical convergence of the simulation results.
A hybrid experimental-numerical technique for determining 3D velocity fields from planar 2D PIV data
NASA Astrophysics Data System (ADS)
Eden, A.; Sigurdson, M.; Mezić, I.; Meinhart, C. D.
2016-09-01
Knowledge of 3D, three component velocity fields is central to the understanding and development of effective microfluidic devices for lab-on-chip mixing applications. In this paper we present a hybrid experimental-numerical method for the generation of 3D flow information from 2D particle image velocimetry (PIV) experimental data and finite element simulations of an alternating current electrothermal (ACET) micromixer. A numerical least-squares optimization algorithm is applied to a theory-based 3D multiphysics simulation in conjunction with 2D PIV data to generate an improved estimation of the steady state velocity field. This 3D velocity field can be used to assess mixing phenomena more accurately than would be possible through simulation alone. Our technique can also be used to estimate uncertain quantities in experimental situations by fitting the gathered field data to a simulated physical model. The optimization algorithm reduced the root-mean-squared difference between the experimental and simulated velocity fields in the target region by more than a factor of 4, resulting in an average error less than 12% of the average velocity magnitude.
A numerical method for interface problems in elastodynamics
NASA Technical Reports Server (NTRS)
Mcghee, D. S.
1984-01-01
The numerical implementation of a formulation for a class of interface problems in elastodynamics is discussed. This formulation combines the use of the finite element and boundary integral methods to represent the interior and the exteriro regions, respectively. In particular, the response of a semicylindrical alluvial valley in a homogeneous halfspace to incident antiplane SH waves is considered to determine the accuracy and convergence of the numerical procedure. Numerical results are obtained from several combinations of the incidence angle, frequency of excitation, and relative stiffness between the inclusion and the surrounding halfspace. The results tend to confirm the theoretical estimates that the convergence is of the order H(2) for the piecewise linear elements used. It was also observed that the accuracy descreases as the frequency of excitation increases or as the relative stiffness of the inclusion decreases.
A hybrid particle-mesh method for incompressible active polar viscous gels
NASA Astrophysics Data System (ADS)
Ramaswamy, Rajesh; Bourantas, George; Jülicher, Frank; Sbalzarini, Ivo F.
2015-06-01
We present a hybrid particle-mesh method for numerically solving the hydrodynamic equations of incompressible active polar viscous gels. These equations model the dynamics of polar active agents, embedded in a viscous medium, in which stresses are induced through constant consumption of energy. The numerical method is based on Lagrangian particles and staggered Cartesian finite-difference meshes. We show that the method is second-order and first-order accurate with respect to grid and time-step sizes, respectively. Using the present method, we simulate the hydrodynamics in rectangular geometries, of a passive liquid crystal, of an active polar film and of active gels with topological defects in polarization. We show the emergence of spontaneous flow due to Fréedericksz transition, and transformation in the nature of topological defects by tuning the activity of the system.
NASA Astrophysics Data System (ADS)
Tong, Ping; Komatitsch, Dimitri; Tseng, Tai-Lin; Hung, Shu-Huei; Chen, Chin-Wu; Basini, Piero; Liu, Qinya
2014-10-01
We present a three-dimensional (3-D) hybrid method that interfaces the spectral-element method (SEM) with the frequency-wave number (FK) technique to model the propagation of teleseismic plane waves beneath seismic arrays. The accuracy of the resulting 3-D SEM-FK hybrid method is benchmarked against semianalytical FK solutions for 1-D models. The accuracy of 2.5-D modeling based on 2-D SEM-FK hybrid method is also investigated through comparisons to this 3-D hybrid method. Synthetic examples for structural models of the Alaska subduction zone and the central Tibet crust show that this method is capable of accurately capturing interactions between incident plane waves and local heterogeneities. This hybrid method presents an essential tool for the receiver function and scattering imaging community to verify and further improve their techniques. These numerical examples also show the promising future of the 3-D SEM-FK hybrid method in high-resolution regional seismic imaging based on waveform inversions of converted/scattered waves recorded by seismic array.
High-resolution seismic array imaging based on an SEM-FK hybrid method
NASA Astrophysics Data System (ADS)
Tong, Ping; Chen, Chin-wu; Komatitsch, Dimitri; Basini, Piero; Liu, Qinya
2014-04-01
We demonstrate the feasibility of high-resolution seismic array imaging based on teleseismic recordings using full numerical wave simulations. We develop a hybrid method that interfaces a frequency-wavenumber (FK) calculation, which provides analytical solutions to 1-D layered background models with a spectral-element (SEM) numerical solver to calculate synthetic responses of local media to plane-wave incidence.This hybrid method accurately deals with local heterogeneities and discontinuity undulations, and represents an efficient tool for the forward modelling of teleseismic coda (including converted and scattered) waves. We benchmark the accuracy of the SEM-FK hybrid method against FK solutions for 1-D media. We then compute sensitivity kernels for teleseismic coda waves by interacting the forward teleseismic waves with an adjoint wavefield, produced by injecting coda waves as adjoint sources, based on adjoint techniques. These sensitivity kernels provide the basis for mapping variations in subsurface discontinuities, density and velocity structures through non-linear conjugate-gradient methods. We illustrate various synthetic imaging experiments, including discontinuity characterization, volumetric structural inversion for the crust or subduction zones. These tests show that using pre-conditioners based upon the scaled product of sensitivity kernels for different phases, combining finite-frequency traveltime and waveform inversion, and/or adopting hierarchical inversions from long- to short-period waveforms could reduce the non-linearity of the seismic inverse problem and speed up its convergence. The encouraging results of these synthetic examples suggest that inversion of teleseismic coda phases based on the SEM-FK hybrid method and adjoint techniques is a promising tool for structural imaging beneath dense seismic arrays.
Numerical method to compute acoustic scattering effect of a moving source.
Song, Hao; Yi, Mingxu; Huang, Jun; Pan, Yalin; Liu, Dawei
2016-01-01
In this paper, the aerodynamic characteristic of a ducted tail rotor in hover has been numerically studied using CFD method. An analytical time domain formulation based on Ffowcs Williams-Hawkings (FW-H) equation is derived for the prediction of the acoustic velocity field and used as Neumann boundary condition on a rigid scattering surface. In order to predict the aerodynamic noise, a hybrid method combing computational aeroacoustics with an acoustic thin-body boundary element method has been proposed. The aerodynamic results and the calculated sound pressure levels (SPLs) are compared with the known method for validation. Simulation results show that the duct can change the value of SPLs and the sound directivity. Compared with the isolate tail rotor, the SPLs of the ducted tail rotor are smaller at certain azimuth. PMID:27610323
Numerical simulation methods for the Rouse model in flow
NASA Astrophysics Data System (ADS)
Howard, Michael P.; Milner, Scott T.
2011-11-01
Simulation of the Rouse model in flow underlies a great variety of numerical investigations of polymer dynamics, in both entangled melts and solutions and in dilute solution. Typically a simple explicit stochastic Euler method is used to evolve the Rouse model. Here we compare this approach to an operator splitting method, which splits the evolution operator into stochastic linear and deterministic nonlinear parts and takes advantage of an analytical solution for the linear Rouse model in terms of the noise history. We show that this splitting method has second-order weak convergence, whereas the Euler method has only first-order weak convergence. Furthermore, the splitting method is unconditionally stable, in contrast to the limited stability range of the Euler method. Similar splitting methods are applicable to a broad class of problems in stochastic dynamics in which noise competes with ordering and flow to determine steady-state order parameter structures.
Hybrid System Reduction Method using Reduced System Regulator
NASA Astrophysics Data System (ADS)
Kuroda, Eisuke; Tsuji, Takao; Oyama, Tsutomu
In order to analyze transient stability of large-scale power systems, it is advantageous to apply system reduction method to external systems. Short-circuit current method is one of the typical engineering reduction techniques. However, the dominant eigenvalues are not necessarily conserved in the reduced system. Therefore, the hybrid reduction method in which controller parameters are adjusted to conserve the dominant eigenvalues was proposed. Automatic voltage regulator (AVR) and power system stabilizer (PSS) have been used for parameters adjustment so far. However, since there are many parameters in AVR and PSS, complicated procedures are required to adjust them. Therefore, in this paper, the reduced system regulator (RSR) is proposed for hybrid system reduction method. The RSR has only two parameters for adjustment. It is easier to adjust the RSR than AVR/PSS. In addition, the initial gains of the RSR are set zero so that dynamic behavior of the system is not influenced before the adjustment. The effect and the accuracy of the hybrid system reduction method with RSR are examined using a typical longitudinal power system, IEEJ WEST 10-machine system model.
Numerical Polynomial Homotopy Continuation Method and String Vacua
Mehta, Dhagash
2011-01-01
Finding vmore » acua for the four-dimensional effective theories for supergravity which descend from flux compactifications and analyzing them according to their stability is one of the central problems in string phenomenology. Except for some simple toy models, it is, however, difficult to find all the vacua analytically. Recently developed algorithmic methods based on symbolic computer algebra can be of great help in the more realistic models. However, they suffer from serious algorithmic complexities and are limited to small system sizes. In this paper, we review a numerical method called the numerical polynomial homotopy continuation (NPHC) method, first used in the areas of lattice field theories, which by construction finds all of the vacua of a given potential that is known to have only isolated solutions. The NPHC method is known to suffer from no major algorithmic complexities and is embarrassingly parallelizable , and hence its applicability goes way beyond the existing symbolic methods. We first solve a simple toy model as a warm-up example to demonstrate the NPHC method at work. We then show that all the vacua of a more complicated model of a compactified M theory model, which has an S U ( 3 ) structure, can be obtained by using a desktop machine in just about an hour, a feat which was reported to be prohibitively difficult by the existing symbolic methods. Finally, we compare the various technicalities between the two methods.« less
Projected discrete ordinates methods for numerical transport problems
Larsen, E.W.
1985-01-01
A class of Projected Discrete-Ordinates (PDO) methods is described for obtaining iterative solutions of discrete-ordinates problems with convergence rates comparable to those observed using Diffusion Synthetic Acceleration (DSA). The spatially discretized PDO solutions are generally not equal to the DSA solutions, but unlike DSA, which requires great care in the use of spatial discretizations to preserve stability, the PDO solutions remain stable and rapidly convergent with essentially arbitrary spatial discretizations. Numerical results are presented which illustrate the rapid convergence and the accuracy of solutions obtained using PDO methods with commonplace differencing methods.
Computational methods for aerodynamic design using numerical optimization
NASA Technical Reports Server (NTRS)
Peeters, M. F.
1983-01-01
Five methods to increase the computational efficiency of aerodynamic design using numerical optimization, by reducing the computer time required to perform gradient calculations, are examined. The most promising method consists of drastically reducing the size of the computational domain on which aerodynamic calculations are made during gradient calculations. Since a gradient calculation requires the solution of the flow about an airfoil whose geometry was slightly perturbed from a base airfoil, the flow about the base airfoil is used to determine boundary conditions on the reduced computational domain. This method worked well in subcritical flow.
Numerical methods for characterization of synchrotron radiation based on the Wigner function method
NASA Astrophysics Data System (ADS)
Tanaka, Takashi
2014-06-01
Numerical characterization of synchrotron radiation based on the Wigner function method is explored in order to accurately evaluate the light source performance. A number of numerical methods to compute the Wigner functions for typical synchrotron radiation sources such as bending magnets, undulators and wigglers, are presented, which significantly improve the computation efficiency and reduce the total computation time. As a practical example of the numerical characterization, optimization of betatron functions to maximize the brilliance of undulator radiation is discussed.
NASA Astrophysics Data System (ADS)
Anis, Fatima; Lou, Yang; Conjusteau, André; Su, Richard; Oruganti, Tanmayi; Ermilov, Sergey A.; Oraevsky, Alexander A.; Anastasio, Mark A.
2014-03-01
In this work, we investigate a novel reconstruction method for laser-induced ultrasound computed tomography (USCT) breast imaging that circumvents limitations of existing methods that rely on ray-tracing. There is currently great interest in developing hybrid imaging systems that combine optoacoustic tomography (OAT) and USCT. There are two primary motivations for this: (1) the speed-of-sound (SOS) distribution reconstructed by USCT can provide complementary diagnostic information; and (2) the reconstructed SOS distribution can be incorporated in the OAT reconstruction algorithm to improve OAT image quality. However, image reconstruction in USCT remains challenging. The majority of existing approaches for USCT breast imaging involve ray-tracing to establish the imaging operator. This process is cumbersome and can lead to inaccuracies in the reconstructed SOS images in the presence of multiple ray-paths and/or shadow zones. To circumvent these problems, we implemented a partial differential equation-based Eulerian approach to USCT that was proposed in the mathematics literature but never investigated for medical imaging applications. This method operates by directly inverting the Eikonal equation without ray-tracing. A numerical implementation of this method was developed and compared to existing reconstruction methods for USCT breast imaging. We demonstrated the ability of the new method to reconstruct SOS maps from TOF data obtained by a hybrid OAT/USCT imager built by our team.
A robust method for handling low density regions in hybrid simulations for collisionless plasmas
Amano, Takanobu Higashimori, Katsuaki; Shirakawa, Keisuke
2014-10-15
A robust method to handle vacuum and near vacuum regions in hybrid simulations for space and astrophysical plasmas is presented. The conventional hybrid simulation model dealing with kinetic ions and a massless charge-neutralizing electron fluid is known to be susceptible to numerical instability due to divergence of the whistler-mode wave dispersion, as well as division-by-density operation in regions of low density. Consequently, a pure vacuum region is not allowed to exist in the simulation domain unless some ad hoc technique is used. To resolve this difficulty, an alternative way to introduce finite electron inertia effect is proposed. Contrary to the conventional method, the proposed one introduces a correction to the electric field rather than the magnetic field. It is shown that the generalized Ohm's law correctly reduces to Laplace's equation in a vacuum which therefore does not involve any numerical problems. In addition, a variable ion-to-electron mass ratio is introduced to reduce the phase velocity of high frequency whistler waves at low density regions so that the stability condition is always satisfied. It is demonstrated that the proposed model is able to handle near vacuum regions generated as a result of nonlinear self-consistent development of the system, as well as pure vacuum regions set up at the initial condition, without losing the advantages of the standard hybrid code.
Fast and stable numerical method for neuronal modelling
NASA Astrophysics Data System (ADS)
Hashemi, Soheil; Abdolali, Ali
2016-11-01
Excitable cell modelling is of a prime interest in predicting and targeting neural activity. Two main limits in solving related equations are speed and stability of numerical method. Since there is a tradeoff between accuracy and speed, most previously presented methods for solving partial differential equations (PDE) are focused on one side. More speed means more accurate simulations and therefore better device designing. By considering the variables in finite differenced equation in proper time and calculating the unknowns in the specific sequence, a fast, stable and accurate method is introduced in this paper for solving neural partial differential equations. Propagation of action potential in giant axon is studied by proposed method and traditional methods. Speed, consistency and stability of the methods are compared and discussed. The proposed method is as fast as forward methods and as stable as backward methods. Forward methods are known as fastest methods and backward methods are stable in any circumstances. Complex structures can be simulated by proposed method due to speed and stability of the method.
Automatic numerical integration methods for Feynman integrals through 3-loop
NASA Astrophysics Data System (ADS)
de Doncker, E.; Yuasa, F.; Kato, K.; Ishikawa, T.; Olagbemi, O.
2015-05-01
We give numerical integration results for Feynman loop diagrams through 3-loop such as those covered by Laporta [1]. The methods are based on automatic adaptive integration, using iterated integration and extrapolation with programs from the QUADPACK package, or multivariate techniques from the ParInt package. The Dqags algorithm from QuadPack accommodates boundary singularities of fairly general types. PARINT is a package for multivariate integration layered over MPI (Message Passing Interface), which runs on clusters and incorporates advanced parallel/distributed techniques such as load balancing among processes that may be distributed over a network of nodes. Results are included for 3-loop self-energy diagrams without IR (infra-red) or UV (ultra-violet) singularities. A procedure based on iterated integration and extrapolation yields a novel method of numerical regularization for integrals with UV terms, and is applied to a set of 2-loop self-energy diagrams with UV singularities.
Numerical method for shear bands in ductile metal with inclusions
Plohr, Jee Yeon N; Plohr, Bradley J
2010-01-01
A numerical method for mesoscale simulation of high strain-rate loading of ductile metal containing inclusions is described. Because of small-scale inhomogeneities, such a composite material is prone to localized shear deformation (adiabatic shear bands). The modeling framework is the Generalized Method of Cells of Paley and Aboudi [Mech. Materials, vol. 14, pp. /27-139, 1992], which ensures that the micromechanical response of the material is reflected in the behavior of the composite at the mesoscale. To calculate the effective plastic strain rate when shear bands are present, the analytic and numerical analysis of shear bands by Glimm, Plohr, and Sharp [Mech. Materials, vol. 24, pp. 31-41, 1996] is adapted and extended.
Numerical Method for the Astronomical Almanac and Orbit Calculations
NASA Astrophysics Data System (ADS)
Kim, Kap-Sung
1993-12-01
We have calculated the astronomical almanac 1994 and simulated the trajectory of a satellite orbit considering all perturbative forces with various initial conditions. In this work, Gauss Jackson multistep integration method has been used to calculate our basic equation of motion with high numerical accuracy. It has been found that our results agree well with the Astronomical Almanac Data distributed by JPL of NASA and the orbit simulations have been carried out with fast speed, stability and excellent round-off error accumulation, comparing with other numerical methods. In order to be carried out our works on almanac and orbit calculations easily by anyone who uses a personal computer, we have made a computer program on graphical user interface to provide various menus for detail works selected by a mouse.
Projected Hybrid Orbitals: A General QM/MM Method
2015-01-01
A projected hybrid orbital (PHO) method was described to model the covalent boundary in a hybrid quantum mechanical and molecular mechanical (QM/MM) system. The PHO approach can be used in ab initio wave function theory and in density functional theory with any basis set without introducing system-dependent parameters. In this method, a secondary basis set on the boundary atom is introduced to formulate a set of hybrid atomic orbtials. The primary basis set on the boundary atom used for the QM subsystem is projected onto the secondary basis to yield a representation that provides a good approximation to the electron-withdrawing power of the primary basis set to balance electronic interactions between QM and MM subsystems. The PHO method has been tested on a range of molecules and properties. Comparison with results obtained from QM calculations on the entire system shows that the present PHO method is a robust and balanced QM/MM scheme that preserves the structural and electronic properties of the QM region. PMID:25317748
Hybrid pathwise sensitivity methods for discrete stochastic models of chemical reaction systems
Wolf, Elizabeth Skubak; Anderson, David F.
2015-01-21
Stochastic models are often used to help understand the behavior of intracellular biochemical processes. The most common such models are continuous time Markov chains (CTMCs). Parametric sensitivities, which are derivatives of expectations of model output quantities with respect to model parameters, are useful in this setting for a variety of applications. In this paper, we introduce a class of hybrid pathwise differentiation methods for the numerical estimation of parametric sensitivities. The new hybrid methods combine elements from the three main classes of procedures for sensitivity estimation and have a number of desirable qualities. First, the new methods are unbiased for a broad class of problems. Second, the methods are applicable to nearly any physically relevant biochemical CTMC model. Third, and as we demonstrate on several numerical examples, the new methods are quite efficient, particularly if one wishes to estimate the full gradient of parametric sensitivities. The methods are rather intuitive and utilize the multilevel Monte Carlo philosophy of splitting an expectation into separate parts and handling each in an efficient manner.
A hybrid method for computing forces on curved dislocations threading to free surfaces
Tang, M; Cai, W; Xu, G; Bulatov, V V
2005-06-06
Dislocations threading to free surfaces present a challenge for numerical implementation of traction-free boundary conditions. The difficulty arises when canonical (singular) solutions of dislocation mechanics are used in combination with the Finite Element or Boundary Element (Green's function) methods. A new hybrid method is developed here in which the singular part and the non-singular (regular) part of the image stress are dealt with separately. A special analytical solution for a semi-infinite straight dislocation intersecting the surface of a half-space is used to account for the singular part of the image stress, while the remaining regular part of the image stress field is treated using the standard Finite Element Method. The numerical advantages of such regularization are demonstrated with examples.
NASA Technical Reports Server (NTRS)
Jin, Jian-Ming; Volakis, John L.
1992-01-01
A hybrid numerical technique is presented for a characterization of the scattering and radiation properties of three-dimensional cavity arrays recessed in a ground plane. The technique combines the finite element and boundary integral methods and invokes Floquet's representation to formulate a system of equations for the fields at the apertures and those inside the cavities. The system is solved via the conjugate gradient method in conjunction with the Fast Fourier Transform (FFT) thus achieving an O(N) storage requirement. By virtue of the finite element method, the proposed technique is applicable to periodic arrays comprised of cavities having arbitrary shape and filled with inhomogeneous dielectrics. Several numerical results are presented, along with new measured data, which demonstrate the validity, efficiency, and capability of the technique.
Asymptotic and Numerical Methods for Rapidly Rotating Buoyant Flow
NASA Astrophysics Data System (ADS)
Grooms, Ian G.
This thesis documents three investigations carried out in pursuance of a doctoral degree in applied mathematics at the University of Colorado (Boulder). The first investigation concerns the properties of rotating Rayleigh-Benard convection -- thermal convection in a rotating infinite plane layer between two constant-temperature boundaries. It is noted that in certain parameter regimes convective Taylor columns appear which dominate the dynamics, and a semi-analytical model of these is presented. Investigation of the columns and of various other properties of the flow is ongoing. The second investigation concerns the interactions between planetary-scale and mesoscale dynamics in the oceans. Using multiple-scale asymptotics the possible connections between planetary geostrophic and quasigeostrophic dynamics are investigated, and three different systems of coupled equations are derived. Possible use of these equations in conjunction with the method of superparameterization, and extension of the asymptotic methods to the interactions between mesoscale and submesoscale dynamics is ongoing. The third investigation concerns the linear stability properties of semi-implicit methods for the numerical integration of ordinary differential equations, focusing in particular on the linear stability of IMEX (Implicit-Explicit) methods and exponential integrators applied to systems of ordinary differential equations arising in the numerical solution of spatially discretized nonlinear partial differential equations containing both dispersive and dissipative linear terms. While these investigations may seem unrelated at first glance, some reflection shows that they are in fact closely linked. The investigation of rotating convection makes use of single-space, multiple-time-scale asymptotics to deal with dynamics strongly constrained by rotation. Although the context of thermal convection in an infinite layer seems somewhat removed from large-scale ocean dynamics, the asymptotic
Numerical simulation of HTPB combustion in a 2D hybrid slab combustor
NASA Astrophysics Data System (ADS)
Gariani, Gabriela; Maggi, Filippo; Galfetti, Luciano
2011-09-01
A code for the numerical simulation of combustion processes in hybrid rockets, developed at the Space Propulsion Laboratory of Politecnico di Milano (SPLab), is presented. The code deals with Navier-Stokes equations solved with RANS approach, blowing effect, combustion kinetics and radiation. The equations are closed with k-epsilon turbulence model and well stirred reactor model. The P1 model, a simplification of the PN radiation model, is adopted. Specific simulation tools were developed using OpenFOAM®open source technology. The computational domain is 2D and split in two subdomains, simulating the reacting gas mixture on one side and the solid fuel grain on the other. The interface between the two regions plays a key role as the solid grain pyrolysis comes from a straight solution of the model without shortcuts. A propellant combination with polybutadiene and gaseous oxygen has been chosen and a reduced kinetic model for combustion of butadiene, considered as the major gaseous constituent coming from polybutadiene pyrolysis, has been developed for reactions occurring in oxygen atmosphere. The computational domain tries to replicate the real experimental setup and is split into three areas: pre-chamber, slab zone and post-chamber. High speed camera visualizations of the combustion processes allow to compare the flame height, obtained by the code and by experimental tests, along the grain for given boundary conditions.
[Numerical methods for multi-fluid flows]. Final progress report
Pozrikidis, C.
1998-07-21
The central objective of this research has been to develop efficient numerical methods for computing multi-fluid flows with large interfacial deformations, and apply these methods to study the rheology of suspensions of deformable particles with viscous and non-Newtonian interfacial behavior. The mathematical formulation employs boundary-integral, immersed-boundary, and related numerical methods. Particles of interest include liquid drops with constant surface tension and capsules whose interfaces exhibit viscoelastic and incompressible characteristics. In one family of problems, the author has considered the shear-driven and pressure-driven flow of a suspension of two-dimensional liquid drops with ordered and random structure. In a second series of investigations, the author carried out dynamic simulations of two-dimensional, unbounded, doubly-periodic shear flows with random structure. Another family of problems addresses the deformation of three-dimensional capsules whose interfaces exhibit isotropic surface tension, viscous, elastic, or incompressible behavior, in simple shear flow. The numerical results extend previous asymptotic theories for small deformations and illuminate the mechanism of membrane rupture.
Numerical method for wave forces acting on partially perforated caisson
NASA Astrophysics Data System (ADS)
Jiang, Feng; Tang, Xiao-cheng; Jin, Zhao; Zhang, Li; Chen, Hong-zhou
2015-04-01
The perforated caisson is widely applied to practical engineering because of its great advantages in effectively wave energy consumption and cost reduction. The attentions of many scientists were paid to the fluid-structure interaction between wave and perforated caisson studies, but until now, most concerns have been put on theoretical analysis and experimental model set up. In this paper, interaction between the wave and the partial perforated caisson in a 2D numerical wave flume is investigated by means of the renewed SPH algorithm, and the mathematical equations are in the form of SPH numerical approximation based on Navier-Stokes equations. The validity of the SPH mathematical method is examined and the simulated results are compared with the results of theoretical models, meanwhile the complex hydrodynamic characteristics when the water particles flow in or out of a wave absorbing chamber are analyzed and the wave pressure distribution of the perforated caisson is also addressed here. The relationship between the ratio of total horizontal force acting on caisson under regular waves and its influence factors is examined. The data show that the numerical calculation of the ratio of total horizontal force meets the empirical regression equation very well. The simulations of SPH about the wave nonlinearity and breaking are briefly depicted in the paper, suggesting that the advantages and great potentiality of the SPH method is significant compared with traditional methods.
An alternative numerical method for the stationary pulsar magnetosphere
NASA Astrophysics Data System (ADS)
Takamori, Yohsuke; Okawa, Hirotada; Takamoto, Makoto; Suwa, Yudai
2014-02-01
Stationary pulsar magnetospheres in the force-free system are governed by the pulsar equation. In 1999, Contopoulos, Kazanas, and Fendt (hereafter CKF) numerically solved the pulsar equation and obtained a pulsar magnetosphere model called the CKF solution that has both closed and open magnetic field lines. The CKF solution is a successful solution, but it contains a poloidal current sheet that flows along the last open field line. This current sheet is artificially added to make the current system closed. In this paper, we suggest an alternative method to solve the pulsar equation and construct pulsar magnetosphere models without a current sheet. In our method, the pulsar equation is decomposed into Ampère's law and the force-free condition. We numerically solve these equations simultaneously with a fixed poloidal current. As a result, we obtain a pulsar magnetosphere model without a current sheet, which is similar to the CKF solution near the neutron star and has a jet-like structure at a distance along the pole. In addition, we discuss physical properties of the model and find that the force-free condition breaks down in a vicinity of the light cylinder due to dissipation that is included implicitly in the numerical method.
Paraxial WKB Method Applied to the Lower Hybrid Wave Propagation
Bertelli, N; Poli, E; Harvey, R; Wright, J C; Bonoli, P T; Phillips, C K; Simov, A P; Valeo, E
2012-07-12
The paraxial WKB (pWKB) approximation, also called beam tracing method, has been employed in order to study the propagation of lower hybrid (LH) waves in a tokamak plasma. Analogous to the well-know ray tracing method, this approach reduces Maxwell's equations to a set of ordinary differential equations, while, in addition, retains the effects of the finite beam cross-section, and, thus, the effects of diffraction. A new code, LHBEAM (Lower Hybrid BEAM tracing), is presented, which solves the pWKB equations in tokamak geometry for arbitrary launching conditions and for analytic and experimental plasma equilibria. In addition, LHBEAM includes linear electron Landau damping for the evaluation of the absorbed power density and the reconstruction of the wave electric field in both the physical and Fourier space. Illustrative LHBEAM calculations are presented along with a comparison with the ray tracing code GENRAY and the full wave solver TORIC-LH.
Hybrid thermal link-wise artificial compressibility method
NASA Astrophysics Data System (ADS)
Obrecht, Christian; Kuznik, Frédéric
2015-10-01
Thermal flow prediction is a subject of interest from a scientific and engineering points of view. Our motivation is to develop an accurate, easy to implement and highly scalable method for convective flows simulation. To this end, we present an extension to the link-wise artificial compressibility method (LW-ACM) for thermal simulation of weakly compressible flows. The novel hybrid formulation uses second-order finite difference operators of the energy equation based on the same stencils as the LW-ACM. For validation purposes, the differentially heated cubic cavity was simulated. The simulations remained stable for Rayleigh numbers up to Ra =108. The Nusselt numbers at isothermal walls and dynamics quantities are in good agreement with reference values from the literature. Our results show that the hybrid thermal LW-ACM is an effective and easy-to-use solution to solve convective flows.
The instanton method and its numerical implementation in fluid mechanics
NASA Astrophysics Data System (ADS)
Grafke, Tobias; Grauer, Rainer; Schäfer, Tobias
2015-08-01
A precise characterization of structures occurring in turbulent fluid flows at high Reynolds numbers is one of the last open problems of classical physics. In this review we discuss recent developments related to the application of instanton methods to turbulence. Instantons are saddle point configurations of the underlying path integrals. They are equivalent to minimizers of the related Freidlin-Wentzell action and known to be able to characterize rare events in such systems. While there is an impressive body of work concerning their analytical description, this review focuses on the question on how to compute these minimizers numerically. In a short introduction we present the relevant mathematical and physical background before we discuss the stochastic Burgers equation in detail. We present algorithms to compute instantons numerically by an efficient solution of the corresponding Euler-Lagrange equations. A second focus is the discussion of a recently developed numerical filtering technique that allows to extract instantons from direct numerical simulations. In the following we present modifications of the algorithms to make them efficient when applied to two- or three-dimensional (2D or 3D) fluid dynamical problems. We illustrate these ideas using the 2D Burgers equation and the 3D Navier-Stokes equations.
Hybrid finite element and Brownian dynamics method for charged particles
NASA Astrophysics Data System (ADS)
Huber, Gary A.; Miao, Yinglong; Zhou, Shenggao; Li, Bo; McCammon, J. Andrew
2016-04-01
Diffusion is often the rate-determining step in many biological processes. Currently, the two main computational methods for studying diffusion are stochastic methods, such as Brownian dynamics, and continuum methods, such as the finite element method. A previous study introduced a new hybrid diffusion method that couples the strengths of each of these two methods, but was limited by the lack of interactions among the particles; the force on each particle had to be from an external field. This study further develops the method to allow charged particles. The method is derived for a general multidimensional system and is presented using a basic test case for a one-dimensional linear system with one charged species and a radially symmetric system with three charged species.
Dielectric Boundary Forces in Numerical Poisson-Boltzmann Methods: Theory and Numerical Strategies.
Cai, Qin; Ye, Xiang; Wang, Jun; Luo, Ray
2011-10-01
Continuum modeling of electrostatic interactions based upon the numerical solutions of the Poisson-Boltzmann equation has been widely adopted in biomolecular applications. To extend their applications to molecular dynamics and energy minimization, robust and efficient methodologies to compute solvation forces must be developed. In this study, we have first reviewed the theory for the computation of dielectric boundary forces based on the definition of the Maxwell stress tensor. This is followed by a new formulation of the dielectric boundary force suitable for the finite-difference Poisson-Boltzmann methods. We have validated the new formulation with idealized analytical systems and realistic molecular systems. PMID:22125339
Dielectric boundary force in numerical Poisson-Boltzmann methods: Theory and numerical strategies
NASA Astrophysics Data System (ADS)
Cai, Qin; Ye, Xiang; Wang, Jun; Luo, Ray
2011-10-01
Continuum modeling of electrostatic interactions based upon the numerical solutions of the Poisson-Boltzmann equation has been widely adopted in biomolecular applications. To extend their applications to molecular dynamics and energy minimization, robust and efficient methodologies to compute solvation forces must be developed. In this study, we have first reviewed the theory for the computation of dielectric boundary force based on the definition of the Maxwell stress tensor. This is followed by a new formulation of the dielectric boundary force suitable for the finite-difference Poisson-Boltzmann methods. We have validated the new formulation with idealized analytical systems and realistic molecular systems.
Battery control system for hybrid vehicle and method for controlling a hybrid vehicle battery
Bockelmann, Thomas R.; Beaty, Kevin D.; Zou, Zhanijang; Kang, Xiaosong
2009-07-21
A battery control system for controlling a state of charge of a hybrid vehicle battery includes a detecting arrangement for determining a vehicle operating state or an intended vehicle operating state and a controller for setting a target state of charge level of the battery based on the vehicle operating state or the intended vehicle operating state. The controller is operable to set a target state of charge level at a first level during a mobile vehicle operating state and at a second level during a stationary vehicle operating state or in anticipation of the vehicle operating in the stationary vehicle operating state. The invention further includes a method for controlling a state of charge of a hybrid vehicle battery.
Taylor, G.; Dong, C.; Sun, S.
2010-03-18
A mathematical model for contaminant species passing through fractured porous media is presented. In the numerical model, we combine two locally conservative methods, i.e. mixed finite element (MFE) and the finite volume methods. Adaptive triangle mesh is used for effective treatment of the fractures. A hybrid MFE method is employed to provide an accurate approximation of velocities field for both the fractures and matrix which are crucial to the convection part of the transport equation. The finite volume method and the standard MFE method are used to approximate the convection and dispersion terms respectively. The model is used to investigate the interaction of adsorption with transport and to extract information on effective adsorption distribution coefficients. Numerical examples in different fractured media illustrate the robustness and efficiency of the proposed numerical model.
Numerical Method for Darcy Flow Derived Using Discrete Exterior Calculus
NASA Astrophysics Data System (ADS)
Hirani, A. N.; Nakshatrala, K. B.; Chaudhry, J. H.
2015-05-01
We derive a numerical method for Darcy flow, and also for Poisson's equation in mixed (first order) form, based on discrete exterior calculus (DEC). Exterior calculus is a generalization of vector calculus to smooth manifolds and DEC is one of its discretizations on simplicial complexes such as triangle and tetrahedral meshes. DEC is a coordinate invariant discretization, in that it does not depend on the embedding of the simplices or the whole mesh. We start by rewriting the governing equations of Darcy flow using the language of exterior calculus. This yields a formulation in terms of flux differential form and pressure. The numerical method is then derived by using the framework provided by DEC for discretizing differential forms and operators that act on forms. We also develop a discretization for a spatially dependent Hodge star that varies with the permeability of the medium. This also allows us to address discontinuous permeability. The matrix representation for our discrete non-homogeneous Hodge star is diagonal, with positive diagonal entries. The resulting linear system of equations for flux and pressure are saddle type, with a diagonal matrix as the top left block. The performance of the proposed numerical method is illustrated on many standard test problems. These include patch tests in two and three dimensions, comparison with analytically known solutions in two dimensions, layered medium with alternating permeability values, and a test with a change in permeability along the flow direction. We also show numerical evidence of convergence of the flux and the pressure. A convergence experiment is included for Darcy flow on a surface. A short introduction to the relevant parts of smooth and discrete exterior calculus is included in this article. We also include a discussion of the boundary condition in terms of exterior calculus.
Teaching Thermal Hydraulics & Numerical Methods: An Introductory Control Volume Primer
Lucas, D.S.
2004-10-03
This paper covers the basics of the implementation of the control volume method in the context of the Homogeneous Equilibrium Model (HEM)(T/H) code using the conservation equations of mass, momentum, and energy. This primer uses the advection equation as a template. The discussion will cover the basic equations of the control volume portion of the course in the primer, which includes the advection equation, numerical methods, along with the implementation of the various equations via FORTRAN into computer programs and the final result for a three equation HEM code and its validation.
Numerical methods for scattering from electrically large objects
NASA Astrophysics Data System (ADS)
Enguist, Bjorn; Murphy, W. D.; Rokhlin, Vladimir; Vassiliou, Marius S.
1991-05-01
A new and computationally very efficient integral equation numerical method for computing electromagnetic scattering and radar cross section (RCS) was developed. A theory of higher order impedance boundary conditions was derived to handle single and multiple dielectric coatings around conductors. The method was tested in two dimensions using a 14,000-line FORTRAN program and was found to be very promising for electrically large objects. Initial ideas for extensions to three dimensions were explored. Treatments of trailing edge and corner singularities were developed.
Integrated numerical methods for hypersonic aircraft cooling systems analysis
NASA Technical Reports Server (NTRS)
Petley, Dennis H.; Jones, Stuart C.; Dziedzic, William M.
1992-01-01
Numerical methods have been developed for the analysis of hypersonic aircraft cooling systems. A general purpose finite difference thermal analysis code is used to determine areas which must be cooled. Complex cooling networks of series and parallel flow can be analyzed using a finite difference computer program. Both internal fluid flow and heat transfer are analyzed, because increased heat flow causes a decrease in the flow of the coolant. The steady state solution is a successive point iterative method. The transient analysis uses implicit forward-backward differencing. Several examples of the use of the program in studies of hypersonic aircraft and rockets are provided.
A comparison of generalized hybrid Monte Carlo methods with and without momentum flip
Akhmatskaya, Elena; Bou-Rabee, Nawaf; Reich, Sebastian
2009-04-01
The generalized hybrid Monte Carlo (GHMC) method combines Metropolis corrected constant energy simulations with a partial random refreshment step in the particle momenta. The standard detailed balance condition requires that momenta are negated upon rejection of a molecular dynamics proposal step. The implication is a trajectory reversal upon rejection, which is undesirable when interpreting GHMC as thermostated molecular dynamics. We show that a modified detailed balance condition can be used to implement GHMC without momentum flips. The same modification can be applied to the generalized shadow hybrid Monte Carlo (GSHMC) method. Numerical results indicate that GHMC/GSHMC implementations with momentum flip display a favorable behavior in terms of sampling efficiency, i.e., the traditional GHMC/GSHMC implementations with momentum flip got the advantage of a higher acceptance rate and faster decorrelation of Monte Carlo samples. The difference is more pronounced for GHMC. We also numerically investigate the behavior of the GHMC method as a Langevin-type thermostat. We find that the GHMC method without momentum flip interferes less with the underlying stochastic molecular dynamics in terms of autocorrelation functions and it to be preferred over the GHMC method with momentum flip. The same finding applies to GSHMC.
Numerical methods for the Poisson-Fermi equation in electrolytes
NASA Astrophysics Data System (ADS)
Liu, Jinn-Liang
2013-08-01
The Poisson-Fermi equation proposed by Bazant, Storey, and Kornyshev [Phys. Rev. Lett. 106 (2011) 046102] for ionic liquids is applied to and numerically studied for electrolytes and biological ion channels in three-dimensional space. This is a fourth-order nonlinear PDE that deals with both steric and correlation effects of all ions and solvent molecules involved in a model system. The Fermi distribution follows from classical lattice models of configurational entropy of finite size ions and solvent molecules and hence prevents the long and outstanding problem of unphysical divergence predicted by the Gouy-Chapman model at large potentials due to the Boltzmann distribution of point charges. The equation reduces to Poisson-Boltzmann if the correlation length vanishes. A simplified matched interface and boundary method exhibiting optimal convergence is first developed for this equation by using a gramicidin A channel model that illustrates challenging issues associated with the geometric singularities of molecular surfaces of channel proteins in realistic 3D simulations. Various numerical methods then follow to tackle a range of numerical problems concerning the fourth-order term, nonlinearity, stability, efficiency, and effectiveness. The most significant feature of the Poisson-Fermi equation, namely, its inclusion of steric and correlation effects, is demonstrated by showing good agreement with Monte Carlo simulation data for a charged wall model and an L type calcium channel model.
Optimization methods and silicon solar cell numerical models
NASA Technical Reports Server (NTRS)
Girardini, K.; Jacobsen, S. E.
1986-01-01
An optimization algorithm for use with numerical silicon solar cell models was developed. By coupling an optimization algorithm with a solar cell model, it is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junction depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm was developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAP1D). SCAP1D uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the performance of a solar cell. A major obstacle is that the numerical methods used in SCAP1D require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the values associated with the maximum efficiency. This problem was alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution.
Numerical integration of population models satisfying conservation laws: NSFD methods.
Mickens, Ronald E
2007-10-01
Population models arising in ecology, epidemiology and mathematical biology may involve a conservation law, i.e. the total population is constant. In addition to these cases, other situations may occur for which the total population, asymptotically in time, approach a constant value. Since it is rarely the situation that the equations of motion can be analytically solved to obtain exact solutions, it follows that numerical techniques are needed to provide solutions. However, numerical procedures are only valid if they can reproduce fundamental properties of the differential equations modeling the phenomena of interest. We show that for population models, involving a dynamical conservation law the use of nonstandard finite difference (NSFD) methods allows the construction of discretization schemes such that they are dynamically consistent (DC) with the original differential equations. The paper will briefly discuss the NSFD methodology, the concept of DC, and illustrate their application to specific problems for population models. PMID:22876826
A Numerical Method for Solving Elasticity Equations with Interfaces
Li, Zhilin; Wang, Liqun; Wang, Wei
2012-01-01
Solving elasticity equations with interfaces is a challenging problem for most existing methods. Nonetheless, it has wide applications in engineering and science. An accurate and efficient method is desired. In this paper, an efficient non-traditional finite element method with non-body-fitting grids is proposed to solve elasticity equations with interfaces. The main idea is to choose the test function basis to be the standard finite element basis independent of the interface and to choose the solution basis to be piecewise linear satisfying the jump conditions across the interface. The resulting linear system of equations is shown to be positive definite under certain assumptions. Numerical experiments show that this method is second order accurate in the L∞ norm for piecewise smooth solutions. More than 1.5th order accuracy is observed for solution with singularity (second derivative blows up) on the sharp-edged interface corner. PMID:22707984
Time-dependent corona models - A numerical method
NASA Astrophysics Data System (ADS)
Korevaar, P.; van Leer, B.
1988-07-01
A time-dependent numerical method for calculating gas flows is described. The method is implicit and especially suitable for finding stationary flow solutions. Although the method is general in its application to ideal compressible fluids, this paper applies it to a stellar atmosphere, heated to coronal temperatures by dissipation of mechanical energy. The integration scheme is based on conservative upwind spatial differencing. The upwind switching is provided by Van Leer's method of differentiable flux-splitting. It is shown that the code can handle large differences in density: up to 14 orders of magnitude. Special attention is paid to the boundary conditions, which are made completely transparent to disturbances. Besides some test-results, converged solutions for various values of the initial mechanical flux are presented which are in good agreement with previous time-independent calculations.
Advanced numerical methods in mesh generation and mesh adaptation
Lipnikov, Konstantine; Danilov, A; Vassilevski, Y; Agonzal, A
2010-01-01
Numerical solution of partial differential equations requires appropriate meshes, efficient solvers and robust and reliable error estimates. Generation of high-quality meshes for complex engineering models is a non-trivial task. This task is made more difficult when the mesh has to be adapted to a problem solution. This article is focused on a synergistic approach to the mesh generation and mesh adaptation, where best properties of various mesh generation methods are combined to build efficiently simplicial meshes. First, the advancing front technique (AFT) is combined with the incremental Delaunay triangulation (DT) to build an initial mesh. Second, the metric-based mesh adaptation (MBA) method is employed to improve quality of the generated mesh and/or to adapt it to a problem solution. We demonstrate with numerical experiments that combination of all three methods is required for robust meshing of complex engineering models. The key to successful mesh generation is the high-quality of the triangles in the initial front. We use a black-box technique to improve surface meshes exported from an unattainable CAD system. The initial surface mesh is refined into a shape-regular triangulation which approximates the boundary with the same accuracy as the CAD mesh. The DT method adds robustness to the AFT. The resulting mesh is topologically correct but may contain a few slivers. The MBA uses seven local operations to modify the mesh topology. It improves significantly the mesh quality. The MBA method is also used to adapt the mesh to a problem solution to minimize computational resources required for solving the problem. The MBA has a solid theoretical background. In the first two experiments, we consider the convection-diffusion and elasticity problems. We demonstrate the optimal reduction rate of the discretization error on a sequence of adaptive strongly anisotropic meshes. The key element of the MBA method is construction of a tensor metric from hierarchical edge
Hybrid star structure with the Field Correlator Method
NASA Astrophysics Data System (ADS)
Burgio, G. F.; Zappalà, D.
2016-03-01
We explore the relevance of the color-flavor locking phase in the equation of state (EoS) built with the Field Correlator Method (FCM) for the description of the quark matter core of hybrid stars. For the hadronic phase, we use the microscopic Brueckner-Hartree-Fock (BHF) many-body theory, and its relativistic counterpart, i.e. the Dirac-Brueckner (DBHF). We find that the main features of the phase transition are directly related to the values of the quark-antiquark potential V1, the gluon condensate G2 and the color-flavor superconducting gap Δ. We confirm that the mapping between the FCM and the CSS (constant speed of sound) parameterization holds true even in the case of paired quark matter. The inclusion of hyperons in the hadronic phase and its effect on the mass-radius relation of hybrid stars is also investigated.
NASA Astrophysics Data System (ADS)
Zi, Bin; Zhou, Bin
2016-07-01
For the prediction of dynamic response field of the luffing system of an automobile crane (LSOAAC) with random and interval parameters, a hybrid uncertain model is introduced. In the hybrid uncertain model, the parameters with certain probability distribution are modeled as random variables, whereas, the parameters with lower and upper bounds are modeled as interval variables instead of given precise values. Based on the hybrid uncertain model, the hybrid uncertain dynamic response equilibrium equation, in which different random and interval parameters are simultaneously included in input and output terms, is constructed. Then a modified hybrid uncertain analysis method (MHUAM) is proposed. In the MHUAM, based on random interval perturbation method, the first-order Taylor series expansion and the first-order Neumann series, the dynamic response expression of the LSOAAC is developed. Moreover, the mathematical characteristics of extrema of bounds of dynamic response are determined by random interval moment method and monotonic analysis technique. Compared with the hybrid Monte Carlo method (HMCM) and interval perturbation method (IPM), numerical results show the feasibility and efficiency of the MHUAM for solving the hybrid LSOAAC problems. The effects of different uncertain models and parameters on the LSOAAC response field are also investigated deeply, and numerical results indicate that the impact made by the randomness in the thrust of the luffing cylinder F is larger than that made by the gravity of the weight in suspension Q . In addition, the impact made by the uncertainty in the displacement between the lower end of the lifting arm and the luffing cylinder a is larger than that made by the length of the lifting arm L .
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-09-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
2D-3D hybrid stabilized finite element method for tsunami runup simulations
NASA Astrophysics Data System (ADS)
Takase, S.; Moriguchi, S.; Terada, K.; Kato, J.; Kyoya, T.; Kashiyama, K.; Kotani, T.
2016-05-01
This paper presents a two-dimensional (2D)-three-dimensional (3D) hybrid stabilized finite element method that enables us to predict a propagation process of tsunami generated in a hypocentral region, which ranges from offshore propagation to runup to urban areas, with high accuracy and relatively low computational costs. To be more specific, the 2D shallow water equation is employed to simulate the propagation of offshore waves, while the 3D Navier-Stokes equation is employed for the runup in urban areas. The stabilized finite element method is utilized for numerical simulations for both of the 2D and 3D domains that are independently discretized with unstructured meshes. The multi-point constraint and transmission methods are applied to satisfy the continuity of flow velocities and pressures at the interface between the resulting 2D and 3D meshes, since neither their spatial dimensions nor node arrangements are consistent. Numerical examples are presented to demonstrate the performance of the proposed hybrid method to simulate tsunami behavior, including offshore propagation and runup to urban areas, with substantially lower computation costs in comparison with full 3D computations.
Meng, Zi-Ming E-mail: lizy@aphy.iphy.ac.cn; Hu, Yi-Hua; Ju, Gui-Fang; Zhong, Xiao-Lan; Ding, Wei; Li, Zhi-Yuan E-mail: lizy@aphy.iphy.ac.cn
2014-07-28
Optical Tamm states (OTSs) in analogy with its electronic counterpart confined at the surface of crystals are optical surface modes at the interfaces between uniform metallic films and distributed Bragg reflectors. In this paper, OTSs are numerically investigated in two-dimensional hybrid plasmonic-photonic crystal nanobeams (HPPCN), which are constructed by inserting a metallic nanoparticle into a photonic crystal nanobeam formed by periodically etching square air holes into dielectric waveguides. The evidences of OTSs can be verified by transmission spectra and the field distribution at resonant frequency. Similar to OTSs in one-dimensional multilayer structures OTSs in HPPCN can be excited by both TE and TM polarization. The physical origin of OTSs in HPPCN is due to the combined contribution of strong reflection imposed by the photonic band gap (PBG) of the photonic crystal (PC) nanobeam and strong backward scattering exerted by the nanoparticle. For TE, incidence OTSs can be obtained at the frequency near the center of the photonic band gap. The transmissivity and the resonant frequency can be finely tuned by the dimension of nanoparticles. While for TM incidence OTSs are observed for relatively larger metallic nanoparticles compared with TE polarization. The differences between TE and TM polarization can be explained by two reasons. For one reason stronger backward scattering of nanoparticles for TE polarization can be achieved by the excitation of localized surface plasmon polariton of nanoparticles. This assumption has been proved by examining the scattering, absorption, and extinction cross section of the metallic nanoparticle. The other can be attributed to the deep and wide PBG available for TE polarization with less number of air holes compared with TM polarization. Our results show great promise in extending the application scope of OTSs from one-dimensional structures to practical integrated photonic devices and circuits.
A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasma
NASA Astrophysics Data System (ADS)
Ku, S.; Hager, R.; Chang, C. S.; Kwon, J. M.; Parker, S. E.
2016-06-01
In order to enable kinetic simulation of non-thermal edge plasmas at a reduced computational cost, a new hybrid-Lagrangian δf scheme has been developed that utilizes the phase space grid in addition to the usual marker particles, taking advantage of the computational strengths from both sides. The new scheme splits the particle distribution function of a kinetic equation into two parts. Marker particles contain the fast space-time varying, δf, part of the distribution function and the coarse-grained phase-space grid contains the slow space-time varying part. The coarse-grained phase-space grid reduces the memory-requirement and the computing cost, while the marker particles provide scalable computing ability for the fine-grained physics. Weights of the marker particles are determined by a direct weight evolution equation instead of the differential form weight evolution equations that the conventional delta-f schemes use. The particle weight can be slowly transferred to the phase space grid, thereby reducing the growth of the particle weights. The non-Lagrangian part of the kinetic equation - e.g., collision operation, ionization, charge exchange, heat-source, radiative cooling, and others - can be operated directly on the phase space grid. Deviation of the particle distribution function on the velocity grid from a Maxwellian distribution function - driven by ionization, charge exchange and wall loss - is allowed to be arbitrarily large. The numerical scheme is implemented in the gyrokinetic particle code XGC1, which specializes in simulating the tokamak edge plasma that crosses the magnetic separatrix and is in contact with the material wall.
Seismic attenuation in fractured porous media: insights from a hybrid numerical and analytical model
NASA Astrophysics Data System (ADS)
Ekanem, A. M.; Li, X. Y.; Chapman, M.; Main, I. G.
2015-04-01
Seismic attenuation in fluid-saturated porous rocks can occur by geometric spreading, wave scattering or the internal dissipation of energy, most likely due to the squirt-flow mechanism. In principle, the pattern of seismic attenuation recorded on an array of sensors contains information about the medium, in terms of material heterogeneity and anisotropy, as well as material properties such as porosity, crack density, and pore-fluid composition and mobility. In practice, this inverse problem is challenging. Here we provide some insights into the effects of internal dissipation by analysing synthetic data produced by a hybrid numerical and analytical model for seismic wave propagation in a fractured medium embedded within a layered geological structure. The model is made up of one anisotropic and three isotropic horizontal layers. The anisotropic layer consists of a porous, fluid-saturated material containing vertically aligned inclusions representing a set of fractures. This combination allows squirt-flow to occur between the pores in the matrix and the model fractures. Our results show that the fluid mobility and the associated relaxation time of the fluid-pressure gradient control the frequency range over which attenuation occurs. This induced attenuation increases with incidence angle and azimuth away from the fracture strike-direction. Azimuthal variations in the induced attenuation are elliptical allowing the fracture orientations to be obtained from the axes of the ellipse. These observations hold out the potential of using seismic attenuation as an additional diagnostic in the characterisation of rock formations for a variety of applications including hydrocarbon exploration and production, subsurface storage of CO2, and geothermal energy extraction.
Numeric Design and Performance Analysis of Solid Oxide Fuel Cell -- Gas Turbine Hybrids on Aircraft
NASA Astrophysics Data System (ADS)
Hovakimyan, Gevorg
The aircraft industry benefits greatly from small improvements in aircraft component design. One possible area of improvement is in the Auxiliary Power Unit (APU). Modern aircraft APUs are gas turbines located in the tail section of the aircraft that generate additional power when needed. Unfortunately the efficiency of modern aircraft APUs is low. Solid Oxide Fuel Cell/Gas Turbine (SOFC/GT) hybrids are one possible alternative for replacing modern gas turbine APUs. This thesis investigates the feasibility of replacing conventional gas turbine APUs with SOFC/GT APUs on aircraft. An SOFC/GT design algorithm was created in order to determine the specifications of an SOFC/GT APU. The design algorithm is comprised of several integrated modules which together model the characteristics of each component of the SOFC/GT system. Given certain overall inputs, through numerical analysis, the algorithm produces an SOFC/GT APU, optimized for specific power and efficiency, capable of performing to the required specifications. The SOFC/GT design is then input into a previously developed quasi-dynamic SOFC/GT model to determine its load following capabilities over an aircraft flight cycle. Finally an aircraft range study is conducted to determine the feasibility of the SOFC/GT APU as a replacement for the conventional gas turbine APU. The design results show that SOFC/GT APUs have lower specific power than GT systems, but have much higher efficiencies. Moreover, the dynamic simulation results show that SOFC/GT APUs are capable of following modern flight loads. Finally, the range study determined that SOFC/GT APUs are more attractive over conventional APUs for longer range aircraft.
Comparison of four stable numerical methods for Abel's integral equation
NASA Technical Reports Server (NTRS)
Murio, Diego A.; Mejia, Carlos E.
1991-01-01
The 3-D image reconstruction from cone-beam projections in computerized tomography leads naturally, in the case of radial symmetry, to the study of Abel-type integral equations. If the experimental information is obtained from measured data, on a discrete set of points, special methods are needed in order to restore continuity with respect to the data. A new combined Regularized-Adjoint-Conjugate Gradient algorithm, together with two different implementations of the Mollification Method (one based on a data filtering technique and the other on the mollification of the kernal function) and a regularization by truncation method (initially proposed for 2-D ray sample schemes and more recently extended to 3-D cone-beam image reconstruction) are extensively tested and compared for accuracy and numerical stability as functions of the level of noise in the data.
Numerical Analysis of a Finite Element/Volume Penalty Method
NASA Astrophysics Data System (ADS)
Maury, Bertrand
The penalty method makes it possible to incorporate a large class of constraints in general purpose Finite Element solvers like freeFEM++. We present here some contributions to the numerical analysis of this method. We propose an abstract framework for this approach, together with some general error estimates based on the discretization parameter ɛ and the space discretization parameter h. As this work is motivated by the possibility to handle constraints like rigid motion for fluid-particle flows, we shall pay a special attention to a model problem of this kind, where the constraint is prescribed over a subdomain. We show how the abstract estimate can be applied to this situation, in the case where a non-body-fitted mesh is used. In addition, we describe how this method provides an approximation of the Lagrange multiplier associated to the constraint.
Dord, Jean-Francois; Farhat, Charbel
2010-08-01
This paper considers the problem of imaging a complex object submerged in shallow waters using a sparse surface sensor array and a hybrid signal processing method. This method is constructed by refining the Kirchhoff migration technique to incorporate a zoning of the sensors and an analysis of multiple reflections, and combining it with the direction of arrival estimation method. Its performance is assessed and analyzed with the shape identification of a mockup submarine by numerical simulation. The obtained numerical results highlight the potential of this approach for identifying underwater intruders. PMID:20707441
Numerical methods for high-dimensional probability density function equations
NASA Astrophysics Data System (ADS)
Cho, H.; Venturi, D.; Karniadakis, G. E.
2016-01-01
In this paper we address the problem of computing the numerical solution to kinetic partial differential equations involving many phase variables. These types of equations arise naturally in many different areas of mathematical physics, e.g., in particle systems (Liouville and Boltzmann equations), stochastic dynamical systems (Fokker-Planck and Dostupov-Pugachev equations), random wave theory (Malakhov-Saichev equations) and coarse-grained stochastic systems (Mori-Zwanzig equations). We propose three different classes of new algorithms addressing high-dimensionality: The first one is based on separated series expansions resulting in a sequence of low-dimensional problems that can be solved recursively and in parallel by using alternating direction methods. The second class of algorithms relies on truncation of interaction in low-orders that resembles the Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) framework of kinetic gas theory and it yields a hierarchy of coupled probability density function equations. The third class of algorithms is based on high-dimensional model representations, e.g., the ANOVA method and probabilistic collocation methods. A common feature of all these approaches is that they are reducible to the problem of computing the solution to high-dimensional equations via a sequence of low-dimensional problems. The effectiveness of the new algorithms is demonstrated in numerical examples involving nonlinear stochastic dynamical systems and partial differential equations, with up to 120 variables.
Higher-order brick-tetrahedron hybrid method for Maxwell's equations in time domain
NASA Astrophysics Data System (ADS)
Winges, Johan; Rylander, Thomas
2016-09-01
We present a higher-order brick-tetrahedron hybrid method for Maxwell's equations in time domain. Brick-shaped elements are used for large homogeneous parts of the computational domain, where we exploit mass-lumping and explicit time-stepping. In regions with complex geometry, we use an unstructured mesh of tetrahedrons that share an interface with the brick-shaped elements and, at the interface, tangential continuity of the electric field is imposed in the weak sense by means of Nitsche's method. Implicit time-stepping is used for the tetrahedrons together with the interface. For cavity resonators, the hybrid method reproduces the lowest non-zero eigenvalues with correct multiplicity and, for geometries without field singularities from sharp corners or edges, the numerical eigenvalues converge towards the analytical result with an error that is approximately proportional to h2p, where h is the cell size and p is the polynomial order of the elements. For a rectangular waveguide, a layer of tetrahedrons embedded in a grid of brick-shaped elements yields a low reflection coefficient that scales approximately as h2p. Finally, we demonstrate hybrid time-stepping for a lossless closed cavity resonator, where the time-domain response is computed for 300,000 time steps without any signs of instabilities.
Improved numerical method for subchannel cross-flow calculations
Kaya, S.; Anghaie, S.
1986-01-01
COBRA-OSU is a fast running computer code for coupled kinetic and thermal-hydraulic analysis of nuclear reactor core subchannels, currently under development at Oregon State University. This code is a modified version of COBRA-IV with two major improved features. First, COBRA-OSU uses the Gaussian elimination method instead of Gauss-Seidel iteration for subchannel cross-flow calculation. Second, COBRA-OSU has an additional model for regionwise point reactor kinetics which includes all major feedback reactivity effects on calculation of the axial power profile during the course of a transient. This paper summarizes the improved numerical features of the COBRA-OSU code.
Hybrid finite element-finite difference method for thermal analysis of blood vessels.
Blanchard, C H; Gutierrez, G; White, J A; Roemer, R B
2000-01-01
A hybrid finite-difference/finite-element technique for the thermal analysis of blood vessels embedded in perfused tissue has been developed and evaluated. This method provides efficient and accurate solutions to the conjugated heat transfer problem of convection by blood coupled to conduction in the tissue. The technique uses a previously developed 3D automatic meshing method for creating a finite element mesh in the tissue surrounding the vessels, coupled iteratively with a 1-D marching finite difference method for the interior of the vessels. This hybrid technique retains the flexibility and ease of automated finite-element meshing techniques for modelling the complex geometry of blood vessels and irregularly shaped tissues, and speeds the solution time by using a simple finite-difference method to calculate the bulk mean temperatures within all blood vessels. The use of the 1D finite-difference technique in the blood vessels also eliminates the large computer memory requirements needed to accurately solve large vessel network problems when fine FE meshes are used in the interior of vessels. The accuracy of the hybrid technique has been verified against previously verified numerical solutions. In summary, the hybrid technique combines the accuracy and flexibility found in automated finite-element techniques, with the speed and reduction of computational memory requirements associated with the 1D finite-difference technique, something which has not been done before. This method, thus, has the potential to provide accurate, flexible and relatively fast solutions for the thermal analysis of coupled perfusion/blood vessel problems, and large vessel network problems. PMID:10949130
A hybrid method for modelling two dimensional non-breaking and breaking waves
NASA Astrophysics Data System (ADS)
Sriram, V.; Ma, Q. W.; Schlurmann, T.
2014-09-01
This is the first paper to present a hybrid method coupling an Improved Meshless Local Petrov Galerkin method with Rankine source solution (IMLPG_R) based on the Navier-Stokes (NS) equations, with a finite element method (FEM) based on the fully nonlinear potential flow theory (FNPT) in order to efficiently simulate the violent waves and their interaction with marine structures. The two models are strongly coupled in space and time domains using a moving overlapping zone, wherein the information from both the solvers is exchanged. In the time domain, the Runge-Kutta 2nd order method is nested with a predictor-corrector scheme. In the space domain, numerical techniques including ‘Feeding Particles’ and two-layer particle interpolation with relaxation coefficients are introduced to achieve the robust coupling of the two models. The properties and behaviours of the new hybrid model are tested by modelling a regular wave, solitary wave and Cnoidal wave including breaking and overtopping. It is validated by comparing the results of the method with analytical solutions, results from other methods and experimental data. The paper demonstrates that the method can produce satisfactory results but uses much less computational time compared with a method based on the full NS model.
Microelectrode arrays fabricated using a novel hybrid microfabrication method.
Merlo, Mark W; Snyder, Russell L; Middlebrooks, John C; Bachman, Mark
2012-02-01
We present novel hybrid microfabrication methods for microelectrode arrays that combine microwire assembly, microelectromechanical systems (MEMS) manufacturing techniques and precision tool-based micromachining. This combination enables hybrid microfabrication to produce complex geometries and structures, increase material selection, and improve integration. A 32-channel shank microelectrode array was fabricated to highlight the hybrid microfabrication techniques. The electrode shank was 130 μm at its narrowest, had a 127 μm thickness and had iridium oxide electrode sites that were 25 μm in diameter with 150 μm spacing. Techniques used to fabricate this electrode include microassembly of insulated gold wires into a micromold, micromolding the microelectrode shank, post molding machining, sacrificial release of the microelectrode and electrodeposition of iridium oxide onto the microelectrode sites. Electrode site position accuracy was shown to have a standard deviation of less than 4 μm. Acute in vivo recordings with the 32-channel shank microelectrode array demonstrated comparable performance to that obtained with commercial microelectrode arrays. This new approach to microelectrode array fabrication will enable new microelectrodes, such as multi-sided arrays, drug eluding electrodes and biodegradable shanks. PMID:21979567
Microelectrode arrays fabricated using a novel hybrid microfabrication method
Merlo, Mark W.; Snyder, Russell L.; Middlebrooks, John C.; Bachman, Mark
2011-01-01
We present novel hybrid microfabrication methods for microelectrode arrays that combine microwire assembly, microelectromechanical systems (MEMS) manufacturing techniques and precision tool-based micromachining. This combination enables hybrid microfabrication to produce complex geometries and structures, increase material selection, and improve integration. A 32-channel shank microelectrode array was fabricated to highlight the hybrid microfabrication techniques. The electrode shank was 130 μm at its narrowest, had a 127 μm thickness and had iridium oxide electrode sites that were 25 μm in diameter with 150 μm spacing. Techniques used to fabricate this electrode include microassembly of insulated gold wires into a micromold, micromolding the microelectrode shank, post molding machining, sacrificial release of the microelectrode and electrodeposition of iridium oxide onto the microelectrode sites. Electrode site position accuracy was shown to have a standard deviation of less than 4 μm. Acute in vivo recordings with the 32-channel shank microelectrode array demonstrated comparable performance to that obtained with commercial microelectrode arrays . This new approach to microelectrode array fabrication will enable new microelectrodes, such as multi-sided arrays, drug eluding electrodes and biodegradable shanks. PMID:21979567
Hybrid method of deterministic and probabilistic approaches for multigroup neutron transport problem
Lee, D.
2012-07-01
A hybrid method of deterministic and probabilistic methods is proposed to solve Boltzmann transport equation. The new method uses a deterministic method, Method of Characteristics (MOC), for the fast and thermal neutron energy ranges and a probabilistic method, Monte Carlo (MC), for the intermediate resonance energy range. The hybrid method, in case of continuous energy problem, will be able to take advantage of fast MOC calculation and accurate resonance self shielding treatment of MC method. As a proof of principle, this paper presents the hybrid methodology applied to a multigroup form of Boltzmann transport equation and confirms that the hybrid method can produce consistent results with MC and MOC methods. (authors)
Optimization methods and silicon solar cell numerical models
NASA Technical Reports Server (NTRS)
Girardini, K.
1986-01-01
The goal of this project is the development of an optimization algorithm for use with a solar cell model. It is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junctions depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm has been developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAPID). SCAPID uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the operation of a solar cell. A major obstacle is that the numerical methods used in SCAPID require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the value associated with the maximum efficiency. This problem has been alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution. Adapting SCAPID so that it could be called iteratively by the optimization code provided another means of reducing the cpu time required to complete an optimization. Instead of calculating the entire I-V curve, as is usually done in SCAPID, only the efficiency is calculated (maximum power voltage and current) and the solution from previous calculations is used to initiate the next solution.
Teaching Thermal Hydraulics & Numerical Methods: An Introductory Control Volume Primer
D. S. Lucas
2004-10-01
A graduate level course for Thermal Hydraulics (T/H) was taught through Idaho State University in the spring of 2004. A numerical approach was taken for the content of this course since the students were employed at the Idaho National Laboratory and had been users of T/H codes. The majority of the students had expressed an interest in learning about the Courant Limit, mass error, semi-implicit and implicit numerical integration schemes in the context of a computer code. Since no introductory text was found the author developed notes taught from his own research and courses taught for Westinghouse on the subject. The course started with a primer on control volume methods and the construction of a Homogeneous Equilibrium Model (HEM) (T/H) code. The primer was valuable for giving the students the basics behind such codes and their evolution to more complex codes for Thermal Hydraulics and Computational Fluid Dynamics (CFD). The course covered additional material including the Finite Element Method and non-equilibrium (T/H). The control volume primer and the construction of a three-equation (mass, momentum and energy) HEM code are the subject of this paper . The Fortran version of the code covered in this paper is elementary compared to its descendants. The steam tables used are less accurate than the available commercial version written in C Coupled to a Graphical User Interface (GUI). The Fortran version and input files can be downloaded at www.microfusionlab.com.
Numerical modeling of spray combustion with an advanced VOF method
NASA Technical Reports Server (NTRS)
Chen, Yen-Sen; Shang, Huan-Min; Shih, Ming-Hsin; Liaw, Paul
1995-01-01
This paper summarizes the technical development and validation of a multiphase computational fluid dynamics (CFD) numerical method using the volume-of-fluid (VOF) model and a Lagrangian tracking model which can be employed to analyze general multiphase flow problems with free surface mechanism. The gas-liquid interface mass, momentum and energy conservation relationships are modeled by continuum surface mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed flow regimes. The objectives of the present study are to develop and verify the fractional volume-of-fluid cell partitioning approach into a predictor-corrector algorithm and to demonstrate the effectiveness of the present approach by simulating benchmark problems including laminar impinging jets, shear coaxial jet atomization and shear coaxial spray combustion flows.