Research article

Hydrothermal carbonization of glucose in saline solution: sequestration of nutrients on carbonaceous materials

  • Received: 17 November 2015 Accepted: 28 January 2016 Published: 02 February 2016
  • In this study, feasibility of selected nutrient sequestration during hydrothermal carbonization (HTC) was tested for three different HTC temperatures (180, 230, and 300 °C). To study the nutrient sequestration in solid from liquid solution, sugar and salt solutions were chosen as HTC feedstock. Glucose was used as carbohydrate source and various salts e.g., ammonium hydrophosphate, potassium chloride, potassium sulfate, and anhydrous ferric chloride were used as source of nitrogen and phosphorus, potassium, and iron, respectively. Solid hydrochar was extensively characterized by means of elemental, ICP-OES, SEM-EDX, surface area, pore volume and size, and ATR-FTIR to determine nutrients’ sequestration as well as hydrochar quality variation with HTC temperatures. The spherical mesoporous hydrochars produced during HTC have low surface area in the range of 1.0–3.5 m2 g−1. Hydrochar yield was increased about 10% with the increase of temperature from 180 °C to 300 °C. Nutrient sequestration was also increased with HTC temperature. In fact, around 71, 31, and 23 wt% nitrogen, iron, and phosphorus were sequestered at 300 °C, respectively. Potassium sequestration was very low throughout the HTC and maximum 5.2% was observed in solid during HTC.

    Citation: M Toufiq Reza, Jessica Nover, Benjamin Wirth, Charles J Coronella. Hydrothermal carbonization of glucose in saline solution: sequestration of nutrients on carbonaceous materials[J]. AIMS Energy, 2016, 4(1): 173-189. doi: 10.3934/energy.2016.1.173

    Related Papers:

    [1] Lei Mao, Yan Zhang . Robust optimal excess-of-loss reinsurance and investment problem with p-thinning dependent risks under CEV model. Quantitative Finance and Economics, 2021, 5(1): 134-162. doi: 10.3934/QFE.2021007
    [2] Sang Phu Nguyen, Toan Luu Duc Huynh . Portfolio optimization from a Copulas-GJR-GARCH-EVT-CVAR model: Empirical evidence from ASEAN stock indexes. Quantitative Finance and Economics, 2019, 3(3): 562-585. doi: 10.3934/QFE.2019.3.562
    [3] Keyue Yan, Ying Li . Machine learning-based analysis of volatility quantitative investment strategies for American financial stocks. Quantitative Finance and Economics, 2024, 8(2): 364-386. doi: 10.3934/QFE.2024014
    [4] James P. Gander . Market Value Volatility and the Volume of Traded Stock for U.S. Industrial Corporations. Quantitative Finance and Economics, 2017, 1(4): 403-410. doi: 10.3934/QFE.2017.4.403
    [5] Yuhang Zheng, Siming Liu . A Dynamic Analysis of the Business Cycle Model with a Fixed-time Delay. Quantitative Finance and Economics, 2017, 1(2): 174-185. doi: 10.3934/QFE.2017.2.174
    [6] Muath Asmar, Susanne Trimbath . Regulatory reform and trade settlement failures in USA equity markets: Does regulatory reform matter?. Quantitative Finance and Economics, 2022, 6(4): 537-552. doi: 10.3934/QFE.2022023
    [7] Lucia Morales, Daniel Rajmil . Investing in virtue and frowning at vice? Lessons from the global economic and financial crisis. Quantitative Finance and Economics, 2023, 7(1): 1-18. doi: 10.3934/QFE.2023001
    [8] Nereida Polovina, Ken Peasnell . The effects of different modes of foreign bank entry in the Turkish banking sector during the 2007–2009 Global financial crisis. Quantitative Finance and Economics, 2023, 7(1): 19-49. doi: 10.3934/QFE.2023002
    [9] Anna Timonina-Farkas . COVID-19: data-driven dynamic asset allocation in times of pandemic. Quantitative Finance and Economics, 2021, 5(2): 198-227. doi: 10.3934/QFE.2021009
    [10] Andrea Ferrario, Massimo Guidolin, Manuela Pedio . Comparing in- and out-of-sample approaches to variance decomposition-based estimates of network connectedness an application to the Italian banking system. Quantitative Finance and Economics, 2018, 2(3): 661-701. doi: 10.3934/QFE.2018.3.661
  • In this study, feasibility of selected nutrient sequestration during hydrothermal carbonization (HTC) was tested for three different HTC temperatures (180, 230, and 300 °C). To study the nutrient sequestration in solid from liquid solution, sugar and salt solutions were chosen as HTC feedstock. Glucose was used as carbohydrate source and various salts e.g., ammonium hydrophosphate, potassium chloride, potassium sulfate, and anhydrous ferric chloride were used as source of nitrogen and phosphorus, potassium, and iron, respectively. Solid hydrochar was extensively characterized by means of elemental, ICP-OES, SEM-EDX, surface area, pore volume and size, and ATR-FTIR to determine nutrients’ sequestration as well as hydrochar quality variation with HTC temperatures. The spherical mesoporous hydrochars produced during HTC have low surface area in the range of 1.0–3.5 m2 g−1. Hydrochar yield was increased about 10% with the increase of temperature from 180 °C to 300 °C. Nutrient sequestration was also increased with HTC temperature. In fact, around 71, 31, and 23 wt% nitrogen, iron, and phosphorus were sequestered at 300 °C, respectively. Potassium sequestration was very low throughout the HTC and maximum 5.2% was observed in solid during HTC.


    1. Introduction

    Investment is an important means of earning profit for insurance companies, reinsurance is in use for controlling insurers' risk exposure. Recently, it became a hot point in optimization and control fields to consider reinsurance and investment synthetically under risk models. There are two main models of reinsurance-investment problem, that is, maximizing the insurer's expected utility of terminal wealth and minimizing the probability of ruin under reinsurance-investment strategies, and each model is based on the quantitative description for insurer's surplus processes. From the known literature, we see that two major types of mathematical model for the surplus of the insurance company are usually exploited, which are the Cramer-Lundberg model and the diffusion model, the former is the so called classical risk process, and the later is an approximation to the former which is feasible as one deals with large insurance portfolios, where an individual claim is relatively small compared to the size of the surplus.

    In our paper, the continuous-time mean-variance reinsurance-investment problem is to be considered, where multiple risky assets for investment exist in the market and diffusion model is exploited to describe the risk process of the insurance company. Under diffusion model of the risk process, related reinsurance-investment problems have attracted a great deal of interest. Promislow et al. (2005) studied the optimal reinsurance-investment strategy to minimize the probability of ruin for insurance companies. Taksar et al. (2003) also obtained the optimal reinsurance-investment strategy to minimize the probability of ruin for insurance companies under diffusion model of surplus process. By solving the corresponding Hamilton-Jacobi-Bellman equations, Cao et al. (2009) obtained the optimal reinsurance-investment strategies for maximizing the utility of terminal wealth with exponential and power utility functions. Luo et al. (2008) studied the similar problem as Cao et al. (2009) in Black-Scholes market subject to some investment constraints. All of the literature (Browne, 1995; Promislow et al., 2005; Taksar et al., 2003; Cao et al., 2009; Luo et al., 2008) researched their optimization problems in Black-Scholes market with single risky asset. Bai et al. (2008) considered the reinsurance-investment problem in Black-Scholes market with multiple risky assets and no-shorting constraint. Zhang et al. (2009) studied the optimal reinsurance-investment problem of maximizing the expected exponential utility of terminal wealth and solve it by using the corresponding Hamilton-Jacobi-Bellman equation. It is noted that the transaction costs when investing in the risky assets and conditional Value-at-Risk to control the whole risk have been considered.

    For the sack of practical need in management, mean-variance portfolio selection model initiated by Markowitz has attracted much attention all the time. There has been few literature up to now which devoted themselves to studying the continuous-time mean-variance reinsurance-investment problem, however, the continuous-time mean-variance portfolio selection problem was researched by Zhou et al. (2000)with the help of stochastic linear-quadratic method. As the subsequent paper, Lim et al. (2002) pointed out: there are many advantages of using the framework of stochastic linear-quadratic control to study dynamic mean-variance problems, meanwhile it was recognized that the solution to the problem could also have been obtained via dynamic programming and the associated Hamilton-Jacobi-Bellman equation when the market coefficients are assumed to be deterministic. It is observed that, when the market coefficients are assumed to be deterministic, the Lagrange method can be applied to solving the continuous-time mean-variance portfolio selection problems, which was expounded in Chow (1997) in a general way and was applied to investment decisions and pricing of contingent claims in Chow (1996) and Chow (1999), respectively. The Lagrange method is more popular than other dynamic optimization methods. For our continuous-time mean-variance reinsurance-investment problem with deterministic market coefficients, the Lagrange method is also a simple one to obtain our results, which is the main reason for this research. It should be recognized that the Lagrange method is limited to complex circumstances such as random market coefficients, which induces more complex backward stochastic differential equation for solving the process of the Lagrange multiplier. For more details on the subject, we refer to the recent work such as Chen et al. (2013), Shen et al. (2015), and so on.

    This paper is organized as follows: In Section 2, we state the mean-variance optimization problem. In Section 3, we use the Lagrange method to analysis the mean-variance problem, necessary and sufficient conditions for the problem are discussed, the optimal strategy and its effective frontier in explicit forms are obtained. In Section 4, we show a numerical example.


    2. Problem Formulation

    Throughout this paper we denote by $M'$ the transpose of any matrix or vector $M=({{m}_{ij}})$, $||M||=\sqrt{\sum\nolimits_{i, j}{m_{ij}^{2}}}$ its norm, ${{R}^{m}}$ the $m$dimensional real space, 1 the vector with unit entries, $[0, T]$ the finite horizon of investment and ${{E}_{t}}$ the expectation at time $t\in [0, T]$.

    Market uncertainty is modeled by a filtered complete probability space $(\Omega, F, P, {{\{{{F}_{t}}\}}_{t\in [0, T]}})$. Assume that the filtration ${{\{{{F}_{t}}\}}_{t\in [0, T]}}$ is generated by an $(m+1)$-dimensional standard Brownian motion $\{(W_{t}^{0}, W_{t}^{1}, \cdots, W_{t}^{m})':t\in [0, T]\}$ for a positive integer $m$. ${{F}_{0}}=\{\Phi, \Omega \}$, ${{F}_{T}}=F$.

    We recall that the space ${{L}_{F, 2}}([0, T];{{R}^{m}})$ of mean square integrable functions $u(t, x)$ from $R\times \Omega $ to ${{R}^{m}}$ possesses a natural Hilbert space structure, inducing the norm

    $ ||u|{{|}_{F,2}}={{\left( {{E}_{0}}\int_{0}^{T}{||u(t,X)|{{|}^{2}}dt} \right)}^{\frac{1}{2}}}<+\infty . $

    Following the framework of Promislow and Young (2005), we model the insurer's claim process $C$ according to a Brownian motion with drift as follows:

    $ dC(t)=adt-bdW_{t}^{0}, $ (1)

    where $a$ and $b$ are positive constants. Assume that the premium be paid continuously at the constant rate ${{c}_{0}}=(1+\theta)a$ with safety loading $\theta >0$, and that the insurer purchase the proportional reinsurance to reduce the underlying risk. Let $q(t)$ represent the proportion reinsured at time $t$. Similar to the premium paid to the insurer, the reinsurance premium is paid continuously at the constant rate ${{c}_{1}}=(1+\eta)aq$ with safety loading $\eta >\theta >0$. Based on Eq. (1), the surplus process $R(t)$ is given by the dynamics

    $ dR(t)={{c}_{0}}dt-(1-q(t))dC(t)-{{c}_{1}}dt=(\theta -\eta q(t))adt+b(1-q(t))dW_{t}^{0}. $ (2)

    In addition to reinsurance, the insurer is allowed to invest its surplus in a financial market in which $(m+1)$ assets are traded continuously over $[0, T]$ in a self-financing fashion (i.e., there is no consumption nor income). The 0th asset is a bond whose price ${{S}_{0}}(t)$ is subject to the following deterministic ordinary differential equation:

    $ d{{S}_{0}}(t)=r(t){{S}_{0}}(t)dt{{,}_{{}}}{{S}_{0}}(0)=1, $ (3)

    where $r(t)\in C([0, T];{{R}^{+}})$ is the interest rate of the risk-free asset. The remaining $m$ assets are risky and their price processes ${{S}_{1}}(t), {{S}_{2}}(t), \cdots, {{S}_{m}}(t)$ satisfy the following stochastic differential equations:

    $ d{{S}_{i}}(t)={{\mu }_{i}}(t){{S}_{i}}(t)dt+{{\sigma }_{i}}(t){{S}_{i}}(t)d{{W}_{t}},{{S}_{i}}(0)={{s}_{i}}>0,i=1,2,\cdots ,m, $ (4)

    where ${{\mu }_{i}}(t)$ and ${{\sigma }_{i}}(t)=[{{\sigma }_{i1}}(t), \cdots, {{\sigma }_{im}}(t)]$ are the appreciation rate and dispersion (or volatility) rate of the $i$th asset, respectively, and $W=(W_{t}^{1}, \cdots, W_{t}^{m})'$. We assume that ${{\mu }_{i}}(t)$ and ${{\sigma }_{i}}(t)$ are scalar-valued deterministic bounded functions. Denoting

    $ \sigma (t)=\left[ σ1(t)σm(t)
    \right]\in {{R}^{m\times m}}, $
    (5)

    we assume throughout that $\sigma (t)$ is uniformly nondegenerate: that is, there exists $\delta >0$ such that

    $ \sigma (t)\sigma (t)'\ge \delta I,\forall t\in [0,T]. $ (6)

    A strategy $s$ is described by a stochastic process $(q(t), \pi (t)')'$, where $\pi (t)$ $:=({{\pi }_{1}}(t), {{\pi }_{2}}(t), \cdots, {{\pi }_{m}}(t))'$, ${{\pi }_{i}}(t)$ is the dollar amount invested in asset $i(i=1, 2, \cdots, m)$ at time $t$.

    Let $X(t)$ denote the resulting surplus process after incorporating strategy $s$ into (2). The dynamics of $X(t)$ can be represented as follows.

    $ dX(t)=[r(t)X(t)+(\theta -\eta q(t))a+\tau (t)'\pi (t)]dt+b(1-q(t))dW_{t}^{0}+\pi (t)'\sigma (t)d{{W}_{t}}, $ (7)

    where $\tau (t)=({{\mu }_{1}}(t)-r(t), \cdots, {{\mu }_{m}}(t)-r(t))'$, $X(0)={{x}_{0}}$.

    The admissible strategy set with initial wealth ${{x}_{0}}$ is defined as

    $ A({{x}_{0}}):=\{s=(q(t),\pi (t)')'|q(t),\pi (t)satisfies(2.7)\} $

    The objective of the insurer is to determine an optimal strategy $s\in A({{x}_{0}})$ such that the variance of the terminal wealth, $Var[X(T)]$, is minimized, subject to that the expected terminal wealth, ${{E}_{0}}[X(T)]$, takes a given level $K>{{x}_{0}}\exp (\int_{0}^{T}{r(t)dt})$. This can be expressed as the optimization problem:

    $ \left\{ minsA(x0)Var[X(T)]=E0[X2(T)K],s.t.E0[X(T)]=K.
    \right. $
    (8)

    We introduce a Lagrange multiplier $2\xi \in R$ and after rearranging terms arrive at the new cost function,

    $ {{E}_{0}}[{{(X{{(T)}^{2}}-K)}^{2}}+2\xi (X(T)-K)]={{E}_{0}}{{[X(T)-(K-\xi )]}^{2}}-{{\xi }^{2}}. $

    Letting $\gamma =K-\xi $ leads to the following optimal stochastic control problem,

    $ \mathop {\max }\limits_{s \in A({x_0})} {\mkern 1mu} \{ {(K - \gamma )^2} - {E_0}{[X{(T)^2} - \gamma ]^2}\} . $ (9)

    Remark 2.1. The link between problem (8) and (9) is provided by the Lagrange duality theorem, see e.g. Luenberger (1968)

    $ \mathop {\min }\limits_{s \in A({x_0})} {\mkern 1mu} Var[X(T)] = \mathop {\max }\limits_{\gamma \in R} {\mkern 1mu} \mathop {\max }\limits_{s \in A({x_0})} {\mkern 1mu} \{ {(K - \gamma )^2} - {E_0}{[X{(T)^2} - \gamma ]^2}\} . $

    3. Mean-Variance Analysis by the Lagrange Method

    For the first step, we consider problem (9) for a fixed $\gamma $. Follow the Lagrange method introduced by Chow (1997), recognizing $dX(t)=X(t+dt)-X(t)$, and let $dt$ be a very small time interval. The Lagrange expression for the optimization problem (9) beginning from $t=0$ is

    $ L={{(K-\gamma )}^{2}}-{{E}_{t}}{{[X{{(T)}^{2}}-\gamma ]}^{2}} $
    $ -\int_{0}^{T}{{{E}_{t}}\lambda (t+dt)[X(t+dt)-X(t)-(r(t)X(t)+(\theta -\eta q(t))}a+\tau (t)'\pi (t))dt $
    $ -b(1-q(t))dW_{t}^{0}-\pi (t)'\sigma (t)d{{W}_{t}}]. $ (10)

    For simplicity, denote by

    $ f(X(t),q(t),{{\pi }_{1}}(t),\cdots ,{{\pi }_{m}}(t))=r(t)X(t)+(\theta -\eta q(t))a+\tau (t)'\pi (t), $
    $ R(q(t))=b(1-q(t)). $

    The Lagrange expression (10) can be rewritten as follows.

    $ L={{(K-\gamma )}^{2}}-{{E}_{t}}{{[X{{(T)}^{2}}-\gamma ]}^{2}} $
    $ -\int_{0}^{T}{{{E}_{t}}\lambda (t+dt)[X(t+dt)-}fdt-R(q)dW_{t}^{0}-\pi (t)'\sigma (t)d{{W}_{t}}]. $ (11)

    Remark 3.1: By theorem 4.2 of Yong and Zhou(1999), problem (9) with dynamic constraint (7) admits a pathwise unique admissible strategy in $A({{x}_{0}})$, which implies, for a given $\gamma \in R$, maximization of the cost function $L$ in (11) is uniquely solvable, henceforth, the Lagrange multiplier $\lambda =\lambda (t, X(t))$ can be uniquely determined.


    3.1. Necessary conditions

    Let the function $\lambda $ be in the form $\lambda =\lambda (t, X(t))$, applying Ito's differentiation rule to evaluate $d\lambda $ yields

    $ d\lambda =\left[ \frac{\partial \lambda }{\partial t}+\frac{\partial \lambda }{\partial X}f+\frac{1}{2}\frac{{{\partial }^{2}}\lambda }{\partial {{X}^{2}}}\left( {{R}^{2}}(q)+||\pi (t)'\sigma (t)|{{|}^{2}} \right) \right]dt $
    $ +\frac{\partial \lambda }{\partial X}R(q)dW_{t}^{0}+\frac{\partial \lambda }{\partial X}\pi (t)'\sigma (t)d{{W}_{t}}. $ (12)

    Differentiating (11) with respect to $X(T)$, we have

    $ \lambda (T)=-2X(T)+2\gamma. $ (13)

    Differentiating (11) with respect to $q(t)$ gives

    $ \frac{\partial L}{\partial q}=\frac{\partial f}{\partial q}{{E}_{t}}\lambda (t+dt)dt+{{E}_{t}}dW_{t}^{0}\frac{\partial R}{\partial q}\lambda (t+dt) $
    $ =\frac{\partial f}{\partial q}(\lambda dt+{{E}_{t}}d\lambda dt)+{{E}_{t}}dW_{t}^{0}\frac{\partial R}{\partial q}\lambda (t+dt) $
    $ =\left[ \frac{\partial f}{\partial q}\lambda (t)+\frac{\partial R}{\partial q}\frac{\partial \lambda }{\partial X}R(q) \right]dt+o(dt) $
    $ =\left[ -a\eta \lambda (t)-bR(q)\frac{\partial \lambda }{\partial X} \right]dt+o(dt)=0. $

    Hence

    $ a\eta \lambda (t)+bR(q)\frac{\partial \lambda }{\partial X}=0. $ (14)

    Differentiating (11) with respect to $\pi (t)$ yields

    $ \frac{\partial L}{\partial \pi }=\left[ \lambda (t)\tau (t)+\frac{\partial \lambda }{\partial X}\sigma (t)\sigma (t)'\pi (t) \right]dt+o(dt)=0. $

    Hence

    $ \lambda (t)\tau (t)+\frac{\partial \lambda }{\partial X}\sigma (t)\sigma (t)'\pi (t)=0. $ (15)

    Noting the term $\lambda (t)[X(t)-\cdots ]$ in the integrand of (11), differentiating (11) with respect to $X(t)(0\le t < T)$ gives

    $ \frac{\partial L}{\partial X}=-\lambda (t)+{{E}_{t}}\lambda (t+dt)+\frac{\partial f}{\partial X}{{E}_{t}}\lambda (t+dt)dt $
    $ =-\lambda (t)+{{E}_{t}}(\lambda (t)+d\lambda (t))+\frac{\partial f}{\partial X}{{E}_{t}}(\lambda (t)+d\lambda (t))dt $
    $ =\left[ \frac{\partial \lambda }{\partial t}+\frac{\partial \lambda }{\partial X}f+\frac{\partial f}{\partial X}\lambda (t)+\frac{1}{2}\frac{{{\partial }^{2}}\lambda }{\partial {{X}^{2}}}\left( {{R}^{2}}(q)+||\pi (t)'\sigma (t)|{{|}^{2}} \right) \right]dt+o(dt) $
    $ =0. $

    Hence

    $ \frac{\partial \lambda }{\partial t}+\frac{\partial \lambda }{\partial X}f+\frac{\partial f}{\partial X}\lambda (t)+\frac{1}{2}\frac{{{\partial }^{2}}\lambda }{\partial {{X}^{2}}}\left( {{R}^{2}}(q)+||\pi (t)'\sigma (t)|{{|}^{2}} \right)=0. $ (16)

    Now we state the necessary conditions for the optimization problem (9) as follows.

    Proposition 3.1. Let $s=(q(t), \pi (t)')'$ be the optimal reinsurance-investment strategy for the optimization problem (9), then the equations (12)-(16) hold.


    3.2. Sufficient conditions

    First, we give the explicit expression of the wealth process $X(t)$. Let

    $ Y(t)=X(t)\exp \left( -\int_{0}^{t}{r(s)ds} \right), $

    which yields

    $ dY(t)=\exp \left( -\int_{0}^{t}{r(s)ds} \right)[dX(t)-r(t)X(t)dt], $

    In view of (7), we have

    $ dY(t)=\exp \left( -\int_{0}^{t}{r(s)ds} \right)\left[ ((\theta -\eta q(t))a+\tau (t)'\pi (t))dt+R(q)dW_{t}^{0}+\pi (t)'\sigma (t)d{{W}_{t}} \right], $

    and thus

    $ Y(t)-Y(0)=\int_{0}^{t}{\exp \left( -\int_{0}^{s}{r(u)du} \right)}\left[ ((\theta -\eta q(s))a+\tau (s)'\pi (s))ds+R(q)dW_{s}^{0}+\pi (s)'\sigma (s)d{{W}_{s}} \right]. $

    So, we have

    $ X(t)={{x}_{0}}\exp \left( \int_{0}^{t}{r(s)ds} \right)+\int_{0}^{t}{\exp \left( \int_{s}^{t}{r(u)du} \right)}[(\theta -\eta q(s))a+\tau (s)'\pi (s))ds. $
    $ +R(q)dW_{s}^{0}+\pi (s)'\sigma (s)dW. $ (17)

    In view of (17), it is not difficulty to find that both ${{E}_{t}}[\partial X{{(T)}^{2}}/\partial {{q}^{2}}]$ and ${{E}_{t}}[\partial X{{(T)}^{2}}/\partial {{\pi }_{i}}^{2}]$ are positive, which guarantees that the necessary conditions stated in proposition 3.1 are also sufficient conditions for the optimization problem (9) according to Chow (1997).


    3.3. Optimal reinsurance-investment strategy

    Noting that $f=r(t)X(t)(\theta -\eta q(t))a+(\theta -\eta q(t))a+\tau (t)'\pi (t)$, from (3.7) and (3.4), we have

    $ \left\{ λt+[r(t)X(t)+(θηq(t))a+τ(t)π(t)]λX    12[R2(q)+||π(t)σ(t)||2]2λX2+r(t)λ(t)=0,λ(T)=2X(T)+2γ.
    \right. $
    (18)

    According to Remark 3.1, problem (18) has an unique solution $\lambda (t, X(t))$. So, in view of (14) and (15), the components of admissible strategy, $q(t)$ and $\pi (t)$, are uniquely determined. For purpose of getting the unique solution of (18), we start by finding $q(t)$ and $\pi (t)$ in linear feedback form as follows.

    $ q(t)={{\phi }_{1}}(t)X(t)+{{\varphi }_{1}}(t),\pi (t)={{\phi }_{2}}(t)X(t)+{{\varphi }_{2}}(t), $ (19)

    where ${{\phi }_{1}}, {{\varphi }_{1}}\in {{C}^{1}}[0, T;R]$ and ${{\phi }_{2}}, {{\varphi }_{2}}\in {{C}^{1}}[0, T;{{R}^{m}}]$ are deterministic functions to be determined.

    From (14) and (19), remarking that $R(q)=b(1-q(t))$, we have

    $ a\eta \lambda (t)+{{b}^{2}}(1-{{\phi }_{1}}(t)X(t)-{{\varphi }_{1}}(t))\frac{\partial \lambda }{\partial X}=0. $ (20)

    Differentiating (20) with respect to $X(t)$ arrives at

    $ (a\eta -{{b}^{2}}{{\phi }_{1}}(t))\frac{\partial \lambda }{\partial X}+bR(q)\frac{{{\partial }^{2}}\lambda }{\partial {{X}^{2}}}=0. $

    Hence

    $ {{R}^{2}}(q)\frac{{{\partial }^{2}}\lambda }{\partial {{X}^{2}}}=({{\phi }_{1}}(t)-a\eta /{{b}^{2}})bR(q)\frac{{{\partial }^{2}}\lambda }{\partial {{X}^{2}}}=-a\eta ({{\phi }_{1}}(t)-a\eta /{{b}^{2}})\lambda, $ (21)

    where (14) is incorporated in the second equality in (21).

    From (15) and (19), we have

    $ \lambda (t)\tau (t)+\frac{\partial \lambda }{\partial X}\sigma (t)\sigma (t)'({{\phi }_{2}}(t)X(t)+{{\varphi }_{2}}(t))=0. $ (22)

    Differentiating (22) with respect to $X(t)$ arrives at

    $ \frac{\partial \lambda }{\partial X}(\tau (t)+\sigma (t)\sigma (t)'{{\phi }_{2}}(t))+\frac{{{\partial }^{2}}\lambda }{\partial {{X}^{2}}}\sigma (t)\sigma (t)'\pi (t)=0. $ (23)

    Hence

    $ ||\pi (t)'\sigma (t)|{{|}^{2}}\frac{{{\partial }^{2}}\lambda }{\partial {{X}^{2}}}=-\frac{\partial \lambda }{\partial X}\pi (t)'(\tau (t)+\sigma (t)\sigma (t)'{{\phi }_{2}}(t)). $
    $ =-\frac{\partial \lambda }{\partial X}\tau (t)'\pi (t)+\lambda (t){{\phi }_{2}}(t)\tau (t), $ (24)

    where (15) is incorporated in the second equality in (23). Substituting from (21) and (24) into (18) gives

    $ \left\{ λt+[r(t)X(t)+(θηq(t))a+12τ(t)π(t)]λX    +12[2r(t)aη(ϕ1(t)aη/b2)+ϕ2(t)τ(t)]λ=0,λ(T)=2X(T)+2γ.
    \right. $
    (25)

    Characteristic equation of the first order partial differential equation in (25) is

    $ \frac{dt}{1}=\frac{dX}{r(t)X(t)+(\theta -\eta q(t))a+\frac{1}{2}\tau (t)'\pi (t)} $
    $ =\frac{d\lambda }{-\frac{1}{2}(2r(t)-a\eta ({{\phi }_{1}}(t)-a\eta /{{b}^{2}})+{{\phi }_{2}}(t)'\tau (t))} $

    In view of

    $ \frac{dt}{1}=\frac{dX}{r(t)X(t)+(\theta -\eta q(t))a+\frac{1}{2}\tau (t)'\pi (t)} $

    we obtain a first integration as follows.

    $ {{C}_{1}}=X(t)\exp \left( -\int_{0}^{t}{(r(s)-a\eta {{\phi }_{1}}(s)+\frac{1}{2}\tau (s)'{{\phi }_{2}}(s)ds} \right) $
    $ -\int_{t}^{T}{(a\theta -a\eta {{\varphi }_{1}}(u)+}\frac{1}{2}\tau (u)'{{\phi }_{2}}(u))\times \exp \left( -\int_{0}^{u}{(r(s)-a\eta {{\phi }_{1}}(s)+\frac{1}{2}\tau (s)'{{\phi }_{2}}(s)ds} \right)du, $

    where ${{C}_{1}}$ is constant. In view of

    $ \frac{dt}{1}=\frac{d\lambda }{-\frac{1}{2}(2r(t)-a\eta ({{\phi }_{1}}(t)-a\eta /{{b}^{2}})+{{\phi }_{2}}(t)'\tau (t))}, $

    we get another first integration:

    $ {{C}_{2}}=\lambda (t)+\frac{1}{2}\int_{t}^{T}{[2r(s)-a\eta ({{\phi }_{1}}(s)-a\eta /{{b}^{2}})+{{\phi }_{2}}(s)'\tau (s)]ds}, $

    where ${{C}_{2}}$ is constant. It is obvious that both the first integrations obtained above are mutual independent. Invoking the boundary condition in (25), that is, $\lambda (T)=-2X(T)+2\gamma $, relationship between ${{C}_{1}}$ and ${{C}_{2}}$ is

    $ {{C}_{2}}=-2\exp \left( \int_{0}^{T}{(r(s)-a\eta {{\phi }_{1}}(s)+\frac{1}{2}\tau (s)'{{\phi }_{2}}(s))ds} \right)\cdot {{C}_{1}}+2\gamma $

    so we get the solution of (25) as following.

    $ \lambda (t)+\frac{1}{2}\int_{t}^{T}{[2r(s)-a\eta ({{\phi }_{1}}(s)-a\eta /{{b}^{2}})+{{\phi }_{2}}(s)'\tau (s)]ds} $
    $ =-2\exp \left( \int_{0}^{T}{(r(s)-a\eta {{\phi }_{1}}(s)+\frac{1}{2}\tau (s)'{{\phi }_{2}}(s))ds} \right) $
    $ \times \left[ X(t)\exp \left( -\int_{0}^{t}{(r(s)-a\eta {{\phi }_{1}}(s)+\frac{1}{2}\tau (s)'{{\phi }_{2}}(s)ds} \right) \right. $
    $ -\int_{t}^{T}{(a\theta -a\eta {{\varphi }_{1}}(u)+\frac{1}{2}\tau (u)'{{\varphi }_{2}}(u))} $
    $ \times \exp \left. \left( -\int_{0}^{u}{(r(s)-a\eta {{\phi }_{1}}(s)+\frac{1}{2}\tau (s)'{{\varphi }_{2}}(s)} \right)du \right]+2\gamma. $ (26)

    From (26), we see that

    $ \left\{ λX=2exp(Tt(r(s)aηϕ1(s)+12τ(s)ϕ2(s)ds)0,2λX2=0.
    \right. $
    (27)

    By (21) and the second equation in (27), we get ${{\phi }_{1}}(t)=a\eta /{{b}^{2}}$. By (23) and (27), we get ${{\phi }_{2}}(t)=-{{(\sigma (t)\sigma (t)')}^{-1}}\tau (t)$. Substituting from the results of ${{\phi }_{1}}(t)$ and ${{\phi }_{2}}(t)$ into (27) shows

    $ \frac{\partial \lambda }{\partial X}=-2\exp \left( \int_{t}^{T}{(r(s)-\frac{{{a}^{2}}{{\eta }^{2}}}{{{b}^{2}}}-\frac{1}{2}||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})ds} \right):=H(t). $ (28)

    Substituting from the results of ${{\phi }_{1}}(t)$ and ${{\phi }_{2}}(t)$ into (26) gives

    $ \lambda (t)=H(t)X(t)+h(t), $ (29)

    where

    $ h(t)=2\int_{t}^{T}{(a\theta -a\eta {{\varphi }_{1}}(u)+\frac{1}{2}\tau (u)'{{\varphi }_{2}}(u))} $
    $ \times \exp \left( \int_{u}^{T}{(r(s)-\frac{{{a}^{2}}{{\eta }^{2}}}{{{b}^{2}}}-\frac{1}{2}||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})ds} \right)du $
    $ -\frac{1}{2}\int_{t}^{T}{[2r(s)-||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}}]ds+2\gamma }. $

    Then $h(T)=2\gamma $, and

    $ \frac{dh}{dt}=\left[ (a\theta -a\eta {{\varphi }_{1}}(t)+\frac{1}{2}\tau (t)'{{\varphi }_{2}}(t)) \right]H(t) $
    $ +\frac{1}{2}[2r(s)-||{{\sigma }^{-1}}(t)\tau (t)|{{|}^{2}}], $ (30)

    where $H(t)$ is given by (28).

    By (22), (29) and ${{\phi }_{2}}(t)=-{{(\sigma (t)\sigma (t)')}^{-1}}\tau (t)$, we have

    $ H(t){{\varphi }_{2}}(t)=-h(t){{(\sigma (t)\sigma (t)')}^{-1}}\tau (t). $ (31)

    By (20), (29) and ${{\phi }_{1}}(t)=a\eta /{{b}^{2}}$, we have

    $ H(t){{\varphi }_{1}}(t)=\frac{a\eta }{{{b}^{2}}}h(t)+H(t). $ (32)

    Substituting from (31) and (32) into (30) gives

    $ \frac{dh}{dt}=-\left( \frac{{{a}^{2}}{{\eta }^{2}}}{{{b}^{2}}}+\frac{1}{2}||{{\sigma }^{-1}}(t)\tau (t)|{{|}^{2}} \right)h(t)+a\theta H(t)+\frac{1}{2}[2r(t)-||{{\sigma }^{-1}}(t)\tau (t)|{{|}^{2}}]. $

    Remarking that $h(T)=2\gamma $, we have

    $ h(t)=2\gamma \exp \left( \frac{{{a}^{2}}{{\eta }^{2}}}{{{b}^{2}}}(T-t)+\frac{1}{2}\int_{t}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}}ds} \right) $
    $ -\int_{t}^{T}{\exp \left( \frac{{{a}^{2}}{{\eta }^{2}}}{{{b}^{2}}}(u-t)+\frac{1}{2}\int_{t}^{u}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}}ds} \right)} $
    $ \times \left[ a\theta H(u)+(r(u)-\frac{1}{2}||{{\sigma }^{-1}}(u)\tau (u)|{{|}^{2}} \right]du. $ (33)

    From (31) and (32), we have

    $ {{\varphi }_{2}}(t)=-\frac{h(t)}{H(t)}{{(\sigma (t)\sigma (t)')}^{-1}}\tau (t),{{\varphi }_{1}}(t)=\frac{a\eta }{{{b}^{2}}}\frac{h(t)}{H(t)}+1, $

    so we have

    $ \left\{ q(t)=ϕ1(t)X(t)+φ1(t)=aηb(X(t)+h(t)H(t))+1;π(t)=ϕ2(t)X(t)+φ2(t)=(σ(t)σ(t))1τ(t)(X(t)+h(t)H(t)),
    \right. $
    (34)

    and

    $ \lambda (t)=H(t)X(t)+h(t), $ (35)

    where $H(t)$ and $h(t)$ are given by (28) and (33), respectively. We state one of our main results as follows.

    Proposition 3.2. Let $\tau (t)=({{\mu }_{1}}(t)-r(t), \cdots, {{\mu }_{m}}(t)-r(t))'$, and $\gamma \in R$ be given. The optimal reinsurance-investment strategy for problem (2.9) are given by (34), and the Lagrange multiplier $\lambda (t, X(t))$ is a linear function of $X(t)$ which is given by (35).


    3.4. The effective frontier

    From (35), we have

    $ \frac{\partial \lambda }{\partial t}=\frac{dH}{dt}X(t)+\frac{dh}{dt},\frac{\partial \lambda }{\partial X}=H(t),\frac{{{\partial }^{2}}\lambda }{\partial {{X}^{2}}}=0. $

    From equations (12) and (13), we have the backward stochastic differential equation (BSDE, for short) with terminal value for the Lagrange multiplier $\lambda $.

    $ \left\{ dλ=[dHdtX(t)+dhdt+H(t)f]dt+H(t)R(q)dW0t+H(t)π(t)σ(t)dWt,λ(T)=2X(T)+2γ.
    \right. $
    (36)

    From (18) we have $\frac{dH}{dt}X(t)+\frac{dh}{dt}+H(t)f=-\lambda r(t)$. Substituting from (34) into (36), one obtains

    $ \left\{ dλ=λ[r(t)dt+aηbdW0t+(σ1(t)τ(t))dWt],λ(T)=2X(T)+2γ.
    \right. $
    (37)

    BSDE (37) has the following solution.

    $ \lambda =-2(X(T)-\gamma )\exp \left[ \int_{t}^{T}{\left( r(s)-\frac{l(s)}{2} \right)ds} \right] $
    $ \times \exp \left[ \frac{a\eta }{b}(W_{T}^{0}-W_{t}^{0})+\int_{t}^{T}{({{\sigma }^{-1}}(s)\tau (s))'}d{{W}_{s}} \right], $ (38)

    where $l={{a}^{2}}{{\eta }^{2}}/{{b}^{2}}+||{{\sigma }^{-1}}\tau |{{|}^{2}}$. By $\lambda $ obtained in (38), we then have

    $ \lambda (0,X(0))=H(0)X(0)+h(0) $
    $ =-2(X(T)-\gamma )\exp \left[ \int_{0}^{T}{\left( r(t)-\frac{l(t)}{2} \right)dt} \right]\times \exp \left[ \frac{a\eta }{b}W_{T}^{0}+\int_{0}^{T}{({{\sigma }^{-1}}(t)\tau (t))'}d{{W}_{t}} \right]. $

    Hence

    $ X(T)=-\frac{1}{2}(H(0){{x}_{0}}+h(0))\exp \left( -\int_{0}^{T}{\left( r(t)-\frac{l(t)}{2} \right)}dt \right) $
    $ \times \exp \left( -\frac{a\eta }{b}W_{T}^{0}-\int_{0}^{T}{({{\sigma }^{-1}}(t)\tau (t))'d{{W}_{t}}} \right)+\gamma, $

    which implies

    $ {{E}_{0}}[X(T)]=-\frac{1}{2}(H(0){{x}_{0}}+h(0))\exp \left( -\int_{0}^{T}{\left( r(t)-l(t) \right)}dt \right)+\gamma. $

    From (28) and (33), we have

    $ H(0)=-2\exp \left( \int_{0}^{T}{(r(s)-\frac{{{a}^{2}}{{\eta }^{2}}}{{{b}^{2}}}-\frac{1}{2}||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})ds} \right), $

    and

    $ h(0)=2\gamma \exp \left( \frac{{{a}^{2}}{{\eta }^{2}}}{{{b}^{2}}}T+\frac{1}{2}\int_{0}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})ds} \right)-\overset{-}{\mathop{h}}\,(0), $

    where

    $ \overset{-}{\mathop{h}}\,(0)=\int_{0}^{T}{\exp \left( \frac{{{a}^{2}}{{\eta }^{2}}}{{{b}^{2}}}u+\frac{1}{2}\int_{0}^{u}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}}ds} \right)} $
    $ \times \left[ a\theta H(u)+(r(u)-\frac{1}{2}||{{\sigma }^{-1}}(u)\tau (u)|{{|}^{2}}) \right]du. $

    Denoting by $A=H(0){{x}_{0}}-\overset{-}{\mathop{h}}\, (0)$, we have

    $ {{E}_{0}}[X(T)]=-\frac{1}{2}A\exp \left( -\int_{0}^{T}{\left( r(t)-l(t) \right)}dt \right) $
    $ +\gamma \left[ 1-\exp \left( -\int_{0}^{T}{(r(s)-2l(s))ds-\frac{1}{2}\int_{0}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})ds}} \right) \right]. $ (39)

    Furthermore, we have

    $ {{X}^{2}}(T)=\frac{1}{4}{{(H(0){{x}_{0}}+h(0))}^{2}}\exp \left( -\int_{0}^{T}{\left( 2r(t)-l(t) \right)}dt \right) $
    $ \times \exp \left( -\frac{2a\eta }{b}W_{T}^{0}-2\int_{0}^{T}{({{\sigma }^{-1}}(t)\tau (t))'d{{W}_{t}}} \right)+{{\gamma }^{2}} $
    $ -\gamma (H(0){{x}_{0}}+h(0))\exp \left( -\int_{0}^{T}{\left( r(t)-\frac{l(t)}{2} \right)}dt \right) $
    $ \times \exp \left( -\frac{a\eta }{b}W_{T}^{0}-\int_{0}^{T}{({{\sigma }^{-1}}(t)\tau (t))'d{{W}_{t}}} \right), $

    which yields

    $ {{E}_{0}}[{{X}^{2}}(T)]=\frac{1}{4}{{A}^{2}}\exp \left( -\int_{0}^{T}{\left( 2r(t)-3l(t) \right)}dt \right) $
    $ +A\left[ \exp \left( -\int_{0}^{T}{(2r(s)-4l(s))ds-\frac{1}{2}\int_{0}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})d}s} \right)-\exp \left( -\int_{0}^{T}{(r(s)-l(s))ds} \right) \right]\gamma $
    $ +\left[ 1+\exp \left( -\int_{0}^{T}{(2r(s)-5l(s))ds-\int_{0}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})ds}} \right) \right] $
    $ -2\exp \left( -\int_{0}^{T}{(r(s)-2l(s))ds-\frac{1}{2}\int_{0}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})d}s} \right){{\gamma }^{2}}. $ (40)

    Using (39) and (40) we can thus give an explicit expression for the optimal cost in problem (8), as a one parameter family in $\gamma $

    $ (Kγ)2E0[X(T)γ]2=K22γK+2γE0[X(T)]E0[X2(T)]=K214A2exp(T0(2r(t)3l(t))dt)  [2K+Aexp(T0(2r(s)4l(s))ds12T0||σ1(s)τ(s)||2)ds)]γ  [exp(T0(2r(s)5l(s))dsT0||σ1(s)τ(s)||2)ds)1]γ2.
    $

    Assume that

    $ \int_{0}^{T}{(5l(s)-2r(s))ds-\int_{0}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})d}s} $
    $ =\int_{0}^{T}{[4||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}}-2r(s)]ds+5{{a}^{2}}{{\eta }^{2}}T/{{b}^{2}}:=\rho >0} $ (41)

    hold, using remark 2.1, we see that the minimum $Var[X(T)]$ is achieved for

    $ {{\gamma }^{*}}=-\frac{K+\frac{A}{2}\exp \left( -\int_{0}^{T}{(2r(s)-4l(s))ds-\frac{1}{2}\int_{0}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}}d}s} \right)}{{{e}^{\rho }}-1} $ (42)

    and

    $\mathop {\min }\limits_{s \in A({x_0})} {\mkern 1mu} Var[X(T)] = \frac{{{e^\rho }}}{{{e^\rho } - 1}} \cdot {\left[ {K + \frac{A}{2}\exp \left( { - \int_0^T {l(s)ds + \frac{1}{2}\int_0^T {||{\sigma ^{ - 1}}(s)\tau (s)|{|^2})d} s} } \right)} \right]^2}. $ (43)

    Now we state the conclusion of the effective frontier of our mean-variance problem as follows.

    Proposition 3.3. Let $K={{E}_{0}}[X(T)]$, $\tau (t)=({{\mu }_{1}}(t)-r(t), \cdots, {{\mu }_{m}}(t)-r(t))'$, $l(t)={{a}^{2}}{{\eta }^{2}}/{{b}^{2}}+||{{\sigma }^{-1}}(t)\tau (t)|{{|}^{2}}$, $A=H(0){{x}_{0}}-\overset{-}{\mathop{h}}\, (0)$. Assume that (41) hold. The effective frontier of the mean-variance reinsurance-investment strategy for problem (8) is given by

    $ Var[X(T)]=\frac{{{e}^{\rho }}}{{{e}^{\rho }}-1}\cdot {{\left[ K+\frac{A}{2}\exp \left( -\int_{0}^{T}{l(s)ds+\frac{1}{2}\int_{0}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})d}s} \right) \right]}^{2}}. $ (44)

    4. Numerical Example

    Let ${{x}_{0}}=1$, $r(t)=0.06$, $m=3$, ${{\mu }_{1}}(t)=0.08$, ${{\mu }_{2}}(t)=-0.10$, ${{\mu }_{3}}(t)=0.15$, $a=1$, $b=0.5$, $\theta =0.2$, $\eta =0.3$, and

    $ \sigma (t)=\left( 0.200.150.060.100.300.250.400.050.50
    \right). $

    Then

    $ {{\sigma }^{-1}}(t)=\left( 2.8701.3570.9481.9712.1570.7652.1740.8701.304
    \right),\tau (t)=(0.02,0.04,0.09)', $

    and

    $ ||{{\sigma }^{-1}}(t)\tau (t)|{{|}^{2}}=0.022865,{{(\sigma (t)\sigma (t)')}^{-1}}=\left( 16.84801.53454.04781.53457.25101.80214.04781.80213.1843
    \right). $

    Hence $l(t)={{a}^{2}}{{\eta }^{2}}/{{b}^{2}}+||{{\sigma }^{-1}}(t)\tau (t)|{{|}^{2}}=0.382865$. In view of (28) and (33), and let $e\approx 2.718$, direct calculation gives $h(t)/H(t)=0.0732-0.9268\times {{1.9794}^{T-t}}$. According to proposition 3.2, the mean-variance optimal reinsurance-investment strategy are given by (34), that is,

    $ \left\{ q(t)=aηb2(X(t)+h(t)H(t))+1=1+1.2(X(t)+0.07320.9268×1.9794Tt)    =1.0878+1.2X(t)1.1122×1.9794Tt;π(t)=(σ(t)σ(t))1τ(t)(X(t)+h(t)H(t))    =(0.63990.09720.2955)(X(t)++0.07320.9268×1.9794Tt)    =(0.04680.00710.0216)(0.63990.09720.2955)X(t)+(0.59310.09010.2739)×1.9794Tt,
    \right. $

    Remark that

    $\rho =\int_{0}^{T}{[4||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}}-2r(s)]ds+5{{a}^{2}}{{\eta }^{2}}T/{{b}^{2}}=1.77146T>0}$, $A=H(0){{x}_{0}}-\overset{-}{\mathop{h}}\, (0)=0.1308-2.5858\times {{1.3653}^{-T}}+0.4550\times {{1.4498}^{T}}$. In view of (44), the effective frontier of the mean-variance reinsurance-investment strategy is given by

    $ Var[X(T)]=\frac{{{e}^{\rho }}}{{{e}^{\rho }}-1}\cdot {{\left[ K+\frac{A}{2}\exp \left( -\int_{0}^{T}{l(s)ds+\frac{1}{2}\int_{0}^{T}{||{{\sigma }^{-1}}(s)\tau (s)|{{|}^{2}})d}s} \right) \right]}^{2}} $
    $ {{\left. =\frac{{{5.8784}^{T}}}{{{5.8784}^{T}}-1}\cdot \left( 0.0654\times {{1.4498}^{T}}+0.2775\times {{2.1019}^{T}} \right.-1.2929\times {{1.0619}^{T}}+{{E}_{0}}[X(T)] \right)}^{2}}. $

    Furthermore, let $T=1$, we have

    $ Var[X(1)]=1.2050\cdot {{\left( -0.6948+{{E}_{0}}[X(1)] \right)}^{2}}. $

    Acknowledgments

    This work is supported by the Foundation of Guangzhou Academy of International Finance (No: 16GFR01A04, 16GFR01A07).


    Conflict of Interest

    All authors declare no conflicts of interest in this paper.


    [1] Reza MT, Lynam JG, Vasquez VR, et al. (2012) Pelletization of Biochar from Hydrothermally Carbonized Wood. Environ Prog Sustain 31: 225-234. doi: 10.1002/ep.11615
    [2] Reza MT, Lynam JG, Uddin MH, et al. (2013) Hydrothermal carbonization: Fate of inorganics. Biomass Bioenerg 49: 86-94. doi: 10.1016/j.biombioe.2012.12.004
    [3] Reza MT, Borrego AG, Wirth B (2014) Optical texture of hydrochar from maize silage and maize silage digestate. Int J of Coal Geology 134–135: 74-79.
    [4] Reza MT, Becker W, Sachsenheimer K, et al. (2014) Hydrothermal carbonization (HTC): Near infrared spectroscopy and partial least-squares regression for determination of selective components in HTC solid and liquid products derived from maize silage. Bioresource Technol 161: 91-101. doi: 10.1016/j.biortech.2014.03.008
    [5] Lynam J, Reza MT, Yan W, et al. (2014) Hydrothermal carbonization of various lignocellulosic biomass. Biomass Conv Bioref : 1-9.
    [6] Coronella C, Lynam J, Reza MT, et al. (2014) Hydrothermal Carbonization of Lignocellulosic Biomass. In: Jin F, editor. Application of Hydrothermal Reactions to Biomass Conversion: Springer Berlin Heidelberg: 275-311.
    [7] Liu ZG, Quek A, Parshetti G, et al. (2013) A study of nitrogen conversion and polycyclic aromatic hydrocarbon (PAH) emissions during hydrochar-lignite co-pyrolysis. Appl Energ 108: 74-81. doi: 10.1016/j.apenergy.2013.03.012
    [8] Reza MT, Andert J, Wirth B, et al. (2014) Hydrothermal Carbonization of Biomass for Energy and Crop Production. Applied Bioenergy 1: 11-29.
    [9] Demir-Cakan R, Makowski P, Antonietti M, et al. (2010) Hydrothermal synthesis of imidazole functionalized carbon spheres and their application in catalysis. Catal Today 150: 115-118. doi: 10.1016/j.cattod.2009.05.003
    [10] Bandura AV, Lvov SN (2006) The ionization constant of water over wide ranges of temperature and density. J Phys Chem Ref Data 35: 15-30. doi: 10.1063/1.1928231
    [11] Reza MT, Yan W, Uddin MH, et al. (2013) Reaction kinetics of hydrothermal carbonization of loblolly pine. Bioresource Technol 139: 161-169. doi: 10.1016/j.biortech.2013.04.028
    [12] Reza MT, Uddin MH, Lynam J, et al. (2014) Hydrothermal carbonization of loblolly pine: reaction chemistry and water balance. Biomass Conv Bioref 4: 311-321. doi: 10.1007/s13399-014-0115-9
    [13] Wiedner K, Naisse C, Rumpel C, et al. (2013) Chemical modification of biomass residues during hydrothermal carbonization - What makes the difference, temperature or feedstock? Org Geochem 54: 91-100. doi: 10.1016/j.orggeochem.2012.10.006
    [14] Funke A (2015) Fate of Plant Available Nutrients during Hydrothermal Carbonization of Digestate. Chemie Ingenieur Technik 87: 1713-1719. doi: 10.1002/cite.201400182
    [15] Heilmann SM, Molde JS, Timler JG, et al. (2014) Phosphorus Reclamation through Hydrothermal Carbonization of Animal Manures. Environ Sci Technol 48: 10323-10329. doi: 10.1021/es501872k
    [16] Funke A, Mumme J, Koon M, et al. (2013) Cascaded production of biogas and hydrochar from wheat straw: Energetic potential and recovery of carbon and plant nutrients. Biomass Bioenerg 58: 229-237.
    [17] Heilmann SM, Davis HT, Jader LR, et al. (2010) Hydrothermal carbonization of microalgae. Biomass Bioenerg 34: 875-882. doi: 10.1016/j.biombioe.2010.01.032
    [18] Sevilla M, Fuertes AB (2009) The production of carbon materials by hydrothermal carbonization of cellulose. Carbon 47: 2281-2289. doi: 10.1016/j.carbon.2009.04.026
    [19] Funke A, Ziegler F (2010) Hydrothermal carbonization of biomass: A summary and discussion of chemical mechanisms for process engineering. Biofuel Bioprod Bior 4: 160-1677. doi: 10.1002/bbb.198
    [20] Kruse A, Funke A, Titirici MM (2013) Hydrothermal conversion of biomass to fuels and energetic materials. Curr Opin Chem Biol 17: 515-521. doi: 10.1016/j.cbpa.2013.05.004
    [21] Zhu XD, Liu YC, Luo G, et al. (2014) Facile Fabrication of Magnetic Carbon Composites from Hydrochar via Simultaneous Activation and Magnetization for Triclosan Adsorption. Environ Sci Technol 48: 5840-5848.
    [22] Kumar S, Loganathan VA, Gupta RB, et al. (2011) An Assessment of U(VI) removal from groundwater using biochar produced from hydrothermal carbonization. J Environ Manage 92: 2504-2512. doi: 10.1016/j.jenvman.2011.05.013
    [23] Regmi P, Moscoso JLG, Kumar S, et al. (2012) Removal of copper and cadmium from aqueous solution using switchgrass biochar produced via hydrothermal carbonization process. J Environ Manage 109: 61-69. doi: 10.1016/j.jenvman.2012.04.047
    [24] Hu B, Wang K, Wu LH, et al. (2010) Engineering Carbon Materials from the Hydrothermal Carbonization Process of Biomass. Adv Mater 22: 813-828. doi: 10.1002/adma.200902812
    [25] Demir-Cakan R, Baccile N, Antonietti M, et al. (2009) Carboxylate-Rich Carbonaceous Materials via One-Step Hydrothermal Carbonization of Glucose in the Presence of Acrylic Acid. Chem Mater 21: 484-490. doi: 10.1021/cm802141h
    [26] Baccile N, Laurent G, Babonneau F, et al. (2009) Structural Characterization of Hydrothermal Carbon Spheres by Advanced Solid-State MAS C-13 NMR Investigations. J Phys Chem C 113: 9644-9654.
    [27] Baccile N, Antonietti M, Titirici MM (2010) One-Step Hydrothermal Synthesis of Nitrogen-Doped Nanocarbons: Albumine Directing the Carbonization of Glucose. Chemsuschem 3: 246-253. doi: 10.1002/cssc.200900124
    [28] Zhao L, Bacsik Z, Hedin N, et al. (2010) Carbon Dioxide Capture on Amine-Rich Carbonaceous Materials Derived from Glucose. Chemsuschem 3: 840-845. doi: 10.1002/cssc.201000044
    [29] Fechler N, Wohlgemuth SA, Jaker P, et al. (2013) Salt and sugar: direct synthesis of high surface area carbon materials at low temperatures via hydrothermal carbonization of glucose under hypersaline conditions. J Mater Chem A 1: 9418-9421. doi: 10.1039/c3ta10674h
    [30] Reza MT, Rottler E, Tolle R, et al. (2015) Production, characterization, and biogas application of magnetic hydrochar from cellulose. Bioresource Technol 186: 34-43. doi: 10.1016/j.biortech.2015.03.044
    [31] Sevilla M, Fuertes AB (2009) Chemical and Structural Properties of Carbonaceous Products Obtained by Hydrothermal Carbonization of Saccharides. Chem-Eur J 15: 4195-4203. doi: 10.1002/chem.200802097
    [32] Titirici MM, Antonietti M (2010) Chemistry and materials options of sustainable carbon materials made by hydrothermal carbonization. Chem Soc Rev 39: 103-116. doi: 10.1039/B819318P
    [33] Sevilla M, Fuertes AB (2009) Chemical and structural properties of carbonaceous products obtained by hydrothermal carbonization of saccharides. Chemistry 15: 4195-4203. doi: 10.1002/chem.200802097
    [34] Cao XD, Ma LN, Gao B, et al. (2009) Dairy-Manure Derived Biochar Effectively Sorbs Lead and Atrazine. Environ Sci Technol 43: 3285-3291. doi: 10.1021/es803092k
    [35] Diakite M, Paul A, Jager C, et al. (2013) Chemical and morphological changes in hydrochars derived from microcrystalline cellulose and investigated by chromatographic, spectroscopic and adsorption techniques. Bioresource Technol 150: 98-105. doi: 10.1016/j.biortech.2013.09.129
    [36] Si Y, Ren T, Li Y, et al. (2012) Fabrication of magnetic polybenzoxazine-based carbon nanofibers with Fe3O4 inclusions with a hierarchical porous structure for water treatment. Carbon 50: 5176-5185.
    [37] Lynam JG, Reza MT, Vasquez VR, et al. (2012) Effect of salt addition on hydrothermal carbonization of lignocellulosic biomass. Fuel 99: 271-273. doi: 10.1016/j.fuel.2012.04.035
    [38] Ryu J, Suh YW, Suh DJ, et al. (2010) Hydrothermal preparation of carbon microspheres from mono-saccharides and phenolic compounds. Carbon 48: 1990-1998. doi: 10.1016/j.carbon.2010.02.006
    [39] Peterson AA, Vogel F, Lachance RP, et al. (2008) Thermochemical biofuel production in hydrothermal media: A review of sub- and supercritical water technologies. Energy Env Sci 1: 32-65.
    [40] Danso-Boateng E, Shama G, Wheatley AD, et al. (2008) Hydrothermal carbonisation of sewage sludge: Effect of process conditions on product characteristics and methane production. Bioresource technology 177: 318-327.
  • This article has been cited by:

    1. Yue Liu, Dehao Ruan, Stability of a Class of Impulsive Neutral Stochastic Functional Partial Differential Equations, 2020, 2020, 1026-0226, 1, 10.1155/2020/9051396
    2. Xia Zhao, Mengjie Li, Qinrui Si, Optimal investment-reinsurance strategy with derivatives trading under the joint interests of an insurer and a reinsurer, 2022, 30, 2688-1594, 4619, 10.3934/era.2022234
  • Reader Comments
  • © 2016 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(7538) PDF downloads(1356) Cited by(13)

Article outline

Figures and Tables

Figures(7)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog