Traveling wave front for the Fisher equation on an infinite band region

Instead of discussing the existence of a one-dimensional traveling wave front solution which connects two constant steady states, the present work deals with the case connecting a constant and a nonhomogeneous steady state on an infinite band region. The corresponding model is the well-known Fisher equation with variational coefficient and Dirichlet boundary condition. (c) 2006 Elsevier Ltd. All rights reserved.

基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

APPLICATION OF THE KHOKHLOV-ZABOLOTSKAYA-KUZNETSOV EQUATION TO MODELING HIGH-INTENSITY FOCUSED ULTRASOUND BEAMS

High-intensity focused ultrasound is a form of therapeutic ultrasound which uses high amplitude acoustic waves to heat and ablate tissue. HIFU employs acoustic amplitudes that are high enough that nonlinear propagation effects are important in the evolution of the sound field. A common model for HIFU beams is the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation which accounts for nonlinearity, diffraction, and absorption. The KZK equation models diffraction using the parabolic or paraxial approximation. Many HIFU sources have an aperture diameter similar to the focal length and the paraxial approximation may not be appropriate. Here, results obtained using the “Texas code,” a time-domain numerical solution to the KZK equation, were used to assess when the KZK equation can be employed. In a linear water case comparison with the O’Neil solution, the KZK equation accurately predicts the pressure field in the focal region. The KZK equation was also compared to simulations of the exact fluid dynamics equations (no paraxial approximation). The exact equations were solved using the Fourier-Continuation (FC) method to approximate derivatives in the equations. Results have been obtained for a focused HIFU source in tissue. For a low focusing gain transducer (focal length 50λ and radius 10λ), the KZK and FC models showed excellent agreement, however, as the source radius was increased to 30λ, discrepancies started to appear. Modeling was extended to the case of tissue with the appropriate power law using a relaxation model. The relaxation model resulted in a higher peak pressure and a shift in the location of the peak pressure, highlighting the importance of employing the correct attenuation model. Simulations from the code that were compared to experimental data in water showed good agreement through the focal plane.

Iterative Linearized Density Matrix Propagation for Modeling Coherent Energy Transfer in Photosynthetic Light Harvesting

We present results of calculations [1] that employ a new mixed quantum classical iterative density matrix propagation approach (ILDM , or so called Is‐Landmap) [2] to explore the survival of coherence in different photo synthetic models. Our model studies confirm the long lived quantum coherence , while conventional theoretical tools (such as Redfield equation) fail to describe these phenomenon [3,4]. Our ILDM method is a numerical exactly propagation scheme and can be served as a bench mark calculation tools[2]. Result get from ILDM and from other recent methods have been compared and show agreement with each other[4,5]. Long lived coherence plateau has been attribute to the shift of harmonic potential due to the system bath interaction, and the harvesting efficiency is a balance between the coherence and dissipation[1]. We use this approach to investigate the excitation energy transfer dynamics in various light harvesting complex include Fenna‐Matthews‐Olsen light harvesting complex[1] and Cryptophyte Phycocyanin 645 [6]. [1] P.Huo and D.F.Coker ,J. Chem. Phys. 133, 184108 (2010) . [2] E.R. Dunkel, S. Bonella, and D.F. Coker, J. Chem. Phys. 129, 114106 (2008). [3] A. Ishizaki and G.R. Fleming, J. Chem. Phys. 130, 234111 (2009). [4] A. Ishizaki and G.R. Fleming, Proc. Natl. Acad. Sci. 106, 17255 (2009). [5] G. Tao and W.H. Miller, J. Phys. Chem. Lett. 1, 891 (2010). [6] P.Huo and D.F.Coker in preparation

PROPAGATION OF SONIC BOOMS THROUGH A REAL, STRATIFIED ATMOSPHERE

Sonic boom propagation in a quiet) stratified) lossy atmosphere is the subject of this dissertation. Two questions are considered in detail: (1) Does waveform freezing occur? (2) Are sonic booms shocks in steady state? Both assumptions have been invoked in the past to predict sonic boom waveforms at the ground. A very general form of the Burgers equation is derived and used as the model for the problem. The derivation begins with the basic conservation equations. The effects of nonlinearity) attenuation and dispersion due to multiple relaxations) viscosity) and heat conduction) geometrical spreading) and stratification of the medium are included. When the absorption and dispersion terms are neglected) an analytical solution is available. The analytical solution is used to answer the first question. Geometrical spreading and stratification of the medium are found to slow down the nonlinear distortion of finite-amplitude waves. In certain cases the distortion reaches an absolute limit) a phenomenon called waveform freezing. Judging by the maturity of the distortion mechanism, sonic booms generated by aircraft at 18 km altitude are not frozen when they reach the ground. On the other hand, judging by the approach of the waveform to its asymptotic shape, N waves generated by aircraft at 18 km altitude are frozen when they reach the ground. To answer the second question we solve the full Burgers equation and for this purpose develop a new computer code, THOR. The code is based on an algorithm by Lee and Hamilton (J. Acoust. Soc. Am. 97, 906-917, 1995) and has the novel feature that all its calculations are done in the time domain, including absorption and dispersion. Results from the code compare very well with analytical solutions. In a NASA exercise to compare sonic boom computer programs, THOR gave results that agree well with those of other participants and ran faster. We show that sonic booms are not steady state waves because they travel through a varying medium, suffer spreading, and fail to approximate step shocks closely enough. Although developed to predict sonic boom propagation, THOR can solve other problems for which the extended Burgers equation is a good propagation model.

Multiscale numerical modelling of crack propagation in two-dimensional metal plate

A novel multiscale model of brittle crack propagation in an Ag plate with macroscopic dimensions has been developed. The model represents crack propagation as stochastic drift-diffusion motion of the crack tip atom through the material, and couples the dynamics across three different length scales. It integrates the nanomechanics of bond rupture at the crack tip with the displacement and stress field equations of continuum based fracture theories. The finite element method is employed to obtain the continuum based displacement and stress fields over the macroscopic plate, and these are then used to drive the crack tip forward at the atomic level using the molecular dynamics simulation method based on many-body interatomic potentials. The linkage from the nanoscopic scale back to the macroscopic scale is established via the Ito stochastic calculus, the stochastic differential equation of which advances the tip to a new position on the macroscopic scale using the crack velocity and diffusion constant obtained on the nanoscale. Well known crack characteristics, such as the roughening transitions of the crack surfaces, crack velocity oscillations, as well as the macroscopic crack trajectories, are obtained.

A multi-scale numerical modelling of crack propagation in a 2D metallic plate

A new multi-scale model of brittle fracture growth in an Ag plate with macroscopic dimensions is proposed in which the crack propagation is identified with the stochastic drift-diffusion motion of the crack-tip atom through the material. The model couples molecular dynamics simulations, based on many-body interatomic potentials, with the continuum-based theories of fracture mechanics. The Ito stochastic differential equation is used to advance the tip position on a macroscopic scale before each nano-scale simulation is performed. Well-known crack characteristics, such as the roughening transitions of the crack surfaces, as well as the macroscopic crack trajectories are obtained.

Computer simulation of crack propagation in power electronics module solder joints

A numerical modelling method for the analysis of solder joint damage and crack propagation has been described in this paper. The method is based on the disturbed state concept. Under cyclic thermal-mechanical loading conditions, the level of damage that occurs in solder joints is assumed to be a simple monotonic scalar function of the accumulated equivalent plastic strain. The increase of damage leads to crack initiation and propagation. By tracking the evolution of the damage level in solder joints, crack propagation path and rate can be simulated using Finite Element Analysis method. The discussions are focused on issues in the implementation of the method. The technique of speeding up the simulation and the mesh dependency issues are analysed. As an example of the application of this method, crack propagation in solder joints of power electronics modules under cyclic thermal-mechanical loading conditions has been analyzed and the predicted cracked area size after 3000 loading cycles is consistent with experimental results.

Human respiration rate estimation using ultra-wideband distributed cognitive radar system

It has been shown that remote monitoring of pulmonary activity can be achieved using ultra-wideband (UWB) systems, which shows promise in home healthcare, rescue, and security applications. In this paper, we first present a multi-ray propagation model for UWB signal, which is traveling through the human thorax and is reflected on the air/dry-skin/fat/muscle interfaces. A geometry-based statistical channel model is then developed for simulating the reception of UWB signals in the indoor propagation environment. This model enables replication of time-varying multipath profiles due to the displacement of a human chest. Subsequently, a UWB distributed cognitive radar system (UWB-DCRS) is developed for the robust detection of chest cavity motion and the accurate estimation of respiration rate. The analytical framework can serve as a basis in the planning and evaluation of future measurement programs. We also provide a case study on how the antenna beamwidth affects the estimation of respiration rate based on the proposed propagation models and system architecture

Effect Of Temperature On The Respiration Rate Of Meiofauna

The effect of temperature on respiration rate has been established, using Cartesian divers, for the meiofaunal sabellid polychaeteManayunkia aestuarina, the free-living nematodeSphaerolaimus hirsutus and the harpacticoid copepodTachidius discipes from a mudflat in the Lynher estuary, Cornwall, U.K. Over the temperature range normally experienced in the field, i.e. 5–20° C the size-compensated respiration rate (R c) was related to the temperature (T) in °C by the equation Log10 R c=-0.635+0.0339T forManayunkia, Log10 R c=0.180+0.0069T forSphaerolaimus and Log10 R c=-0.428+0.0337T forTachidius, being equivalent toQ 10 values of 2.19, 1.17 and 2.17 respectively. In order to derive the temperature response forManayunkia a relationship was first established between respiration rate and body size: Log10 R=0.05+0.75 Log10 V whereR=respiration in nl·O2·ind-1·h-1 andV=body volume in nl. TheQ 10 values are compared with values for other species derived from the literature. From these limited data a dichotomy emerges: species with aQ 10≏2 which apparently feed on diatoms and bacteria, the abundance of which are subject to large short term variability, and species withQ 10≏1 apparently dependent on more stable food sources.

Reduced estimated glomerular filtration rate in Alzheimer's disease

OBJECTIVES:

Renal disease is increasingly regarded as an independent risk factor for vascular disease which in itself is believed to influence risk of AD. Alterations in amyloid homeostasis via reduced renal clearance of peripheral beta-amyloid (A|*beta*|) may represent another potential role for variation in renal function leading to increased risk of AD. We sought to examine estimates of glomerular filtration rate in AD and control groups.
METHODS:

AD patients were randomly recruited from the Memory Clinic of the Belfast City Hospital (n = 83). Genomic DNA was extracted from peripheral leucocytes and was genotyped for Apolipoprotein E using standard methods. Using creatinine values, age and gender, estimated Glomerular Filtration Rates (eGFR) were calculated using the isotope dilution mass spectrometry (IDMS)-traceable Modification of Diet in Renal Disease (MDRD) Study equation (using the United Kingdom National External Quality Assessment Scheme (UKNEQAS) correction factor). IDMS eGFR values were then compared between AD and control groups.

RESULTS:

Significant baseline differences in age, diastolic blood pressure, education level attained and APOE |*epsilon*|4 carriage were noted between cases and controls. The AD group had a significantly lower eGFR versus controls (69 vs 77 ml/min) which persisted after adjustment for possible confounders (p = 0.045).

CONCLUSIONS:

This case-control analysis suggests that using a relatively accurate estimate of renal function, patients with AD have greater renal impairment than cognitively normal controls. This may reflect impaired renal clearance of peripheral A|*beta*| or be a marker of shared vascular processes altering cerebral and renal functioning.

The finding of reduced estimated glomerular filtration rate is associated with increased mortality in a large UK population

BACKGROUND: CKD as defined by KDIGO/KDOQI has been shown to affect ~ 8.5% of the UK population. The prevalence of CKD in the UK is similar to that in the USA, yet incident dialysis rates are dramatically different. This retrospective cohort study investigates the association between reduced kidney function and mortality in a large UK population. METHODS: All serum creatinine results covering Northern Ireland's 1.7 million population were collected between 1 January 2001 and 31 December 2002. Estimated glomerular filtration rates (eGFR) were calculated for all serum creatinine measurements using four-variable MDRD equation (IDMS aligned). Patients were followed up for both all-cause and cardiovascular mortality data until the end of December 2006. Patients on renal replacement therapy were excluded. Subgroup analysis in the 75 345 subjects enrolled within a parallel primary care study permitted additional survival analysis with adjustment for traditional cardiovascular risk factors. RESULTS: A total of 1 967 827 serum creatinine results from 533 798 patients were collected. During the period of follow-up, 59 980 deaths occurred. In multivariate survival analysis, using eGFR as a time-varying covariate, a graded association between CKD (defined by eGFR) and all-cause mortality was identified. Compared with participants with an eGFR of > 60 mL/min/1.73 m(2), the adjusted hazard ratios (and 95% confidence intervals) for participants with an eGFR of 45-59 mL/min/1.73 m(2) was 1.02 (0.99-1.04), an eGFR of 30-44 mL/min/1.73 m(2) was 1.44 (1.40-1.47), an eGFR of 15-29 mL/min/1.73 m(2) was 2.12 (2.05-2.20) and an eGFR of

Oblique amplitude modulation of dust-acoustic plasma waves

Theoretical and numerical studies are presented of the nonlinear amplitude modulation of dust-acoustic (DA) waves propagating in an unmagnetized three component, weakly-coupled, fully ionized plasma consisting of electrons, positive ions and charged dust particles, considering perturbations oblique to the carrier wave propagation direction. The stability analysis, based on a nonlinear Schrodinger-type equation (NLSE), shows that the wave may become unstable; the stability criteria depend on the angle theta between the modulation and propagation directions. Explicit expressions for the instability rate and threshold have been obtained in terms of the dispersion laws of the system. The possibility and conditions for the existence of different types of localized excitations have also been discussed.

Nonlinear propagation of electromagnetic waves in negative-refraction-index composite materials

We investigate the nonlinear propagation of electromagnetic waves in left-handed materials. For this purpose, we consider a set of coupled nonlinear Schrodinger (CNLS) equations, which govern the dynamics of coupled electric and magnetic field envelopes. The CNLS equations are used to obtain a nonlinear dispersion, which depicts the modulational stability profile of the coupled plane-wave solutions in left-handed materials. An exact (in)stability criterion for modulational interactions is derived, and analytical expressions for the instability growth rate are obtained.

Nonlinear propagation of modulated ion-acoustic plasma waves in the presence of an electron beam

Theoretical and numerical studies are presented of the amplitude modulation of ion-acoustic waves (IAWs) in a plasma consisting of warm ions, Maxwellian electrons, and a cold electron beam. Perturbations parallel to the carrier IAW propagation direction have been investigated. The existence of four distinct linear ion acoustic modes is shown, each of which possesses a different behavior from the modulational stability point of view. The stability analysis, based on a nonlinear Schrodinger equation (NLSE) reveals that the IAW may become unstable. The stability criteria depend on the IAW carrier wave number, and also on the ion temperature, the beam velocity and the beam electron density. Furthermore, the occurrence of localized envelope structures (solitons) is investigated, from first principles. The numerical analysis shows that the two first modes (essentially IAWs, modified due to the beam) present a complex behavior, essentially characterized by modulational stability for large wavelengths and instability for shorter ones. Dark-type envelope excitations (voids, holes) occur in the former case, while bright-type ones (pulses) appear in the latter. The latter two modes are characterized by an intrinsic instability, as the frequency develops a finite imaginary part for small ionic temperature values. At intermediate temperatures, both bright- and dark-type excitations may exist, although the numerical landscape is intertwined between stability and instability regions.(c) 2006 American Institute of Physics.