首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
A new probe design is presented for obtaining homonuclear, heteronuclear, and inverse detected NMR spectra from more than one sample in the same total data acquisition time as for a single sample, thus increasing data acquisition efficiency. Specifically, a two-coil system, with each solenoidal coil impedance matched to 50 Omega at both proton and nitrogen frequencies, has been designed for operation at 11.7 T with an observe volume of 15 microL for each coil. Isolation between the two frequencies for each individual coil, and at each frequency between coils, was greater than 30 dB. Two-dimensional COSY and HMQC spectra were obtained with negligible NMR cross-talk between the two coils.  相似文献   

2.
We demonstrate the feasibility of using a non-conforming, piecewise harmonic finite element method on an unstructured grid in solving a magnetospheric physics problem. We use this approach to construct a global discrete model of the magnetic field of the magnetosphere that includes the effects of shielding currents at the outer boundary (the magnetopause). As in the approach of F. R. Toffolettoet al.(1994,Geophys. Res. Lett.21, 7) the internal magnetospheric field model is that of R. V. Hilmer and G.-H. Voigt (1995,J. Geophys. Res.) while the magnetopause shape is based on an empirically determined approximation (1997, J. Shueet al.,J. Geophys. Res.102, 9497). The results is a magnetic field model whose field lines are completely confined within the magnetosphere. The presented numerical results indicate that the discrete non-conforming finite element model is well-suited for magnetospheric field modeling.  相似文献   

3.
An r-adaptive finite-element method based on moving-mesh partial differential equations (PDEs) and an error indicator is presented. The error indicator is obtained by applying a technique developed by Bank and Weiser to elliptic equations which result in this case from temporal discretization of the underlying physical PDEs on moving meshes. The construction of the monitor function based on the error indicator is discussed. Numerical results obtained with the current method and the commonly used method based on solution gradients are presented and analyzed for several examples.  相似文献   

4.
We have implemented a scheme, SPECMON, for monitoring various parameters of a spectrometer, such as nitrogen pressure and sample temperature, and taking corrective action. The scheme is based on considerations of protection management which are of general application. Evaluation of the spectrometer state is incorporated in macros of the application software (VNMR) and is therefore very flexible. In contrast, corrective action is limited to the single one which is deemed fully safe: complete shutdown of the spectrometer and logging. Shutdown is implemented by a minor hardware modification of the spectrometer: the introduction of a second input to a relay already present for protection of the spectrometer power supply. Monitoring is handled by the host computer, and the shutdown command is transmitted via control lines of its series port, independent of the standard connection between the host computer and the NMR system console. The monitoring system (software and hardware) is unobtrusive in normal conditions, and it can be tested without affecting the operation of the spectrometer.  相似文献   

5.
A new finite volume method is presented for discretizing general linear or nonlinear elliptic second-order partial-differential equations with mixed boundary conditions. The advantage of this method is that arbitrary distorted meshes can be used without the numerical results being altered. The resulting algorithm has more unknowns than standard methods like finite difference or finite element methods. However, the matrices that need to be inverted are positive definite, so the most powerful linear solvers can be applied. The method has been tested on a few elliptic and parabolic equations, either linear, as in the case of the standard heat diffusion equation, or nonlinear, as in the case of the radiation diffusion equation and the resistive diffusion equation with Hall term.  相似文献   

6.
A new method is presented for the prediction of unsteady axisymmetric inviscid flows. By combining a triangulated vortex approach with a novel evaluation technique for the Biot–Savart integrals, a Lagrangian vortex method is developed which eliminates the singularities usually present in axisymmetric methods, without recourse to normalizations or other approximations. Furthermore, the computational effort scales as the number of control points N and, in the large N limit, depends only on the order of quadrature employed. The accuracy and computational effort are assessed by comparison with the velocity field of a Gaussian core vortex ring and the use of the technique is illustrated by computation of the motion of Norbury rings and of vortex ring pairing.  相似文献   

7.
We propose a new model and a solution method for two-phase compressible flows. The model involves six equations obtained from conservation principles applied to each phase, completed by a seventh equation for the evolution of the volume fraction. This equation is necessary to close the overall system. The model is valid for fluid mixtures, as well as for pure fluids. The system of partial differential equations is hyperbolic. Hyperbolicity is obtained because each phase is considered to be compressible. Two difficulties arise for the solution: one of the equations is written in non-conservative form; non-conservative terms exist in the momentum and energy equations. We propose robust and accurate discretisation of these terms. The method solves the same system at each mesh point with the same algorithm. It allows the simulation of interface problems between pure fluids as well as multiphase mixtures. Several test cases where fluids have compressible behavior are shown as well as some other test problems where one of the phases is incompressible. The method provides reliable results, is able to compute strong shock waves, and deals with complex equations of state.  相似文献   

8.
An extension of the alias sampling technique for distribution functions depending on a number of parameters was developed. It takes advantage of modern computer architectures with large amounts of cheap memory, by using discrete representations of probability distribution functions. The sampling is done by fast interpolation techniques involving only elementary logical and arithmetical operations, allowing one to keep a higher degree of accuracy as the grids spacing is controlled by the user. By this method it is possible to obtain the value of interest by direct interpolation between the sampled values obtained with the same set of random numbers for the grid values of the parameters adjacent to the values of interest. Sampling tests carried for the case of Molière electron multi-scatter angle distribution show that this method can be successfully used in Monte Carlo codes for sampling complex probability distributions.  相似文献   

9.
We present a new version of the fast multipole method (FMM) for screened Coulomb interactions in three dimensions. Existing schemes can compute such interactions in O(N) time, where N denotes the number of particles. The constant implicit in the O(N) notation, however, is dominated by the expense of translating far-field spherical harmonic expansions to local ones. For each box in the FMM data structure, this requires 189p4 operations per box, where p is the order of the expansions used. The new formulation relies on an expansion in evanescent plane waves, with which the amount of work can be reduced to 40p2+6p3 operations per box.  相似文献   

10.
11.
An algorithm is presented for the solution of the time dependent reaction-diffusion systems which arise in non-equilibrium radiation diffusion applications. This system of nonlinear equations is solved by coupling three numerical methods, Jacobian-free Newton–Krylov, operator splitting, and multigrid linear solvers. An inexact Newton's method is used to solve the system of nonlinear equations. Since building the Jacobian matrix for problems of interest can be challenging, we employ a Jacobian–free implementation of Newton's method, where the action of the Jacobian matrix on a vector is approximated by a first order Taylor series expansion. Preconditioned generalized minimal residual (PGMRES) is the Krylov method used to solve the linear systems that come from the iterations of Newton's method. The preconditioner in this solution method is constructed using a physics-based divide and conquer approach, often referred to as operator splitting. This solution procedure inverts the scalar elliptic systems that make up the preconditioner using simple multigrid methods. The preconditioner also addresses the strong coupling between equations with local 2×2 block solves. The intra-cell coupling is applied after the inter-cell coupling has already been addressed by the elliptic solves. Results are presented using this solution procedure that demonstrate its efficiency while incurring minimal memory requirements.  相似文献   

12.
Three-dimensional multiphase flow and flow with phase change are simulated using a simplified method of tracking and reconstructing the phase interface. The new level contour reconstruction technique presented here enables front tracking methods to naturally, automatically, and robustly model the merging and breakup of interfaces in three-dimensional flows. The method is designed so that the phase surface is treated as a collection of physically linked but not logically connected surface elements. Eliminating the need to bookkeep logical connections between neighboring surface elements greatly simplifies the Lagrangian tracking of interfaces, particularly for 3D flows exhibiting topology change. The motivation for this new method is the modeling of complex three-dimensional boiling flows where repeated merging and breakup are inherent features of the interface dynamics. Results of 3D film boiling simulations with multiple interacting bubbles are presented. The capabilities of the new interface reconstruction method are also tested in a variety of two-phase flows without phase change. Three-dimensional simulations of bubble merging and droplet collision, coalescence, and breakup demonstrate the new method's ability to easily handle topology change by film rupture or filamentary breakup. Validation tests are conducted for drop oscillation and bubble rise. The susceptibility of the numerical method to parasitic currents is also thoroughly assessed.  相似文献   

13.
In this paper we outline a new particle-mesh method for rapidly rotating shallow water flows based on a set of regularized equations of motion. The time-stepping method uses an operator splitting of the equations into an Eulerian gravity wave part and a Lagrangian advection part. An essential ingredient is the advection of absolute vorticity by means of translated radial basis functions. We show that this implies exact conservation of enstrophy. The method is tested on two model problems based on the qualitative features of the solutions obtained (i.e., dispersion or smoothness of potential vorticity contours) as well as on the increase in mean divergence level.  相似文献   

14.
In this paper we introduce a high-order discontinuous Galerkin method for two-dimensional incompressible flow in the vorticity stream-function formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method. The stream function is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability. The method is efficient for inviscid or high Reynolds number flows. Optimal error estimates are proved and verified by numerical experiments.  相似文献   

15.
This paper presents a numerical method directed towards the simulation of flows with mass transfer due to changes of phase. We use a volume of fluid (VOF) based interface tracking method in conjunction with a mass transfer model and a model for surface tension. The bulk fluids are viscous, conducting, and incompressible. A one-dimensional test problem is developed with the feature that a thin thermal layer propagates with the moving phase interface. This test problem isolates the ability of a method to accurately calculate the thermal layers responsible for driving the mass transfer in boiling flows. The numerical method is tested on this problem and then is used in simulations of horizontal film boiling.  相似文献   

16.
We have derived approximate analytic solutions to the master equation describing the evolution of the spin I=3/2 density operator in the presence of a radio-frequency (RF) field and both static and fluctuating quadrupolar interactions. Spectra resulting from Fourier transformation of the evolutions of the on-resonance spin-locked magnetization into the various coherences display two satellite pairs and, in some cases, a central line. The central line is generally trimodal, consisting of a narrow component related to a slowly relaxing mode and two broad components pertaining to two faster relaxing modes. The rates of the fast modes are sensitive to slow molecular motion. Neither the amplitude nor the width of the narrow component is affected by the magnitude of the static coupling, whereas the corresponding features of the broad components depend in a rather complicated manner on the spin-lock field strength and static quadrupolar interaction. Under certain experimental conditions, the dependencies of the amplitudes on the dynamics are seen to vanish and the relaxation rates reduce to relatively simple expressions. One of the promising emerging features is the fact that the evolutions into the selectively detected quadrupolar spin polarization order and the rank-two double-quantum coherence do not exhibit a slowly relaxing mode and are particularly sensitive to slow molecular motion. Furthermore, these coherences can only be excited in the presence of a static coupling and this makes it possible to discern nuclei in anisotropic from those in isotropic environment. The feasibility of the spin-lock pulse sequences with limited RF power and a nonvanishing average electric field gradient has been demonstrated through experiments on sodium in a dense lyotropic DNA liquid crystal.  相似文献   

17.
New NMR broadband inversion pulses that compensate both for resonance offset and radiofrequency (RF) inhomogeneity are described. The approach described is a straightforward computer optimization of an initial digitized waveform generated from either a constant-amplitude frequency sweep or from an existing composite inversion pulse. Problems with convergence to local minima are alleviated by the way the optimization is carried out. For a given duration and maximum allowable RF field strength B1 (but not necessarily given RMS power deposition), the resultant broadband inversion pulse (BIP) shows superior inversion compared to inversion pulses obtained from previous methods, including adiabatic inversion pulses. Any existing BIP can be systematically elaborated to build up longer inversion pulses that perform over larger and larger bandwidths. The resulting pulse need not be adiabatic throughout its duration or across the entire operational bandwidth.  相似文献   

18.
The particle-based lattice solid model developed to study the physics of rocks and the nonlinear dynamics of earthquakes is refined by incorporating intrinsic friction between particles. The model provides a means for studying the causes of seismic wave attenuation, as well as frictional heat generation, fault zone evolution, and localisation phenomena. A modified velocity–Verlat scheme that allows friction to be precisely modelled is developed. This is a difficult computational problem given that a discontinuity must be accurately simulated by the numerical approach (i.e., the transition from static to dynamical frictional behaviour). This is achieved using a half time step integration scheme. At each half time step, a nonlinear system is solved to compute the static frictional forces and states of touching particle-pairs. Improved efficiency is achieved by adaptively adjusting the time step increment, depending on the particle velocities in the system. The total energy is calculated and verified to remain constant to a high precision during simulations. Numerical experiments show that the model can be applied to the study of earthquake dynamics, the stick–slip instability, heat generation, and fault zone evolution. Such experiments may lead to a conclusive resolution of the heat flow paradox and improved understanding of earthquake precursory phenomena and dynamics.  相似文献   

19.
In many realistic fluid-dynamical simulations the specification of the boundary conditions, the error sources, and the number of time steps to reach a steady state are important practical considerations. In this paper we study these issues in the case of the lattice-BGK model. The objective is to present a comprehensive overview of some pitfalls and shortcomings of the lattice-BGK method and to introduce some new ideas useful in practical simulations. We begin with an evaluation of the widely used bounce-back boundary condition in staircase geometries by simulating flow in an inclined tube. It is shown that the bounce-back scheme is first-order accurate in space when the location of the non-slip wall is assumed to be at the boundary nodes. Moreover, for a specific inclination angle of 45 degrees, the scheme is found to be second-order accurate when the location of the non-slip velocity is fitted halfway between the last fluid nodes and the first solid nodes. The error as a function of the relaxation parameter is in that case qualitatively similar to that of flat walls. Next, a comparison of simulations of fluid flow by means of pressure boundaries and by means of body force is presented. A good agreement between these two boundary conditions has been found in the creeping-flow regime. For higher Reynolds numbers differences have been found that are probably caused by problems associated with the pressure boundaries. Furthermore, two widely used 3D models, namelyD3Q15andD3Q19, are analysed. It is shown that theD3Q15model may induce artificial checkerboard invariants due to the connectivity of the lattice. Finally, a new iterative method, which significantly reduces the saturation time, is presented and validated on different benchmark problems.  相似文献   

20.
Recovery of the magnetic dipolar interaction between nuclei bearing the same gyromagnetic ratio in rotating solids can be promoted by synchronous rf irradiation. Determination of the dipolar interaction strength can serve as a tool for structural elucidation in polycrystalline powders. Spinning frequency dependent narrow-band (nb) RFDR and SEDRA experiments are utilized as simple techniques for the determination of dipolar interactions between the nuclei in coupled homonuclear spin pairs. The magnetization exchange and coherence dephasing due to a fixed number of rotor-synchronously applied pi-pulses is monitored at spinning frequencies in the vicinity of the rotational resonance (R(2)) conditions. The powder nbRFDR and nbSEDRA decay curves of spin magnetizations and coherences, respectively, as a function of the spinning frequency can be measured and analyzed using simple rate equations providing a quantitative measure of the dipolar coupling. The effects of the phenomenological relaxation parameters in these rate equations are discussed and an improved methodology is suggested for analyzing nbRFDR data for small dipolar couplings. The distance between the labeled nuclei in the 1,3-(13)C(2)-hydroxybutyric acid molecule is rederived using existing nbRFDR results and the new simulation procedure. A nbSEDRA experiment has been performed successfully on a powder sample of singly labeled 1-(13)C-L-leucine measuring the dipolar interaction between the labeled carboxyl carbon and the natural abundant beta-carbon. Both narrowband techniques are employed for the determination of the nuclear distances between the side-chain carbons of leucine and its carbonyl carbon in a tripeptide Leu-Gly-Phe that is singly (13)C-labeled at the leucine carbonyl carbon position.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号